Single-cell RNA sequencing reveals dynamics of gene expression for 2D elongation and 3D growth in Physcomitrium patens.

The transition from two-dimensional (2D) to 3D growth likely facilitated plants to colonize land, but its heterogeneity is not well understood. In this study, we utilized single-cell RNA sequencing to analyze the moss Physcomitrium patens, whose morphogenesis involves a transition from 2D to 3D growth. We profiled over 17,000 single cells covering all major vegetative tissues, including 2D filaments (chloronema and caulonema) and 3D structures (bud and gametophore). Pseudotime analyses revealed larger numbers of candidate genes that determine cell fates for 2D tip elongation or 3D bud differentiation. Using weighted gene co-expression network analysis, we identified a module that connects ß-type carbonic anhydrases (ßCAs) with auxin. We further validated the cellular expression patterns of ßCAs and demonstrated their roles in 3D gametophore development. Overall, our study provides insights into cellular heterogeneity in a moss and identifies molecular signatures that underpin the 2D-to-3D growth transition at single-cell resolution.

A translation proofreader of archaeal origin imparts multi-aldehyde stress tolerance to land plants.

Aldehydes, being an integral part of carbon metabolism, energy generation, and signalling pathways, are ingrained in plant physiology. Land plants have developed intricate metabolic pathways which involve production of reactive aldehydes and its detoxification to survive harsh terrestrial environments. Here, we show that physiologically produced aldehydes, i.e., formaldehyde and methylglyoxal in addition to acetaldehyde, generate adducts with aminoacyl-tRNAs, a substrate for protein synthesis. Plants are unique in possessing two distinct chiral proofreading systems, D-aminoacyl-tRNA deacylase1 (DTD1) and DTD2, of bacterial and archaeal origins, respectively. Extensive biochemical analysis revealed that only archaeal DTD2 can remove the stable D-aminoacyl adducts on tRNA thereby shielding archaea and plants from these system-generated aldehydes. Using Arabidopsis as a model system, we have shown that the loss of DTD2 gene renders plants susceptible to these toxic aldehydes as they generate stable alkyl modification on D-aminoacyl-tRNAs, which are recycled only by DTD2. Bioinformatic analysis identifies the expansion of aldehyde metabolising repertoire in land plant ancestors which strongly correlates with the recruitment of archaeal DTD2. Finally, we demonstrate that the overexpression of DTD2 offers better protection against aldehydes than in wild type Arabidopsis highlighting its role as a multi-aldehyde detoxifier that can be explored as a transgenic crop development strategy.

Distinct localization of chiral proofreaders resolves organellar translation conflict in plants.

Plants have two endosymbiotic organelles originated from two bacterial ancestors. The transition from an independent bacterium to a successful organelle would have required extensive rewiring of biochemical networks for its integration with archaeal host. Here, using Arabidopsis as a model system, we show that plant D-aminoacyl-tRNA deacylase 1 (DTD1), of bacterial origin, is detrimental to organellar protein synthesis owing to its changed tRNA recognition code. Plants survive this conflict by spatially restricting the conflicted DTD1 to the cytosol. In addition, plants have targeted archaeal DTD2 to both the organelles as it is compatible with their translation machinery due to its strict D-chiral specificity and lack of tRNA determinants. Intriguingly, plants have confined bacterial-derived DTD1 to work in archaeal-derived cytosolic compartment whereas archaeal DTD2 is targeted to bacterial-derived organelles. Overall, the study provides a remarkable example of the criticality of optimization of biochemical networks for survival and evolution of plant mitochondria and chloroplast.

Dirigent protein subfamily function and structure in terrestrial plant phenol metabolism.

Aquatic plant transition to land, and subsequent terrestrial plant species diversification, was accompanied by the emergence and massive elaboration of plant phenol chemo-diversity. Concomitantly, dirigent protein (DP) and dirigent-like protein subfamilies, derived from large multigene families, emerged and became extensively diversified. DP biochemical functions as gateway entry points into new and diverse plant phenol skeletal types then markedly expanded. DPs have at least eight non-uniformly distributed subfamilies, with different DP subfamily members of known biochemical/physiological function now implicated as gateway entries to lignan, lignin, aromatic diterpenoid, pterocarpan and isoflavene pathways. While some other DP subfamily members have jacalin domains, both these and indeed the majority of DPs throughout the plant kingdom await discovery of their biochemical roles. Methods and approaches were developed to discover DP biochemical function as gateway entry points to distinct plant phenol skeletal types in land plants. Various DP 3D X-ray structural determinations enabled structure-based comparative sequence analysis and modeling to understand similarities and differences among the different DP subfamilies. We consider that the core DP ß-barrel fold and associated characteristics are likely common to all DPs, with several residues conserved and nearly invariant. There is also considerable variation in residue composition and topography of the putative substrate binding pockets, as well as substantial differences in several loops, such as the ß1-ß2 loop. All DPs likely bind and stabilize quinone methide intermediates, while guiding distinctive regio- and/or stereo-chemical entry into Nature's chemo-diverse land plant phenol metabolic classes.

PIN-FORMED is required for shoot phototropism/gravitropism and facilitates meristem formation in Marchantia polymorpha.

PIN-FORMED auxin efflux transporters, a subclass of which is plasma membrane-localised, mediate a variety of land-plant developmental processes via their polar localisation and subsequent directional auxin transport. We provide the first characterisation of PIN proteins in liverworts using Marchantia polymorpha as a model system. Marchantia polymorpha possesses a single PIN-FORMED gene, whose protein product is predicted to be plasma membrane-localised, MpPIN1. To characterise MpPIN1, we created loss-of-function alleles and produced complementation lines in both M. polymorpha and Arabidopsis. In M. polymorpha, gene expression and protein localisation were tracked using an MpPIN1 transgene encoding a translationally fused fluorescent protein. Overexpression of MpPIN1 can partially complement loss of an orthologous gene, PIN-FORMED1, in Arabidopsis. In M. polymorpha, MpPIN1 influences development in numerous ways throughout its life cycle. Most notably, MpPIN1 is required to establish gemmaling dorsiventral polarity and for orthotropic growth of gametangiophore stalks, where MpPIN1 is basally polarised. PIN activity is largely conserved within land plants, with PIN-mediated auxin flow providing a flexible mechanism to organise growth. Specifically, PIN is fundamentally linked to orthotropism and to the establishment of de novo meristems, the latter potentially involving the formation of both auxin biosynthesis maxima and auxin-signalling minima.

Peptide signaling through leucine-rich repeat receptor kinases: insight into land plant evolution.

Multicellular organisms need mechanisms for communication between cells so that they can fulfill their purpose in the organism as a whole. Over the last two decades, several small post-translationally modified peptides (PTMPs) have been identified as components of cell-to-cell signaling modules in flowering plants. Such peptides most often influence growth and development of organs not universally conserved among land plants. PTMPs have been matched to subfamily XI leucine-rich repeat receptor-like kinases with > 20 repeats. Phylogenetic analyses, facilitated by recently published genomic sequences of non-flowering plants, have identified seven clades of such receptors with a history back to the common ancestor of bryophytes and vascular plants. This raises a number of questions: When did peptide signaling arise during land plant evolution? Have orthologous peptide-receptor pairs preserved their biological functions? Has peptide signaling contributed to major innovations, such as stomata, vasculature, roots, seeds, and flowers? Using genomic, genetic, biochemical, and structural data and non-angiosperm model species, it is now possible to address these questions. The vast number of peptides that have not yet found their partners suggests furthermore that we have far more to learn about peptide signaling in the coming decades.

An optimized transformation protocol for Anthoceros agrestis and three more hornwort species.

Land plants comprise two large monophyletic lineages, the vascular plants and the bryophytes, which diverged from their most recent common ancestor approximately 480 million years ago. Of the three lineages of bryophytes, only the mosses and the liverworts are systematically investigated, while the hornworts are understudied. Despite their importance for understanding fundamental questions of land plant evolution, they only recently became amenable to experimental investigation, with Anthoceros agrestis being developed as a hornwort model system. Availability of a high-quality genome assembly and a recently developed genetic transformation technique makes A. agrestis an attractive model species for hornworts. Here we describe an updated and optimized transformation protocol for A. agrestis, which can be successfully used to genetically modify one more strain of A. agrestis and three more hornwort species, Anthoceros punctatus, Leiosporoceros dussii, and Phaeoceros carolinianus. The new transformation method is less laborious, faster, and results in the generation of greatly increased numbers of transformants compared with the previous method. We have also developed a new selection marker for transformation. Finally, we report the development of a set of different cellular localization signal peptides for hornworts providing new tools to better understand the hornwort cell biology.

Physcomitrium patens PpRIC, an ancestral CRIB-domain ROP effector, inhibits auxin-induced differentiation of apical initial cells.

RHO guanosine triphosphatases are important eukaryotic regulators of cell differentiation and behavior. Plant ROP (RHO of plant) family members activate specific, incompletely characterized downstream signaling. The structurally simple land plant Physcomitrium patens is missing homologs of key animal and flowering plant RHO effectors but contains a single CRIB (CDC42/RAC interactive binding)-domain-containing RIC (ROP-interacting CRIB-containing) protein (PpRIC). Protonemal P. patens filaments elongate based on regular division and PpROP-dependent tip growth of apical initial cells, which upon stimulation by the hormone auxin differentiate caulonemal characteristics. PpRIC interacts with active PpROP1, co-localizes with this protein at the plasma membrane at the tip of apical initial cells, and accumulates in the nucleus. Remarkably, PpRIC is not required for tip growth but is targeted to the nucleus to block caulonema differentiation downstream of auxin-controlled gene expression. These observations establish functions of PpRIC in mediating crosstalk between ROP and auxin signaling, which contributes to the maintenance of apical initial cell identity.

A nuclear Pandora's box: functions of nuclear envelope proteins in cell division.

The nucleus is characteristic of eukaryotic cells and nuclear envelope proteins are conserved across the kingdoms. Over the years, the function of these proteins was studied in the intact nuclear envelope. Knowledge regarding the localization and function of nuclear envelope proteins during mitosis, after the nuclear envelope breaks down, is limited. Until recently, the localization of nuclear envelope proteins during mitosis has been observed with the mitotic apparatus. In this context, research in plant cell biology is more advanced compared to non-plant model systems. Although current studies shed light on the localization of nuclear envelope proteins, further experiments are required to determine what, if any, functional role different nuclear envelope proteins play during mitosis. This review will highlight our current knowledge about the role of nuclear envelope proteins and point out the unanswered questions as future direction.

BAHD Company: The Ever-Expanding Roles of the BAHD Acyltransferase Gene Family in Plants.

Plants' ability to chemically modify core structures of specialized metabolites is the main reason why the plant kingdom contains such a wide and rich array of diverse compounds. One of the most important types of chemical modifications of small molecules is the addition of an acyl moiety to produce esters and amides. Large-scale phylogenomics analyses have shown that the enzymes that perform acyl transfer reactions on the myriad small molecules synthesized by plants belong to only a few gene families. This review is focused on describing the biochemistry, evolutionary origins, and chemical ecology implications of one of these families-the BAHD acyltransferases. The growth of advanced metabolomic studies coupled with next-generation sequencing of diverse plant species has confirmed that the BAHD family plays critical roles in modifying nearly all known classes of specialized metabolites. The current and future outlook for research on BAHDs includes expanding their roles in synthetic biology and metabolic engineering.

Establishment and optimization of a new model organism to study early land plant evolution: Germination, cultivation and oospore variation of Chara braunii Gmelin, 1826.

For studying land plant evolution, the establishment and optimization of model organisms representing streptophytic algae, sister to land plants, is essential. Long-term cultivation experiments with Chara braunii S276 were performed over 8 years, since 4 years (Nov. 2018) under constant conditions. Additionally, short-term experiments for optimization of culture conditions were performed with three strains of C. braunii (S276, NIES-1604 and Lausiger Teiche, LaT-2708). Germination success after application of sterilization agents, addition of gibberellic acid and under different incubation conditions with respect to pre-treatment, irradiance regime and substrate was investigated in order to develop protocols for generative cultivation of at least unialgal cultures. The resulting cultivation protocols for C. braunii S276, allowing maintenance of vegetative as well as generative cultures are presented in detail, including protocols for germination induction and growth of sterilized and unsterilized oospores.

Genome-wide characterization of C2H2 zinc-finger gene family provides insight into the mechanisms and evolution of the dehydration-rehydration responses in Physcomitrium and Arabidopsis.

Dehydration tolerance is a vital factor for land plant evolution and world agricultural production. Numerous studies enlightened that the plant-specific C2H2-type zinc-finger proteins (C2H2-ZFPs) as master regulators played pivotal roles in the abiotic stress responses of plants. However, a comprehensive understanding of the evolution of C2H2-ZFPs in terrestrial plants and its regulatory mechanism in dehydration and rehydration response remains a mystery. In this study, the genome-wide identification of C2H2-ZFP genes revealed 549 homologs in the representatives of terrestrial plant lineages from liverwort to angiosperms. Based on the characteristics of the conserved C2H2-ZF domains, four major C2H2-ZF types (M-, Z-, Q-, and D-type) were identified in the C2H2-ZFPs, with the dominants of M-type in all selected species and followed by Z-type in non-seed plants and Q-type in seed plants, respectively. Phylogenetic analyses of the identified C2H2-ZFPs supported four major groups in the land plant representatives, among which the members from the desiccation-tolerant Physcomitrium patens and the dehydration-sensitive Arabidopsis thaliana displayed different topological relationships in the phylogenies reconstructed for a single species. C2H2-ZFPs clustered in the same subclades shared similar features in their conserved domains and gene structures. Approximately, 81% of the C2H2-ZFP promoters of all 549 identified C2H2-ZFPs harbored the conserved ABA-responsive elements (ABREs) and/or dehydration-responsive elements (DREs). Comparative transcriptomic analyses showed that 50 PpZFPs and 56 AtZFPs significantly changed their transcripts abundance. Interestingly, most of the dehydration- and rehydration-responsive PpZPFs and AtZFPs had been predicted to contain the ABRE and DRE elements in their promoter regions and with over half of which phylogenetically belonging to group III. The differences in the expression patterns of C2H2-ZFPs in responses to dehydration and rehydration between P. patens and A. thaliana reflected their different strategies to adapt to dehydration. The identified candidate PpZFPs were specifically induced by moderate dehydration and reached the peak transcript abundance in severe dehydration. Our study lays the foundations for further functional investigation of C2H2-ZFPs in dehydration responses from an evolutionary perspective in land plants. The findings will provide us with genetic resources and potential targets for drought tolerance breeding in crops and beyond.

KANADI promotes thallus differentiation and FR-induced gametangiophore formation in the liverwort Marchantia.

In angiosperms, KANADI transcription factors have roles in the sporophyte generation regulating tissue polarity, organogenesis and shade avoidance responses, but are not required during the gametophyte generation. Whether these roles are conserved in the gametophyte-dominant bryophyte lineages is unknown, which we examined by characterising the sole KANADI ortholog, MpKAN, in the liverwort Marchantia polymorpha. In contrast to angiosperm orthologs, MpKAN functions in the gametophyte generation in Marchantia, where it regulates apical branching and tissue differentiation, but does not influence tissue polarity in either generation. MpKAN can partially rescue the sporophyte polarity defects of kanadi mutants in Arabidopsis, indicating that MpKAN has conserved biochemical activity to its angiosperm counterparts. Mpkan loss-of-function plants display defects in far-red (FR) light responses. Mpkan plants have reduced FR-induced growth tropisms, have a delayed transition to sexual reproduction and fail to correctly form gametangiophores. Our results indicate that MpKAN is a modulator of FR responses, which may reflect a conserved role for KANADI across land plants. Under FR, MpKAN negatively regulates MpDELLA expression, suggesting that MpKAN and MpDELLA act in a pathway regulating FR responses, placing MpKAN in a gene regulatory network exhibiting similarities with those of angiosperms.

The RGF/GLV/CLEL Family of Short Peptides Evolved Through Lineage-Specific Losses and Diversification and Yet Conserves Its Signaling Role Between Vascular Plants and Bryophytes.

Short secreted plant peptides act as key signaling molecules and control a plethora of developmental and physiological processes. The ROOT GROWTH FACTOR (RGF)/GOLVEN (GLV)/CLE-Like (CLEL) family of peptides was discovered to be involved in root development in Arabidopsis thaliana. In contrast to active research efforts, which have been revealing receptors and downstream signaling components, little attention has been paid to evolutionary processes that shaped the RGF signaling system as we know it in angiosperms today. As a first step toward understanding how RGF signaling emerged and evolved, this study aimed to elucidate the phylogenetic distribution and functional conservation of RGF-like sequences. Using publicly available, genome and transcriptome data, RGF-like sequences were searched in 27 liverworts, 22 mosses, 8 hornworts, 23 lycophytes, 23 ferns, 38 gymnosperms, and 8 angiosperms. This led to the identification of more than four hundreds of RGF-like sequences in all major extant land plant lineages except for hornworts. Sequence comparisons within and between taxonomic groups identified lineage-specific characters. Notably, one of the two major RGF subgroups, represented by A. thaliana RGF6/GLV1/CLEL6, was found only in vascular plants. This subgroup, therefore, likely emerged in a common ancestor of vascular plants after its divergence from bryophytes. In bryophytes, our results infer independent losses of RGF-like sequences in mosses and hornworts. On the other hand, a single, highly similar RGF-like sequence is conserved in liverworts, including Marchantia polymorpha, a genetically tractable model species. When constitutively expressed, the M. polymorpha RGF-like sequence (MpRGF) affected plant development and growth both in A. thaliana and M. polymorpha. This suggests that MpRGF can exert known RGF-like effects and that MpRGF is under transcriptional control so that its potent activities are precisely controlled. These data suggest that RGFs are conserved as signaling molecules in both vascular plants and bryophytes and that lineage-specific diversification has increased sequence variations of RGFs. All together, our findings form a basis for further studies into RGF peptides and their receptors, which will contribute to our understandings of how peptide signaling pathways evolve.

Plasticity in plastid redox networks: evolution of glutathione-dependent redox cascades and glutathionylation sites.

BACKGROUND: Flexibility of plant metabolism is supported by redox regulation of enzymes via posttranslational modification of cysteine residues, especially in plastids. Here, the redox states of cysteine residues are partly coupled to the thioredoxin system and partly to the glutathione pool for reduction. Moreover, several plastid enzymes involved in reactive oxygen species (ROS) scavenging and damage repair draw electrons from glutathione. In addition, cysteine residues can be post-translationally modified by forming a mixed disulfide with glutathione (S-glutathionylation), which protects thiol groups from further oxidation and can influence protein activity. However, the evolution of the plastid glutathione-dependent redox network in land plants and the conservation of cysteine residues undergoing S-glutathionylation is largely unclear. RESULTS: We analysed the genomes of nine representative model species from streptophyte algae to angiosperms and found that the antioxidant enzymes and redox proteins belonging to the plastid glutathione-dependent redox network are largely conserved, except for lambda- and the closely related iota-glutathione S-transferases. Focussing on glutathione-dependent redox modifications, we screened the literature for target thiols of S-glutathionylation, and found that 151 plastid proteins have been identified as glutathionylation targets, while the exact cysteine residue is only known for 17% (26 proteins), with one or multiple sites per protein, resulting in 37 known S-glutathionylation sites for plastids. However, 38% (14) of the known sites were completely conserved in model species from green algae to flowering plants, with 22% (8) on non-catalytic cysteines. Variable conservation of the remaining sites indicates independent gains and losses of cysteines at the same position during land plant evolution. CONCLUSIONS: We conclude that the glutathione-dependent redox network in plastids is highly conserved in streptophytes with some variability in scavenging and damage repair enzymes. Our analysis of cysteine conservation suggests that S-glutathionylation in plastids plays an important and yet under-investigated role in redox regulation and stress response.

PpGRAS12 acts as a positive regulator of meristem formation in Physcomitrium patens.

KEY MESSAGE: This study focused on the key regulatory function of Physcomitrium patens GRAS12 gene underlying an increasing plant complexity, an important step in plant terrestrialization and the evolutionary history of life. The miR171-GRAS module has been identified as a key player in meristem maintenance in angiosperms. PpGRAS12 is a member of the GRAS family and a validated target for miR171 in Physcomitrium (Physcomitrella) patens. Here we show a regulatory function of miR171 at the gametophytic vegetative growth stage and targeted deletion of the PpGRAS12 gene adversely affects sporophyte production since fewer sporophytes were produced in ΔPpGRAS12 knockout lines compared to wild type moss. Furthermore, highly specific and distinct growth arrests were observed in inducible PpGRAS12 overexpression lines at the protonema stage. Prominent phenotypic aberrations including the formation of multiple apical meristems at the gametophytic vegetative stage in response to elevated PpGRAS12 transcript levels were discovered via scanning electron microscopy. The production of multiple buds in the PpGRAS12 overexpression lines similar to ΔPpCLV1a/1b disruption mutants is accompanied by an upregulation of PpCLE and downregulation of PpCLV1, PpAPB, PpNOG1, PpDEK1, PpRPK2 suggesting that PpGRAS12 acts upstream of these genes and negatively regulates the proposed pathway to specify simplex meristem formation. As CLV signaling pathway components are not present in the chlorophytic or charophytic algae and arose with the earliest land plants, we identified a key regulatory function of PpGRAS12 underlying an increasing plant complexity, an important step in plant terrestrialization and the evolutionary history of life.

NO GAMETOPHORES 2 Is a Novel Regulator of the 2D to 3D Growth Transition in the Moss Physcomitrella patens.

The colonization of land by plants was one of the most transformative events in the history of life on Earth. The transition from water, which coincided with and was likely facilitated by the evolution of three-dimensional (3D) growth, enabled the generation of morphological diversity on land. In many plants, the transition from two-dimensional (2D) to 3D growth occurs during embryo development. However, in the early divergent moss Physcomitrella patens, 3D growth is preceded by an extended filamentous phase that can be maintained indefinitely. Here, we describe the identification of the cytokinin-responsive NO GAMETOPHORES 2 (PpNOG2) gene, which encodes a shikimate o-hydroxycinnamoyltransferase. In mutants lacking PpNOG2 function, transcript levels of CLAVATA and SCARECROW genes are significantly reduced, excessive gametophore initial cells are produced, and buds undergo premature developmental arrest. Mutants also exhibit misregulation of auxin-responsive genes. Our results suggest that PpNOG2 functions in the ascorbic acid pathway leading to cuticle formation and that NOG2-related genes were co-opted into the lignin biosynthesis pathway after the divergence of bryophytes and vascular plants. We present a revised model of 3D growth in which PpNOG2 comprises part of a feedback mechanism that is required for the modulation of gametophore initial cell frequency. We also propose that the 2D to 3D growth transition in P. patens is underpinned by complex auxin-cytokinin crosstalk that is regulated, at least in part, by changes in flavonoid metabolism.

Transcriptional and Morpho-Physiological Responses of Marchantia polymorpha upon Phosphate Starvation.

Phosphate (Pi) is a pivotal nutrient that constraints plant development and productivity in natural ecosystems. Land colonization by plants, more than 470 million years ago, evolved adaptive mechanisms to conquer Pi-scarce environments. However, little is known about the molecular basis underlying such adaptations at early branches of plant phylogeny. To shed light on how early divergent plants respond to Pi limitation, we analyzed the morpho-physiological and transcriptional dynamics of Marchantia polymorpha upon Pi starvation. Our phylogenomic analysis highlights some gene networks present since the Chlorophytes and others established in the Streptophytes (e.g., PHR1-SPX1 and STOP1-ALMT1, respectively). At the morpho-physiological level, the response is characterized by the induction of phosphatase activity, media acidification, accumulation of auronidins, reduction of internal Pi concentration, and developmental modifications of rhizoids. The transcriptional response involves the induction of MpPHR1, Pi transporters, lipid turnover enzymes, and MpMYB14, which is an essential transcription factor for auronidins biosynthesis. MpSTOP2 up-regulation correlates with expression changes in genes related to organic acid biosynthesis and transport, suggesting a preference for citrate exudation. An analysis of MpPHR1 binding sequences (P1BS) shows an enrichment of this cis regulatory element in differentially expressed genes. Our study unravels the strategies, at diverse levels of organization, exerted by M. polymorpha to cope with low Pi availability.

With Over 60 Independent Losses, Stomata Are Expendable in Mosses.

Because stomata in bryophytes are uniquely located on sporangia, the physiological and evolutionary constraints placed on bryophyte stomata are fundamentally different from those on leaves of tracheophytes. Although losses of stomata have been documented in mosses, the extent to which this evolutionary process occurred remains relatively unexplored. We initiated this study by plotting the known occurrences of stomata loss and numbers per capsule on the most recent moss phylogeny. From this, we identified 40 families and 74 genera that lack stomata, of which at least 63 are independent losses. No trends in stomata losses or numbers are evident in any direction across moss diversity. Extant taxa in early divergent moss lineages either lack stomata or produce pseudostomata that do not form pores. The earliest land plant macrofossils from 400 ma exhibit similar sporangial morphologies and stomatal distribution to extant mosses, suggesting that the earliest mosses may have possessed and lost stomata as is common in the group. To understand why stomata are expendable in mosses, we conducted comparative anatomical studies on a range of mosses with and without stomata. We compared the anatomy of stomate and astomate taxa and the development of intercellular spaces, including substomatal cavities, across mosses. Two types of intercellular spaces that develop differently are seen in peristomate mosses, those associated with stomata and those that surround the spore sac. Capsule architecture in astomate mosses ranges from solid in the taxa in early divergent lineages to containing an internal space that is directly connected to the conducing tissue and is involved in capsule expansion and the nourishment, hydration and development of spores. This anatomy reveals there are different architectural arrangements of tissues within moss capsules that are equally effective in accomplishing the essential processes of sporogenesis and spore dispersal. Stomata are not foundational to these processes.