{'fraction_non_alphanumeric': 0.04216825061598029, 'fraction_numerical': 0.031608588525167194, 'mean_word_length': 4.947456562696253, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 12, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Normal oil body formation in Marchantia polymorpha requires functional coat protein complex I proteins

Eukaryotic cells possess endomembrane organelles equipped with specific sets of proteins, lipids, and polysaccharides that are fundamental for realizing each organelle’s specific function and shape. A tightly regulated membrane trafficking system mediates the transportation and localization of these substances. Generally, the secretory/exocytic pathway is responsible for transporting cargo to the plasma membrane and/or the extracellular space. However, in the case of oil body cells in the liverwort Marchantia polymorpha, the oil body, a liverwort-unique organelle, is thought to be formed by secretory vesicle fusion through redirection of the secretory pathway inside the cell. Although their formation mechanism remains largely unclear, oil bodies exhibit a complex and bumpy surface structure. In this study, we isolated a mutant with spherical oil bodies through visual screening of mutants with abnormally shaped oil bodies. This mutant harbored a mutation in a coat protein complex I (COPI) subunit MpSEC28, and a similar effect on oil body morphology was also detected in knockdown mutants of other COPI subunits. Fluorescently tagged MpSEC28 was localized to the periphery of the Golgi apparatus together with other subunits, suggesting that it is involved in retrograde transport from and/or in the Golgi apparatus as a component of the COPI coat. The Mpsec28 mutants also exhibited weakened stiffness of the thalli, suggesting impaired cell–cell adhesion and cell wall integrity. These findings suggest that the mechanism of cell wall biosynthesis is also involved in shaping the oil body in M. polymorpha, supporting the redirection of the secretory pathway inward the cell during oil body formation.


Introduction
The trafficking system between endomembranous organelles, including the plasma membrane (PM) in eukaryotic cells, is responsible for the precise transportation and localization of proteins, lipids, and polysaccharides, and is referred to as membrane traffic. The membrane trafficking system mainly consists of the secretory, endocytic, and lysosomal/ vacuolar transport pathways. The secretory, endocytic, and lysosomal/vacuolar transport pathways mediate transport from intracellular organelles to the extracellular space or PM, the extracellular space or PM to intracellular organelles, and to the lysosomes/vacuole, respectively (Kanazawa and Ueda, 2017;Aniento et al., 2022). A single round of trafficking between two compartments is initiated by sorting cargo and forming transport vesicles on a donor organelle membrane, generally mediated by SAR1/ARF GTPases and coat protein complexes, such as coat protein complexes I and II (COPI and COPII) and the clathrin coat (Dacks and Robinson, 2017;Bethune and Wieland, 2018;Dell'Angelica and Bonifacino, 2019). RAB GTPases mediate the targeting and tethering of transport vesicles to target organelle membranes. SNARE proteins execute membrane fusion between these membranes to unload cargo molecules (Kanazawa and Ueda, 2017;Aniento et al., 2022;Ito and Uemura, 2022). Among the coat protein complexes, COPI acts in retrograde transport from the Golgi apparatus to the endoplasmic reticulum (ER) and in intra-Golgi trafficking. The COPI coat is composed of hetero-heptameric subunits α, β, β′, γ, δ, ε, and ζ, which are highly conserved among eukaryotes (Gao et al., 2014;Ito et al., 2014;Dodonova et al., 2015;Dacks and Robinson, 2017).
Although the fundamental molecular framework and major membrane trafficking pathways are highly conserved among eukaryotic lineages, some lineages have also acquired lineagespecific organelles, which should be accompanied by the development of new trafficking pathways to or from these organelles. Comparative genomics analyses have inferred the molecular mechanisms underlying the emergence of new endomembranous organelles and membrane trafficking pathways and detected lineagespecific acquisition of key machinery components, such as coat protein complexes, RAB GTPases, and SNARE proteins. Given that practically all extant membrane trafficking pathways occur through the actions of these machinery components, acquiring a novel membrane trafficking pathway should be accompanied by obtaining a new set of these components. Attaining these components with differentiated functions is thought to be accomplished through gene duplication followed by neo-and/or sub-functionalization by accumulated mutations (Dacks and Field, 2007;Schlacht et al., 2014). Co-opting evolutionarily conserved machinery components with distinct membrane trafficking pathways would be another possibility to explore new trafficking pathways. These theoretical frameworks have been proposed for many years, based on translineage phylogenetic analyses. However, empirical support for these frameworks is limited.
The liverwort oil body is a striking example of a lineage-unique organelle specifically acquired by distinctive eukaryotic lineages during evolution. Notably, the liverwort oil body is an organelle independent from the oil body in seed plants (also known as the oleosome or lipid body) in regard to its origin or function; the liverwort oil body is surrounded by the phospholipid bilayer membrane similar to endomembrane organelles and probably acts in chemical protection against microbes and herbivores, whereas the seed plant oil body is surrounded by a lipid monolayer derived from the ER and accumulates neutral lipids as energy sources (Huang, 1992;He et al., 2013;Huang et al., 2013;Kanazawa et al., 2020;Romani et al., 2020Romani et al., , 2022. The morphology and distribution patterns of the liverwort oil body, such as the shape, color, and number in a cell, vary depending on the liverwort species and have thus been used as an index for taxonomically classifying liverworts (Müller, 1939;He et al., 2013;Romani et al., 2022). In the liverwort Marchantia polymorpha, one oil body is formed in a specific idioblastic cell, referred to as the oil body cell. Recent analyses have revealed that at least two key transcription factors, MpC1HDZ and MpERF13, mediate oil body cell differentiation and oil body development in M. polymorpha Romani et al., 2020Romani et al., , 2022. The Qa-SNARE protein MpSYP12B, which is homologous to the PM-resident MpSYP13B protein, is targeted to the oil-body membrane (Kanazawa et al., 2016). Artificial secretory cargo, sec-mRFP, is transported to the oil body lumen when expressed under MpSYP12B promoter regulation in oil body cells, whereas sec-mRFP, driven by the MpSYP13B promoter, is secreted to the extracellular space . Alongside other supporting evidence, it was proposed that cellular state-specific redirection of the secretory pathway forms the oil body.
Genome information and a series of versatile molecular genetic tools are available for M. polymorpha studies (Ishizaki et al., 2016;Bowman et al., 2017;Montgomery et al., 2020;Kohchi et al., 2021;Bowman, 2022;Bowman et al., in press), including the high-efficiency Agrobacterium-mediated transformation of sporelings (Ishizaki et al., 2008). Using this transformation method can generate several transgenic plants in which transfer DNA (T-DNA) is randomly inserted into the genome, enabling mutants defective in gene functions involved in various biological processes, such as air chamber formation , rhizoid development (Honkanen et al., , 2018Proust et al., 2016), lateral organ development (Naramoto et al., 2019), and oil body formation , to be screened. In this study, to obtain further insights into oil body biogenesis, we screened a T-DNA-mutagenized population of M. polymorpha for mutants with abnormally shaped oil bodies and isolated a mutant with oil bodies with increased circularity. By analyzing the gene responsible for this mutation and other related genes, we demonstrated that normal oil body formation in M. polymorpha requires COPI-mediated secretory activity.
The T-DNA flanking sequences were amplified by TAIL-PCR and sequenced, and the obtained sequences were subjected to a BLAST search using MarpolBase. 2

PCR-based genotyping
Sample preparation and PCR for genotyping were performed as described by Kanazawa et al., 2020. Supplementary (Edgar, 2004a(Edgar, , 2004b. Alignment gaps were removed manually or with the Gblocks program version 0.91b (Castresana, 2000;Talavera and Castresana, 2007). The maximum likelihood phylogenetic analyses were performed using PhyML 3.0 under the LG model (Guindon et al., 2010). Bootstrap analyses were performed by resampling 1,000 sets.

Vector construction
Supplementary Tables 2, 3 list the oligonucleotides and plasmids used in this study, respectively. Open reading frames (ORFs) and genomic sequences were amplified by PCR from complementary DNA (cDNA) and genomic DNA prepared from M. polymorpha accession Tak-1, respectively, for entry vector construction. The amplified products were subcloned into pENTR/D-TOPO (Thermo Fisher) according to the manufacturer's instructions. mCherry cDNA with a linker sequence was inserted into the AscI site of the pENTR vector containing the cDNA for MpSEC28 using the In-Fusion HD cloning kit (Clontech), according to the manufacturer's instructions. The PCR-amplified DNA fragment for mCitrine was inserted into the SacI site of pMpGWB301 using an In-Fusion cloning kit to construct the pMpGWB301-mCitrine destination vector harboring the Gateway cassette. Target sequences for genome editing were selected from the CRISPRdirect website 6 (Naito et al., 2015). Double-stranded DNA fragments of the target sequences were subcloned into pMpGE_En04 to construct dual-gRNA vectors (Koide et al., 2020). The dual-gRNA system was used to generate mutants with a large genomic sequence deletion (Sugano et al., 2018). However, gRNA destined for the 3′-end of the MpSEC28 gene did not work in this study, resulting in the mutations indicated in Figure 1F. For the amiRNA experiments, MpMIR160 was used as the backbone according to Flores-Sandoval et al., 2016. The amiRNA sequences were amplified by PCR and subcloned into the pENTR/D-TOPO vector.
Frontiers in Plant Science 04 frontiersin.org

Confocal laser scanning microscopy and image analyses
Five-day-old thalli were used for observation. For BODIPY staining, thalli were incubated in 200 nM BODIPY 493/503, dissolved in water for 30 min, washed twice in water, and used for microscopic observation. For β-estradiol treatment, gemmae were grown on a 1/2 × Gamborg's B5 medium plate for 3 days, then soaked in a liquid medium containing 20 μM β-estradiol or 0.1% (v/v) DMSO for 48 h. The samples were mounted in a 1/2 × Gamborg's B5 liquid medium and observed using an LSM780 confocal microscope (Carl Zeiss) equipped with an oil immersion lens (63 ×, numerical aperture = 1.4) and lambda and Airyscan detectors. For spectral imaging, the samples were excited at 488 nm (Argon 488) and 561 nm (DPSS 561-10), and emissions between 482 and 659 nm were collected. For high-resolution imaging using the Airyscan detector, samples were excited at 488 and 561 nm, and the emission was separated using a BP495-550 + BP570-630 filter. The images were acquired using line scanning. Spectral unmixing and Airy processing of the obtained images were performed using Abnormal oil body development in Mpsec28 mutants. (A-D) Confocal images merged with bright-field images of oil bodies in wild-type Tak-1 (A), Mpsec28-1 (B), Mpsec28-2 ge (C), and Mpsec28-3 ge (D). Green and blue pseudo-colors indicate the fluorescence from BODIPY and chlorophyll, respectively. Bars = 10 μm. (E) The circularity of the oil bodies in each genotype. Data are presented as the means ± s.d. Differences between wildtype and mutant plants were statistically significant based on a two-tailed Welch's t-test. ****p < 0.001. The sample numbers were 40, 47, 44, and 42 oil bodies for Tak-1, Mpsec28-1, Mpsec28-2 ge , and Mpsec28-3 ge , respectively. p = 8.23 × 10 −32 for Mpsec28-1, p = 2.73 × 10 −24 for Mpsec28-2 ge , and p = 4.40 × 10 −28 for Mpsec28-3 ge . (F) Schematic representation of the MpSEC28 gene structure and mutations generated in this study. Boxes indicate the UTRs (gray) and coding sequences (black). (G,H) PCR-based genotyping of Tak-1 and Mpsec28-1 (G), Tak-1, Mpsec28-2 ge , and Mpsec28-3 ge (H). (F) Shows the primer (A-C) annealing sites. (I) Relative MpSEC28 mRNA accumulation in Tak-1 and Mpsec28-1 measured by qRT-PCR. MpAPT was used as an internal reference. Data are presented as the means ± s.d. The difference between Tak-1 and Mpsec28-1 was statistically analyzed by a two-tailed Welch's t-test. ****p < 0.001. Three biological replicates were prepared, and experiments were performed three times for each gene. Each biological replicate is indicated by a circle (replicate 1), triangle (replicate 2), or square (replicate 3). (J) Genes encoding COPI complex subunits in Arabidopsis thaliana, Marchantia polymorpha, and Saccharomyces cerevisiae.
Frontiers in Plant Science 05 frontiersin.org ZEN2.3 SP1 software (Carl Zeiss). We calculated the circularity of the oil bodies and analyzed colocalization using ImageJ and the Coloc_2 plugin in ImageJ Fiji (Schindelin et al., 2012), respectively.

Quantitative reverse transcription (qRT)-PCR
The β-estradiol treatment was performed as described in the microscopy section. Total RNA was extracted as described by Kanazawa et al., 2013 andKanazawa et al., 2020. cDNA was synthesized from total RNA using oligo dT(18) primers and SuperScript III Reverse Transcriptase (Thermo Fisher) according to the manufacturer's instructions. qRT-PCR was performed with a LightCycler 96 (Roche) using the KAPA SYBR Fast qPCR Master Mix Kit (KAPA Biosystems) according to the manufacturer's protocol. MpAPT was used as the housekeeping reference for normalization (Saint-Marcoux et al., 2015). Three biological replicates were prepared, and each experiment was repeated two or three times. Semi-quantitative RT-PCR was performed using KOD FX Neo (Toyobo), according to the manufacturer's instructions.

Statistics and reproducibility
We present representative confocal microscopy images of biological and technical replicates. The Source data list the results of quantitative image analyses and qRT-PCR, the number of generated transgenic lines, and the number of observed samples. Statistical analyses were performed using the Excel 365 software (ver. 2,202) or GraphPad Prism 9.3.1. (GraphPad Software), and graphs were drawn using GraphPad Prism 9.3.1.

Semi-quantitative chopping assay
Scheme of the assay is shown in Supplementary Figure 9. Twelve-day-old thalli were placed on 5 mm-gapped blocks. The blocks were placed on an electronic force balance, and the thalli were pressed using tweezers to record the force that severed the thalli.

Isolating a T-DNA-inserted mutant with abnormally shaped oil bodies
The oil bodies of M. polymorpha can be visualized by staining with the lipophilic dye BODIPY 493/503 . This allowed us to screen mutants with aberrantly shaped oil bodies. One of these mutants, Mpsec28-1 (identification of the responsible gene is described below), exhibited rounded oil body morphology. Oil bodies in wild-type plants exhibited a bumpy and complex morphology, with a mean circularity of 0.445 ± 0.0937 (mean ± s.d.; Figures 1A,E). The circularity of the oil body in the Mpsec28-1 mutant was 0.821 ± 0.0740, which was significantly higher than that of the wild-type (Figures 1B,E). This phenotype was heritable. The abnormal oil body phenotype was detected in 109 of 224 F1 plants obtained by crossing the Mpsec28-1 mutant with wild-type plants (p = 0.6885 in χ-squared test). This oil body phenotype was linked to the hygromycin resistance conferred by the HPT gene in the T-DNA, indicating that the phenotype was caused by a single nuclear mutation resulting from T-DNA insertion.
We performed a thermal asymmetric interlaced (TAIL)polymerase chain reaction (PCR), sequenced the obtained DNA fragments, and conducted a BLAST search against the M. polymorpha genome sequence to identify the T-DNA insertion site. T-DNA was inserted into the 5′ untranslated region (UTR) of Mp1g28220 ( Figures 1F,G). Mp1g28220 is annotated as the gene encoding the epsilon COPI subunit orthologous to Sec28 and ε-COP in Saccharomyces cerevisiae and mammals, respectively, which was supported by our maximum likelihood phylogenetic analysis ( Figure 1J; Supplementary Figure 1; Hara-Kuge et al., 1994;Duden et al., 1998). Thus, we named this gene MpSEC28. The Mpsec28-1 allele was likely a knockdown allele, because mRNA accumulation was detected by qRT-PCR at a significantly reduced level compared to that in wild-type ( Figure 1I). We generated other mutant alleles (Mpsec28-2 ge and Mpsec28-3 ge ) using CRISPR-based genome editing, which should be knockout mutants because the start codon was deleted in these mutants ( Figures 1F,H), to confirm that a mutation in MpSEC28 caused the abnormally shaped oil body phenotype. The mean circularities of the oil bodies in Mpsec28-2 ge and Mpsec28-3 ge were 0.785 ± 0.1173 and 0.790 ± 0.0875, respectively, similar to that of Mpsec28-1 ( Figure 1E). These results indicated that normal oil body formation requires MpSEC28.

Defects in COPI function result in abnormally shaped oil bodies
The COPI coat consists of α, β, β′, γ, δ, ε, and ζ subunits and mediates intra-Golgi transport and retrograde transport from the cis-Golgi to the ER ( Figure 1J; Glick and Nakano, 2009;Ito et al., 2014;Aniento et al., 2022). We then examined whether the oil body phenotype was specifically caused by defective MpSEC28 function or whether defective functions of other COPI subunits exerted a similar effect. In S. cerevisiae, six genes encoding COPI subunits other than SEC28 are essential for viability (Hosobuchi et al., 1992;Duden et al., 1994Duden et al., , 1998Letourneur et al., 1994;Cosson et al., 1996), which could also be the case for M. polymorpha COPI subunits. Therefore, we generated β-estradiol-inducible knockdown mutants of the αand γ-subunits of the COPI-coat complex using artificial microRNA (amiRNA) for MpRET1 and MpSEC21, respectively (Supplementary Figures 2, 3). The constructs XVE > > amiR-MpRET1-1, XVE> > amiR-MpRET1-2, XVE> > amiR-MpSEC21-1, or XVE> > amiR-MpSEC21-2 were introduced into plants expressing mGFP-MpSYP12B, an oil-body  Figure 4). These results strongly suggest that the abnormal oil body phenotype can be attributed to defective COPI function. We then examined the effects of impairments in the function of COPI components on thallus growth. The thalli of Mpsec28 mutants were viable but were smaller than those of the wild-type ( Figures 3A-E). MpRET1 and MpSEC21 knockdown conferred severe effects. Five-day-old thalli of XVE > > amiR-MpRET1 and XVE > > amiR-MpSEC21 plants grown on the normal medium exhibited no detectable abnormalities in thallus development. However, following 7 days of incubation on the medium containing β-estradiol, these plants showed severely impaired growth, whereas mock treatment did not affect thallus growth ( Figures 3F-R and Supplementary Figure 5). The more severe effects of impairments in MpRET1 and MpSEC21 functions compared to the loss of function of MpSEC28 could reflect a difference in essentiality between these genes, as shown for orthologs in S. cerevisiae (Hosobuchi et al., 1992;Letourneur et al., 1994;Gerich et al., 1995;Duden et al., 1998;Ishii et al., 2016).

Effects of Mpsec28 mutations on thalli stiffness and transport to the oil body
MpSEC28-mCherry and MpSEC21-mCitrine were expressed in almost all cells in five-day-old thalli, including oil body cells (Supplementary Figure 6). The mRNA levels were not significantly changed in the Mperf13-1 ge mutant, which possessed no oil body cells (Supplementary Figure 6C), indicating that MpSEC28 also acts in non-oil body cells. We noted that the thalli of Mpsec28 mutants were fragile during routine handling, suggesting that Mpsec28 mutations altered thalli stiffness. We performed a semi-quantitative chopping assay using 12-dayold thalli (Supplementary Figure 9). The mean force required to sever wild-type thalli was higher than that for Mpsec28 mutant thalli ( Figure 6A), implying that Mpsec28 mutations affect cell wall integrity and cell-cell adhesion, probably due to the faulty secretion of substances required for cell wall integrity. Given that redirection of the secretory pathway is thought to form the oil body, transport to the oil body could also be compromised in Mpsec28 mutants. However, we did not detect markedly affected transport of a general secretion marker sec-mRFP driven by the MpSYP12B promoter or the oil-body-resident membrane protein mGFP-MpSYP12B to the oil body in Mpsec28 mutants ( Figures 6B-G). These results suggest that partially defective COPI coat function confers selective effects on secretory activity, which might be required for cell wall integrity and oil body morphogenesis in M. polymorpha.

Discussion
The secretory pathway begins with transport from the ER to the Golgi apparatus. Most de novo synthesized secretory and vacuolar proteins begin their journey from the ER, mediated by COPII-coated vesicles (Barlowe et al., 1994;Bethune and Wieland, 2018). Subpopulations of ER-and Golgi-resident proteins cycle between the ER and Golgi, and COPI-coated vesicles mediate Frontiers in Plant Science 08 frontiersin.org Data are presented as the means ± s.d. Significant differences between mGFP-MpSYP12B and each of the amiR plants were statistically analyzed using a two-tailed Welch's t-test. n.s., not significant; ****p < 0.001. Twenty (5-day-old thalli) or ten (mock and β-estradiol) thalli were analyzed.
Frontiers in Plant Science 10 frontiersin.org body morphology similar to that in Mpsec28 mutants. Intriguingly, transport of sec-mRFP and mGFP-MpSYP12B was not markedly affected by Mpsec28 mutations. This limited effect of mutations in MpSEC28 on transport of these cargos could reflect the limited requirement of MpSEC28 protein for the total secretory activity in M. polymorpha, which is consistent to what is known in yeast; Sec28 is an only subunit that is dispensable for viability among COPI subunits (Duden et al., 1998). Consistently, knockdown mutants of Mpret1 and Mpsec21 exhibited severer developmental inhibition than Mpsec28 mutants. Whereas the shared effect on the oil body morphology among mutations in COPI subunit genes suggests that the COPI activity is critical for normal oil body formation, we do not rule out the possibility that the severe inhibition of plant growth indirectly affected the oil body development when MpRET1 or MpSEC21 was knocked down.
It is interesting to ask why the mutants defective in COPI function exhibited abnormal oil body morphology. Mechanical support is required to shape and maintain the textured surface of the oil body in M. polymorpha. Given that the oil body lumen should be topologically equivalent to the extracellular space , a cell wall-like structure might line the luminal side of the oil-body membrane, providing mechanical support to the oil body surface. Mutations in MpSEC28 resulted in thalli fragility suggesting that cell wall integrity requires COPImediated exocytic trafficking. The rounder shape of the oil body in mutants defective in COPI function could be due to the compromised integrity of the lining of the oil body lumen.
Arf GTPase activation triggers COPI coat assembly in yeast and animal cells and is catalyzed by Arf guanine nucleotide exchange factors (Arf GEFs; Stearns et al., 1990;Serafini et al., 1991;Waters et al., 1991;Donaldson et al., 1992;Hara-Kuge et al., 1994;Peyroche et al., 1996;Donaldson and Jackson, 2011;Nielsen, 2020). In Arabidopsis, two Arf GEFs with partially redundant functions, GNOM and GNL1, have been shown to act in the early secretory pathway by regulating COPI-coated vesicle formation (Richter et al., 2007;Teh and Moore, 2007;Naramoto et al., 2014;Singh et al., 2018). Calli of a mutant GNOM allele (emb30) exhibit defective cellcell adhesion and increased fragility (Shevell et al., 2000), suggesting that cell-cell adhesion requires COPI-mediated trafficking in Arabidopsis. This is not surprising, given that the Golgi apparatus in plants is the site for pectic polysaccharide biosynthesis, which are crucial in cell-cell adhesion and the mechanical stiffness of plant tissues (Atmodjo et al., 2013;Hoffmann et al., 2021;Codjoe et al., 2022). In M. polymorpha, these Golgi apparatus functions should be conserved, which could account for the fragile nature of Mpsec28 mutants. Notably, major secretory activity, including the transport of PM proteins and general secretion, is redirected toward the oil body during oil body formation in oil body cells . In such situations, the transport of cell wall-related enzymes and polysaccharides could also be redirected toward the oil body, accommodating the polysaccharide lining with the oil-body membrane to confer mechanical strength for the formation and/or maintenance of the complex and bumpy shape of the oil body. Verifying this possibility and unraveling the precise role of the COPI coat in oil body formation would be an interesting future project. SEC28/ε-COP knockdown causes mislocalization of Golgiresident Arabidopsis endomembrane protein 12 (AtEMP12), which harbors the KXD/E motif that mediates the interaction with the COPI coat at the C-terminal cytosolic tail (Woo et al., 2015). In this plant, the number of cisternae in the Golgi stack is significantly decreased compared to the Golgi in wild-type plants, suggesting that COPI is crucial in maintaining Golgi functions in plant cells. In our previous study, we proposed that the oil body is formed through redirection of the secretory pathway inward of the cell , but the involvement of the Golgi apparatus and machinery components of membrane traffic in oil body formation remains elusive. In this study, we showed that the COPI coat is pivotal in normal oil body formation, strongly suggesting that the Golgi-dependent transport pathway substantially contributes to oil body formation. The results obtained in this study further support the proposal that the redirection of the secretory pathway forms the oil body in M. polymorpha. However, it is still possible that other subcellular activities dependent on the ER-Golgi trafficking are also involved in oil body formation. Further studies are needed to elucidate precise molecular and regulatory mechanisms of oil body formation in liverworts.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions
TK performed the experiments and data analyses. RN designed and constructed vectors. TU supervised the study. TK and TU wrote the manuscript. All authors contributed to the article and approved the submitted version. Effects of Mpsec28 mutations on thalli stiffness and transport to the oil body. (A) Semi-quantitative chopping assay. The force required to sever 12-day-old thalli was measured. Data are presented as the means ± s.d. Significant differences between Tak-1 and each genotype were statistically analyzed by a two-tailed Welch's t-test. ***p < 0.005, ****p < 0.001. Ten thalli were analyzed for each genotype. p = 1.33 × 10 −03 for Mpsec28-1, p = 3.21 × 10 −05 for Mpsec28-2 ge , and p = 1.79 × 10 −04 for Mpsec28-3 ge . (B-E) Single confocal images of M. polymorpha thallus cells expressing sec-mRFP in Tak-1 (B), Mpsec28-1 (C), Mpsec28-2 ge (D), and Mpsec28-3 ge (E) under the control of the MpSYP12B promoter. (F,G) Single confocal images of M. polymorpha thallus cells expressing mGFP-MpSYP12B under the control of its regulatory elements in Tak-1 (F) and Mpsec28-1 (G). Magenta, green, and blue pseudo-colors indicate fluorescence from mRFP, mGFP, and chlorophyll, respectively. Bars = 10 μm. Representative images are shown among replicates. The Source data lists the number of observed cells.
Frontiers in Plant Science 12 frontiersin.org

Funding
This work was financially supported by Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology of Japan [19H05670, 19H05675, 21H02515 (TU), and 22K06272 (TK)].