Systems acclimation to osmotic stress in zygnematophyte cells

Systems acclimation to osmotic stress in zygnematophyte cells | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Systems acclimation to osmotic stress in zygnematophyte cells View ORCID Profile Jaccoline M.S. Zegers , View ORCID Profile Lukas Pfeifer , View ORCID Profile Tatyana Darienko , View ORCID Profile Kerstin Schmitt , Craig A. Dorsch , View ORCID Profile Cäcilia F. Kunz , View ORCID Profile Gerhard H. Braus , View ORCID Profile Birgit Classen , View ORCID Profile Oliver Valerius , View ORCID Profile Jan de Vries doi: https://doi.org/10.1101/2025.02.06.636795 Jaccoline M.S. Zegers 1 —University of Göttingen, Institute for Microbiology and Genetics, Department of Applied Bioinformatics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jaccoline M.S. Zegers Lukas Pfeifer 2 —Christian-Albrechts-University of Kiel, Pharmaceutical Institute, Department of Pharmaceutical Biology , 24118 Kiel, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Lukas Pfeifer Tatyana Darienko 1 —University of Göttingen, Institute for Microbiology and Genetics, Department of Applied Bioinformatics , 37077 Göttingen, Germany 3 —University of Göttingen, Albrecht-von-Haller-Institute for Plant Sciences, Experimental Phycology and Culture Collection of Algae at Göttingen , 37073 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Tatyana Darienko Kerstin Schmitt 4 —University of Göttingen, Institute of Microbiology and Genetics, Department of Molecular Microbiology and Genetics , 37077 Göttingen, Germany 5 —University of Göttingen, Göttingen Center for Molecular Biosciences (GZMB), Service Unit LCMS Protein Analytics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Kerstin Schmitt Craig A. Dorsch 1 —University of Göttingen, Institute for Microbiology and Genetics, Department of Applied Bioinformatics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Cäcilia F. Kunz 1 —University of Göttingen, Institute for Microbiology and Genetics, Department of Applied Bioinformatics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Cäcilia F. Kunz Gerhard H. Braus 4 —University of Göttingen, Institute of Microbiology and Genetics, Department of Molecular Microbiology and Genetics , 37077 Göttingen, Germany 5 —University of Göttingen, Göttingen Center for Molecular Biosciences (GZMB), Service Unit LCMS Protein Analytics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Gerhard H. Braus Birgit Classen 2 —Christian-Albrechts-University of Kiel, Pharmaceutical Institute, Department of Pharmaceutical Biology , 24118 Kiel, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Birgit Classen Oliver Valerius 4 —University of Göttingen, Institute of Microbiology and Genetics, Department of Molecular Microbiology and Genetics , 37077 Göttingen, Germany 5 —University of Göttingen, Göttingen Center for Molecular Biosciences (GZMB), Service Unit LCMS Protein Analytics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Oliver Valerius Jan de Vries 1 —University of Göttingen, Institute for Microbiology and Genetics, Department of Applied Bioinformatics , 37077 Göttingen, Germany 6 —University of Göttingen, Göttingen Center for Molecular Biosciences (GZMB), Department of Applied Bioinformatics , 37077 Göttingen, Germany 7 —University of Göttingen, Campus Institute Data Science (CIDAS), Department of Applied Bioinformatics , 37077 Göttingen, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jan de Vries For correspondence: devries.jan{at}uni-goettingen.de Abstract Full Text Info/History Metrics Supplementary material Preview PDF ABSTRACT Zygnematophytes are the closest algal relatives of land plants. They hold key information to infer how the earliest land plants overcame the barrage of terrestrial stressors, prime of which is osmotic stress. Here, we applied two osmotic stressors on a unicellular and a multicellular representative of zygnematophytes and studied their response over a 25-hour time course generating 130, 60, and 30 of transcriptomic, proteomic, and metabolomic samples combined with photophysiology, sugar analysis, immunocytochemical glycoprotein analysis, and microscopy. Our data highlight a shared protein chassis that shows divergent responses with the same outcome: successful acclimation to osmotic challenges. We establish a model of how the algal sisters of land plants can overcome a prime stressor in the terrestrial habitat and highlight components of the plant terrestrialization toolkit. INTRODUCTION Phylogenomic research over the last decade has consistently recovered the Zygnematophyceae as the closest algal relatives of embryophytes (land plants) 1 – 5 . Zygnematophyceae diversified into a species-rich class of unicellular and filamentous algae that live in diverse habitats 6 , 7 . Studying their genomes 8 – 11 revealed a whole set of genes that are homologous to the molecular chassis known from environmental adaptation in land plants. Here the focus was particularly on signaling cascades 12 – 17 . Yet, our knowledge on the proteins that constitute the stress physiology of zygnematophyte cells is very limited. Here, we focused on osmotic stress, the prime challenge on land. Land plants exhibit diverse strategies for managing changes in water availability. Vascular plants are homoiohydric, actively maintaining water homeostasis, while non-vascular plants and algae are generally poikilohydric, passively equilibrating their water status with their surroundings. However, even poikilohydric organisms regulate water fluxes to some extent, reducing their dependence on environmental conditions. For example, bryophytes use waxy cuticles to minimize water loss, while flagellated freshwater algae expel excess water through contractile vacuoles. We can assume that the earliest land plants were poikilohydric 13 . Our study delves into poikilohydry in the closest algal relatives of land plants. At the cellular level, water regulation is vital for maintaining structural integrity and function. Plant and algal cells possess cell walls, onto which intracellular fluid mediates an outward force known as turgor pressure, critical for maintaining cell structure and function. Under osmotic stress—such as when osmolyte concentrations in the surrounding environment increase due to water loss—turgor pressure decreases as water movement into the cell diminishes or reverses. This pressure drop triggers a cascade of physiological and biochemical responses, including structural alterations, membrane leakage, and other forms of cellular damage 18 , 19 . To mitigate the effects of osmotic stress, plants have evolved several adaptive mechanisms. A palpable response well-known from tracheophytes is the closure of stomata to minimize water loss 20 , 21 . On the cellular level, plants and algae initiate membrane reorganization and synthesize osmotically active metabolites, such as proline, sugars, and sugar alcohols, to stabilize the intracellular osmotic potential 22 . When dehydration results from osmotic stress, some algae can down-regulate the water-intensive process that is photosynthesis 23 . Reactive oxygen species (ROS) are also commonly produced during drought-related stresses, being both signals and culprits emerging from an imbalance in photosystems 21 , 24 – 26 . In response, plants activate ROS-scavenging systems to prevent oxidative damage. Plants further produce a variety of stress-responsive proteins that play essential roles in protecting cellular membranes, maintaining ionic homeostasis, and degrading damaged proteins to facilitate repair 18 , 27 . Conserved structural modifications to the cell wall and cytoskeleton 28 bolster cellular strength and resilience to the mechanical stresses associated with turgor pressure loss. Collectively, these mechanisms represent a multifaceted strategy that allows plants to survive and adapt to osmotic stress conditions. Some Zygnematophyceae exhibit exceptional adaptability to osmotic stress, thriving in habitats subject to frequent desiccation 29 – 31 . Previous studies, often with a focus on one specific technical aspect, have identified key coping strategies and osmotic stress-responsive genes in these algae; these include upregulated glycosyltransferase family 2 genes in Zygnema circumcarinatum 32 and microscopic observations of cell wall and endomembrane adaptations in Penium margaritaceum 33 . Here, we employ a systems biology approach of temporal transcriptomics and proteomics to uncover the core osmotic stress response in Zygnematophyceae. From our data, we infer key components of the plant terrestrialization toolkit. RESULTS Quantitative impacts of osmotic stress in Mesotaenium and Zygnema We worked with two representatives of the closest algal relatives of land plants: the unicellular Mesotaenium endlicherianum SAG12.97 and the filamentous Zygnema circumcarinatum SAG698-1, hereafter referred to as Mesotaenium and Zygnema . Next to representing the two predominant body plans in zygnematophytes, these species were selected due to (i) their natural exposure to osmotic stress in their habitats: Mesotaenium originates from a lake subject to annual drying cycles, while Zygnema is found in a meadow ditch; (ii) the availability of fully sequenced genomes and updated gene models 8 , 10 , 16 ; and (iii) the presence of all protein components involved in the canonical abscisic acid (ABA) signaling cascade 10 , 34 , critical in osmotic stress response of embryophytes 35 , 36 —albeit these proteins likely act independent of ABA in zygnematophytes 14 . Mesotaenium and Zygnema were cultivated on cellophane disks placed on agarised medium, allowing transfer to other agar plates without significantly disturbing the algae ( Fig. 1a ). Unlike cultivation in liquid medium, this approach allowed for stress experiments that more closely mimic terrestrial habitats while enabling easy and rapid sampling. The osmotic stress treatments for this study were initially based on Kilian et al. 37 : 300 mM mannitol and 150 mM NaCl for 24 hours. However, during preliminary trials, no clear morphological response was observed under 300 mM mannitol within 24 hours. Consequently, the mannitol concentration was increased to 800 mM. In contrast, treatment with 150 mM NaCl (both osmotic and ionic stress) led to a clear and rapid drop in the maximum quantum yield of photosystem II, Fv/Fm, that is a non-invasive and comparable measure of gross physiology 38 . Download figure Open in new tab Figure 1. Physiological impact of osmotic stress in two zygnematophytes A) Cladogram of the Streptophyta, highlighting the two zygnematophytes (Zyg.) used and schematic overview of experimental setup, depicting with a dot when during the stress experiment samples were taken for which analysis. “M” stands for Mesotaenium only. B) Weight per plate and water content measured after 25 hours of osmotic stress. C) Relative photosynthetic efficiency (Fv/Fm) estimated from chlorophyll fluorescence measurements. D+E) Microscopic features observed under light microscopy after 25 hours of osmotic stress in Mesotaenium (D) and Zygnema (E). The highest Fv/Fm values were observed 25 hours after transferring the algae to a new growth plate without osmolytes, likely reflecting optimal conditions provided by the fresh medium. In this treatment, Fv/Fm values ranged from 0.6 to 0.73 for Mesotaenium and 0.7 to 0.76 for Zygnema ( Fig. 1b ). Following the onset of osmotic stress, Fv/Fm values decreased in both algae ( P <0.01 for both treatments after 3h, Fig. 1b ). However, recovery patterns differed: Zygnema showed a rapid rebound in the NaCl treatment, whereas Mesotaenium exhibited a slower recovery. In contrast, Fv/Fm values remained 10% below their starting values under mannitol treatment for both species, though they never approached zero, indicating that the algae retained some potential for photosynthesis. Water moves from regions of higher concentration to regions of lower concentration; upon osmotic stress, we expected an efflux and/or a reduced influx into the cell. The water content in both species was significantly ( P < 0.001) reduced upon NaCl and mannitol treatments ( Fig. 1a,c ). Notably, the dry weight of the mannitol-treated samples was substantially higher compared to the NaCl-treated and control samples but this is explained by an accumulation mannitol in the samples as later confirmed through metabolite fingerprinting (see below). A rapid decrease of water content triggers plasmolysis in plants: the movement of water out of the cell decreases the volume of the protoplast, which could detach from the rigid structure of the cell wall as also previously well-described for Zygnema circumcarinatum 39 . Although plasmolysis was not complete in all cells, with the plasma membrane remaining connected to the cell wall in most areas, it was the most significant feature observed during the morphological analysis of both Mesotaenium and Zygnema under high osmotic stress induced by mannitol treatment ( Figs. 1c,d ). In addition to plasmolysis, other structural changes included distortions of the chloroplast, bending and breaking of cells/filaments, changes in vacuolar structures and other inclusions ( Figs. 1c,d ). In Mesotaenium “cups” were formed, which are non-green and often reddish cell-wall-enclosed compartments that can detach from the mother cell. Dynamics in stress-responsive proteins and transcripts To investigate the molecular response of the algae, we applied RNA sequencing at t0, 3h, 6h, 9h, and 25h and shotgun proteomics at t0, 6h, and 25h. Accurate transcript levels were determined for 21,651 transcripts in Mesotaenium and 9,631 transcripts in Zygnema . Protein levels were reliably quantified for 2,171 and 1,713 proteins, respectively, 88% and 92% of which were detected in 100% of samples for at least one treatment, enabling their inclusion in differential abundance analyses. Further, a de novo transcriptome was assembled as a complementary database for proteomics analysis, yielding 90 and 200 additional protein groups in Mesotaenium and Zygnema , respectively. In Mesotaenium , many of these newly identified proteins corresponded to isoforms of incomplete genes located at the ends of contigs, while in Zygnema , they were predominantly new genes already annotated in Zygnema circumcarinatum UTEX 1559 (Supplementary Fig. I). Principal component analysis (PCA, Fig. 2a ) revealed that proteins from different treatments were better separated at the 25-hour timepoint compared to the 6-hour timepoint. The PCA of the transcript data provided further insights, highlighting distinct differences between Mesotaenium and Zygnema : In Mesotaenium , controls showed minimal variation across timepoints, whereas Zygnema exhibited clear shifts along the second principal component over time. Each treatment displayed unique temporal patterns: in Mesotaenium , early NaCl-treated samples were farthest from controls but later converged, indicating a rapid and strong response. In contrast, the transcript PCA for Mesotaenium revealed greater variability among biological replicates in the latest timepoints of the mannitol treatment, with some samples clustering near the controls while others diverged substantially, suggesting a complex acclimation process. Download figure Open in new tab Figure 2. Differential proteome and transcriptome analysis in Mesotaenium and Zygnema A) Principal component analysis (PCA) plots of log2-transformed label-free quantification (LFQ) values from the proteomics dataset and transcripts per million (TPM) values from the RNA-seq dataset. B) Log2 fold changes in protein abundance for the 0.8 M mannitol treatment relative to the control, plotted against those for the 0.15 M NaCl treatment relative to the control. The proteins in bold are: AKR = NADPH-dependent Aldo-Keto Reductase, EXL = Exordium-like protein, FAS1 = protein with at least one Fasciclin-like 1 domain, SVR3 = putative chloroplast TypA translation elongation GTPase, XTH = Xyloglucan endotransglucosylase/hydrolase. The full names of the proteins can be found in supplementary data I. C) Log2 fold changes in transcript abundance for the 0.8 M mannitol treatment relative to the control, plotted against those for the 0.15 M NaCl treatment relative to the control. We performed a differential abundance analysis, comparing each treatment to the control for all timepoints ( Figs. 2b,c ). Consistent with the PCA results, protein fold changes were more pronounced at 25h than at 6h while transcript changes were more extreme at 6h than at 25h, indicating a faster transcriptomic response compared to the proteomic response. An intriguing difference was observed between Zygnema and Mesotaenium : in Mesotaenium , the fold changes for mannitol and NaCl treatments were similar, whereas in Zygnema , the mannitol treatment exhibited more pronounced fold changes than NaCl. These trends were consistent across both transcriptomic and proteomic levels. Several proteins changed in differential abundance (Supplementary Data I). Focusing on proteins identified in both Mesotaenium and Zygnema , five types were found. Among these, SVR3, a putative chloroplast TypA translation elongation GTPase, was the only protein downregulated. SVR3 has been implicated in chloroplast development under cold stress in plants 40 , but its broader functional roles remain largely unexplored. Among the upregulated proteins, one of the most differentially abundant during osmotic stress was an NADPH-dependent aldo-keto reductase (AKR). AKR catalyzes the reduction of carbonyl groups to primary and secondary alcohols. It is frequently associated with abiotic stress responses, functioning in detoxification, osmolyte production, or the synthesis of secondary metabolites 41 . The remaining three upregulated proteins are predicted to localize in the cell wall: an exordium-like protein (EXL), xyloglucan endotransglucosylase/hydrolase (XTH), and proteins containing at least one fasciclin-like (FAS1) domain. XTH likely contributes to cell wall remodeling by modifying xyloglucan, whereas the functions of the other cell wall-associated proteins are less well understood. The number of transcripts showing significant changes is substantial (>4,000 with P adj < 0.01), making it impractical to analyze all of them individually. Focusing on the top five differentially expressed transcripts, we observe that both Zygnema and Mesotaenium share a common upregulated gene encoding a LEA (Late Embryogenesis Abundant) protein. LEA proteins play a crucial protective role during environmental stress in plants and are frequently upregulated in response to drought- and osmotic-related stressors across various species 42 . Integrated analysis of transcript and protein correlations Our parallelized experimental setup—with both transcriptome and proteome samples originating from the same plate—enabled us to perform robust comparison between transcript and protein levels. When examining how osmotic-stress-induced changes in specific transcript–protein pairs correlate ( Figs. 3a , c, Extended Data Fig. 1, and Supplementary Figure II), we observe a strong correlation in Mesotaenium (R=0.63) and a moderate correlation in Zygnema (R=0.39). The correlations improve on average by 1.8x when using 25h protein samples compared to 6h samples, a difference that exceeds the variability observed across transcript timepoints (maximum 1.2x). In this study, conserved stress responders such as AKR (both species), ELIP ( Mesotaenium ), and heat shock proteins (HSP, Zygnema ) exhibited the strongest correlations. For both species, several FAS1 proteins displayed a reversed pattern: higher protein abundance under osmotic stress conditions but lower transcript levels. This pattern was also observed for the photosystem II subunit S (PsbS, Mesotaenium), ribosomal protein S5B (RPS5B, Zygnema ) and expansin ( Zygnema ), suggesting the involvement of post-transcriptional regulatory mechanisms. Download figure Open in new tab Figure 3. Relationship between RNA and protein abundance during osmotic stress A) Relative protein abundance at 25 hours compared to relative RNA abundance (integrated from 0 to 25 hours) for proteins showing a significant fold change (FDR < 0.05) in either salt or mannitol treatments relative to the control. B) Absolute protein abundance in Mesotaenium at 0, 6, and 25 hours plotted against absolute RNA levels at corresponding time points and replicates, with data points colored according to the size of the contig from which the protein-coding gene originates. C) Pearson correlation coefficients between relative (log2 fold change) protein abundance at 6 or 25 hours and corresponding relative RNA abundance. D) Pearson correlation coefficients between absolute protein and RNA abundance at 6 or 25 hours. Examining the same data from a different perspective, we compared absolute transcript and protein levels ( Figs. 3b , 3d , and Supplementary Figure II), finding a moderate correlation for Zygnema (R≥0.34), except under mannitol treatment at 25h, which suggests the potential involvement of post-transcriptional regulatory mechanisms. For Mesotaenium , the initial correlation between absolute protein and transcript levels was overall poor (R≤0.28). Upon further examination, we found that many proteins with high abundance but unexpectedly low transcript levels are encoded on extremely small contigs (<4096 bp) or the plastome. Excluding genes located on small contigs significantly improves the correlation, surpassing that observed in Zygnema . Temporal clustering unveils conserved osmotic stress responses To delve deeper into the biological significance of the proteomic and transcriptomic datasets, we first divided the proteomic dataset into quadrants ( Fig. 4a , Methods). Functional enrichment analysis yielded results consistent with the manual annotation performed for the statistically relevant genes. To cluster the substantially larger transcriptomic dataset we grouped temporal patterns for each transcript while retaining the distinctions between the treatment conditions. Temporal patterns were first grouped using Gaussian Mixture Model (GMM) clustering; groups were normalized and further clustered using the same algorithm to reduce complexity. This approach resulted in 12 and 9 temporal modules for Mesotaenium and Zygnema , respectively, which were then subjected to functional enrichment analysis ( Fig. 4b ). In Mesotaenium , most modules exhibited a control pattern that remained consistently flat. Conversely, in Zygnema , transcript levels in the control treatment fluctuated markedly over time. Interestingly, in Zygnema , the temporal patterns of the salt treatment often appeared intermediate between those of the control and mannitol treatments. For Mesotaenium , the largest modules displayed similar response patterns for both mannitol and salt treatments, albeit with slightly faster responses observed in the salt treatment. Yet, the smaller modules also revealed treatment-specific responses unique to the mannitol and salt conditions. Download figure Open in new tab Figure 4. Functions enriched in two zyngematophytes challenged with osmotic stress Functional Enrichment of GO Terms and InterPro Domains During Osmotic Stress (FDR < 0.01). Top) Proteins were grouped into quadrants based on log2 fold changes (≥1 or ≤-1), highlighting significant changes in response to 0.15 M NaCl and 0.8 M mannitol treatments. Functional enrichment analysis was performed on each group to identify enriched GO terms and InterPro domains. Bottom) Transcript temporal changes were clustered into patterns using a Gaussian Mixture Model, with normalized patterns further organized into modules. Functional enrichment analysis identified enriched terms for each module, with the top nine displayed. A total of 670 and 320 enriched terms ( p < 0.01) were identified for Mesotaenium and Zygnema , respectively. Intersection network analysis ( Fig. 5 ) revealed several highly connected terms, including oxidoreductase activity, membrane, and chlorophyll binding. These terms are functionally broad and are represented by both upregulated and downregulated proteins under osmotic stress. For instance, “chlorophyll binding” includes both Early Light-Induced Proteins (ELIPs), commonly stress-inducible 43 , and classical Light-Harvesting Complexes (LHCs). Download figure Open in new tab Figure 5. Integrated network of comparative proteomic and transcriptomic analysis in two zygnematophytes This network illustrates the overlap of enriched Gene Ontology (GO) terms and protein domains across transcript modules and/or protein quadrants in Zygnema and Mesotaenium . Long names have been abbreviated for clarity. Module 7 of Zygnema and Module 9 of Mesotaenium show the highest overlap, both exhibiting upregulation under osmotic stress. Uniquely shared terms are highlighted in green, with representative proteins common to both species schematically depicted in the algae models (AHA = Plasma Membrane H+-ATPase, SuSy = Sucrose Synthase, TIP = Tonoplast Intrinsic Protein, SYT = Synaptotagmin). Terms and corresponding proteins showing opposing transcriptional regulation under osmotic stress are highlighted in red (upregulated in Mesotaenium ) and blue (upregulated in Zygnema ) within the schematic algae model (HSP20 = Small Heat Shock Protein 20, PetC = Component of the Cytochrome b6-f Complex). Proteins with significantly increased or decreased abundance (see Figure 2 ) in both algae are depicted as well, in yellow and purple, respectively. Locations and functions depicted in the algae model are hypothetical. The two most interconnected modules are Module 9 of Mesotaenium and Module 7 of Zygnema . Both modules exhibit an upregulated pattern under osmotic stress ( Fig. 4b ), suggesting a conserved osmotic stress response within the Zygnematophyceae. The uniquely shared terms within these modules are represented in both Mesotaenium and Zygnema by autoinhibited plasma membrane proton (H⁺) ATPases (AHA), sucrose synthase (SuSy), tonoplast intrinsic protein (TIP), and synaptotagmin (SYT). The AHA proteins are responsible for generating an electrochemical proton gradient across the plasma membrane, a crucial process for pH regulation and transport across the membrane 44 . Sucrose synthase catalyzes the reversible reaction converting sucrose into fructose and UDP-glucose, playing a vital role in energy metabolism, osmolyte production, and complex carbohydrate biosynthesis 45 . TIPs are aquaporins located in the tonoplast—the vacuolar membrane—that facilitate the transport of water and small solutes across this membrane 46 . SYTs are membrane-associated proteins with an N-terminal domain located in a membrane, e.g. the At SYT1 is located in the ER-membrane; their C-terminal domain binds to other membranes, enabling them to bring membranes together, organize membrane structure, and facilitate lipid transfer 47 . Several terms merit closer examination, particularly two that exhibit opposing regulation under osmotic stress in Mesotaenium and Zygnema : the Hsp20 domain and 2Fe-2S-binding. These terms encompass multiple proteins, but in both Mesotaenium and Zygnema , they include Heat Shock Protein Family 20 (Hsp20) and the Rieske iron-sulfur protein PetC, respectively. Hsp20 is a small heat shock protein with a protective role, while PetC, part of the cytochrome b6/f complex, is crucial for electron transport in photosynthesis. Both proteins are involved in stress response mechanisms: Hsp20 prevents the aggregation of proteins and lipids 48 , while PetC contributes to photosynthetic redox control 49 . In our experimental setup, Zygnema predominantly relies on the Hsp20 family, whereas Mesotaenium utilizes PetC, reflecting distinct strategies to cope with osmotic stress. We attempted to predict a Gene Regulatory Networks (GRNs) from temporal transcriptional data using Sliding Window Inference for Network Generation (SWING; supplementary Fig. III). While the GRNs were not conclusive (likely due to lack of additional and earlier time points), examination of differentially expressed transcripts revealed various signaling components. For instance, the ABA-responsive Element Binding Factor (ABF) was strongly upregulated in both Zygnema and Mesotaenium (Supplementary Fig. IV). Conserved feedback ah2t the cell wall under osmotic and ionic stress The cell wall is the first barrier of a plant cell. Plants sense the integrity of their cell wall during growth and stress through mechanisms concertedly maintain cell wall integrity (CWI) 50 . In previous work, we found that homologs that could constitute a CWI sensing and signaling loop co-express in Mesotaenium 16 . Our combined proteomic and transcriptomic data point to several layers of such biophysical feedback, integrating homeostasis of the cell wall into the overall response ( Fig. 5 ). Cell wall remodeling is key during acclimation to osmotic challenges 51 . Several cell wall-associated proteins were highlighted in the overlap between both species, including expansins, xyloglucan metabolism and more ( Fig. 5 ). Zygnematophytes are rich in glycosyl hydrolases (GH) akin to land plants 10 . We screened the proteome (FDR 3) for increases in GHs. While a cellulase (GH5) was identified in the proteome, its increase was relatively modest. The most pronounced increases were observed in the XTHs and enzymes involved in starch degradation present in the NaCl treatment and the mannitol treatment. To follow up on this, we used LC-MS-based metabolite fingerprinting on Mesotaenium and identified a small set of robust features that significantly differed between treatments ( Fig. 6a , Supplementary Fig. V). These features were classified into three main categories: 1) Mannitol, which was expectedly abundant in the mannitol treatment and not detected in the other samples. 2) Sodium formate, a known artifact resulting from the addition of NaCl to biological samples 52 . 3) Carbohydrates build from aldo/keto-hexoses, predominantly detected with sodium adducts in the NaCl treatment. Among the carbohydrate-associated features, several were identified as corresponding to a disaccharide, a trisaccharide, a tetrasaccharide, and a pentasaccharide eluted in the same retention time range. When scanning specifically for glycosides in negative ion mode, several features showed a marked increase in intensity following osmotic treatment. This suggests that Mesotaenium may utilize glycosides, such as sucrose or maltose, as osmolytes. Interestingly, the feature corresponding best to proline, often seen as osmolyte in plants 53 , did not exhibit any increase upon osmotic treatment. Download figure Open in new tab Figure 6. Effect of osmotic and ionic stress on sugar and cell wall biology of zygnematophytes A) Log2 fold changes of metabolites from metabolite fingerprinting analysis of methanol-extracted samples after 9 hours of stress treatment with 0.8M mannitol versus 0.15M NaCl in positive ion mode. Identified metabolites are labeled. B) Long-term (26 days) effects of 0.15M NaCl treatment on algae morphology. Left: Mesotaenium . Right, top: Zygnema after prolonged treatment. Right, bottom: Zygnema morphology after recovery for 14 days after 14 months of 0.15M NaCl treatment. C) Domain architecture, names, and IDs of proteins significantly more abundant (FDR < 0.05) during 0.15M NaCl and/or 0.8M mannitol treatments that have predicted arabinogalactan glycosylation sites. Species abbreviations: Me = Mesotaenium , Zc = Zygnema . D) Arabinogalactan protein (AGP)-focused analysis of aqueous-extracted samples after 25 hours of 0.15M NaCl or 0.8M mannitol treatment. Top: β-Yariv radial diffusion assay. Bottom: Sugar epitope ELISA. Right: Schematic representation of an AGP chain (not representative for Mesotaenium and Zygnema ). E) Staining of Mesotaenium with β-Yariv reagent after 25 hours of 0.15M NaCl treatment, with quantification of extracellular β-Yariv-reactive material across different treatments. F) Phylogenetic tree of glycosyl hydrolase family 16. The expanded clade represents subfamily 20, which includes xyloglucan endotransglucosylase/hydrolases (XTHs). The gene IDs in bold black indicate transcripts upregulated (ΔFC > 3) during osmotic stress, while those highlighted in orange correspond to proteins that are significantly more abundant (FDR < 0.05) under osmotic stress conditions. G) Schematic representation of hypothesized major changes in cell wall components during osmotic and ionic stress in Zygnema and Mesotaenium . The respective enzymes showing significant changes at both protein and transcript levels are indicated. The depicted cell wall structure is hypothetical. In the differential abundance analysis, we noted an upregulation of cell wall-situated proteins, including the aforementioned XTHs and cell wall-associated (predicted) proteins containing a fasciclin-like (FAS1) domain. XTH is a well-characterized enzyme belonging to the GH family 16, which remodels xyloglucan chains. XTH has been frequently reported to be strongly upregulated in response to osmotic stress in land plants 54 , assumingly contributing to cell wall strengthening and thereby indirectly preventing water loss 55 . We constructed a phylogenetic framework encompassing all glycosyl hydrolases from GH family 16 ( Fig. 6f ) and recovered the XTHs (GH16, subfamily 20) exclusively for Embryophyta and Zygnematophyceae—not in any other streptophyte or chlorophyte algae; in Coleochaetophyceae, some proteins were initially annotated as GH16 subfamily 20 using dbCAN3 but were phylogenetically closer to those of subfamily 27. The XTHs from Embryophyta and Zygnematophyceae each form distinct monophyletic groups, in case of the latter further divided into two deep clades ( Fig. 6f ). Notably, all XTHs responsive to osmotic stress clustered within the same clade. The FAS1 domain is a conserved structural motif found in extracellular proteins across diverse life forms, with a notable prevalence in plants. This domain provides multiple interaction surfaces, facilitating binding to a variety of ligands. In plants, FAS1 domains occur almost exclusively in the form of fasciclin-like arabinogalactan-proteins (FLA) 56 . Their precise biochemical mechanism is obscure but likely includes interaction and signaling processes within the extracellular matrix 57 . FLAs are a subgroup of the highly O -glycosylated plant glycoproteins called arabinogalactan-proteins (AGPs) 58 . They contain one or two (or more) FAS1 domains as well as glycosylated and non-glycosylated intrinsically disordered regions 59 . We examined our significantly differentially abundant proteins for potential AG glycosylation sites and identified one FLA in Zygnema , two FLAs in Mesotaenium , as well as one pectin methylesterase with potential AG glycosylation sites ( Fig. 6c ). All can be defined as chimeric AGPs, consisting of a functional, structured protein domain and disordered regions 60 . To investigate the influence of osmotic stress on AGP glycosylation, we analyzed cellular extracts from different treatments using the β-Yariv reagent and various AGP-specific antibodies. In Mesotaenium , radial diffusion assays indicated an increase in the galactose backbone of AGPs under salt and mannitol treatments compared to the control ( Fig. 6d ). Conversely, the two AGP-specific antibodies 61 , 62 KM1 and JIM13 showed signal decrease in both treatments ( Fig. 6d ). Both treatments showed significantly elevated levels of hydroxyproline (Hyp) per glycan (Supplementary Fig. VI) corresponding to a more condensed glycan structure with reduced side-chain length while the overall galactose content remains relatively stable. The β-Yariv reagent did not stain Zygnema extracts ( Fig. 6b ), possibly due to differences in AGP glycosylation structure with higher levels of uronic acids ( Fig. 6d , Supplementary Fig. VI). Nonetheless, trends in the KM1 antibody results in Zygnema were similar to those in Mesotaenium . The discrepancy in signals from the backbone and side chains under osmotic stress suggests alterations at the post-translational level. We further scanned upregulated transcripts for glycosyl hydrolases that might hydrolyze AGP side chains. This analysis proved challenging due to the limited functional characterization of many glycosyl hydrolases and the incomplete understanding of AGP glycosylation structure in the investigated algae 63 . Among the identified candidates, GH30— a galactosidase 64 —was upregulated in both species. However, these proteins lacked signal peptides, complicating their functional interpretation. We stained live Mesotaenium with β-Yariv reagent after 25h of stress treatment ( Fig. 6e , Supplementary Fig. VII). Substantially more staining was observed inside the cells subjected to mannitol treatment. However, when staining with an unrelated dye, fuchsin red, we also noted an increase in intracellular staining, suggesting higher cell membrane permeability in general (Supplementary Fig. VII). This could be attributed to a hypotonic shock following the stress treatment. In contrast, during salt treatment, we frequently observed staining of extracellular material, indicating the release of AGPs into the extracellular space ( Fig. 6e ). This phenomenon of AGP release has previously been documented in tobacco 65 , 66 . It was hypothesized that phospholipase C could cleave the glycosylphosphatidylinositol (GPI) anchor of FLAs, thereby releasing them into the extracellular space. However, in our experiments, especially the smaller FLA isoform lacking a GPI anchor was significantly upregulated in protein abundance under salt treatment and simultaneously downregulated at the transcriptional level. Additionally, none of the predicted phospholipase C genes were upregulated, making this mechanism unlikely. To further investigate, we analyzed two modules of transcripts that were strongly upregulated exclusively during salt treatment, Module 3 and Module 10 ( Fig. 4 ), likely reflecting the ionic-specific stress response in Mesotaenium . The modules contained some genes related to photosynthesis, many uncharacterized proteins, and pectin methylesterases (PMEs) alongside a PME inhibitor. PMEs are enzymes that demethylate galacturonic acid residues, a key component of the pectin backbone in plant cell walls. Random demethylation loosens the pectin structure, enhancing enzymatic degradation while block-wise demethylation exposes free carboxyl groups on galacturonic acid residues, which can then crosslink with calcium ions to form a gel-like, compact pectin structure 67 . The PME in Mesotaenium predicted to have AG sites shares its closest Arabidopsis homolog with PME53, and thus likely removes methyl groups in a random manner 68 . Both pectin and AG groups exhibit a strong affinity for binding calcium ions. However, sodium ions, introduced in large quantities during NaCl treatment, can compete with calcium for these binding sites 69 . This competition can disrupt the calcium-mediated crosslinks between pectin molecules, potentially resulting in cell wall softening and destabilization 70 . During long-term (>26 days) NaCl treatment, we observed significant deterioration of the cell walls in Zygnema ( Fig. 6b ). While the algae survived, no growth was recorded under these conditions. In contrast, Mesotaenium showed remarkable resilience: it remained healthy, maintained active growth, and preserved chloroplasts with a proper shape and vibrant green pigmentation. Notably, some differences were observed compared to the control treatment, including the formation of small filaments up to eight cells in length and an increased accumulation of cell wall remnants (see Figure 6b and Supplementary Figure VIII). Overall, our data pinpoint deeply conserved proteins that zygnematophyte cells use in homeostatic mechanisms—that include cell wall remodeling—to successfully acclimate to osmotic stress. Discussion We pinpoint systems biological responses to osmotic stress in Mesotaenium and Zygnema . Overall, the responses differed markedly between the two algae, yet led to acclimation to osmotic challenges in both systems. Further, in both cases the responses built on a deeply conserved protein chassis. In Mesotaenium , the control treatment—serving as a mock treatment conducted alongside the tested treatments for comparison—exhibited minimal changes over time. In contrast, Zygnema showed a pronounced and temporally dynamic response under the control conditions over the course of 25h. This highlights the resilient internal and diurnal growth program of Zygnema . It seems that Zygnema prioritizes growth and rapid adaptation to improved conditions rather than immediate responses to worsening stress. This is consistent with other observations highlighting the resilience of Zygnema 30 , 71 , 72 . Another striking difference between the species was their response to salt and mannitol treatments. In Zygnema , the response to salt generally fell between the control and mannitol treatments in magnitude, consistent with the varying degrees of osmotic stress. By contrast, in Mesotaenium , the salt-induced response was often comparable to, but faster than, the response to mannitol. Moreover, Mesotaenium exhibited a distinct reaction exclusive to the salt treatment, suggesting that ionic components not only accelerate the osmotic stress response but also trigger a unique ionic-specific pathway. In Mesotaenium , many genes specifically upregulated under salt stress remain unidentified, one notable exception being PME, which modify and thereby regulate pectin structure in the cell wall. Another significant observation is the release of AGPs into the medium under salt stress, a phenomenon also reported previously in tobacco 65 , 66 . A physicochemical property of both pectin and AGPs is their potential to bind to positively charged ions. Pectin can form complexes with calcium to crosslink its own polymers 73 and AGPs, particularly their arabinogalactan (AG) domains, are known to function as calcium reservoirs due to their strong calcium-binding properties 69 , 74 . Under NaCl treatment, sodium ions could compete with calcium ions for binding sites in both calcium-pectin complexes and calcium-AG interactions ( Fig. 6g ). This competition likely drives an increase in block-wise PME activity to stabilize the cell wall structure. Furthermore, the AG regions of Mesotaenium ’s PME53 may serve as ion sensors or calcium reservoirs. The increased tendency of Mesotaenium to form short filaments under prolonged salt treatment remains unclear but may also be related to alterations in pectin-modifying enzymes, as pectin plays a crucial role in cell adhesion and the formation of multicellular complexes 67 . AGPs, particularly the FLAs, are influenced not only by ionic stress but also by other types of stress like osmotic stress. Just recently, some FLAs from Physcomitrium paten s were shown to be involved in the cellular response to rapid dehydration stress 75 . FLAs might play a structural role in the cell wall, and alterations in glycosylation patterns could potentially modify cell wall properties 76 . Alternatively, FLAs have been proposed to play significant roles in signaling 56 or to. be transported via vesicles to the vacuole under osmotic stress to foster ion transport 66 . If this mechanism occurs, FLAs might be degraded within the vacuole, thereby increasing the osmotic potential inside it and helping to restore the cell’s turgor pressure. While this is an intriguing hypothesis, it warrants further detailed research to elucidate the underlying mechanisms. The vacuole is key in osmotic stress responses. Tonoplast intrinsic proteins (TIPs), which facilitate the transport of water and small solutes, are upregulated under such conditions. Mannitol treatment appears to visibly alter the tonoplast vacuolar system in Mesotaenium . Similar changes, characterized by a transition from a simple globular vacuolar shape to a more network-like structure, have been reported in tobacco under osmotic stress 77 . In Zygnema , this shift in vacuolar structure was less evident. However, changes in chloroplast morphology were observed, indirectly suggesting vacuolar shape alterations. The star-shaped chloroplasts, seemingly anchored to the cell wall, lose their distinctive shape leading up to and during plasmolysis. In the zygnematophyte Penium , similar vacuolar shape changes were observed under extreme osmotic stress 33 . Such alterations in the vacuolar system have been hypothesized to signal the onset of programmed cell death or to serve as a means of membrane storage until conditions become more favorable 33 , 77 . We propose that these vacuolar changes may also play a role in modifying cellular structures to either prevent or prepare for plasmolysis. We identified numerous proteins that are strongly upregulated under osmotic stress, including several well-established abiotic-stress-responsive genes in plants, such as AKR, AHA, LEA, XTH, and ELIP 27 , 55 , 71 , 78 , 79 . Additionally, SuSy and enzymes from the starch degradation pathway were upregulated upon osmotic stress, potentially playing a role in osmolyte production. Consistent with findings in other streptophyte algae, previous studies have shown that sucrose synthase (SuSy) activity is significantly elevated in Klebsormidium crenulatum under desiccation stress 80 . Similarly, sucrose levels were found to be markedly increased in Chara braunii in response to salinity-induced stress 81 . Species-specific patterns include sugar-degrading enzymes, such as pullanase in Mesotaenium . Our data point to various biophysical feedback mechanisms enriched in the osmotic stress response, including SYTs, which have only infrequently been associated with osmotic stress 82 . A further link to the cell wall and growth exists through EXL homologs, which were consistently enriched and facilitate growth under challenging conditions 83 , 84 ; indeed, EXL consistently appears as a stress hub shared between land plants and zygnematophytes 16 , 72 . This goes hand in hand with the changes in OSCA, one the key mechanosensitive ion channels 85 , and a major hub in a predicted gene regulatory network shared by zygnematophyceae and land plants 72 . And the cell wall-loosening proteins expansin 86 . Overall our comparative approach describes multiple layers of biophysical feedback that concertedly acclimate cell wall and cell interior based on mechanisms shared by land plants and their algal sisters upon a prime stressor on land. Methods Cultivation of the algae Mesotaenium endlicherianum SAG 12.97 and Zygnema circumcarinatum SAG 698-1b were sourced from the Culture Collection of Algae at Göttingen University (SAG). They were cultivated on cellophane disks as previously described 71 , with the modification that Mesotaenium was grown on 1.5% agar-solidified C medium 87 for a total of 8 days at 60-80 μmol photons m −2 s −1 and Zygnema was grown on 1% agar-solidified WHM medium (MZCH recipe) for a total of 14 days at 15-20 μmol photons m −2 s −1 . Stress treatments The stress treatments were performed under the same conditions as the algae cultivation. After two hours of light exposure within their 16/8h light-dark cycle, the algae cultures were transferred to a sterile workbench. Within a maximum of 15 minutes—during which minor fluctuations in light, temperature, and humidity could occur—the cellophane discs with algae were placed onto fresh plates containing the same medium used for cultivation, supplemented with either no additive (control), 0.15 M NaCl, or 0.8 M D-mannitol. The treated algae were then returned to the cultivation chamber to resume growth. Physiological measurements At 0, 3, 6, 9, and 25 hours during the stress treatments, algae samples were incubated in the dark for 20 ± 5 minutes. Each sample was then placed on the sensor of the Walz PAM-210 Chlorophyll Fluorometer, and the maximum quantum yield of PSII (Fv/Fm) was measured. The following instrument settings were used: ML = 6, AL = 3, FR = 8, and SP = 10. For each biological replicate, two technical replicates were measured, and the mean Fv/Fm value was recorded. Algal samples used for PAM measurements were excluded from further stress experiments to prevent potential measurement interference. Microscopy At 25 hours, algae samples were transferred to a liquid medium, and microscopic images were immediately captured using a Zeiss Axioscope 7 light microscope at 400× magnification. For Mesotaenium , individual cells were counted, while for Zygnema , individual filaments were enumerated. Additionally, the presence or absence of specific traits was recorded. For each trait, the percentage of cells or filaments exhibiting the trait was calculated per biological replicate. Detailed definitions of all traits, along with the number of cells or filaments analyzed, are provided in the Supplementary Information (Supplementary Table I). For microscopic analysis using the β-Yariv reagent, treated algae (control, 0.15 M NaCl, or 0.8 M D-mannitol) were returned to the growth chamber and incubated for 25 hours. Staining was performed with β-Yariv reagent and 0.05%-1% Fuchsin stain, the latter serving as a control for nonspecific staining. Algae were then incubated with 100–200 µM β-Yariv dissolved in C-medium for 10 minutes. Imaging was conducted as described above. To ensure direct comparison, images of each biological replicate were acquired in a staggered manner over approximately 10 minutes. For quantification, additional imaging was carried out for up to 2 hours, with treatments observed for about 3 minutes each. A total of 176 photos were captured, and extracellular staining was quantified by assessing the presence or absence of red-brown staining of extracellular material in each image (Supplementary Table I). Measurements of AGPs, sugar, and weight After 25 hours of stress treatment, eight technical replicates (equals eight petri dishes) were collected into Eppendorf tubes (making one biological bulk replicate) and immediately frozen in liquid nitrogen. The samples were then lyophilized using a VaCo 2 Zirbus lyophilizer at −80°C. The weight of each Eppendorf tube was recorded at three stages: empty, with the frozen sample, and with the lyophilized sample. Two biological bulk replicates were used for analyses. The lyophilized samples were ground with a beat-mill (VWR Star-Beater, VWR International, UK) and pre-extracted two times (2h and 21h respectively) with acetone:water (70:30, v/v) in a ratio of 1:100 (v/v) at 4°C on a shaker. Afterwards, the samples were extracted with double-distilled water for 21h at 4°C on a shaker. In between all extraction steps, the supernatant was separated from the pellet by centrifugation at 19.000 g, 20 min, 4°C (Heraeus Multifuge X3R, Thermo Fisher Scientific Corp., USA). The resulting aqueous extract was lyophilized (Christ Alpha 1-4 LSC, Martin Christ GmbH, Germany) and afterwards dissolved at a concentration of 10 mg/ml under homogenization in the beat-mill (VWR Star-Beater, VWR International, UK, 10 Hz, 1 h with a metal beat). All following experiments were performed with aliquots of this solution. Total sugar and uronic acid contents were determined colorimetrically by the protocols of Dubois et al. 88 and Blumenkrantz and Asboe-Hansen 89 , respectively. Radial gel diffusion assay for interaction of samples with β-D-glucosyl-Yariv (β-D-GlcY) reagent was performed by using the protocol of Castilleux et al. 90 with modifications. In brief, 20 µl of the aqueous extracts were loaded into cavities of an agarose gel containing 0.25 mol/l NaCl and 10 µg/ml β-D-GlcY. Gum arabic solutions in different concentrations (1.0, 0.5, 0.25, 0.125, 0.063 mg/ml in double-distilled water) were used as controls. Hydroxyproline content was determined by colorimetric quantification using the protocol of Stegemann and Stalder 91 . Reactivity of aqueous extracts with the monoclonal antibodies KM1, LM2, LM6 and JIM13 was determined by enzyme-linked immunosorbent assays following the workflow of Mueller et al. 92 . Neutral monosaccharide composition was determined the derivatization of Blakeney et al. 93 with modifications (see Mueller et al. 92 ). To determine mannitol content, the hydrolysis and reduction steps of the derivatization were omitted. All experiments were performed in three technical replicates per each biological bulk replicate. Transcriptomics At 0, 3, 6, 9, and 25 hours, algae samples were collected and immediately frozen in liquid nitrogen, followed by storage at −70°C. For each time point and biological replicate (5 replicates in total), two halves of a culture plate were sampled using a spatula, pooled into an Eppendorf tube, and frozen in liquid nitrogen. The remaining plate halves were reserved for proteomics or metabolomics analyses. RNA was extracted using the Spectrum™ Plant Total RNA Kit following the supplemented protocol A, with modifications made to the cell lysis step: Zygnema samples were lyophilized for 48 hours, then ground using a Retch mill at 20 Hz for 1 minute. One milliliter of lysis buffer was added to the ground, dried samples, which were vortexed for 40 seconds and heated at 56°C for 5 minutes. Mesotaenium samples were placed in 500 µL of lysis buffer, vortexed for 30 seconds, sonicated in a DK SONIC™ ultrasonic cleaner, and heated at 56°C for 5 minutes. Following lysis, DNase treatment with DNase I was performed. For Zygnema, this step was carried out directly on the binding column, while for Mesotaenium , it was performed post-elution of the binding column. RNA sequencing, as well as quality control and raw data trimming, were conducted according to the methods detailed in Dadras et al. 16 . The sample “salt_9hours_BR2” from Zygnema was excluded from further analysis due to a low read count (<0.1 million reads). Trimmed reads were mapped to the transcriptome (see section Transcript-Based Proteome Assembly for details) using Kallisto with bootstrap-samples set to 100. Quantified reads were adjusted in Sleuth to account for technical variability in transcript quantification, and the transcripts per million (TPM) values were subsequently exported. Within the Sleuth visualization environment, the distribution of TPM values was examined. Based on this analysis, the sample “mannitol_3hours_BR3” from Mesotaenium was excluded from further analysis due to an abnormal TPM distribution. Transcripts were filtered based on abundance thresholds to ensure data quality. Transcripts where a single sample had a TPM value below 2 −6 and the average abundance across all samples was less than 2 TPM were removed. Additionally, transcripts with more than five values below 2 −6 and an average abundance across all samples below 4 TPM were excluded from the dataset. Remaining TPM values below 2 −6 were set to 2 −6 . All TPM values were log-transformed prior to further analysis to normalize the data. Proteomics Proteomics analysis was performed on the transcriptomics-paired samples from Mesotaenium (0, 6, and 25 hours) and Zygnema (6 and 25 hours). Protein extraction followed an adapted version of the methods of Niu et al. 94 , with species-specific protocols. For extraction, each sample was homogenized in extraction buffer ( Mesotaenium buffer: 1 mL/sample, 50 mM Tris-HCl, pH 7.5, 10 mM KCl, 0.2 mM PMSF; Zygnema buffer: 1.5 mL/sample, 1% SDS, 0.1 M Tris-HCl, pH 6.8, 2 mM EDTA-Na₂, 20 mM DTT, 2 mM PMSF). Samples were ground in a Ten Broek grinder on ice for 1 minute ( Mesotaenium ) or 4 minutes ( Zygnema ), transferred to clean 1.5 mL eppendorf tubes, and centrifuged ( Mesotaenium : 1,000 × g for 10 seconds; Zygnema : 10,400 × g at 4°C for 7 minutes). The supernatant was then processed further. For precipitation, 100 µL ( Mesotaenium ) or 2 × 200 µL ( Zygnema ) of supernatant was added to 1 mL of pre-cooled ethanol ( Mesotaenium ) or acetone with 5 mM DTT ( Zygnema ), briefly vortexed, and incubated overnight at 4°C ( Mesotaenium ) or −20°C ( Zygnema ). After precipitation, samples were centrifuged ( Mesotaenium : 2,000 × g at 4°C for 15 minutes; Zygnema : 5,000 × g at 4°C for 10 minutes) and the supernatant discarded. The pellets were washed twice with pre-cooled ethanol ( Mesotaenium ) or acetone with 5 mM DTT ( Zygnema ). Mesotaenium samples were dried in a speed-vac at 50°C for 30 minutes, while Zygnema samples were air-dried for 5–8 minutes. The pellets were dissolved in either 30 µL of 5% SDS with 6 M urea ( Mesotaenium ) or 100 µL of the extraction buffer ( Zygnema ). Protein concentration was determined using a 20× dilution of the protein solution in water and the Thermo Scientific™ Pierce™ BCA Protein Assay Kit. Retrospective analysis revealed that the Zygnema protocol, particularly the acetone-based precipitation, also performed slightly better for Mesotaenium compared to ethanol precipitation (data not shown). Protein samples (50 µg) were loaded onto 12% polyacrylamide SDS-PAGE gels (1 mm thickness) and electrophoresed at 100 V using a Bio-Rad™ Mini-PROTEAN® 3 Cell system. The gels were run until the dye front had migrated approximately 1.8 mm into the separation gel. Following electrophoresis, gels were briefly rinsed with distilled water and incubated for 60 minutes in a fixation solution containing 40% (v/v) ethanol and 10% (v/v) acetic acid. For staining, gels were immersed in Coomassie Brilliant Blue stain (0.1% (w/v) Coomassie Brilliant Blue G-250, 5% (w/v) aluminum sulfate-(14–18)-hydrate, 10% (v/v) ethanol, and 2% (v/v) ortho-phosphoric acid) for 30 minutes. Excess stain was removed by washing the gels three times with distilled water for 10 minutes per wash. All staining and washing steps were performed on a rocking platform. Protein-containing gel pieces were excised using a scalpel and subjected to an overnight wash in fixation solution, followed by three additional washes with distilled water (10 minutes each) on a rotating platform. Proteins were subsequently digested and peptides extracted in-gel according to the method described by Shevchenko et al. 95 . Peptides were then desalted and LC-MS/MS analysis was performed as described in the methods of Niemeyer et al. 96 . For initial quality control of the acquired LC-MS/MS data, the raw files were converted to mzML format using MSconvert with the PeakPicking algorithm. Two-dimensional LC-MS chromatograms were plotted in MSnbase to identify abnormalities in the spectra. Based on this assessment, the sample “mannitol_25hours_BR3” from Zygnema and the samples “mannitol_25hours_BR1” and “0hours_BR2” from Mesotaenium were excluded from further analysis. Subsequently, the raw files were analyzed using MaxQuant with the monoisotopic option on a high-performance computing cluster. Quality reports were generated with the R package PTXQC in two consecutive runs. The first MaxQuant analysis aimed to identify batch effects and discover novel protein groups (details provided in the “Transcript-Based Proteome Assembly” section). During this analysis, it was observed that the variation of Mesotaenium samples in the PCA was more strongly influenced by batch effects than by the experimental treatment. As a result, samples from Mesotaenium originating from two SDS-PAGE gels were excluded from the second run. For the second MaxQuant run, the updated FASTA database was employed, and LFQ as well as iBAQ values were calculated. The LFQ values were imported into Perseus for further processing. These values were divided into three distinct matrices: 1) A matrix for the 6-hour timepoint (used for 6-hour differential expression analysis). 2) A matrix for the 25-hour timepoint (used for 25-hour differential expression analysis). 3) A comprehensive matrix containing all samples (used for PCA analysis). Within each matrix, protein groups in which all conditions contained ≥1 missing value were excluded. For the remaining missing values, imputation was performed using random sampling from a normal distribution (downshift = 1.9, width = 0.3). The resulting imputed data were utilized for differential abundance analysis, relative protein-RNA comparison, and functional enrichment analysis. The iBAQ values, in contrast, were not processed in Perseus. Instead, they were directly used for absolute protein-RNA comparisons. Metabolomics The transcriptomics-paired samples from Mesotaenium (3 and 9 hours) were lyophilized using a VaCo 2 Zirbus lyophilizer at −80°C. The dry mass of each sample was adjusted to 9–11 mg per technical replicate. Two types of metabolite extracts were prepared for subsequent LC-MS analysis: methanol-based extracts and lipid extracts. For the methanol-based extracts, 1 mL of 70% (v/v) methanol (pre-cooled to 8°C) was added to the lyophilized, ground samples. The samples were ultrasonicated twice, each for 15 minutes. Following ultrasonication, the samples were centrifuged at 16,000 × g for 20 minutes at 4°C. From the resulting supernatant, 200 µL was transferred to a vial containing an insert. The extracts were dried using a SpeedVac concentrator and reconstituted sequentially: first in 20 µL of 100% methanol, followed by 40 µL of HLJO, with vortexing after each addition. The reconstituted extracts were stored at 4°C until LC-MS analysis. For the lipid extractions, 600 µL of tert-butyl methyl ether (MTBE):methanol (2:1, v/v) was added per 10 mg of lyophilized, ground sample. The mixture was vortexed until phase separation and then sonicated for 15 minutes. After adding 100 µL of HLJO, the samples were vortexed again and centrifuged at 14,000 rpm for 20 minutes at 4°C. From the upper (MTBE) phase, 200 µL was transferred to a vial containing an insert and dried using a nitrogen evaporator. The dried extracts were reconstituted by adding 40 µL of methanol, followed by 10 µL of HLJO, with vortexing after each addition. Extracts were analysed by LC QTOF MS as described in Kasper et al. 97 . The raw Agilent .d data files were converted to .mzML format using MSConvert with the PeakPicking algorithm applied. The resulting .mzML files were processed in MZmine 3 through the following sequential steps: 1) Mass Detection: Initial detection of mass-to-charge (m/z) values, 2) Chromatogram Builder: Construction of chromatograms based on m/z values and retention time. 3) Local Minimum Feature Resolver: Refinement of chromatographic peaks to resolve overlapping features. 4) Join Aligner: Alignment of features across samples to account for retention time shifts. 5)Feature Filtering: Retention of features based on minimal intensity thresholds and presence in a minimum number of samples. 6) Gap Filling: Use of the Peak Finder algorithm to fill gaps in feature matrices where peaks were not initially detected. The identified features were exported and imported into Perseus for statistical analysis. Within Perseus, the data were log-transformed, missing values were imputed from a normal distribution (width = 0.3, downshift = 1.8) or with the constant value of 1000 when combining features from both ion modes, and an ANOVA test with FDR-correction was conducted to identify statistically significant differences between groups. Additionally, autoMSMS files were available were integrated into the MZmine analysis to match fragmentation spectra with their corresponding features. However, their quantification data were excluded from the statistical analysis in Perseus. Transcript-based proteome assembly A transcriptome for both Zygnema and Mesotaenium was assembled using Trinity, incorporating the samples “mannitol_9hours_BR1,” “salt_3hours_BR4,” and “salt_25hours_BR5.” The quality-trimming parameters, adapter sequences, and settings applied during trimming were consistent with those used in Trimmomatic during the regular transcriptomics analysis. Subsequently, Kallisto was employed to quantify the read counts for the newly assembled transcriptome across the aforementioned samples. Transcripts with read counts below 2 TPM for Mesotaenium and below 3 TPM for Zygnema were excluded from the transcriptome. Coding sequences were then identified using TransDecoder with default settings. For Mesotaenium , the --complete_orfs_only flag was applied, ensuring only complete open reading frames were included. This setting was unintentionally omitted for Zygnema , leading to the grouping of incomplete translations with their corresponding complete transcripts for gene ontology (GO) enrichment analysis. The inferred proteome was analyzed using MaxQuant, incorporating .raw files from test runs and samples disqualified due to batch effects. This analysis was performed with the following existing protein databases: Mesotaenium assembly version 1 with annotation version 2 gene models 16 and Mesotaenium chloroplast-encoded proteins 98 for Mesotaenium , and Zygnema circumcarinatum SAG698-1b (the exact strain used for experimentation) gene model 10 and UTEX1559 (a strain relative) plastome 99 for Zygnema . Only protein groups with at least one unique peptide and a Q-value of 0 were retained. Protein groups exclusively originating from the newly transcriptome-based proteome were incorporated into the standard proteome to generate an updated database. This database was subsequently assessed for completeness and quality using BUSCO with the -l viridiplantae_odb10 dataset. Downstream Analysis of Protein and Transcript Data: Differential abundance analysis The log2 fold changes of the processed LFQ values were calculated for mannitol-treated samples versus control samples and NaCl-treated samples versus control samples. These calculations were performed for both the 6-hour and 25-hour time points. The data were visualized using the ggplot package in R, with points color-coded by p-values. P-values were computed per matrix in Perseus using an ANOVA test with parameters S0 = 0.15 and permutation-based FDR = 0.05 to determine significance. Significant protein groups were manually annotated using blastp (default settings, databases nr and SwissProt), domain annotations from SMART, and, where necessary, results from TargetP, dbCAN3, and NetGPI. For the filtered TPM values, log2 fold changes were computed for each combination of NaCl-treated and mannitol-treated time points relative to their respective controls. Statistical analysis was conducted using Python’s scipy.stats and statsmodels packages, where one-way ANOVA tests were performed, followed by FDR correction using a multitest adjustment. These data were visualized similarly to the LFQ results, but with the Seaborn package in Python. Downstream Analysis of Protein and Transcript Data: Comparative Analysis Transcript and corresponding protein levels were compared using both relative and absolute approaches. For the relative comparison, the previously computed log2 fold changes for LFQ and TPM values from the differential abundance analysis were directly plotted against each other without further processing. Additionally, log2 fold changes of the integrated TPM levels for all transcripts were calculated for each treatment-control pair and plotted against the log2 fold changes for LFQ values as well. For the absolute comparison, transcript and protein group levels were compared on a sample-by-sample basis. This was feasible as the RNA sequencing and shotgun proteomics were performed on the same biological samples. Biological replicates were analyzed individually without grouping. Protein levels were quantified using iBAQ values, which were ranked within each sample, with a rank of 1 assigned to the most abundant protein group. The TPM levels were ranked in the same way. The ranks of the protein groups were then matched with the ranks of their corresponding transcripts, and the relationships were visualized using the R package ggplot2. In cases where a protein group corresponded to multiple transcripts, the transcript with the highest abundance was selected for comparison. For the Mesotaenium dataset, each protein-transcript pair was further annotated by coloring points according to contig length from the genome 8 . Downstream Analysis of Protein and Transcript Data: Functional enrichment InterProScan 100 results from previous work 72 were utilized for Gene Ontology (GO) term and InterPro domain annotation. Newly identified transcripts were also annotated using InterProScan. For the functional enrichment analysis at the protein level, the data from the 6-hour and 25-hour time points were merged. Proteins were grouped into nine quadrants based on their log2 fold change (log2FC), using thresholds of 1 and −1. These protein groups were then mapped back to the gene level, and functional enrichment was performed for all quadrants except the central one (where |log2FC| < 1). The analysis was conducted using FuncE with an e-value cutoff of <0.01. To analyze transcript functional enrichment, a time-dependent clustering approach was employed. Log2 fold changes were computed for each treatment-timepoint combination relative to timepoint 0. Transcripts were filtered based on criteria: max(tpm) > 5, max(log2 fold change) > 1, and adjusted p-value < 0.01. The filtered log2 fold changes were clustered using a Gaussian Mixture Model (GMM) implemented in Python’s Scikit-learn library (sklearn.mixture.GaussianMixture). To determine the optimal number of clusters, the Bayes Information criterion (BIC) were calculated for models with 1 to 100 clusters, and the configuration with the lowest BIC score was selected. Transcripts were then grouped into temporal patterns based on the selected number of clusters. Next, these patterns were normalized, and the clustering process was refined with GMM. Starting with 15 clusters, the number was progressively reduced until two seemingly distinct patterns were merged, indicating the final grouping. These final transcript groups, referred to as modules, were mapped back to their corresponding genes. Functional enrichment analysis of these gene modules was performed using FUNC-E with an e-value cutoff of <0.01. An overview was created of enriched terms showing overlap between transcript modules and/or protein quadrants from Zygnema and/or Mesotaenium . Terms associated with only a single gene were excluded, as were terms with 100% gene overlap. In cases of identical overlap, the term with the best p-value was retained; if p-values were equal, one term was selected at random. The resulting list of terms was visualized as a network using the igraph package in R. Phylogenetic analysis XTH The protein data for 16 embryophytes, 11 streptophyte algae, and 8 chlorophytes, as described 101 , were combined into a single database, incorporating the updated Mesotaenium gene model 16 . CAZyme families were identified using dbCAN3 with default settings. Based on results from the HMMer search, proteins identified and annotated as Glycosyl Hydrolase family 16 (GH16) or any of its subfamilies were selected for phylogenetic analysis. The selected protein sequences were aligned using MAFFT 102 with the L-INS-i algorithm. From these aligned sequences, a maximum likelihood phylogenetic tree was constructed using IQ-Tree 103 with 1,000 ultrafast bootstrap replicates 104 . Automatic model selection identified WAG+R9 as the best-fit model. The resulting phylogeny was visualized using iTOL 105 . Downstream Analysis of Transcript Data: Gene Regulatory Network (GRN) From the time-series transcriptomics data, a Gene Network Reconstruction (GNR) was constructed for each species and treatment individually using Sliding Window Inference for Network Generation (SWING). Log-transformed TPM values were utilized, and the data was filtered for transcripts with significant temporal variation (ANOVA, unadjusted p-value < 0.0001) across all timepoints (0, 3, 6, 9, and 25 hours). Edges were computed in SWING using the following parameters: k_min = 1, k_max = 3, window size (w) = 3, method = “RandomForest,” and trees = 500. To construct the GNRs, the top 0.1% of scoring edges were retained, consistent with the cut-off previously applied 72 . For the osmotic stress GNR of each species, common edges between the salt and mannitol time-series networks were retained, except for edges that were also present in the control time-series network. Unfortunately, no edges remained for constructing a network for either Zygnema or Mesotaenium . Long term salt stress For one Zygnema and one Mesotaenium plate, cellophane discs containing algae were transferred onto fresh plates with the same cultivation medium, supplemented with 0.15 M NaCl. These plates were then placed in a long-term storage cabinet under constant light conditions (light intensity ≈ 10 µE) and a temperature of 14°C. After 26 days, we examined the algal morphology using a Zeiss Axioscope 7 light microscope. In Zygnema , we observed deterioration of the cell wall, while in Mesotaenium , we noted the formation of short filaments (maximum observed length = 10 cells/filament). This prompted us to set up a growth experiment: The strains SAG 12.97 Mesotaenium endlicherianum and SAG 698-1b Zygnema circumcarinatum were cultivated on solid Warris-H medium under standard culture conditions (20°C, 30 µE, 14:10 L/D cycle). In the middle of log-phase after 14 days of growth, the algae were resuspended in 2 ml of liquid Warris-H medium containing 0.15 M NaCl. 300 µl of this suspension were placed on Petri dishes (Ø 9 cm) with solid Warris-H medium containing 0.15 M NaCl with and without sterile cellophane disks. The suspension was carefully spreaded out on the top of the agar surface with sterile Delta™ Disposable Cell Spreader (Merck, Darmstadt, Germany). The algae were grown for eight weeks. The microscopical control was performed every seven days. The Olympus BX-60 microscope (Olympus, Tokyo, Japan) equipped with camera Gryphax Prokyon and the software Gryphax ver. 2.2.0.1234 (both from Jenoptik, Jena, Germany) was used for taking microphotographs. The experiments were performed in three technical replicates. Zygnema bleached out after 48 hours of exposure to the salt stress and was therefore excluded from further experiments. In contrast, Mesotaenium showed satisfied growth on medium containing salt. The number of short filaments were measured for all taken images; chains of loosely connected cells were counted as individual cells, while a continuous cell wall spanning multiple cells was required for a structure to be considered as filament. Data availability All RNAseq data have been deposited under Bioproject ID PRJNA1217932 Micrographs used for quantifications have been uploaded on 10.5281/zenodo.14777872 The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE 106 partner repository with the dataset identifier PXD060327 Extended Data Figure 1. Relative protein abundance at 25 hours is compared to relative RNA abundance (integrated from 0 to 25 hours) for proteins exhibiting a significant fold change (FDR < 0.05) under either salt or mannitol treatments relative to the control. Protein abbreviations are labeled when the summed absolute values of log2 fold changes (LFQ and TPM) exceed 2.5. Full protein names are provided in Supplementary Data I. Software list View this table: View inline View popup Acknowledgements J.d.V. thanks (a) the European Research Council for funding under the European Union’s Horizon 2020 research and innovation program (Grant Agreement No. 852725; ERC-StG “TerreStriAL”) and (b) the German Research Foundation (DFG) for support through project “SHOAL” (514060973), project “Klebsome” (509535047) and the framework of the Priority Program “MAdLand – Molecular Adaptation to Land: Plant Evolution to Change” (SPP 2237; 440231723 and 528076711), in which C.F.K. and J.M.S.Z. partake as associate members. C.F.K. and J.M.S.Z. are grateful for support through the IMPRS Genome Science. B.C. and L.P. thank the German Research Foundation (DFG) for support through the framework of the Priority Program “MAdLand – Molecular Adaptation to Land: Plant Evolution to Change” (SPP 2237; 440046237). K.S. is financed by DFG INST 186/1465-1. The HF LC-MS system ultimate 3000-Q Exactive was granted by DFG INST 186/1230-1 FUGG. K.S. is financed by DFG INST 186/1465-1. We extend our heartfelt gratitude to Dr. Ilka Nacif Abreu, Prof. Dr. Ivo Feußner, and Dr. Kirstin Feußner for their expert guidance in metabolomics analysis and their invaluable contribution to performing the metabolomic measurements. We also wish to thank René Heisse for providing exceptional technical assistance during the stress experiments. REFERENCES 1. ↵ Wickett , N. J. et al. Phylotranscriptomic analysis of the origin and early diversification of land plants . Proc Natl Acad Sci U S A 111 , E4859 – 4868 ( 2014 ). doi: 10.1073/pnas.1323926111 OpenUrl Abstract / FREE Full Text 2. One Thousand Plant Transcriptomes, I . One thousand plant transcriptomes and the phylogenomics of green plants . Nature 574 , 679 - 685 ( 2019 ). doi: 10.1038/s41586-019-1693-2 OpenUrl CrossRef PubMed 3. Puttick , M. N. et al. The Interrelationships of Land Plants and the Nature of the Ancestral Embryophyte . Curr Biol 28 , 733 – 745 e732 ( 2018 ). doi: 10.1016/j.cub.2018.01.063 OpenUrl CrossRef PubMed 4. Bierenbroodspot , M. et al. Phylogeny and evolution of streptophyte algae . Ann Bot ( 2024 ). doi: 10.1093/aob/mcae091 OpenUrl CrossRef PubMed 5. ↵ Harris , B. J. et al. Divergent evolutionary trajectories of bryophytes and tracheophytes from a complex common ancestor of land plants . Nat Ecol Evol 6 , 1634 – 1643 ( 2022 ). doi: 10.1038/s41559-022-01885-x OpenUrl CrossRef 6. ↵ Gontcharov , A. A. , Marin , B. A. & Melkonian , M. A . Molecular phylogeny of conjugating green algae (Zygnemophyceae, Streptophyta) inferred from SSU rDNA sequence comparisons . J Mol Evol 56 , 89 – 104 ( 2003 ). doi: 10.1007/s00239-002-2383-4 OpenUrl CrossRef PubMed Web of Science 7. ↵ Hess , S. et al. A phylogenomically informed five-order system for the closest relatives of land plants . Current Biology 32 , 4473 – 4482 . e4477 ( 2022 ). doi: 10.1016/j.cub.2022.08.022 OpenUrl CrossRef PubMed 8. ↵ Cheng , S. et al. Genomes of Subaerial Zygnematophyceae Provide Insights into Land Plant Evolution . Cell 179 , 1057 – 1067 e1014 ( 2019 ). doi: 10.1016/j.cell.2019.10.019 OpenUrl CrossRef PubMed 9. Jiao , C. et al. The Penium margaritaceum Genome: Hallmarks of the Origins of Land Plants . Cell 181 , 1097 – 1111 e1012 ( 2020 ). doi: 10.1016/j.cell.2020.04.019 OpenUrl CrossRef PubMed 10. ↵ Feng , X. et al. Genomes of multicellular algal sisters to land plants illuminate signaling network evolution . Nature Genetics 56 , 1018 – 1031 ( 2024 ). doi: 10.1038/s41588-024-01737-3 OpenUrl CrossRef PubMed 11. ↵ Sekimoto , H. et al. A divergent RWP-RK transcription factor determines mating type in heterothallic Closterium . New Phytologist 237 , 1636 – 1651 ( 2023 ). OpenUrl CrossRef PubMed 12. ↵ Ju , C. et al. Conservation of ethylene as a plant hormone over 450 million years of evolution . Nature plants 1 , 14004 ( 2015 ). 13. ↵ de Vries , J. & Archibald , J. M . Plant evolution: landmarks on the path to terrestrial life . New Phytol 217 , 1428 – 1434 ( 2018 ). doi: 10.1111/nph.14975 OpenUrl CrossRef PubMed 14. ↵ Sun , Y. et al. A ligand-independent origin of abscisic acid perception . Proceedings of the National Academy of Sciences 116 , 24892 – 24899 ( 2019 ). OpenUrl Abstract / FREE Full Text 15. Kuhn , A. et al. RAF-like protein kinases mediate a deeply conserved, rapid auxin response . Cell 187 , 130 – 148 e117 ( 2024 ). doi: 10.1016/j.cell.2023.11.021 OpenUrl CrossRef PubMed 16. ↵ Dadras , A. et al. Environmental gradients reveal stress hubs pre-dating plant terrestrialization . Nat Plants 9 , 1419 – 1438 ( 2023 ). doi: 10.1038/s41477-023-01491-0 OpenUrl CrossRef 17. ↵ Zegers , J. M. S. , Irisarri , I. , de Vries , S. & de Vries , J . Evolving circuitries in plant signaling cascades . J Cell Sci 137 ( 2024 ). doi: 10.1242/jcs.261712 OpenUrl CrossRef 18. ↵ Yu , B. , Chao , D. Y. & Zhao , Y . How plants sense and respond to osmotic stress . J Integr Plant Biol 66 , 394 – 423 ( 2024 ). doi: 10.1111/jipb.13622 OpenUrl CrossRef PubMed 19. ↵ Abid , M. et al. Physiological and biochemical changes during drought and recovery periods at tillering and jointing stages in wheat (Triticum aestivum L .). Sci Rep 8 , 4615 ( 2018 ). doi: 10.1038/s41598-018-21441-7 OpenUrl CrossRef PubMed 20. ↵ McAdam , S. A. & Brodribb , T. J . Stomatal innovation and the rise of seed plants . Ecol Lett 15 , 1 – 8 ( 2012 ). doi: 10.1111/j.1461-0248.2011.01700.x OpenUrl CrossRef PubMed 21. ↵ Cruz de Carvalho , M. H. Drought stress and reactive oxygen species: Production, scavenging and signaling . Plant Signal Behav 3 , 156 – 165 ( 2008 ). doi: 10.4161/psb.3.3.5536 OpenUrl CrossRef PubMed 22. ↵ Nagao , M. , Matsui , K. & Uemura , M . Klebsormidium flaccidum, a charophycean green alga, exhibits cold acclimation that is closely associated with compatible solute accumulation and ultrastructural changes . Plant Cell Environ 31 , 872 – 885 ( 2008 ). doi: 10.1111/j.1365-3040.2008.01804.x OpenUrl CrossRef PubMed Web of Science 23. ↵ Holzinger , A. & Karsten , U . Desiccation stress and tolerance in green algae: consequences for ultrastructure, physiological and molecular mechanisms . Front Plant Sci 4 , 327 ( 2013 ). doi: 10.3389/fpls.2013.00327 OpenUrl CrossRef 24. ↵ Suzuki , N. et al. Temporal-spatial interaction between reactive oxygen species and abscisic acid regulates rapid systemic acclimation in plants . Plant Cell 25 , 3553 – 3569 ( 2013 ). doi: 10.1105/tpc.113.114595 OpenUrl Abstract / FREE Full Text 25. Mittler , R. & Blumwald , E . The roles of ROS and ABA in systemic acquired acclimation . Plant Cell 27 , 64 – 70 ( 2015 ). doi: 10.1105/tpc.114.133090 OpenUrl Abstract / FREE Full Text 26. ↵ Mittler , R. ROS Are Good . Trends Plant Sci 22 , 11 – 19 ( 2017 ). doi: 10.1016/j.tplants.2016.08.002 OpenUrl CrossRef PubMed 27. ↵ Xiong , L. & Zhu , J. K . Molecular and genetic aspects of plant responses to osmotic stress . Plant Cell Environ 25 , 131 – 139 ( 2002 ). doi: 10.1046/j.1365-3040.2002.00782.x OpenUrl CrossRef PubMed Web of Science 28. ↵ Vilarrasa-Blasi , J. et al. Multi-omics analysis of green lineage osmotic stress pathways unveils crucial roles of different cellular compartments . Nat Commun 15 , 5988 ( 2024 ). doi: 10.1038/s41467-024-49844-3 OpenUrl CrossRef 29. ↵ Rippin , M. , Becker , B. & Holzinger , A . Enhanced Desiccation Tolerance in Mature Cultures of the Streptophytic Green Alga Zygnema circumcarinatum Revealed by Transcriptomics . Plant Cell Physiol 58 , 2067 – 2084 ( 2017 ). doi: 10.1093/pcp/pcx136 OpenUrl CrossRef PubMed 30. ↵ Rippin , M. et al. Metatranscriptomic and metabolite profiling reveals vertical heterogeneity within a Zygnema green algal mat from Svalbard (High Arctic) . Environmental Microbiology 21 , 4283 – 4299 ( 2019 ). OpenUrl CrossRef 31. ↵ Herburger , K. , Lewis , L. A. & Holzinger , A . Photosynthetic efficiency, desiccation tolerance and ultrastructure in two phylogenetically distinct strains of alpine Zygnema sp. (Zygnematophyceae, Streptophyta): role of pre-akinete formation . Protoplasma 252 , 571 – 589 ( 2015 ). doi: 10.1007/s00709-014-0703-3 OpenUrl CrossRef PubMed 32. ↵ Fitzek , E. et al. Cell wall enzymes in Zygnema circumcarinatum UTEX 1559 respond to osmotic stress in a plant-like fashion . Frontiers in plant science 10 , 732 ( 2019 ). 33. ↵ Domozych , D. S. , Kozel , L. & Palacio-Lopez , K . The effects of osmotic stress on the cell wall-plasma membrane domains of the unicellular streptophyte, Penium margaritaceum . Protoplasma 258 , 1231 – 1249 ( 2021 ). doi: 10.1007/s00709-021-01644-y OpenUrl CrossRef PubMed 34. ↵ de Vries , J. , Curtis , B. A. , Gould , S. B. & Archibald , J. M . Embryophyte stress signaling evolved in the algal progenitors of land plants . Proceedings of the National Academy of Sciences 115 , E3471 – E3480 ( 2018 ). OpenUrl Abstract / FREE Full Text 35. ↵ Fabregas , N. , Yoshida , T. & Fernie , A. R . Role of Raf-like kinases in SnRK2 activation and osmotic stress response in plants . Nat Commun 11 , 6184 ( 2020 ). doi: 10.1038/s41467-020-19977-2 OpenUrl CrossRef PubMed 36. ↵ Chen , K. et al. Abscisic acid dynamics, signaling, and functions in plants . J Integr Plant Biol 62 , 25 – 54 ( 2020 ). doi: 10.1111/jipb.12899 OpenUrl CrossRef PubMed 37. ↵ Kilian , J. et al. The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses . Plant J 50 , 347 – 363 ( 2007 ). doi: 10.1111/j.1365-313X.2007.03052.x OpenUrl CrossRef PubMed Web of Science 38. ↵ Baker , N. R . Chlorophyll fluorescence: a probe of photosynthesis in vivo . Annu Rev Plant Biol 59 , 89 – 113 ( 2008 ). doi: 10.1146/annurev.arplant.59.032607.092759 OpenUrl CrossRef PubMed Web of Science 39. ↵ Schlösser , U. G. & Renkert , U. (IWF (Göttingen) , 1992 ). 40. ↵ Liu , X. , Rodermel , S. R. & Yu , F . A var2 leaf variegation suppressor locus, SUPPRESSOR OF VARIEGATION3, encodes a putative chloroplast translation elongation factor that is important for chloroplast development in the cold . BMC Plant Biol 10 , 287 ( 2010 ). doi: 10.1186/1471-2229-10-287 OpenUrl CrossRef PubMed 41. ↵ Sengupta , D. , Naik , D. & Reddy , A. R . Plant aldo-keto reductases (AKRs) as multi-tasking soldiers involved in diverse plant metabolic processes and stress defense: A structure-function update . J Plant Physiol 179 , 40 – 55 ( 2015 ). doi: 10.1016/j.jplph.2015.03.004 OpenUrl CrossRef PubMed 42. ↵ Hand , S. C. , Menze , M. A. , Toner , M. , Boswell , L. & Moore , D . LEA proteins during water stress: not just for plants anymore . Annu Rev Physiol 73 , 115 – 134 ( 2011 ). doi: 10.1146/annurev-physiol-012110-142203 OpenUrl CrossRef PubMed 43. ↵ Hutin , C. et al. Early light-induced proteins protect Arabidopsis from photooxidative stress . Proceedings of the National Academy of Sciences 100 , 4921 – 4926 ( 2003 ). OpenUrl Abstract / FREE Full Text 44. ↵ Li , J. , Guo , Y. & Yang , Y . The molecular mechanism of plasma membrane H(+)-ATPases in plant responses to abiotic stress . J Genet Genomics 49 , 715 – 725 ( 2022 ). doi: 10.1016/j.jgg.2022.05.007 OpenUrl CrossRef PubMed 45. ↵ Stein , O. & Granot , D . An Overview of Sucrose Synthases in Plants . Front Plant Sci 10 , 95 ( 2019 ). doi: 10.3389/fpls.2019.00095 OpenUrl CrossRef PubMed 46. ↵ Sudhakaran , S. et al. Significance of solute specificity, expression, and gating mechanism of tonoplast intrinsic protein during development and stress response in plants . Physiol Plant 172 , 258 – 274 ( 2021 ). doi: 10.1111/ppl.13386 OpenUrl CrossRef 47. ↵ Benitez-Fuente , F. & Botella , M. A . Biological roles of plant synaptotagmins . Eur J Cell Biol 102 , 151335 ( 2023 ). doi: 10.1016/j.ejcb.2023.151335 OpenUrl CrossRef PubMed 48. ↵ Guihur , A. , Rebeaud , M. E. & Goloubinoff , P . How do plants feel the heat and survive? Trends Biochem Sci 47 , 824 – 838 ( 2022 ). doi: 10.1016/j.tibs.2022.05.004 OpenUrl CrossRef 49. ↵ Hura , T. et al. Rieske iron-sulfur protein of cytochrome-b6f is involved in plant recovery after drought stress . Environmental and Experimental Botany 156 , 228 – 239 ( 2018 ). doi: 10.1016/j.envexpbot.2018.09.003 OpenUrl CrossRef 50. ↵ Gigli-Bisceglia , N. , Engelsdorf , T. & Hamann , T . Plant cell wall integrity maintenance in model plants and crop species-relevant cell wall components and underlying guiding principles . Cell Mol Life Sci 77 , 2049 – 2077 ( 2020 ). doi: 10.1007/s00018-019-03388-8 OpenUrl CrossRef PubMed 51. ↵ Colin , L. et al. The cell biology of primary cell walls during salt stress . Plant Cell 35 , 201 – 217 ( 2023 ). doi: 10.1093/plcell/koac292 OpenUrl CrossRef 52. ↵ Carnevale Neto , F. , Pascua , V. & Raftery , D. Formation of sodium cluster ions complicates liquid chromatography-mass spectrometry metabolomics analyses . Rapid Commun Mass Spectrom 35 , e9175 ( 2021 ). doi: 10.1002/rcm.9175 OpenUrl CrossRef PubMed 53. ↵ Delauney , A. J. & Verma , D. P. S . Proline biosynthesis and osmoregulation in plants . The Plant Journal 4 , 215 – 223 ( 1993 ). doi: 10.1046/j.1365-313X.1993.04020215.x OpenUrl CrossRef Web of Science 54. ↵ Tenhaken , R . Cell wall remodeling under abiotic stress . Front Plant Sci 5 , 771 ( 2014 ). doi: 10.3389/fpls.2014.00771 OpenUrl CrossRef PubMed 55. ↵ Choi , J. Y. , Seo , Y. S. , Kim , S. J. , Kim , W. T. & Shin , J. S . Constitutive expression of CaXTH3, a hot pepper xyloglucan endotransglucosylase/hydrolase, enhanced tolerance to salt and drought stresses without phenotypic defects in tomato plants (Solanum lycopersicum cv . Dotaerang). Plant Cell Rep 30 , 867 – 877 ( 2011 ). doi: 10.1007/s00299-010-0989-3 OpenUrl CrossRef PubMed Web of Science 56. ↵ Ma , Y. et al. Distinct functions of FASCILIN-LIKE ARABINOGALACTAN PROTEINS relate to domain structure . Plant Physiol 192 , 119 – 132 ( 2023 ). doi: 10.1093/plphys/kiad097 OpenUrl CrossRef PubMed 57. ↵ Johnson , K. L. , Jones , B. J. , Bacic , A. & Schultz , C. J . The fasciclin-like arabinogalactan proteins of Arabidopsis. A multigene family of putative cell adhesion molecules . Plant Physiol 133 , 1911 – 1925 ( 2003 ). doi: 10.1104/pp.103.031237 OpenUrl Abstract / FREE Full Text 58. ↵ Seifert , G. J . Fascinating Fasciclins: A Surprisingly Widespread Family of Proteins that Mediate Interactions between the Cell Exterior and the Cell Surface . Int J Mol Sci 19 ( 2018 ). doi: 10.3390/ijms19061628 OpenUrl CrossRef PubMed 59. ↵ Shafee , T. , Bacic , A. & Johnson , K . Evolution of Sequence-Diverse Disordered Regions in a Protein Family: Order within the Chaos . Mol Biol Evol 37 , 2155 – 2172 ( 2020 ). doi: 10.1093/molbev/msaa096 OpenUrl CrossRef PubMed 60. ↵ Ma , Y. & Johnson , K . Arabinogalactan proteins - Multifunctional glycoproteins of the plant cell wall . Cell Surf 9 , 100102 ( 2023 ). doi: 10.1016/j.tcsw.2023.100102 OpenUrl CrossRef PubMed 61. ↵ Classen , B. , Csavas , M. , Borbas , A. , Dingermann , T. & Zundorf , I . Monoclonal antibodies against an arabinogalactan-protein from pressed juice of Echinacea purpurea . Planta Med 70 , 861 – 865 ( 2004 ). doi: 10.1055/s-2004-827236 OpenUrl CrossRef PubMed 62. ↵ Yates , E. A. et al. Characterization of carbohydrate structural features recognized by anti-arabinogalactan-protein monoclonal antibodies . Glycobiology 6 , 131 – 139 ( 1996 ). doi: 10.1093/glycob/6.2.131 OpenUrl CrossRef PubMed Web of Science 63. ↵ Pfeifer , L. et al. Search for evolutionary roots of land plant arabinogalactan-proteins in charophytes: presence of a rhamnogalactan-protein in Spirogyra pratensis (Zygnematophyceae) . Plant J 109 , 568 – 584 ( 2022 ). doi: 10.1111/tpj.15577 OpenUrl CrossRef PubMed 64. ↵ Knoch , E. , Dilokpimol , A. & Geshi , N . Arabinogalactan proteins: focus on carbohydrate active enzymes . Front Plant Sci 5 , 198 ( 2014 ). doi: 10.3389/fpls.2014.00198 OpenUrl CrossRef PubMed 65. ↵ Lamport , D. T. , Kieliszewski , M. J. & Showalter , A. M . Salt stress upregulates periplasmic arabinogalactan proteins: using salt stress to analyse AGP function . New Phytol 169 , 479 – 492 ( 2006 ). doi: 10.1111/j.1469-8137.2005.01591.x OpenUrl CrossRef PubMed Web of Science 66. ↵ Olmos , E. , Garcia De La Garma , J. , Gomez-Jimenez , M. C. & Fernandez-Garcia , N. Arabinogalactan Proteins Are Involved in Salt-Adaptation and Vesicle Trafficking in Tobacco by-2 Cell Cultures . Front Plant Sci 8 , 1092 ( 2017 ). doi: 10.3389/fpls.2017.01092 OpenUrl CrossRef PubMed 67. ↵ Domozych , D. S. et al. Pectin metabolism and assembly in the cell wall of the charophyte green alga Penium margaritaceum . Plant Physiol 165 , 105 – 118 ( 2014 ). doi: 10.1104/pp.114.236257 OpenUrl Abstract / FREE Full Text 68. ↵ Wu , H. C. , Yu , S. Y. , Wang , Y. D. & Jinn , T. L . Guard Cell-Specific Pectin METHYLESTERASE53 Is Required for Abscisic Acid-Mediated Stomatal Function and Heat Response in Arabidopsis . Front Plant Sci 13 , 836151 ( 2022 ). doi: 10.3389/fpls.2022.836151 OpenUrl CrossRef 69. ↵ Lamport , D. T. , Varnai , P. & Seal , C. E . Back to the future with the AGP-Ca2+ flux capacitor . Ann Bot 114 , 1069 – 1085 ( 2014 ). doi: 10.1093/aob/mcu161 OpenUrl CrossRef PubMed 70. ↵ Shin , Y. , Chane , A. , Jung , M. & Lee , Y . Recent Advances in Understanding the Roles of Pectin as an Active Participant in Plant Signaling Networks . Plants (Basel ) 10 ( 2021 ). doi: 10.3390/plants10081712 OpenUrl CrossRef PubMed 71. ↵ Rieseberg , T. P. et al. Divergent responses in desiccation experiments in two ecophysiologically different Zygnematophyceae . Physiol Plant 175 , e14056 ( 2023 ). doi: 10.1111/ppl.14056 OpenUrl CrossRef 72. ↵ Rieseberg , T. P. et al. Time-resolved oxidative signal convergence across the algae– embryophyte divide . bioRxiv ( 2024 ). doi: 10.1101/2024.03.11.584470 OpenUrl Abstract / FREE Full Text 73. ↵ Braccini , I. & Perez , S . Molecular basis of C(2+)-induced gelation in alginates and pectins: the egg-box model revisited . Biomacromolecules 2 , 1089 – 1096 ( 2001 ). doi: 10.1021/bm010008g OpenUrl CrossRef PubMed Web of Science 74. ↵ Zhou , K . The regulation of the cell wall by glycosylphosphatidylinositol-anchored proteins in Arabidopsis . Front Cell Dev Biol 10 , 904714 ( 2022 ). doi: 10.3389/fcell.2022.904714 OpenUrl CrossRef 75. ↵ Qi , Z. et al. The dehydration-responsive protein PpFAS1.3 in moss Physcomitrium patens plays a regulatory role in lipid metabolism . J Plant Physiol 297 , 154253 ( 2024 ). doi: 10.1016/j.jplph.2024.154253 OpenUrl CrossRef 76. ↵ Basu , D. et al. Glycosylation of a Fasciclin-Like Arabinogalactan-Protein (SOS5) Mediates Root Growth and Seed Mucilage Adherence via a Cell Wall Receptor-Like Kinase (FEI1/FEI2) Pathway in Arabidopsis . PLoS One 11 , e0145092 ( 2016 ). doi: 10.1371/journal.pone.0145092 OpenUrl CrossRef PubMed 77. ↵ Reisen , D. , Marty , F. & Leborgne-Castel , N . New insights into the tonoplast architecture of plant vacuoles and vacuolar dynamics during osmotic stress . BMC Plant Biol 5 , 13 ( 2005 ). doi: 10.1186/1471-2229-5-13 OpenUrl CrossRef PubMed 78. ↵ Kumar , D. et al. The Xerophyta viscosa aldose reductase (ALDRXV4) confers enhanced drought and salinity tolerance to transgenic tobacco plants by scavenging methylglyoxal and reducing the membrane damage . Mol Biotechnol 54 , 292 – 303 ( 2013 ). doi: 10.1007/s12033-012-9567-y OpenUrl CrossRef PubMed 79. ↵ Kanayama , Y. et al. Characterization of an uncharacterized aldo-keto reductase gene from peach and its role in abiotic stress tolerance . Phytochemistry 104 , 30 – 36 ( 2014 ). doi: 10.1016/j.phytochem.2014.04.008 OpenUrl CrossRef PubMed 80. ↵ Holzinger , A. et al. Transcriptomics of desiccation tolerance in the streptophyte green alga Klebsormidium reveal a land plant-like defense reaction . PLoS One 9 , e110630 ( 2014 ). doi: 10.1371/journal.pone.0110630 OpenUrl CrossRef PubMed 81. ↵ Hess , D. , Heise , C. M. , Schubert , H. , Hess , W. R. & Hagemann , M . The impact of salt stress on the physiology and the transcriptome of the model streptophyte green alga Chara braunii . Physiol Plant 175 , e14123 ( 2023 ). doi: 10.1111/ppl.14123 OpenUrl CrossRef 82. ↵ Lee , E. et al. Ionic stress enhances ER-PM connectivity via phosphoinositide-associated SYT1 contact site expansion in Arabidopsis . Proc Natl Acad Sci U S A 116 , 1420 – 1429 ( 2019 ). doi: 10.1073/pnas.1818099116 OpenUrl Abstract / FREE Full Text 83. ↵ Schroder , F. , Lisso , J. , Lange , P. & Mussig , C . The extracellular EXO protein mediates cell expansion in Arabidopsis leaves . BMC Plant Biol 9 , 20 ( 2009 ). doi: 10.1186/1471-2229-9-20 OpenUrl CrossRef PubMed 84. ↵ Schroder , F. , Lisso , J. & Mussig , C . EXORDIUM-LIKE1 promotes growth during low carbon availability in Arabidopsis . Plant Physiol 156 , 1620 – 1630 ( 2011 ). doi: 10.1104/pp.111.177204 OpenUrl Abstract / FREE Full Text 85. ↵ Yuan , F. et al. OSCA1 mediates osmotic-stress-evoked Ca2+ increases vital for osmosensing in Arabidopsis . Nature 514 , 367 – 371 ( 2014 ). doi: 10.1038/nature13593 OpenUrl CrossRef PubMed Web of Science 86. ↵ Cosgrove , D. J . Plant expansins: diversity and interactions with plant cell walls . Curr Opin Plant Biol 25 , 162 – 172 ( 2015 ). doi: 10.1016/j.pbi.2015.05.014 OpenUrl CrossRef PubMed 87. ↵ Ichimura , T. 88. ↵ DuBois , M. , Gilles , K. A. , Hamilton , J. K. , Rebers , P. A. & Smith , F . Colorimetric Method for Determination of Sugars and Related Substances . Analytical Chemistry 28 , 350 – 356 ( 1956 ). doi: 10.1021/ac60111a017 OpenUrl CrossRef Web of Science 89. ↵ Blumenkrantz , N. & Asboe-Hansen , G . New method for quantitative determination of uronic acids . Anal Biochem 54 , 484 – 489 ( 1973 ). doi: 10.1016/0003-2697(73)90377-1 OpenUrl CrossRef PubMed Web of Science 90. ↵ Zoë A. Castilleux , R. et al. in The Plant Cell Wall: Methods and Protocols (ed Zoë A. Popper ) 383 – 402 ( Springer New York , 2020 ). 91. ↵ Stegemann , H. & Stalder , K. Determination of hydroxyproline . Clin Chim Acta 18 , 267 – 273 ( 1967 ). doi: 10.1016/0009-8981(67)90167-2 OpenUrl CrossRef PubMed Web of Science 92. ↵ Mueller , K. K. et al. Fern cell walls and the evolution of arabinogalactan proteins in streptophytes . Plant J 114 , 875 – 894 ( 2023 ). doi: 10.1111/tpj.16178 OpenUrl CrossRef 93. ↵ Blakeney , A. B. , Harris , P. J. , Henry , R. J. & Stone , B. A . A simple and rapid preparation of alditol acetates for monosaccharide analysis . Carbohydrate Research 113 , 291 – 299 ( 1983 ). doi: 10.1016/0008-6215(83)88244-5 OpenUrl CrossRef Web of Science 94. ↵ Niu , L. et al. Modified TCA/acetone precipitation of plant proteins for proteomic analysis . PLoS One 13 , e0202238 ( 2018 ). doi: 10.1371/journal.pone.0202238 OpenUrl CrossRef PubMed 95. ↵ Shevchenko , A. et al. Linking genome and proteome by mass spectrometry: large-scale identification of yeast proteins from two dimensional gels . Proc Natl Acad Sci U S A 93 , 14440 – 14445 ( 1996 ). doi: 10.1073/pnas.93.25.14440 OpenUrl Abstract / FREE Full Text 96. ↵ Niemeyer , P. W. et al. A seed-like proteome in oil-rich tubers . Plant J 112 , 518 – 534 ( 2022 ). doi: 10.1111/tpj.15964 OpenUrl CrossRef PubMed 97. ↵ Kasper , K. et al. Multi-omics analysis of xylem sap uncovers dynamic modulation of poplar defenses by ammonium and nitrate . Plant J 111 , 282 – 303 ( 2022 ). doi: 10.1111/tpj.15802 OpenUrl CrossRef PubMed 98. ↵ Civan , P. , Foster , P. G. , Embley , M. T. , Seneca , A. & Cox , C. J . Analyses of charophyte chloroplast genomes help characterize the ancestral chloroplast genome of land plants . Genome Biol Evol 6 , 897 – 911 ( 2014 ). doi: 10.1093/gbe/evu061 OpenUrl CrossRef PubMed 99. ↵ Orton , L. M. et al. Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution . Journal of experimental botany 71 , 3361 – 3373 ( 2020 ). OpenUrl CrossRef PubMed 100. ↵ Jones , P. et al. InterProScan 5: genome-scale protein function classification . Bioinformatics 30 , 1236 – 1240 ( 2014 ). doi: 10.1093/bioinformatics/btu031 OpenUrl CrossRef PubMed Web of Science 101. ↵ de Vries , S. et al. The evolution of the phenylpropanoid pathway entailed pronounced radiations and divergences of enzyme families . Plant J 107 , 975 – 1002 ( 2021 ). doi: 10.1111/tpj.15387 OpenUrl CrossRef PubMed 102. ↵ Katoh , K. & Standley , D. M . MAFFT multiple sequence alignment software version 7: improvements in performance and usability . Molecular Biology and Evolution 30 , 772 – 780 ( 2013 ). OpenUrl CrossRef PubMed Web of Science 103. ↵ Minh , B. Q. et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era . Molecular biology and evolution 37 , 1530 – 1534 ( 2020 ). OpenUrl CrossRef PubMed 104. ↵ Hoang , D. T. , Chernomor , O. , Von Haeseler , A. , Minh , B. Q. & Vinh , L. S . UFBoot2: improving the ultrafast bootstrap approximation . Molecular biology and evolution 35 , 518 – 522 ( 2018 ). OpenUrl CrossRef PubMed 105. ↵ Letunic , I. & Bork , P . Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation . Nucleic acids research 49 , W293 – W296 ( 2021 ). OpenUrl CrossRef PubMed 106. ↵ Perez-Riverol , Y. et al. The PRIDE database at 20 years: 2025 update . Nucleic Acids Res 53 , D543 – D553 ( 2025 ). doi: 10.1093/nar/gkae1011 OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted February 07, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Systems acclimation to osmotic stress in zygnematophyte cells Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Systems acclimation to osmotic stress in zygnematophyte cells Jaccoline M.S. Zegers , Lukas Pfeifer , Tatyana Darienko , Kerstin Schmitt , Craig A. Dorsch , Cäcilia F. Kunz , Gerhard H. Braus , Birgit Classen , Oliver Valerius , Jan de Vries bioRxiv 2025.02.06.636795; doi: https://doi.org/10.1101/2025.02.06.636795 Share This Article: Copy Citation Tools Systems acclimation to osmotic stress in zygnematophyte cells Jaccoline M.S. Zegers , Lukas Pfeifer , Tatyana Darienko , Kerstin Schmitt , Craig A. Dorsch , Cäcilia F. Kunz , Gerhard H. Braus , Birgit Classen , Oliver Valerius , Jan de Vries bioRxiv 2025.02.06.636795; doi: https://doi.org/10.1101/2025.02.06.636795 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Evolutionary Biology Subject Areas All Articles Animal Behavior and Cognition (7622) Biochemistry (17650) Bioengineering (13871) Bioinformatics (41881) Biophysics (21424) Cancer Biology (18566) Cell Biology (25461) Clinical Trials (138) Developmental Biology (13365) Ecology (19866) Epidemiology (2067) Evolutionary Biology (24290) Genetics (15590) Genomics (22476) Immunology (17713) Microbiology (40331) Molecular Biology (17148) Neuroscience (88473) Paleontology (666) Pathology (2827) Pharmacology and Toxicology (4816) Physiology (7635) Plant Biology (15114) Scientific Communication and Education (2044) Synthetic Biology (4286) Systems Biology (9815) Zoology (2268)