Eukaryotic life cycles alternate between haploid and diploid phases and in phylogenetically diverse unicellular eukaryotes, expression of paralogous homeodomain genes in gametes primes the haploid-to-diploid transition. In the unicellular chlorophyte alga Chlamydomonas, KNOX and BELL TALE-homeodomain genes mediate this transition. We demonstrate that in the liverwort Marchantia polymorpha, paternal (sperm) expression of three of five phylogenetically diverse BELL genes, MpBELL234, and maternal (egg) expression of both MpKNOX1 and MpBELL34 mediate the haploid-to-diploid transition. Loss-of-function alleles of MpKNOX1 result in zygotic arrest, whereas a loss of either maternal or paternal MpBELL234 results in variable zygotic and early embryonic arrest. Expression of MpKNOX1 and MpBELL34 during diploid sporophyte development is consistent with a later role for these genes in patterning the sporophyte. These results indicate that the ancestral mechanism to activate diploid gene expression was retained in early diverging land plants and subsequently co-opted during evolution of the diploid sporophyte body.
Eukaryotic life cycles alternate between haploid and diploid phases and in phylogenetically diverse unicellular eukaryotes, expression of paralogous homeodomain genes in gametes primes the haploid-to-diploid transition. In the unicellular chlorophyte alga Chlamydomonas, KNOX and BELL TALE-homeodomain genes mediate this transition. We demonstrate that in the liverwort Marchantia polymorpha, paternal (sperm) expression of three of five phylogenetically diverse BELL genes, MpBELL234, and maternal (egg) expression of both MpKNOX1 and MpBELL34 mediate the haploid-to-diploid transition. Loss-of-function alleles of MpKNOX1 result in zygotic arrest, whereas a loss of either maternal or paternal MpBELL234 results in variable zygotic and early embryonic arrest. Expression of MpKNOX1 and MpBELL34 during diploid sporophyte development is consistent with a later role for these genes in patterning the sporophyte. These results indicate that the ancestral mechanism to activate diploid gene expression was retained in early diverging land plants and subsequently co-opted during evolution of the diploid sporophyte body.
Life cycles of eukaryotes alternate between haploid and diploid phases, initiated by meiosis and gamete fusion, respectively. Expression of paralogous homeodomain (HD) genes in the two gametes and the subsequent heterodimerization of the respective proteins in the zygote direct the haploid-to-diploid transition in gene expression in phylogenetically diverse eukaryotes, including the ascomycete fungus Saccharomyces cerevisiae (Goutte and Johnson, 1988; Herskowitz, 1989), the basidiomycete fungi Coprinopsis cinerea and Ustilago maydis (Gillissen et al., 1992; Hull et al., 2005; Kues et al., 1992; Spit et al., 1998, Urban, 1996), the Amoebozoa Dictyostelium discoideum (Hedgethorne et al., 2017), the brown alga Ectocarpus (Arun et al., 2019), the red alga Pyropia yezoensis (Mikami et al., 2019), and the unicellular chlorophyte alga Chlamydomonas reinhardtii (Ferris and Goodenough, 1987; Lee et al., 2008; Nishimura et al., 2012; Zhao et al., 2001). This broad phylogenetic distribution suggests this was an ancestral function of HD genes (reviewed by Bowman et al., 2016b). In Viridiplantae, the paralogs are two subclasses, KNOX and BELL, of TALE class HD genes. In Chlamydomonas, the minus (-) gamete expresses a KNOX protein (GSM1) and the plus (+) gamete expresses a BELL protein (GSP1), and upon gamete fusion, the two proteins heterodimerize and translocate to the nucleus, activating zygotic gene expression (Lee et al., 2008). GSM1 and GSP1 are necessary for diploid gene expression and, when ectopically expressed together in vegetative haploid cells, are sufficient to induce the diploid genetic program (Ferris and Goodenough, 1987; Lee et al., 2008; Nishimura et al., 2012; Zhao et al., 2001). Biologically, the expression of a unique paralog in each type of gamete, coupled with the requirement for heterodimerization for functionality, is a mechanism to preload gametes such that immediately following gamete fusion/fertilization a distinct diploid genetic program is initiated.Heterodimerization of TALE-HD paralogs is mediated by subclass-specific protein domains N-terminal of the homeodomain conserved in phylogenetically diverse eukaryotes (Bellaoui et al., 2001; Bürglin, 1997). For example, Viridiplantae KNOX proteins and metazoan MEIS proteins share a homologous heterodimerization domain, thus defining a TALE-HD subclass, named MEINOX, present in the ancestral eukaryote (Bürglin, 1997). In contrast, heterodimerization domains of MEINOX partners are not as well conserved (Joo et al., 2018). Some MEINOX partners in animals have a characteristic PBC-C heterodimerization domain (Bürglin, 1997), and PBC-related domains are found in some algal BELL-related TALE-HD proteins, suggesting its presence in the ancestral eukaryote (Joo et al., 2018). In contrast, previously characterized land plant BELL proteins possess conserved SKY and BELL domains, collectively termed POX (pre-homeobox), that mediate heterodimerization with KNOX partners (Bellaoui et al., 2001; Hackbusch et al., 2005; Smith et al., 2002). POX domains are either highly divergent from, or unrelated to, PBC domains. Finally, in other algal BELL-related TALE-HD proteins, neither heterodimerization domain is evident (Joo et al., 2018). Regardless, when tested for the dimerization potential, Archaeplastida BELL-related TALE-HD proteins interact only with MEINOX TALE-HD partners and not with other BELL-related TALE-HD proteins (Joo et al., 2018).Land plants are characterized by an alternation of generations whereby complex multicellular bodies develop in both haploid (gametophyte) and diploid (sporophyte) phases of the life cycle (Hofmeister, 1862), and elements of the KNOX/BELL system have been implicated in regulating the land plant haploid-to-diploid transition. In the moss Physcomitrium patens, one subclass of KNOX genes, KNOX1, is required for proper proliferation and differentiation in the diploid body (Sakakibara et al., 2008; Singer and Ashton, 2007), while another subclass, KNOX2, acts to suppress the haploid genetic program during diploid development (Sakakibara et al., 2013). In addition to being expressed during sporophyte development, both classes of KNOX genes are expressed in the egg cell during gamete formation (Sakakibara et al., 2013; Sakakibara et al., 2008). The expression of KNOX genes in egg cells suggests that this gamete may correspond to the (-) gamete in Chlamydomonas, which also expresses a KNOX protein, implying that the male gamete in land plants might express BELL proteins, as does the (+) gamete in Chlamydomonas. Ectopic expression of a Physcomitrium BELL paralog, PpBELL1, has been noted to induce the diploid genetic program in specific cell types of the gametophyte, but PpBELL1 expression was detected in the egg cell (Horst et al., 2016), rather than the male gamete, as might have been expected if it is functionally similar to Chlamydomonas BELL. However, expression of PpBELL1 has also been reported to be induced in the sperm by activation of glutamate channels, with PpBELL1 subsequently functionally active in the sporophyte (Ortiz-Ramírez et al., 2017). These observations prompted our investigation of the homologous genetic programs in the liverwort Marchantia polymorpha.
Results
M. polymorpha possesses phylogenetically diverse TALE homeodomain proteins
The M. polymorpha genome encodes nine TALE-HD-related family members: four KNOX genes and five BELL genes (Bowman et al., 2017; Figure 1), all of which are expressed in either the haploid sexual organs or the diploid sporophyte or both, with minimal or no expression detected by reverse transcription-polymerase chain reaction (RT-PCR) in the haploid vegetative thallus (Figure 1—figure supplement 1). Of the four KNOX genes, three are of KNOX1 subclass, but with only one encoding a HD, and one of KNOX2 subclass (Figure 1—figure supplement 1). The two subclasses arose via gene duplication in an ancestral charophycean alga (Sakakibara, 2016; Joo et al., 2018; Frangedakis et al., 2017). Expression of MpKNOX1 (Mp5g01600), which encodes a HD, is predominantly detected in archegoniophores and young sporophytes. In contrast, the two other KNOX1 genes (MpKNOX1A, Mp4g12450; MpKNOX1B, Mp2g11140) lacking a HD are predominantly expressed in antheridiophores. Phylogenetic analysis indicates that these HD-less paralogs are liverwort specific and are likely evolved in the ancestral liverwort (Figure 1—figure supplement 1). In contrast, expression of MpKNOX2 (Mp7g05320) is not detected in unfertilized reproductive organs and is expressed primarily during sporophyte development.
Figure 1.
Schematic representations of the nine TALE-HD loci of Marchantia polymorpha.
White, 5' and 3' untranslated regions (UTRs); black, coding exons; thin lines, introns; green triangles, start codon; black triangles, stop codon; red triangles, guide RNA-targeted positions. The mutant alleles generated via CRISPR-Cas9 are described in more detail in table 1; molecular lesions of some alleles can be found in Figure 1—source data 1. All protein annotation models are based on the Marchantia genome assembly of v5.1 except for MpKNOX2. The MpKNOX2 model is based on sequences derived from reverse transcription-polymerase chain reaction (RT-PCR). In genes with multiple homeodomains, they are denoted a, b, and c. Gene models were assembled using wormweb (http://wormweb.org/exonintron).
RNA extracted from various tissues of Marchantia polymorpha ssp ruderalis, ecotype BoGa, were subjected to semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) (30 cycles, except MpKNOX1 and MpKNOX2 at 32 cycles; 55°C anneal; 1 min extension, except MpKNOX1, MpKNOX2, and MpBELL1 at 2 min). The expected size band from complementary DNA (cDNA) is indicated by a red arrow at the left of each panel. We were unable to amplify a product from MpBELL2. See (Supplementary file 1) for primer sequences.
Tree constructed using nucleotide alignment of the homeodomain and MEINOX and ELK domains. The four Marchantia polymorpha KNOX genes are highlighted (purple arrows). The tree topology suggests a single KNOX gene in the common ancestor of the Viridiplantae, and possibly the Streptophyta, with a gene duplication in an ancestral charophycean alga after the divergence of Mesostigma, producing ancestral KNOX1 and KNOX2 genes (Sakakibara, 2016; Joo et al., 2018; Frangedakis et al., 2017). The KNOX2 genes, including charophyte sequences, form a well-supported monophyletic clade (posterior probability = 1). The KNOX1 clade is less well-supported with regard to the early diverging charophyte sequences; however, land plant KNOX1 genes, including MpKNOX1, share a conserved MEINOX domain (see also Figure 3—figure supplement 3). The ancestral land plant likely had single KNOX1 and KNOX2 orthologs, a condition retained in M. polymorpha. MpKNOX1A and MpKNOX1B orthologs are found throughout liverworts, suggesting their origins in the ancestral liverwort. Chlorophyte, light blue; charophyte, dark blue; liverwort, purple; moss, green; hornwort, dark green; lycophyte, brown; fern, orange; seed plant, red. Numbers indicate posterior probability values >50%.
Figure 1—figure supplement 1.
Expression of Marchantia polymorpha TALE-HD genes.
RNA extracted from various tissues of Marchantia polymorpha ssp ruderalis, ecotype BoGa, were subjected to semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) (30 cycles, except MpKNOX1 and MpKNOX2 at 32 cycles; 55°C anneal; 1 min extension, except MpKNOX1, MpKNOX2, and MpBELL1 at 2 min). The expected size band from complementary DNA (cDNA) is indicated by a red arrow at the left of each panel. We were unable to amplify a product from MpBELL2. See (Supplementary file 1) for primer sequences.
Schematic representations of the nine TALE-HD loci of Marchantia polymorpha.
White, 5' and 3' untranslated regions (UTRs); black, coding exons; thin lines, introns; green triangles, start codon; black triangles, stop codon; red triangles, guide RNA-targeted positions. The mutant alleles generated via CRISPR-Cas9 are described in more detail in table 1; molecular lesions of some alleles can be found in Figure 1—source data 1. All protein annotation models are based on the Marchantia genome assembly of v5.1 except for MpKNOX2. The MpKNOX2 model is based on sequences derived from reverse transcription-polymerase chain reaction (RT-PCR). In genes with multiple homeodomains, they are denoted a, b, and c. Gene models were assembled using wormweb (http://wormweb.org/exonintron).
Expression of Marchantia polymorpha TALE-HD genes.
RNA extracted from various tissues of Marchantia polymorpha ssp ruderalis, ecotype BoGa, were subjected to semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) (30 cycles, except MpKNOX1 and MpKNOX2 at 32 cycles; 55°C anneal; 1 min extension, except MpKNOX1, MpKNOX2, and MpBELL1 at 2 min). The expected size band from complementary DNA (cDNA) is indicated by a red arrow at the left of each panel. We were unable to amplify a product from MpBELL2. See (Supplementary file 1) for primer sequences.
Unrooted Bayesian phylogram of Viridiplantae KNOX-TALE class homeodomain genes.
Tree constructed using nucleotide alignment of the homeodomain and MEINOX and ELK domains. The four Marchantia polymorpha KNOX genes are highlighted (purple arrows). The tree topology suggests a single KNOX gene in the common ancestor of the Viridiplantae, and possibly the Streptophyta, with a gene duplication in an ancestral charophycean alga after the divergence of Mesostigma, producing ancestral KNOX1 and KNOX2 genes (Sakakibara, 2016; Joo et al., 2018; Frangedakis et al., 2017). The KNOX2 genes, including charophyte sequences, form a well-supported monophyletic clade (posterior probability = 1). The KNOX1 clade is less well-supported with regard to the early diverging charophyte sequences; however, land plant KNOX1 genes, including MpKNOX1, share a conserved MEINOX domain (see also Figure 3—figure supplement 3). The ancestral land plant likely had single KNOX1 and KNOX2 orthologs, a condition retained in M. polymorpha. MpKNOX1A and MpKNOX1B orthologs are found throughout liverworts, suggesting their origins in the ancestral liverwort. Chlorophyte, light blue; charophyte, dark blue; liverwort, purple; moss, green; hornwort, dark green; lycophyte, brown; fern, orange; seed plant, red. Numbers indicate posterior probability values >50%.
Figure 3—figure supplement 3.
MpKNOX1A and MpKNOX1B lack the homeodomain and are both primarily expressed in the antheridiophores.
(A) Alignment of MpKNOX1 MEINOX domains with land plant KNOX1 proteins. A consensus sequence of the MEINOX domain (Bürglin, 1997; Crooks et al., 2004) of land plant KNOX1 proteins was computed from sequences in Figure 1—figure supplement 1, with the alignment of the three Marchantia polymorpha KNOX1 proteins displayed below. The lines above the sequence highlight the previously defined KNOX subdomains of the MEINOX domain. Gray shading indicates amino acids conserved with the land plant consensus and asterisks highlight amino acids conserved in all three M. polymorpha KNOX1 sequences. (B) Panels as described in Figure 3. Expression of both genes is predominantly detected in antheridiophores. (C–E) Translational fusion reporter line expression containing 1 kb (MpKNOX1A) and 2.3 kb (MpKNOX1B) of the putative 5' regulatory sequences suggests that antheridiophore expression identified by RNA-sequencing (RNAseq) is primarily in the antheridia. Expression of a translational fusion of MpKNOX1A with β-glucuronidase (GUS) was detectable in the sperm (D) and in mucilage papillae of the archegoniophore (E), but not in mucilage papillae in other parts of the thallus (not shown). Scale bar = 200 µm (antheridiophores and archegoniophores), 10 µm (sperm), and 50 µm (mucilage cells).
The M. polymorpha BELL genes reside in three phylogenetically distinct clades, each with its own characteristic HD (Figure 2; Figure 2—figure supplement 1). Like some algal BELL-related proteins (Joo et al., 2018), the presence of a PBC domain is not easily discernible in any of the predicted M. polymorpha BELL protein sequences. Only one M. polymorpha BELL gene (MpBELL1, Mp8g18310) is phylogenetically related to previously described, ‘canonical’, land plant BELL genes. MpBELL1 harbors a conserved canonical land plant BELL homeodomain sequence and a discernible, albeit divergent, BELL domain characteristic of other land plant BELL genes. Similar to MpKNOX2, MpBELL1 expression is not detected in reproductive organs prior to fertilization and is expressed predominantly during sporophyte development (Figure 1—figure supplement 1). The other four M. polymorpha BELL genes fall into two distinct clades and are more closely related to algal TALE-HD genes than to the canonical land plant BELL genes (Figure 2). Three of the algal-like BELL genes, MpBELL2 (Mp4g09650), MpBELL3 (Mp8g02970), and MpBELL4 (Mp8g07680), each encode two or three homeodomains (Figure 1). The carboxyl HD is most conserved with those upstream becoming progressively less conserved, consistent with their origins being via successive intragenic duplications. Notably, MpBELL3 and MpBELL4 both encode large proteins, while MpBELL2 encodes a much smaller protein consisting of only two HDs. These three genes, named here BELL-gamete, arose via gene duplications within liverworts, and related genes are found throughout liverworts and in some moss lineages that are sister groups to the vast majority of moss diversity (Figure 2). The fourth algal-like sequence, MpBELL5 (Mp5g11060), is phylogenetically distinct from the other four M. polymorpha BELL genes, and resides in the previously defined ‘GLX-basal’ clade (Joo et al., 2018). Orthologs of MpBELL5 exist in other liverworts (Figure 2). The most parsimonious interpretation is that a diversity of BELL paralogs arose in a charophyte algal ancestor, with the ancestral land plant likely possessing three paralogs representing the canonical land plant, GLX-basal and BELL-gamete lineages (Figure 2). This diversity has persisted in M. polymorpha and other liverworts and, to a lesser extent, in mosses.
Figure 2.
Unrooted Bayesian phylogram of Viridiplantae BELL-TALE class homeodomain genes.
Tree was constructed using a nucleotide alignment of the homeodomain. The five Marchantia polymorpha BELL genes are highlighted (purple arrows). A clade of genes representing previously identified ‘canonical’ land plant BELL genes is recovered with high support. This clade includes MpBELL1 and charophycean algal sequences dating as far back as prior to the divergence of Klebsormidium. These canonical land plant BELL genes also share significant sequence similarity outside the homeodomain, including the SKY and BELL domains. The remaining BELL genes, including all chlorophyte sequences, sequences from charophycean algae, and the other four M. polymorpha sequences, reside in either other well-supported clades or in polytomies. Several subclades (labeled) previously identified are evident (Joo et al., 2018). MpBELL2/3/4 resides in a polytomy composed of liverwort and charophyte sequences, with other liverwort, moss, and charophyte sequences residing in a second clade; these are labeled ‘BELL-gamete’. The multiple homeodomains of MpBELL2, MpBELL3 and MpBELL4 are designated a, b, and c (see Figure 1). MpBELL5 resides, with other liverwort and charophyte sequences, in the GLX-basal clade. Chlorophyte, light blue; charophyte, dark blue; liverwort, purple; moss, green; hornwort, dark green; lycophyte, brown; fern, orange; seed plant, red. Numbers at branches indicate posterior probability values >50%; branches explicitly shown have probability values >50%, whereas polytomies represent nodes with probability values <50%; values within subclades are omitted for clarity.
The alignments used in the construction of phylogenetic trees in Figure 2 and Figure 1—figure supplement 1 were used to generate consensus sequences for nine different classes of Viridiplantae TALE-homeodomain (HD) proteins. Homeodomain is flanked by red arrows. Some TALE-HD classes have insertions relative to others and some subclasses of BELL protein have conserved motifs carboxyl to the HD. The invariant asparagine (N) at position 50 (61 in this figure) within the HD is critical for specific binding to AT-rich sequences, with the amino acids at positions 47, 50, and 54 (57, 60, and 64 in this figure highlighted by blue arrows and gray shading) influencing sequence specificity adjacent to the AT-rich core (Berger et al., 2008; Noyes et al., 2008). Of note is that, based on these critical residues, the three different subclasses of land plant BELL HD proteins are predicted to have distinct DNA-recognition specificities.
Figure 2—figure supplement 1.
Comparison of Viridiplantae TALE homeodomains.
The alignments used in the construction of phylogenetic trees in Figure 2 and Figure 1—figure supplement 1 were used to generate consensus sequences for nine different classes of Viridiplantae TALE-homeodomain (HD) proteins. Homeodomain is flanked by red arrows. Some TALE-HD classes have insertions relative to others and some subclasses of BELL protein have conserved motifs carboxyl to the HD. The invariant asparagine (N) at position 50 (61 in this figure) within the HD is critical for specific binding to AT-rich sequences, with the amino acids at positions 47, 50, and 54 (57, 60, and 64 in this figure highlighted by blue arrows and gray shading) influencing sequence specificity adjacent to the AT-rich core (Berger et al., 2008; Noyes et al., 2008). Of note is that, based on these critical residues, the three different subclasses of land plant BELL HD proteins are predicted to have distinct DNA-recognition specificities.
Unrooted Bayesian phylogram of Viridiplantae BELL-TALE class homeodomain genes.
Tree was constructed using a nucleotide alignment of the homeodomain. The five Marchantia polymorpha BELL genes are highlighted (purple arrows). A clade of genes representing previously identified ‘canonical’ land plant BELL genes is recovered with high support. This clade includes MpBELL1 and charophycean algal sequences dating as far back as prior to the divergence of Klebsormidium. These canonical land plant BELL genes also share significant sequence similarity outside the homeodomain, including the SKY and BELL domains. The remaining BELL genes, including all chlorophyte sequences, sequences from charophycean algae, and the other four M. polymorpha sequences, reside in either other well-supported clades or in polytomies. Several subclades (labeled) previously identified are evident (Joo et al., 2018). MpBELL2/3/4 resides in a polytomy composed of liverwort and charophyte sequences, with other liverwort, moss, and charophyte sequences residing in a second clade; these are labeled ‘BELL-gamete’. The multiple homeodomains of MpBELL2, MpBELL3 and MpBELL4 are designated a, b, and c (see Figure 1). MpBELL5 resides, with other liverwort and charophyte sequences, in the GLX-basal clade. Chlorophyte, light blue; charophyte, dark blue; liverwort, purple; moss, green; hornwort, dark green; lycophyte, brown; fern, orange; seed plant, red. Numbers at branches indicate posterior probability values >50%; branches explicitly shown have probability values >50%, whereas polytomies represent nodes with probability values <50%; values within subclades are omitted for clarity.
Comparison of Viridiplantae TALE homeodomains.
The alignments used in the construction of phylogenetic trees in Figure 2 and Figure 1—figure supplement 1 were used to generate consensus sequences for nine different classes of Viridiplantae TALE-homeodomain (HD) proteins. Homeodomain is flanked by red arrows. Some TALE-HD classes have insertions relative to others and some subclasses of BELL protein have conserved motifs carboxyl to the HD. The invariant asparagine (N) at position 50 (61 in this figure) within the HD is critical for specific binding to AT-rich sequences, with the amino acids at positions 47, 50, and 54 (57, 60, and 64 in this figure highlighted by blue arrows and gray shading) influencing sequence specificity adjacent to the AT-rich core (Berger et al., 2008; Noyes et al., 2008). Of note is that, based on these critical residues, the three different subclasses of land plant BELL HD proteins are predicted to have distinct DNA-recognition specificities.
Maternal MpKNOX1 is required for post-zygotic embryo development
During gametophyte generation, MpKNOX1 is expressed specifically in the egg cell but is not detected at the stage prior, when the venter canal cell is present (Figure 3A–B). During differentiation of the egg cell, MpKNOX1 appears to be expressed in a pulse early, with the signal diminishing in older archegonia (Figure 3B). Following fertilization, MpKNOX1 is expressed throughout the developing sporophyte up until the time when future sporogenous cells become distinct (Figure 3C–E), with expression thereafter becoming undetectable during meiotic stages (Figure 3F). A female harboring a loss-of-function allele consisting of a 2.3 kb deletion spanning the MpKNOX1 locus (Figure 1; Table 1) was crossed with wild-type males. In contrast to control crosses using wild-type males and females where sporophyte production was observed (Figure 3G–I), mature sporophytes never developed on gametangiophores from >100 crosses between wild-type males and Mpknox1-6ge females (Figure 3J–L). However, on closer inspection, on rare occasions (3/129 fertilization events) did we encounter arrested multicellular embryos within senescing archegoniophores (Figure 3—figure supplement 1). A closer examination of female Mpknox1-6 mutants revealed a developmental arrest of the zygote following fertilization, with the zygote failing to undergo cytokinesis (Figure 3K). However, two gametophytic tissues, the calyptra and the pseudoperianth, whose growth is dependent upon fertilization (Hofmeister, 1862), commence development normally. In Mpknox1-6 mutants, as in wild type, the calyptra undergoes periclinal divisions indicating successful fertilization and implying an intergenerational (sporophyte to gametophyte) signal inducing its development (Figure 3H–I and K–L). Likewise, the pseudoperianth, a ring of tissue surrounding the archegonium that only develops post-fertilization, initially develops normally in Mpknox1-6 mutants (Figure 3H–I and K–L). Subsequently, both the calyptra and pseudoperianth are developmentally arrested between 1 and 2 weeks post fertilization (wpf) in Mpknox1-6 mutants. In addition, Mpknox1-6 mutants exhibit a senescent archegoniophore phenotype. Unfertilized wild-type archegoniophores remain green for multiple weeks before they senesce (Figure 3M). In contrast, Mpknox1-6 archegoniophores begin to senesce around 2 weeks post maturation (Figure 3N), suggesting that an MpKNOX1-dependent signal is required for archegoniophore maintenance.
Figure 3.
Maternal MpKNOX1 is required for post-zygotic embryo development.
(A–F) MpKNOX1 expression pattern. MpKNOX1 transcripts could be detected via in situ hybridization in the egg (e) cells of unfertilized archegonia (B), but not in the venter canal cell (A; vcc). After fertilization, expression is detected in young developing sporophytes (approximately 1 week post fertilization (wpf); C, D), with expression continuing until foot and seta begin to differentiate (E). MpKNOX1 is not expressed in older sporophytes where sporogenous tissue has differentiated (F). (G–L) Comparison of wild-type and Mpknox1 development; unfertilized (G, J), 1 wpf (H, K), and 3 wpf (I, L). Mpknox1-6 egg cells appear wild-type-like (unfertilized; G, J). In contrast to wild-type (H, I), the embryo of Mpknox1-6 mutants does not develop post fertilization (K, L). The zygote fails to undergo cytokinesis, with the nucleus disappearing after 3 weeks, leaving an empty space (approximately 3 wpf; L). The initial outgrowth of the pseudoperianth (black arrows) and calyptra (red arrows) post-fertilization is not affected in Mpknox1-6 mutants, but their development is arrested as well approximately 1 wpf (K, L). (M, N) Archegoniophores of Mpknox1-6 mutants begin to senesce approximately 2 weeks after maturation (N), while wild-type archegoniophores of the same age remain green (M). (O) Based on RNA-sequencing (RNA-seq) data and associated gene models, the 3' end of MpKNOX1 overlaps with the 3' untranslated region (UTR) of MpSUK1, which is transcribed from the opposite strand. Predicted full-length MpKNOX1 transcripts, that is, those including exons 2 and 3, are present primarily in the archegoniophore, consistent with the in situ data and semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) data in Figure 1—figure supplement 1. FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) scales are set to a standard within each tissue to allow potential weakly expressed transcripts to be visualized. (P) Signal of MpSUK1 in situ hybridization is detected in antheridia of the antheridiophore. (Q, R) Expression of MpKNOX1:GUS:SUK1 in an archegoniophore (R), an archegonium (S), and an antheridiophore (Q). Scale bar = 50 µm (A, E, G, J, K, L, S); 100 µm (B, C, D, H, I); 200 µm (P); 500 µm (F, R); 1000 µm (Q).
(A–D) Cross-sections of Mpknox1-6 female mutants crossed with wild-type males and observed 2 weeks post fertilization (wpf). Note that some cells are deformed and collapsed, a characteristic of senescing tissue. (B) is the close up of (A). (E) Percentages of developmental stages observed from this cross. Plants were crossed once, and then examined after 1 or 2 wpf. Observed embryos were grouped into the following developmental stages: aborted zygotes, sporophytes that were arrested at a stage younger than the maximum stage possible, and sporophytes that had reached the maximum stage expected. n = total number of observed fertilization events. See Figure 3—figure supplement 1—source data 1 for details. Scale bar = 50 µm (B–D); 200 µm (A).
Sequences in purple are from the endogenous locus, while sequences in other colors are from the heterologous locus. (A) Endogenous MpKNOX1/MpSUK1 locus. (B) A construct, MpKNOX1:GUS, including a 4.6 kb sequence upstream (5') of the MpKNOX1 translational start site fused with the β-glucuronidase (GUS) reporter gene followed by a nopaline synthase (NOS) transcriptional terminator (250 bp), drives expression in both the egg cell (left) and antheridia (right). (C) Inclusion of the MpSUK1 locus, MpKNOX1:GUSterm:SUK1, including 1.5 kb upstream (5') of the MpSUK1 translational start site, downstream of the NOS terminator, results in the same expression pattern, that is, in both the egg cell and antheridia. (D) Only in constructs featuring an assembly in which MpKNOX1:GUS and MpSUK1:MpSUK1 are opposing and not interrupted by a transcriptional terminator (MpKNOX1:GUS:SUK1) is expression of MpKNOX1 restricted to the female egg cell (Figure 3) and the developing sporophyte, shown here at 1 week post fertilization (wpf), 2 wpf, and 4 wpf. No expression is observed in males (Figure 3). Scale bar = 50 µm (archegonia in B, C and 1-week-old sporophyte in D), 250 µm (2- and 4-week-old sporophytes in D), 500 µm (archegoniophores in B, C), and 1000 µm (antheridiophores in B, C and archegoniophores in D).
(A) Alignment of MpKNOX1 MEINOX domains with land plant KNOX1 proteins. A consensus sequence of the MEINOX domain (Bürglin, 1997; Crooks et al., 2004) of land plant KNOX1 proteins was computed from sequences in Figure 1—figure supplement 1, with the alignment of the three Marchantia polymorpha KNOX1 proteins displayed below. The lines above the sequence highlight the previously defined KNOX subdomains of the MEINOX domain. Gray shading indicates amino acids conserved with the land plant consensus and asterisks highlight amino acids conserved in all three M. polymorpha KNOX1 sequences. (B) Panels as described in Figure 3. Expression of both genes is predominantly detected in antheridiophores. (C–E) Translational fusion reporter line expression containing 1 kb (MpKNOX1A) and 2.3 kb (MpKNOX1B) of the putative 5' regulatory sequences suggests that antheridiophore expression identified by RNA-sequencing (RNAseq) is primarily in the antheridia. Expression of a translational fusion of MpKNOX1A with β-glucuronidase (GUS) was detectable in the sperm (D) and in mucilage papillae of the archegoniophore (E), but not in mucilage papillae in other parts of the thallus (not shown). Scale bar = 200 µm (antheridiophores and archegoniophores), 10 µm (sperm), and 50 µm (mucilage cells).
Figure 3—figure supplement 1.
Mpknox1-6 mutants rarely form embryos on senescing archegoniophores.
(A–D) Cross-sections of Mpknox1-6 female mutants crossed with wild-type males and observed 2 weeks post fertilization (wpf). Note that some cells are deformed and collapsed, a characteristic of senescing tissue. (B) is the close up of (A). (E) Percentages of developmental stages observed from this cross. Plants were crossed once, and then examined after 1 or 2 wpf. Observed embryos were grouped into the following developmental stages: aborted zygotes, sporophytes that were arrested at a stage younger than the maximum stage possible, and sporophytes that had reached the maximum stage expected. n = total number of observed fertilization events. See Figure 3—figure supplement 1—source data 1 for details. Scale bar = 50 µm (B–D); 200 µm (A).
Maternal MpKNOX1 is required for post-zygotic embryo development.
(A–F) MpKNOX1 expression pattern. MpKNOX1 transcripts could be detected via in situ hybridization in the egg (e) cells of unfertilized archegonia (B), but not in the venter canal cell (A; vcc). After fertilization, expression is detected in young developing sporophytes (approximately 1 week post fertilization (wpf); C, D), with expression continuing until foot and seta begin to differentiate (E). MpKNOX1 is not expressed in older sporophytes where sporogenous tissue has differentiated (F). (G–L) Comparison of wild-type and Mpknox1 development; unfertilized (G, J), 1 wpf (H, K), and 3 wpf (I, L). Mpknox1-6 egg cells appear wild-type-like (unfertilized; G, J). In contrast to wild-type (H, I), the embryo of Mpknox1-6 mutants does not develop post fertilization (K, L). The zygote fails to undergo cytokinesis, with the nucleus disappearing after 3 weeks, leaving an empty space (approximately 3 wpf; L). The initial outgrowth of the pseudoperianth (black arrows) and calyptra (red arrows) post-fertilization is not affected in Mpknox1-6 mutants, but their development is arrested as well approximately 1 wpf (K, L). (M, N) Archegoniophores of Mpknox1-6 mutants begin to senesce approximately 2 weeks after maturation (N), while wild-type archegoniophores of the same age remain green (M). (O) Based on RNA-sequencing (RNA-seq) data and associated gene models, the 3' end of MpKNOX1 overlaps with the 3' untranslated region (UTR) of MpSUK1, which is transcribed from the opposite strand. Predicted full-length MpKNOX1 transcripts, that is, those including exons 2 and 3, are present primarily in the archegoniophore, consistent with the in situ data and semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) data in Figure 1—figure supplement 1. FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) scales are set to a standard within each tissue to allow potential weakly expressed transcripts to be visualized. (P) Signal of MpSUK1 in situ hybridization is detected in antheridia of the antheridiophore. (Q, R) Expression of MpKNOX1:GUS:SUK1 in an archegoniophore (R), an archegonium (S), and an antheridiophore (Q). Scale bar = 50 µm (A, E, G, J, K, L, S); 100 µm (B, C, D, H, I); 200 µm (P); 500 µm (F, R); 1000 µm (Q).
Mpknox1-6 mutants rarely form embryos on senescing archegoniophores.
(A–D) Cross-sections of Mpknox1-6 female mutants crossed with wild-type males and observed 2 weeks post fertilization (wpf). Note that some cells are deformed and collapsed, a characteristic of senescing tissue. (B) is the close up of (A). (E) Percentages of developmental stages observed from this cross. Plants were crossed once, and then examined after 1 or 2 wpf. Observed embryos were grouped into the following developmental stages: aborted zygotes, sporophytes that were arrested at a stage younger than the maximum stage possible, and sporophytes that had reached the maximum stage expected. n = total number of observed fertilization events. See Figure 3—figure supplement 1—source data 1 for details. Scale bar = 50 µm (B–D); 200 µm (A).
Expression profile of MpKNOX1 transcriptional β-glucuronidase reporter gene fusions unravels regulation by antisense transcript MpSUK1.
Sequences in purple are from the endogenous locus, while sequences in other colors are from the heterologous locus. (A) Endogenous MpKNOX1/MpSUK1 locus. (B) A construct, MpKNOX1:GUS, including a 4.6 kb sequence upstream (5') of the MpKNOX1 translational start site fused with the β-glucuronidase (GUS) reporter gene followed by a nopaline synthase (NOS) transcriptional terminator (250 bp), drives expression in both the egg cell (left) and antheridia (right). (C) Inclusion of the MpSUK1 locus, MpKNOX1:GUSterm:SUK1, including 1.5 kb upstream (5') of the MpSUK1 translational start site, downstream of the NOS terminator, results in the same expression pattern, that is, in both the egg cell and antheridia. (D) Only in constructs featuring an assembly in which MpKNOX1:GUS and MpSUK1:MpSUK1 are opposing and not interrupted by a transcriptional terminator (MpKNOX1:GUS:SUK1) is expression of MpKNOX1 restricted to the female egg cell (Figure 3) and the developing sporophyte, shown here at 1 week post fertilization (wpf), 2 wpf, and 4 wpf. No expression is observed in males (Figure 3). Scale bar = 50 µm (archegonia in B, C and 1-week-old sporophyte in D), 250 µm (2- and 4-week-old sporophytes in D), 500 µm (archegoniophores in B, C), and 1000 µm (antheridiophores in B, C and archegoniophores in D).
MpKNOX1A and MpKNOX1B lack the homeodomain and are both primarily expressed in the antheridiophores.
(A) Alignment of MpKNOX1 MEINOX domains with land plant KNOX1 proteins. A consensus sequence of the MEINOX domain (Bürglin, 1997; Crooks et al., 2004) of land plant KNOX1 proteins was computed from sequences in Figure 1—figure supplement 1, with the alignment of the three Marchantia polymorpha KNOX1 proteins displayed below. The lines above the sequence highlight the previously defined KNOX subdomains of the MEINOX domain. Gray shading indicates amino acids conserved with the land plant consensus and asterisks highlight amino acids conserved in all three M. polymorpha KNOX1 sequences. (B) Panels as described in Figure 3. Expression of both genes is predominantly detected in antheridiophores. (C–E) Translational fusion reporter line expression containing 1 kb (MpKNOX1A) and 2.3 kb (MpKNOX1B) of the putative 5' regulatory sequences suggests that antheridiophore expression identified by RNA-sequencing (RNAseq) is primarily in the antheridia. Expression of a translational fusion of MpKNOX1A with β-glucuronidase (GUS) was detectable in the sperm (D) and in mucilage papillae of the archegoniophore (E), but not in mucilage papillae in other parts of the thallus (not shown). Scale bar = 200 µm (antheridiophores and archegoniophores), 10 µm (sperm), and 50 µm (mucilage cells).
Alleles used in this study.
Nomenclature is as outlined previously (Bowman et al., 2016c; Montgomery et al., 2020).At the MpKNOX1 locus, a convergently transcribed gene, SUPPRESSOR OF KNOX1 (MpSUK1, Mp5g01590), is expressed predominantly in antheridia (Figure 3O–P), suggesting that MpKNOX1 may be regulated by an antisense transcript as has been described for MpFGMYB (Mp1g17210) (Hisanaga, 2019). Consistent with this hypothesis, a transcriptional fusion of 4.6 kb 5' of MpKNOX1 coding sequence with a β-glucuronidase (GUS) reporter coding sequence and 3' terminator results in expression in both egg cells and antheridia (Figure 3—figure supplement 2). In contrast, incorporation of 2.7 kb of the MpSUK1 antisense transcriptional unit immediately downstream of the GUS coding sequence (MpKNOX1:GUS:SUK1) results in female-specific expression (Figure 3; Figure 3—figure supplement 2). However, if a 3' terminator sequence is inserted between the GUS coding sequence and the MpSUK1 sequences, the antisense repression is lost, suggesting a requirement of transcriptional overlap but in a non-sequence-specific mechanism (Figure 3—figure supplement 2). MpSUK1 also encodes a 225 aa protein, including a 55 aa domain conserved across most liverworts but not found in other land plant lineages, that does not overlap with the MpKNOX1 coding sequence (Figure 3O).
Figure 3—figure supplement 2.
Expression profile of MpKNOX1 transcriptional β-glucuronidase reporter gene fusions unravels regulation by antisense transcript MpSUK1.
Sequences in purple are from the endogenous locus, while sequences in other colors are from the heterologous locus. (A) Endogenous MpKNOX1/MpSUK1 locus. (B) A construct, MpKNOX1:GUS, including a 4.6 kb sequence upstream (5') of the MpKNOX1 translational start site fused with the β-glucuronidase (GUS) reporter gene followed by a nopaline synthase (NOS) transcriptional terminator (250 bp), drives expression in both the egg cell (left) and antheridia (right). (C) Inclusion of the MpSUK1 locus, MpKNOX1:GUSterm:SUK1, including 1.5 kb upstream (5') of the MpSUK1 translational start site, downstream of the NOS terminator, results in the same expression pattern, that is, in both the egg cell and antheridia. (D) Only in constructs featuring an assembly in which MpKNOX1:GUS and MpSUK1:MpSUK1 are opposing and not interrupted by a transcriptional terminator (MpKNOX1:GUS:SUK1) is expression of MpKNOX1 restricted to the female egg cell (Figure 3) and the developing sporophyte, shown here at 1 week post fertilization (wpf), 2 wpf, and 4 wpf. No expression is observed in males (Figure 3). Scale bar = 50 µm (archegonia in B, C and 1-week-old sporophyte in D), 250 µm (2- and 4-week-old sporophytes in D), 500 µm (archegoniophores in B, C), and 1000 µm (antheridiophores in B, C and archegoniophores in D).
In addition to the full-length MpKNOX1 gene, two additional KNOX1-related genes, MpKNOX1A and MpKNOX1B, are encoded in the M. polymorpha genome. Both genes encode a conserved KNOX1 MEINOX domain (Figure 3—figure supplement 3A) but lack a HD and are predominantly expressed in the antheridia and sperm cells (Figure 3—figure supplement 3B-E).
MpBELL234 are expressed in antheridia, archegonia, and young embryos
The acropetal development of antheridia within antheridiophores results in the youngest being peripheral. MpBELL2, MpBELL3, and MpBELL4 expression was detected in stage 2–4 antheridia (Higo et al., 2016) within antheridiophores as assayed by in situ hybridization (Figure 4A, C and E). Likewise, translational MpBELL2 and transcriptional MpBELL3 and MpBELL4 GUS reporter gene fusion lines using 3.5 kb (MpBELL2), 5.6 kb (MpBELL3), and 5.8 kb (MpBELL4) of sequence 5’ to a transcriptional start site exhibit signals in antheridia of developing antheridiophores (Figure 4B, D, F and G). In these reporter lines, no signal was detected in unfertilized archegoniophores or developing sporophytes, nor was any expression observed in the vegetative gametophyte. However, in contrast to the above reporter lines, strand-specific transcriptomic analysis suggests that shorter transcripts, likely produced from an alternative promoter, are generated at the MpBELL3 and MpBELL4 loci in archegoniophores and sporophytes (Figure 4—figure supplement 1A, B). Indeed, expression of both genes was detected in egg cells of unfertilized archegonia, gametophytic tissues surrounding fertilized archegonia at 1 wpf, as well as sporophytes up to at least 3 wpf via in situ hybridization (Figure 4I; Figure 4—figure supplement 1C). An MpBELL4 reporter construct using a 4.2 kb sequence upstream of a CAGE (Cap Analysis of Gene Expression) signal (Montgomery et al., 2020) corresponding to the presumed alternative transcriptional start site of the archegonial/sporophytic transcript results in an egg-cell and post-fertilization gametophytic signal; however, no expression was observed in the older sporophytes (Figure 4G–H and J). Directed RT-PCR experiments also confirmed the existence of a shorter transcript encoding homeodomain sequences produced from MpBELL3 in archegoniophores and during sporophyte development (Figure 4—figure supplement 2).
Figure 4.
Expression patterns of MpBELL2, MpBELL3, and MpBELL4.
In situ localization of all three MpBELL mRNAs in antheridia as seen in antheridiophore cross-sections; younger antheridia are toward the left (A, C, E). Antheridiophores with translational (B) or transcriptional (D, F) β-glucuronidase reporter gene fusions of MpBELL2 (B), MpBELL3 (D), and MpBELL4 (F). These reporter genes all harbor sequences 5’ of the longest predicted transcript at each of the loci; for example, MpBELL4 in (G). Signal appears strongest in young- to medium-aged antheridia (stage 3 and stage 4; Higo et al., 2016), with the signal being lost in older antheridia toward the center of the antheridiophore (B, F), possibly due to draining of spermatogenous tissue. Staining was not observed in unfertilized archegonia (B, D, F). (I) In situ hybridization using the full-length MpBELL4 coding sequence exhibits signals in egg cells of unfertilized archegonia, gametophytic tissues surrounding fertilized archegonia at 1 week post fertilization (wpf), as well as sporophytes up to at least 3 wpf. (H, J) Reporter genes constructed using an alternative promoter internal of the MpBELL4 locus (MpBELL4) marked with a Cap Analysis of Gene Expression (CAGE) signal (G) exhibit a signal in unfertilized archegonia, but not in antheridia (H). Signal remains visible in developing unfertilized archegonia within the center of the archegoniophore, but no signal was observed in sporophytes at 2 wpf (J). Scale bar = 200 μm for left panels in A, C, E and the right panel in C, panels B, D, F, H, and the right panel in J; 100 μm for right panels in A, E, I and the left panel in J; 50 μm for the left and middle panel in I.
Panels as described in Figure 3. Vertical pink lines highlight exons encoding the homeodomains, each spanning three exons, in MpBELL3, MpBELL4, and MpBELL5. (A) Expression in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017). FPKM scales are set to a standard within each outlined box (if less than 10, rounded to.1). (B) Strand-specific expression at each MpBELL locus (marchantia.info). Potential alternative shorter transcripts (highlighted in red boxes) are detected for MpBELL3 and MpBELL4 in archegoniophores and sporophytes. Scales in the strand-specific expression panels vary according to tissue to allow potential weakly expressed transcripts to be visualized.
(A) Schematics of MpBELL3 transcripts predicted to be produced in each of the male gametophytes, female gametophytes, and sporophytes are depicted, showing Cap Analysis of Gene Expression (CAGE) signals of transcription start sites (TSS), intron/exon boundaries (from https://marchantia.info), position of homeodomains, and primers used in subsequent experiments. The (B) reverse transcription polymerase chain reaction (RT-PCR) using different primer combinations (as shown in A) suggests a continuous long transcript in antheridia. (C) RT-PCR indicates that transcripts including homeodomains HDb and HDc are only expressed in antheridiophores (yellow), while a shorter 3’ transcript (blue) including HDa is also expressed in archegoniophores. (D) MpBELL3 nested 5’ and 3’ RACE (Rapid Amplification of cDNA Ends) on complementary DNA (cDNA) derived from archegoniophores, archegoniophores harboring 1- + 2-week-old sporophytes, and antheridiophores. Red arrows show the expected amplicon length, while red asterisks indicate that a correct amplicon could be verified via direct sequencing. Both panels show that a truncated 3’ transcript indicated by a CAGE signal around the exon 9/10 boundary including homeodomain HDa is present in all three tissues. (E) Strand-specific expression profile of MpBELL3 in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017) as in Figure 4—figure supplement 1 (B), highlighting the positions of 5’ and 3’ RACE primers in relation to the transcript signal. The transcriptome signal suggests that some introns may experience alternative splicing. (F) In situ experiments using the truncated MpBELL3 3’ transcript spanning exons 11–16 (blue transcript in A) give signal in archegonia, archegonia harboring 1- and 3-week-old sporophytes as well as antheridia. Scale bar = 50 μm for the two left panels in (F) and 200 μm for the two right panels in (F).
(A) Schematics of MpBELL4 transcripts predicted to be produced in each of the male gametophytes, female gametophytes, and sporophytes are depicted, showing Cap Analysis of Gene Expression (CAGE) signals, intron/exon boundaries (from https://marchantia.info), position of homeodomains, and primers used in subsequent experiments. (B) Reverse transcription-polymerase chain reaction (RT-PCR) using different primer combinations (as shown in A) suggests that transcripts harboring HDa are only expressed in antheridia, whereas a truncated 3’ transcript is expressed in archegonia harboring sporophytes at 3 weeks post fertilization (wpf). Tissues from which RNA was extracted: 1, wild-type antheridiophores; 2, wild-type archegoniophores (not fertilized); 3, wild-type 1 wpf; 4, wild-type 3 wpf; 5, wild-type female x Mpbell34-4 (1 wpf); 6, wild-type female x Mpbell2/3/4-15 (1 wpf); 7, genomic DNA; 8, no template control. (C) MpBELL4 nested 5’ and 3’ RACE on complementary DNA (cDNA) derived from archegoniophores, archegoniophores harboring 1- + 2-week-old sporophytes, and antheridiophores. Red arrows show the expected amplicon length, while red asterisks indicate that a matching amplicon could be verified via direct sequencing. 5’ RACE and 3’ RACE on untranslated regions (UTRs) (panels 1 and 3) only give signal in antheridia, suggesting that a full-length transcript harboring the carboxy-terminal homeodomain (HDa) is only present in antheridia. A truncated 3’ transcript could be detected by 5’ RACE from exon 11 in all tissue examined (panel 2), although only antheridia-derived transcripts were spliced. (D) Strand-specific expression profile of MpBELL4 in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017) as used in Figure 4—figure supplement 1 (B), highlighting the positions of 5’ and 3’ RACE primers in relation to the transcript signal. The transcriptome signal suggests that some introns may experience alternative splicing.
Figure 4—figure supplement 1.
Expression profiles of MpBELL loci.
Panels as described in Figure 3. Vertical pink lines highlight exons encoding the homeodomains, each spanning three exons, in MpBELL3, MpBELL4, and MpBELL5. (A) Expression in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017). FPKM scales are set to a standard within each outlined box (if less than 10, rounded to.1). (B) Strand-specific expression at each MpBELL locus (marchantia.info). Potential alternative shorter transcripts (highlighted in red boxes) are detected for MpBELL3 and MpBELL4 in archegoniophores and sporophytes. Scales in the strand-specific expression panels vary according to tissue to allow potential weakly expressed transcripts to be visualized.
Figure 4—figure supplement 2.
MpBELL3 locus is a source of multiple transcripts in antheridia versus archegonia/sporophytes.
(A) Schematics of MpBELL3 transcripts predicted to be produced in each of the male gametophytes, female gametophytes, and sporophytes are depicted, showing Cap Analysis of Gene Expression (CAGE) signals of transcription start sites (TSS), intron/exon boundaries (from https://marchantia.info), position of homeodomains, and primers used in subsequent experiments. The (B) reverse transcription polymerase chain reaction (RT-PCR) using different primer combinations (as shown in A) suggests a continuous long transcript in antheridia. (C) RT-PCR indicates that transcripts including homeodomains HDb and HDc are only expressed in antheridiophores (yellow), while a shorter 3’ transcript (blue) including HDa is also expressed in archegoniophores. (D) MpBELL3 nested 5’ and 3’ RACE (Rapid Amplification of cDNA Ends) on complementary DNA (cDNA) derived from archegoniophores, archegoniophores harboring 1- + 2-week-old sporophytes, and antheridiophores. Red arrows show the expected amplicon length, while red asterisks indicate that a correct amplicon could be verified via direct sequencing. Both panels show that a truncated 3’ transcript indicated by a CAGE signal around the exon 9/10 boundary including homeodomain HDa is present in all three tissues. (E) Strand-specific expression profile of MpBELL3 in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017) as in Figure 4—figure supplement 1 (B), highlighting the positions of 5’ and 3’ RACE primers in relation to the transcript signal. The transcriptome signal suggests that some introns may experience alternative splicing. (F) In situ experiments using the truncated MpBELL3 3’ transcript spanning exons 11–16 (blue transcript in A) give signal in archegonia, archegonia harboring 1- and 3-week-old sporophytes as well as antheridia. Scale bar = 50 μm for the two left panels in (F) and 200 μm for the two right panels in (F).
Expression patterns of MpBELL2, MpBELL3, and MpBELL4.
In situ localization of all three MpBELL mRNAs in antheridia as seen in antheridiophore cross-sections; younger antheridia are toward the left (A, C, E). Antheridiophores with translational (B) or transcriptional (D, F) β-glucuronidase reporter gene fusions of MpBELL2 (B), MpBELL3 (D), and MpBELL4 (F). These reporter genes all harbor sequences 5’ of the longest predicted transcript at each of the loci; for example, MpBELL4 in (G). Signal appears strongest in young- to medium-aged antheridia (stage 3 and stage 4; Higo et al., 2016), with the signal being lost in older antheridia toward the center of the antheridiophore (B, F), possibly due to draining of spermatogenous tissue. Staining was not observed in unfertilized archegonia (B, D, F). (I) In situ hybridization using the full-length MpBELL4 coding sequence exhibits signals in egg cells of unfertilized archegonia, gametophytic tissues surrounding fertilized archegonia at 1 week post fertilization (wpf), as well as sporophytes up to at least 3 wpf. (H, J) Reporter genes constructed using an alternative promoter internal of the MpBELL4 locus (MpBELL4) marked with a Cap Analysis of Gene Expression (CAGE) signal (G) exhibit a signal in unfertilized archegonia, but not in antheridia (H). Signal remains visible in developing unfertilized archegonia within the center of the archegoniophore, but no signal was observed in sporophytes at 2 wpf (J). Scale bar = 200 μm for left panels in A, C, E and the right panel in C, panels B, D, F, H, and the right panel in J; 100 μm for right panels in A, E, I and the left panel in J; 50 μm for the left and middle panel in I.
Expression profiles of MpBELL loci.
Panels as described in Figure 3. Vertical pink lines highlight exons encoding the homeodomains, each spanning three exons, in MpBELL3, MpBELL4, and MpBELL5. (A) Expression in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017). FPKM scales are set to a standard within each outlined box (if less than 10, rounded to.1). (B) Strand-specific expression at each MpBELL locus (marchantia.info). Potential alternative shorter transcripts (highlighted in red boxes) are detected for MpBELL3 and MpBELL4 in archegoniophores and sporophytes. Scales in the strand-specific expression panels vary according to tissue to allow potential weakly expressed transcripts to be visualized.
MpBELL3 locus is a source of multiple transcripts in antheridia versus archegonia/sporophytes.
(A) Schematics of MpBELL3 transcripts predicted to be produced in each of the male gametophytes, female gametophytes, and sporophytes are depicted, showing Cap Analysis of Gene Expression (CAGE) signals of transcription start sites (TSS), intron/exon boundaries (from https://marchantia.info), position of homeodomains, and primers used in subsequent experiments. The (B) reverse transcription polymerase chain reaction (RT-PCR) using different primer combinations (as shown in A) suggests a continuous long transcript in antheridia. (C) RT-PCR indicates that transcripts including homeodomains HDb and HDc are only expressed in antheridiophores (yellow), while a shorter 3’ transcript (blue) including HDa is also expressed in archegoniophores. (D) MpBELL3 nested 5’ and 3’ RACE (Rapid Amplification of cDNA Ends) on complementary DNA (cDNA) derived from archegoniophores, archegoniophores harboring 1- + 2-week-old sporophytes, and antheridiophores. Red arrows show the expected amplicon length, while red asterisks indicate that a correct amplicon could be verified via direct sequencing. Both panels show that a truncated 3’ transcript indicated by a CAGE signal around the exon 9/10 boundary including homeodomain HDa is present in all three tissues. (E) Strand-specific expression profile of MpBELL3 in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017) as in Figure 4—figure supplement 1 (B), highlighting the positions of 5’ and 3’ RACE primers in relation to the transcript signal. The transcriptome signal suggests that some introns may experience alternative splicing. (F) In situ experiments using the truncated MpBELL3 3’ transcript spanning exons 11–16 (blue transcript in A) give signal in archegonia, archegonia harboring 1- and 3-week-old sporophytes as well as antheridia. Scale bar = 50 μm for the two left panels in (F) and 200 μm for the two right panels in (F).
MpBELL4 transcripts in archegonia/sporophytes lack homeodomains.
(A) Schematics of MpBELL4 transcripts predicted to be produced in each of the male gametophytes, female gametophytes, and sporophytes are depicted, showing Cap Analysis of Gene Expression (CAGE) signals, intron/exon boundaries (from https://marchantia.info), position of homeodomains, and primers used in subsequent experiments. (B) Reverse transcription-polymerase chain reaction (RT-PCR) using different primer combinations (as shown in A) suggests that transcripts harboring HDa are only expressed in antheridia, whereas a truncated 3’ transcript is expressed in archegonia harboring sporophytes at 3 weeks post fertilization (wpf). Tissues from which RNA was extracted: 1, wild-type antheridiophores; 2, wild-type archegoniophores (not fertilized); 3, wild-type 1 wpf; 4, wild-type 3 wpf; 5, wild-type female x Mpbell34-4 (1 wpf); 6, wild-type female x Mpbell2/3/4-15 (1 wpf); 7, genomic DNA; 8, no template control. (C) MpBELL4 nested 5’ and 3’ RACE on complementary DNA (cDNA) derived from archegoniophores, archegoniophores harboring 1- + 2-week-old sporophytes, and antheridiophores. Red arrows show the expected amplicon length, while red asterisks indicate that a matching amplicon could be verified via direct sequencing. 5’ RACE and 3’ RACE on untranslated regions (UTRs) (panels 1 and 3) only give signal in antheridia, suggesting that a full-length transcript harboring the carboxy-terminal homeodomain (HDa) is only present in antheridia. A truncated 3’ transcript could be detected by 5’ RACE from exon 11 in all tissue examined (panel 2), although only antheridia-derived transcripts were spliced. (D) Strand-specific expression profile of MpBELL4 in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017) as used in Figure 4—figure supplement 1 (B), highlighting the positions of 5’ and 3’ RACE primers in relation to the transcript signal. The transcriptome signal suggests that some introns may experience alternative splicing.
Paternal and maternal MpBELL are required for proper embryo development
To examine whether MpBELL234 could provide the male counterpart to the female MpKNOX1, we created loss-of-function alleles for each gene (Table 1). When used as male parents in crosses with wild-type females, Mpbell2, Mpbell3, and Mpbell4 single mutants did not exhibit an aberrant phenotype; however, when Mpbell3 Mpbell4 (Mpbell34) double mutant or Mpbell2 Mpbell3 Mpbell4 (Mpbell234) triple mutant males were used in crosses with wild-type females, mature sporophytes were formed at a reduced frequency compared to crosses between wild-type parents (Figure 5A–B). To quantify the reduction in fertility, we first examined crosses between wild-type parents in greater detail. A majority of fertilization events in crosses between wild-type males and wild-type females produced mature sporophytes (Figure 5A–B). However, nearly 40% of wild-type crosses were either aborted at the zygote stage or arrested at a globular multicellular stage, characteristic of about 1 wpf. As we did not observe sporophytes arrested at later stages, sporophytes that surpass the 1 wpf stage must generally progress to maturity. In contrast, most fertilization events derived from crosses between wild-type females and Mpbell(2)34 males resulted in sporophyte development that only progressed to a globular stage characteristic of the first week of wild-type development, followed by sporophyte arrest. A smaller number of fertilization events were aborted at zygote formation, resembling the phenotype observed in female Mpknox1-6 mutants (Figure 5A–B). However, mature sporophytes with viable spores were produced in about 17% of fertilization events. In individual crosses in which mature sporophytes form, the corresponding archegoniophores remain viable until after the sporophytes have matured, while unfertilized archegoniophores undergo senescence. The archegoniophores with only arrested sporophytes senesced in a manner similar to unfertilized archegoniophores.
Table 1.
Alleles used in this study.
Nomenclature is as outlined previously (Bowman et al., 2016c; Montgomery et al., 2020).
Sporophyte development entails BELL activity in both male and female gametes.
(A) Wild-type females crossed with wild-type males (left column), reciprocal crosses between wild-type and Mpbell mutant alleles (middle two columns), and crosses between Mpbell234 homozygotes (right column) observed after 1, 2, and 3 weeks post fertilization (wpf). All crosses produced aborted embryos or arrested sporophytes (left pictures in each panel). However, crosses between Mpbell234-18 and Mpbell234-83 only produced embryos aborted at the zygote stage or arrested at approximately the 1 wpf stage. See Table 1 for details on mutant alleles. (B) Percentages of developmental stages observed for each of the crosses. Plants were crossed once and then examined after 1, 2, or 3 wpf. Observed embryos were grouped into the following developmental stages: aborted zygotes, sporophytes that were arrested at a stage younger than the maximum stage possible, and sporophytes that had reached the maximum stage expected. N = the total number of observed fertilization events. See Figure 5—source data 1 for details. (C–F) BiFC (Bimolecular fluorescence complementation) assay of protein-protein interaction of MpKNOX1 and MpBELL4 in Nicotiana benthamiana leaves. Homodimers of MpKNOX1 show cytoplasmic localization (C), while MpBELL4 homodimers localize to the nucleus (D). Co-expression of both proteins resulted in nuclear localization of heterodimers (E, F). (G–J) BiFC assays of the protein derived from the shorter 3’ transcript of MpBELL3 (Figure 4—figure supplement 2) with MpKNOX1 show weak cytoplasmic interaction (G, H). However, MpBELL4 interacts strongly with shortMpBELL3 and results in nuclear localization (I, J). Scale bar = 50 µm, for all left pictures in (A) panels, all double mutant pictures, and all 1 wpf pictures; 100 µm for right pictures of the 2 wpf panel; 200 µm for all right pictures of the 3 wpf panel.
(a) Developmental staging data of crosses between wild-type females and wild-type males. (b) Developmental staging data of crosses between wild-type females and Mpbell234-83 males. (c) Developmental staging data of crosses between Mpbell234-18 females and wild-type males. (d) Developmental staging data of crosses between Mpbell234-18 females and Mpbell234-83 males. (e) Summary and calculations of developmental staging data.
BiFC assay of protein-protein interaction in Nicotiana benthamiana leaves. MpKNOX1 homodimers showed cytoplasmic localization (A), but no interaction could be observed when MpKNOX1 was co-expressed with MpBELL1 (B, C). MpBELL4 homodimers localized to the nucleus (D), but co-expression of both MpBELL4 and MpKNOX2 showed no interaction (E, F). MpKNOX2 homodimers alone showed cytoplasmic localization (G); however, when MpKNOX2 and MpBELL1 were co-expressed, the interaction signal was located to the nucleus (H, I). MpBELL1 did not homodimerize (J–L). Note that (A) and (D) are taken from Figure 5C and D for reference.
(A) Plants were doubly transformed using one construct in cis and one in trans: proEF1:MpKNOX1 constitutively expressing MpKNOX1 CDS and MpBELL3 expressing a 4.8 kb fragment of the putative 9 kb coding sequence after induction. Note that the pair of constructs features two different selection markers (nptII for selection on G418 and hptII for hygromycin selection) for selection of doube transformants. (B) Semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) of mock and 5 μM 17-β-estradiol-treated plants after 72 hr of induction. MpKNOX1 transcription was induced after 72 hr in MpBELL3 lines (red arrowhead). Antheridia-specific MpBELL3 expression was induced in MpBELL3 lines, in case of the double transformants also in the un-induced line (blue arrowheads). MpKNOX2 and MpBELL1, both normally expressed primarily in the sporophyte, were de-repressed after induction via MpBELL3 expression (green arrowheads). Th = thallus; arch = archegoniophore; 3 wk sporo = archegoniophore including sporophytes at 3 weeks post fertilization (wpf); anth = antheridiophores; gDNA = genomic DNA; nc = negative control.
Sporophyte development entails BELL activity in both male and female gametes.
(A) Wild-type females crossed with wild-type males (left column), reciprocal crosses between wild-type and Mpbell mutant alleles (middle two columns), and crosses between Mpbell234 homozygotes (right column) observed after 1, 2, and 3 weeks post fertilization (wpf). All crosses produced aborted embryos or arrested sporophytes (left pictures in each panel). However, crosses between Mpbell234-18 and Mpbell234-83 only produced embryos aborted at the zygote stage or arrested at approximately the 1 wpf stage. See Table 1 for details on mutant alleles. (B) Percentages of developmental stages observed for each of the crosses. Plants were crossed once and then examined after 1, 2, or 3 wpf. Observed embryos were grouped into the following developmental stages: aborted zygotes, sporophytes that were arrested at a stage younger than the maximum stage possible, and sporophytes that had reached the maximum stage expected. N = the total number of observed fertilization events. See Figure 5—source data 1 for details. (C–F) BiFC (Bimolecular fluorescence complementation) assay of protein-protein interaction of MpKNOX1 and MpBELL4 in Nicotiana benthamiana leaves. Homodimers of MpKNOX1 show cytoplasmic localization (C), while MpBELL4 homodimers localize to the nucleus (D). Co-expression of both proteins resulted in nuclear localization of heterodimers (E, F). (G–J) BiFC assays of the protein derived from the shorter 3’ transcript of MpBELL3 (Figure 4—figure supplement 2) with MpKNOX1 show weak cytoplasmic interaction (G, H). However, MpBELL4 interacts strongly with shortMpBELL3 and results in nuclear localization (I, J). Scale bar = 50 µm, for all left pictures in (A) panels, all double mutant pictures, and all 1 wpf pictures; 100 µm for right pictures of the 2 wpf panel; 200 µm for all right pictures of the 3 wpf panel.
Developmental staging data of crosses between wild-type and Mpbell234 mutants.
(a) Developmental staging data of crosses between wild-type females and wild-type males. (b) Developmental staging data of crosses between wild-type females and Mpbell234-83 males. (c) Developmental staging data of crosses between Mpbell234-18 females and wild-type males. (d) Developmental staging data of crosses between Mpbell234-18 females and Mpbell234-83 males. (e) Summary and calculations of developmental staging data.
MpKNOX1/MpBELL4 and MpKNOX2/MpBELL1 exhibit specificity in their interactions.
BiFC assay of protein-protein interaction in Nicotiana benthamiana leaves. MpKNOX1 homodimers showed cytoplasmic localization (A), but no interaction could be observed when MpKNOX1 was co-expressed with MpBELL1 (B, C). MpBELL4 homodimers localized to the nucleus (D), but co-expression of both MpBELL4 and MpKNOX2 showed no interaction (E, F). MpKNOX2 homodimers alone showed cytoplasmic localization (G); however, when MpKNOX2 and MpBELL1 were co-expressed, the interaction signal was located to the nucleus (H, I). MpBELL1 did not homodimerize (J–L). Note that (A) and (D) are taken from Figure 5C and D for reference.
Co-expression strategy to express MpKNOX1 and MpBELL3 simultaneously in the vegetative gametophyte.
(A) Plants were doubly transformed using one construct in cis and one in trans: proEF1:MpKNOX1 constitutively expressing MpKNOX1 CDS and MpBELL3 expressing a 4.8 kb fragment of the putative 9 kb coding sequence after induction. Note that the pair of constructs features two different selection markers (nptII for selection on G418 and hptII for hygromycin selection) for selection of doube transformants. (B) Semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) of mock and 5 μM 17-β-estradiol-treated plants after 72 hr of induction. MpKNOX1 transcription was induced after 72 hr in MpBELL3 lines (red arrowhead). Antheridia-specific MpBELL3 expression was induced in MpBELL3 lines, in case of the double transformants also in the un-induced line (blue arrowheads). MpKNOX2 and MpBELL1, both normally expressed primarily in the sporophyte, were de-repressed after induction via MpBELL3 expression (green arrowheads). Th = thallus; arch = archegoniophore; 3 wk sporo = archegoniophore including sporophytes at 3 weeks post fertilization (wpf); anth = antheridiophores; gDNA = genomic DNA; nc = negative control.Since both MpBELL3 and MpBELL4 expression was detected in the egg cell, we performed reciprocal crosses between Mpbell(2)34 females and wild-type males and we observed a similar distribution of sporophyte phenotypes compared to wild-type females crossed with Mpbell234 males: aborted zygotes, aborted multicellular embryos, and phenotypically normal sporophytes (Figure 5A–B). However, in crosses where both the male and female harbored mutations in both MpBELL3 and MpBELL4, no mature sporophytes were formed, with a majority of embryos arrested at the zygote stage (Figure 5A–B), a phenotype similar to, albeit less severe than, that observed for Mpknox1 mutants. Similar sporophyte phenotype ratios in reciprocal crosses between Mpknox1 and Mpbell34 have been described independently (Hisanaga, 2021).We next examined whether MpKNOX1 and MpBELL proteins could interact. We chose antheridial-expressed full-length MpBELL4 for analysis due to the truncated nature of MpBELL2 and the extreme length of MpBELL3. In a split Yellow Fluorescent Protein (YFP) BiFC assay in Nicotiana benthamia leaf epidermal cells, MpKNOX1 on its own was cytoplasmically localized, but when co-expressed with nuclear localized MpBELL4, MpKNOX1 signal became nuclear (Figure 5C–F). A similar interaction was observed with MpKNOX2 and MpBELL1 (Figure 5—figure supplement 1). These interactions are selective, since neither interaction between MpKNOX1 and MpBELL1 nor between MpKNOX2 and MpBELL4 was observed (Figure 5—figure supplement 1). In contrast to full-length MpBELL4, we did not observe any efficient interaction between the shorter sporophyte-expressed MpBELL3 and MpKNOX1 in this heterologous system (Figure 5G–H). Surprisingly, however, in this system, full-length MpBELL4 and sporophyte-expressed MpBELL3 interacted (Figure 5I–J).
Figure 5—figure supplement 1.
MpKNOX1/MpBELL4 and MpKNOX2/MpBELL1 exhibit specificity in their interactions.
BiFC assay of protein-protein interaction in Nicotiana benthamiana leaves. MpKNOX1 homodimers showed cytoplasmic localization (A), but no interaction could be observed when MpKNOX1 was co-expressed with MpBELL1 (B, C). MpBELL4 homodimers localized to the nucleus (D), but co-expression of both MpBELL4 and MpKNOX2 showed no interaction (E, F). MpKNOX2 homodimers alone showed cytoplasmic localization (G); however, when MpKNOX2 and MpBELL1 were co-expressed, the interaction signal was located to the nucleus (H, I). MpBELL1 did not homodimerize (J–L). Note that (A) and (D) are taken from Figure 5C and D for reference.
Expression of the fifth M. polymorpha MpBELL gene, MpBELL5, may be restricted to the archegoniophore among tissues examined (Figure 4—figure supplement 1A, B). As with MpKNOX1 (Figure 3) and MpFGMYB (Hisanaga, 2019), the MpBELL5 transcript may also be opposed by a convergently transcribed locus in antheridia (Figure 4—figure supplement 1A, B).
Activation of sporophyte gene expression in the vegetative gametophyte
The expression patterns and mutant phenotypes of MpBELL2/3/4 and MpKNOX1 genes are consistent with a role in activating diploid gene expression, and thus we examined whether ectopic (co-)expression of these genes in vegetative gametophyte is sufficient to activate diploid gene expression (Figure 5—figure supplement 2A). Ectopic expression of MpKNOX1 alone does not alter expression of MpKNOX2, MpBELL1, or MpBELL3 in the vegetative gametophyte (Figure 5—figure supplement 2B). However, ectopic expression of either MpBELL3 or co-expression of MpKNOX1 and MpBELL3 in the vegetative gametophyte for 72 hr is sufficient to activate both MpKNOX2 and MpBELL1, whose expression is normally limited to sporophyte development, or, in the case of MpBELL1, also after continuous far-red light induction (Inoue et al., 2019; Figure 5—figure supplement 2B). Thus, in this context, MpBELL3 alone can activate MpKNOX1 expression, a scenario perhaps reminiscent of post-zygotic activation of MpKNOX1 after fertilization (Figure 3C, Figure 3—figure supplement 2D).
Figure 5—figure supplement 2.
Co-expression strategy to express MpKNOX1 and MpBELL3 simultaneously in the vegetative gametophyte.
(A) Plants were doubly transformed using one construct in cis and one in trans: proEF1:MpKNOX1 constitutively expressing MpKNOX1 CDS and MpBELL3 expressing a 4.8 kb fragment of the putative 9 kb coding sequence after induction. Note that the pair of constructs features two different selection markers (nptII for selection on G418 and hptII for hygromycin selection) for selection of doube transformants. (B) Semiquantitative reverse transcription-polymerase chain reaction (sqRT-PCR) of mock and 5 μM 17-β-estradiol-treated plants after 72 hr of induction. MpKNOX1 transcription was induced after 72 hr in MpBELL3 lines (red arrowhead). Antheridia-specific MpBELL3 expression was induced in MpBELL3 lines, in case of the double transformants also in the un-induced line (blue arrowheads). MpKNOX2 and MpBELL1, both normally expressed primarily in the sporophyte, were de-repressed after induction via MpBELL3 expression (green arrowheads). Th = thallus; arch = archegoniophore; 3 wk sporo = archegoniophore including sporophytes at 3 weeks post fertilization (wpf); anth = antheridiophores; gDNA = genomic DNA; nc = negative control.
Discussion
An ancestral function for TALE-HD genes in Viridiplantae
The eukaryotic life cycle alternates between haploid and diploid phases, initiated by meiosis and gamete fusion, respectively. Organisms spanning the phylogenetic diversity of eukaryotes, including ascomycete and basiomycete fungi, Amoebozoa, brown algae, and chlorophyte algae, utilize paralogous homeodomain proteins that heterodimerize following gamete fusion to initiate the diploid genetic program, lending support for the idea that this may have been an ancestral function of homeodomain proteins. Our observations in M. polymorpha indicate that TALE-HD proteins, specifically MpKNOX1 and MpBELL2/3/4, are initially supplied in gametes, and that this gametic expression is required for diploid sporophyte development (Figure 6). MpKNOX1 is absolutely required in the egg cell, as is MpBELL3/4 in either the sperm or egg. This scenario is reminiscent of that characterized in Chlamydomonas, wherein KNOX is supplied by the (-) gamete and BELL via the (+) gamete, and that once together in the zygote, the diploid genetic program can be activated. Thus, the basic tenets of the genetic regulation of the haploid-to-diploid transition are conserved in two widely divergent Viridiplantae lineages, consistent with the notion that such a system was present in the common ancestor of extant Viridiplantae.
Figure 6.
Model of KNOX and BELL function in Marchantia polymorpha.
Three MpBELL genes (MpBELL2/3/4) are expressed in antheridia in the cells that will differentiate into sperm cells. MpKNOX1 is expressed in the egg cell of the archegonia, along with three MpBELL genes (MpBELL5 and MpBELL3/4). In contrast to the antheridia, where full-length transcripts are found, alternative short transcripts of MpBELL3/4 (MpBELL3s and MpBELL4s; see Figure 4—figure supplements 2 and 3) are produced in the egg cell. After fertilization, MpBELL3/4 derived from both sperm and egg cells and MpKNOX1 derived from the egg cell act together to activate the zygotic program. MpBELL1 and MpKNOX2 form a distinct heterodimer that acts at later stages of sporophyte development. The homeodomain lacking MpKNOX1A and MpKNOX1B genes expressed in the antheridia along with MpBELL2/3/4 potentially prevents functionality by sequestering MpBELL2/3/4 proteins. The function of MpBELL5 is unknown. Images are obtained from Marchant, 1713; Mirbel, 1835; Thuret, 1851; and Unger, 1837.
Unrooted phylogenetic tree adapted from Adl et al., 2012 and Keeling, 2014. Triangles at the end of branches represent monophyletic groups; triangles extending to the node represent paraphyletic groups. Taxa for which functional data are available are listed, as is the mode of homeodomain (HD) involvement if known. While the ancestral eukaryote possessed both TALE- and non-TALE-HD genes, the predicted composition of the heterodimer in the ancestral eukaryote is equivocal. T = TALE-HD; nt = non-TALE-HD; dv = divergent HD.
Model of KNOX and BELL function in Marchantia polymorpha.
Three MpBELL genes (MpBELL2/3/4) are expressed in antheridia in the cells that will differentiate into sperm cells. MpKNOX1 is expressed in the egg cell of the archegonia, along with three MpBELL genes (MpBELL5 and MpBELL3/4). In contrast to the antheridia, where full-length transcripts are found, alternative short transcripts of MpBELL3/4 (MpBELL3s and MpBELL4s; see Figure 4—figure supplements 2 and 3) are produced in the egg cell. After fertilization, MpBELL3/4 derived from both sperm and egg cells and MpKNOX1 derived from the egg cell act together to activate the zygotic program. MpBELL1 and MpKNOX2 form a distinct heterodimer that acts at later stages of sporophyte development. The homeodomain lacking MpKNOX1A and MpKNOX1B genes expressed in the antheridia along with MpBELL2/3/4 potentially prevents functionality by sequestering MpBELL2/3/4 proteins. The function of MpBELL5 is unknown. Images are obtained from Marchant, 1713; Mirbel, 1835; Thuret, 1851; and Unger, 1837.
Figure 4—figure supplement 3.
MpBELL4 transcripts in archegonia/sporophytes lack homeodomains.
(A) Schematics of MpBELL4 transcripts predicted to be produced in each of the male gametophytes, female gametophytes, and sporophytes are depicted, showing Cap Analysis of Gene Expression (CAGE) signals, intron/exon boundaries (from https://marchantia.info), position of homeodomains, and primers used in subsequent experiments. (B) Reverse transcription-polymerase chain reaction (RT-PCR) using different primer combinations (as shown in A) suggests that transcripts harboring HDa are only expressed in antheridia, whereas a truncated 3’ transcript is expressed in archegonia harboring sporophytes at 3 weeks post fertilization (wpf). Tissues from which RNA was extracted: 1, wild-type antheridiophores; 2, wild-type archegoniophores (not fertilized); 3, wild-type 1 wpf; 4, wild-type 3 wpf; 5, wild-type female x Mpbell34-4 (1 wpf); 6, wild-type female x Mpbell2/3/4-15 (1 wpf); 7, genomic DNA; 8, no template control. (C) MpBELL4 nested 5’ and 3’ RACE on complementary DNA (cDNA) derived from archegoniophores, archegoniophores harboring 1- + 2-week-old sporophytes, and antheridiophores. Red arrows show the expected amplicon length, while red asterisks indicate that a matching amplicon could be verified via direct sequencing. 5’ RACE and 3’ RACE on untranslated regions (UTRs) (panels 1 and 3) only give signal in antheridia, suggesting that a full-length transcript harboring the carboxy-terminal homeodomain (HDa) is only present in antheridia. A truncated 3’ transcript could be detected by 5’ RACE from exon 11 in all tissue examined (panel 2), although only antheridia-derived transcripts were spliced. (D) Strand-specific expression profile of MpBELL4 in different tissues of wild-type Marchantia polymorpha based on RNA-sequencing (RNA-seq) experiments (Bowman et al., 2017) as used in Figure 4—figure supplement 1 (B), highlighting the positions of 5’ and 3’ RACE primers in relation to the transcript signal. The transcriptome signal suggests that some introns may experience alternative splicing.
Phylogenetic perspective on HD gene-mediated haploid-to-diploid transitions in eukaryotes.
Unrooted phylogenetic tree adapted from Adl et al., 2012 and Keeling, 2014. Triangles at the end of branches represent monophyletic groups; triangles extending to the node represent paraphyletic groups. Taxa for which functional data are available are listed, as is the mode of homeodomain (HD) involvement if known. While the ancestral eukaryote possessed both TALE- and non-TALE-HD genes, the predicted composition of the heterodimer in the ancestral eukaryote is equivocal. T = TALE-HD; nt = non-TALE-HD; dv = divergent HD.The obvious major difference between the roles of KNOX and BELL in Marchantia when compared to Chlamydomonas is the activity of MpBELL34 in both the egg and sperm, rather than being confined to the sperm. Thus, questions on how the system operates mechanistically and why it might have evolved are paramount. The functions of maternal MpKNOX1 and paternal MpBELL34 could be presumed to represent the ancestral condition, with the role of maternal MpBELL34 being a derived condition. The absolute requirement for maternal MpKNOX1 (Figure 3) and its ability to heterodimerize with full-length MpBELL3/4 (Figure 5), which is expressed only paternally (Figure 4), suggest that this ancestral system functions in Marchantia, albeit with a reduced efficacy. In Marchantia, although neither maternal nor paternal MpBELL is absolutely required, both are necessary for efficient sporophyte development. One plausible hypothesis is that maternally supplied MpBELL3/4 evolved as a backup system to ensure diploid development following fertilization, concomitant with the evolution of sperm with reduced cytoplasmic content and highly condensed nuclei, anatomical features that could reduce the efficiency of MpBELL protein delivery by sperm. That both MpKNOX1 and shortMpBELL3 are present in the egg cell prior to fertilization and that they do not efficiently interact, at least in a split YFP BiFC assay in Nicotiana leaf cells (Figure 5), suggests that the act of fertilization may trigger a biochemical change required for their interaction. Although shortMpBELL3 can interact with full-length MpBELL4 (Figure 5), it is likely not critical in planta, given the absolute requirement of MpKNOX1 for sporophyte development. The requirement of both maternal and paternal MpBELL for efficient sporophyte production may also provide a mechanism to ensure activation of sporophyte development only occurs following fertilization. Finally, it is of note that a similar system, wherein two TALE-HD genes are both expressed in male and female gametes, has been described for the brown alga Ectocarpus (Arun et al., 2019). As Ectocarpus also possesses multicellular haploid and diploid generations and is anisogamous, perhaps similar evolutionary pressures contributed to the independent evolution of similar systems.Given the phylogenetic affinity of mosses as a sister group to liverworts (Renzaglia et al., 2018), we might expect a similar KNOX/BELL program to regulate the haploid-to-diploid transition in Physcomitrium. The phenotype of Mpknox1 alleles is more extreme than triple-null loss-of-function alleles of all PpKNOX1 loci where a sporophyte with viable spores forms (Sakakibara et al., 2008), suggesting that PpKNOX2 may also play a role in the transition in this species (Sakakibara et al., 2013). In Physcomitrium, PpBELL1 expression has been observed in sperm (Ortiz-Ramírez et al., 2017) via transcriptome analysis and in the egg via reporter gene (Horst et al., 2016). While the latter study reported a lack of male expression, this could be due to a loss of expression in this Physcomitrium accession (Meyberg et al., 2020), or alternatively, if the Physcomitrium PpBELL1 genomic architecture is similar to MpBELL3/4, the reporter gene insertion may have excluded regulatory sequences required for alternative transcription start sites. Hence, aspects of the ancestral system regulating the haploid-to-diploid transition could be present in multiple early diverging land plant lineages.In contrast to the situation in bryophytes, in angiosperms, paternally supplied pluripotency factors related to BABYBOOM, an AP2/ERF transcription factor, can activate the zygotic genetic program (Khanday et al., 2019; Conner et al., 2017; Conner et al., 2015). Furthermore, despite the extensive literature on KNOX/BELL function in angiosperms, the only potential known vestige of this ancestral system existing in angiosperms is the presence of KNOX2 in the Arabidopsis female gametophyte, possibly in the egg cell (Pagnussat et al., 2007). Thus, it seems that the ancestral function has been lost entirely in some derived land plants, as it has presumably in the metazoan lineage (Bowman et al., 2016b).Since liverworts are placental organisms, with the diploid sporophyte nourished by the maternal haploid gametophyte, there is opportunity for extensive intergenerational communication. For example, that the maternal gametophytic calyptra and pseudoperianth initiate their development in Mpknox1 mutants implicates a MpKNOX1-independent non-cell autonomous zygote-derived signal activated post-fertilization. In contrast, the premature senescence of Mpknox1 archegoniophores suggests a MpKNOX1-dependent intra-gametophytic signal for the maintenance of archegonial viability. Further, the premature senescence of wild-type archegoniophores bearing only aborted embryos suggests a signal emanating from older sporophytes to maintain the viability of the maternal archegoniophore tissues until sporophyte development is completed. Finally, that aborted embryos are found even in crosses involving only wild-type parents suggests possible maternal control over allocation of resources across multiple embryos, a phenomenon similar to maternal control over fruit set in angiosperms (Stephenson, 1981).
Roles of MpKNOX1 in Marchantia sporophyte
A key land plant innovation was the evolution of a multicellular diploid generation, the embryo, via mitoses, interpolated between gamete fusion and meiosis (Bower, 1908). In seed plants, ferns, and mosses, KNOX1 activity is associated with continued sporophyte cell proliferation, including sporophyte apical meristem activity (Hay and Tsiantis, 2010; Sakakibara et al., 2008; Sano et al., 2005). While M. polymorpha and other liverworts lack apical meristems during sporophyte development (Kienitz-gerloff, 1874), MpKNOX1 expression is detected throughout developing sporophytes until the inception of sporogenous cell differentiation (Figure 3), consistent with MpKNOX1 having a role in cell proliferation during M. polymorpha sporophyte development. Thus, while MpKNOX1 retains the ancestral function in the haploid-to-diploid transition, the lack of zygotic cell division in Mpknox1 mutants and the MpKNOX1 sporophytic expression pattern suggest that neofunctionalization of KNOX1 contributed to the evolution of the embryo via stimulation of cell proliferation. That MpBELL3/4 are also expressed during these stages of sporophyte development (Figure 4; Figure 4—figure supplements 1 and 2, 3) suggests these genes encode the heterodimeric partners of MpKNOX1 during sporophyte development. This is supported by the subfunctionalization of MpBELL genes with respect to heterodimerization partners—MpBELL3/4 specifically heterodimerize with MpKNOX1 while MpBELL1 interacts only with MpKNOX2 (Figure 5—figure supplement 2). This subfunctionalization resembles that observed for KNOX/BELL interactions in angiosperms (Furumizu et al., 2015) but represents an independent evolutionary event since all angiosperm BELL genes are MpBELL1 orthologs.In eukaryotic lineages in which multicellularity evolved in the diploid phase of the life cycle, genes involved in the haploid-to-diploid transition, or paralogs thereof, have been repeatedly co-opted into roles directing development in the diploid generation. Most conspicuously, despite the loss of their role in zygotic gene activation (see below), both non-TALE- (eg, Hox) and TALE-HD genes act to pattern the metazoan body (Merabet and Mann, 2016; Pearson et al., 2005; Lewis, 1978). Likewise, in the basidiomycete fungus Coprinopsis, the same HD heterodimer that initiates the diploid genetic program also directs early developmental stages in the multicellular diploid (Kamada, 2002). In the syncytial chlorophyte alga Caulerpa lentillifera, differential TALE-HD expression in the diploid body was speculated to influence differentiation of fronds and stolons, but functional data are lacking (Arimoto, 2019). Finally, in land plants, as the multicellular sporophyte evolved increasing complexity, KNOX/BELL genetic modules were co-opted to direct development of novel organs and tissues via regulation of cell proliferation (KNOX1) or differentiation (KNOX2) in conjunction with other gene regulatory networks (Furumizu et al., 2015; Hay and Tsiantis, 2010).
The remaining TALE-HD genes of M. polymorpha
The two KNOX1-related sequences lacking homeodomains (MpKNOX1A and MpKNOX1B) are reminiscent of similar, albeit independently evolved, proteins in angiosperms that act as inhibitors of KNOX/BELL function by forming inactive heterodimers primarily with BELL partners (Magnani and Hake, 2008; Kimura et al., 2008). The antheridial expression of MpKNOX1A and MpKNOX1B prompts the hypothesis that the encoded proteins may act as a safeguard against the inappropriate activity of MpBELL in the male gametophyte. This is particularly pertinent since ectopic expression of MpBELL3 alone in the vegetative gametophyte can activate MpKNOX1, MpKNOX2, and MpBELL1 (Figure 5—figure supplement 2). The lack of an aberrant phenotype observed in lines constitutively expressing MpKNOX1 in male gametophytes would be consistent with this hypothesis.In contrast to the other MpBELL genes, MpBELL5 transcripts are only detected in archegoniophores. Given the strong heterodimerization affinity preferences of BELL proteins for KNOX proteins (Bellaoui et al., 2001; Joo et al., 2018), the only conspicuous potential partner in M. polymorpha would be MpKNOX1. Hence, MpBELL5 might play an as-yet-undefined role in the activation of zygotic gene expression. At the MpBELL5 locus, a convergently transcribed locus is expressed predominantly in antheridia (Figure 4—figure supplement 1), suggesting MpBELL5 may be repressed by an antisense transcript in a manner similar to MpKNOX1 (Figure 3) and MpFGMYB (Hisanaga, 2019). Such male-expressed antisense transcript-mediated repression apparently provides a general mechanism for female-specific expression of autosomal genes, with sex-specific regulation ultimately being linked to the feminizing locus on the female sex chromosome (Bowman, 2016a; Knapp, 1935; Lorbeer, 1936).
Marchantia BELL gene diversity resembles that of charophyte algae
The phylogenetic diversity of BELL genes in M. polymorpha, and liverworts more broadly, more closely resembles the diversity observed in charophyte algae than that of other land plants (Joo et al., 2018; Lee et al., 2008). All previously described land plant BELL genes form a single clade that evolved within the charophycean algae (Figure 2). The single M. polymorpha gene in this clade, MpBELL1, is the only MpBELL gene predominantly expressed during sporophyte development. The other MpBELL genes reside in two other phylogenetically distinct clades, GLX-basal and BELL-gamete, both of which include charophyte sequences (Figure 2). This phylogenetic pattern suggests an ancient proliferation of BELL paralogs within the charophycean algal ancestor, with M. polymorpha (and other liverworts and some mosses) retaining this diversity that was subsequently lost in most other land plant lineages. The structural diversity (two or three homeodomains) of the MpBELL2/3/4 paralogs and the extensive sequence diversity of BELL-gamete clade genes are consistent with rapid evolution of genes involved in reproductive processes (Swanson and Vacquier, 2002). Finally, the expression pattern of MpBELL5, a member of the GLX-basal clade that evolved early in the charophytes, also suggests an as-yet-unresolved function in reproduction. Similarly, the Chlamydomonas genome encodes a second BELL paralog, HDG1, that is expressed in both (+) and (-) gametes, suggesting a role in reproduction, but whose function is unknown.
An ancestral function for TALE-HD genes in eukaryotes
In Viridiplantae (Chlamydomonas and Marchantia), the haploid-to-diploid transition is mediated by two TALE-HD genes. In the red alga Pyropia yezoensis, KNOX gene expression is detected in the diploid conchosporangium, but not in haploid thalli (Mikami et al., 2019). The phylogenetic distribution of KNOX and BELL subfamilies and their heterodimerization affinities (Joo et al., 2018) suggest that the Archaeplastida common ancestor utilized KNOX/BELL TALE-HD genes to mediate the haploid-to-diploid transition. Likewise, in the brown alga Ectocarpus, two TALE-HD proteins, OUROBOROS and SAMSARA, mediate the transition (Arun et al., 2019). In contrast, in both ascomycete and basidiomycete fungi, the haploid-to-diploid transition is mediated by heterodimerization of a TALE-HD and a non-TALE-HD protein (Gillissen et al., 1992; Goutte and Johnson, 1988; Herskowitz, 1989; Hull et al., 2005; Kues et al., 1992; Spit et al., 1998; Urban, 1996). In the Amoebozoa Dictyostelium, the two homeodomain-like proteins controlling the haploid-to-diploid transition are highly divergent, rendering the phylogenetic affinities enigmatic (Hedgethorne et al., 2017). As these taxa span eukaryotic phylogenetic diversity, one role of homeodomain genes in the ancestral eukaryote was to regulate the haploid-to-diploid transition—the evolution of the homeodomain in the ancestral eukaryote was associated with evolution of a novel life cycle (Figure 6—figure supplement 1). In arguably the most intensively studied taxon, the Metazoa, zygotic gene activation has been replaced by maternally derived pluripotency factors (Schulz and Harrison, 2019), but homeodomain genes have been retained.
Figure 6—figure supplement 1.
Phylogenetic perspective on HD gene-mediated haploid-to-diploid transitions in eukaryotes.
Unrooted phylogenetic tree adapted from Adl et al., 2012 and Keeling, 2014. Triangles at the end of branches represent monophyletic groups; triangles extending to the node represent paraphyletic groups. Taxa for which functional data are available are listed, as is the mode of homeodomain (HD) involvement if known. While the ancestral eukaryote possessed both TALE- and non-TALE-HD genes, the predicted composition of the heterodimer in the ancestral eukaryote is equivocal. T = TALE-HD; nt = non-TALE-HD; dv = divergent HD.
Given the ancestral eukaryote possessed a minimum of two HD genes, including one TALE-HD and one non-TALE-HD (Bharathan et al., 1997; Bürglin, 1997; Derelle et al., 2007), it is an intriguing question as to whether, ancestrally, both TALE and non-TALE paralogs acted in the diploid-to-haploid transition. Either possible ancestral scenario (TALE + non-TALE or TALE + TALE) requires either a co-option of a pre-existing paralog (eg, non-TALE), or alternatively, a TALE gene duplication and transference of function to the new paralog, respectively. To resolve the ancestral condition, broader phylogenetic sampling and functional analyses across additional unicellular Bikont lineages, particularly in the Excavata and the paraphyletic sister groups of the Metazoa (choanoflagellates, Filasterea, and Ichthyosporea), might be informative (Figure 6—figure supplement 1). Resolution of the ancestral condition could inform whether the primary ancestral function of both TALE-HD and non-TALE-HD was regulating the haploid-to-diploid transition or whether non-TALE-HD genes had another fundamental function in ancestral eukaryotes.
Materials and methods
Ecotypes, plant growth, transformation, and induction
M. polymorpha ssp ruderalis, ecotype BoGa, obtained from the Botanical Garden of Osnabrueck, Germany, was used in SZ lab. Plant cultivation, transformation, and induction were performed according to Ishizaki et al., 2008 and Althoff, 2014. M. polymorpha ssp ruderalis, ecotype MEL, was used in JLB lab and grown, transformed, and induced according to Ishizaki et al., 2008 and Flores-Sandoval et al., 2015. In case of double transformations, sporelings were co-transformed using two constructs featuring different selectable markers. Spores were drained and plated on selection media for approximately 2 weeks and then subjected to a second round of selection before being transferred to ½ B5 plates. For induction of reproductive organs, plants were transferred to white light supplemented with far red light (735 nm; 45 µmol/m2/s) on ½ B5 media supplemented with 1% glucose.
DNA was extracted using a modified protocol from Edwards et al., 1991. Instead of vacuum drying, the pelleted DNA was air-dried. Amplicons were directly sequenced. RNA was extracted from ≈100 mg wet weight tissue using Rneasy kit from Qiagen following the manufacturer’s instructions, including an on-column Dnase I digestion. The same amount of total RNA was used for complementary DNA (cDNA) synthesis using either Superscript II (Invitrogen) (SZ lab) or Bioscript (Bioline) (JLB lab) reverse transcriptase according to the manufacturer’s instructions with oligo-dT15 primer. 5’ and 3’ RACE-ready cDNA synthesis was performed using SMARTScribe reverse transcriptase (Takara) according to the manufacturer’s instructions but with the addition of random primer mix (NEB) to allow reverse transcription of long mRNA templates. See Supplementary file 1, Table 1a for primer sequences.
In situ hybridization, sectioning, staining, and microscopy
M. polymorpha ssp ruderalis, ecotype BoGa, tissue fixation, embedding, sectioning, and hybridization with digoxigenin (DIG)-labeled antisense RNA probes were performed according to Zachgo, 2002. Sections were either stained in toluidine blue O (TBO) alone or counterstained with ruthenium red (RR) according to Retamales and Scharaschkin, 2014. Staining time for RR and TBO was 90 s and 60 s, respectively. Microscopic slides were observed using a Leica MZ16 FA microscope and pictures were taken with a Leica DFC490 camera. Plants were observed using either a Lumar dissecting microscope (Zeiss) and photographed with AxioCam HRc and AxioVision software (both Zeiss) or a Leica stereomicroscope (Leica M165 FC) and photographed with an integrated Leica DFC490 camera. See Supplementary file 1, Table 1b for probe sequences.
Bimolecular fluorescence complementation
Coding sequences of KNOX and BELL genes were seamlessly cloned into the pGREEN derivatives (Hellens et al., 2000) pAMON/pSUR for N-terminal split VENUS (I152L) fusion and into pURIL/pDOX for C-terminal split VENUS (I152L) fusion (Lampugnani et al., 2016). All constructs feature the 35 S promoter driving the N- (pAMON/pURIL) or C-terminal (pSUR/pDOX) part of the split VENUS (I152L) fluorophore for translational fusion of the gene-of-interest in frame with the split VENUS (I152L). All constructs were transformed into Agrobacterium tumefaciens strain GV3101 for transfection of plant leaves (Bracha-Drori et al., 2004). For fluorescence microscopy, usually three, approximately 50 × 50 mm, abaxial leaf tissue fragments were examined under an Axioskop two mot plus (Zeiss) microscope and photographed using AxioCam HRc and AxioVision software. For YFP filter set 46, excitation BP 500/20; beam splitter FT 515; emission BP 535/30 was used. See Supplementary file 1, Table 1c for primer sequences.
Overexpression and transcriptional and translational fusion constructs
Primers used for overexpression with the endogenous, constitutive MpEF1α promoter (Althoff, 2014), inducible overexpression, and transcriptional and translational fusion constructs are shown in Supplementary file 1, Table 1c. The complete 1.8 kb MpKNOX1 coding sequence and a 4.8 kb fragment of the MpBELL3 coding sequence were amplified from cDNA and seamlessly cloned (NEBuilder, NEB) into pENTR-D (Invitrogen) via NotI/AscI sites. The MpKNOX1 coding sequence was subsequently recombined into pMpGWB403 (Addgene entry #68668) for constitutive expression, whereas the MpBELL3 coding sequence was recombined into the estrogen-inducible binary pHART XVE (Flores-Sandoval et al., 2016). Regulatory sequences upstream of the transcriptional start site of MpKNOX1 (4.6 kb), MpBELL3 (5.6 kb), and MpBELL4 (5.8 kb) and internal of the MpBELL4 locus (4.1 kb) were amplified and cloned into pENTR-D and subsequently recombined using LR-clonase II (Invitrogen) into pMpGWB104 (Addgene entry #68558) featuring the GUS reporter gene. To elucidate the transcriptional regulation at the MpKNOX1/MpSUK1 locus, the 4.6 kb regulatory region of MpKNOX1 was seamlessly cloned (NEBuilder, NEB) into the HindIII/XbaI site 5’ of the Gateway cassette of pMpGWB401 (Addgene entry #68666). The reversely transcribed MpSUK1 locus (2.7 kb) was seamlessly cloned into either the SacI site or AscI site (5’ or 3’ of the NOS terminator). Subsequently, the GUS reporter gene was recombined using pENTR-GUS (Invitrogen). For translational fusions of MpBELL2, MpKNOX1A, and MpKNOX1B, the regulatory sequence 5’ upstream of the transcriptional start site including the genomic gene locus and excluding the stop codon was amplified and cloned into pRITA, a shuttle vector featuring the GUS reporter gene and NOS terminator. MpBELL2 and MpKNOX1A were cloned into the KpnI or SalI/KpnI site, respectively. MpKNOX1B was seamlessly cloned using Infusion cloning (Clontech). All constructs were subsequently transferred to the binary vector pHART using NotI.
CRISPR/Cas9-mediated gene knock-out
Potential guide RNAs (gRNAs) were obtained and screened for off targets in the M. polymorpha genome using CRISPOR website (http://crispor.tefor.net/). Quadruple gRNAs were each cloned into pMpGE_En04, pBC-14, pBC-23, or pBC-34 using BsaI sites. Double gRNAs were cloned into pMpGE_En04 and pBC-14 (Hisanaga, 2019). These guides were subsequently assembled into the pMpGE_En04 backbone using BglI sites. pMpGE_En04 were gateway cloned into pMpGE010 (Sugano et al., 2018; Sugano et al., 2014; Hisanaga, 2019; Addgene entry #71,536), featuring Cas9 nuclease. In case of six gRNAs, the remaining two were cloned into pMpGWB401 (Ishizaki et al., 2015; Addgene entry #68,666) and plants doubly transformed (Ishizaki et al., 2008). See Supplementary file 1, Table 1d for gRNA sequences.
GUS staining
Generally, detection of GUS activity in a minimum of three independent lines was performed overnight at 37°C in GUS staining solution (0.05 M NaPO4, pH 7.2, 2 mM K3[Fe(CN)6], 2 mM K4[Fe(CN)6], 0.3% (v/v) Triton X-100) supplemented with 0.6 mg/ml b-D-glucopyranosiduronic acid (X-gluc) dissolved in N,N-dimethylformamide (DMF). Only detection of MpKNOX1:GUS activity in archegoniophores was performed for 4 hr in staining solution containing 5 mM K3[Fe(CN)6] and 5 mM K4[Fe(CN)6]. Antheridiophores were cleared in 70% (w/v) chloral hydrate in 10% (v/v) glycerol solution.
Phylogenetics
Predicted TALE-HD sequences were assembled from land plants, charophytes, and chlorophytes, via Genbank or additional sources as detailed below. Gymnosperm sequences were obtained from Congenie (congenie.org), fern sequences from Equisetum and Azolla transcriptomes (de Vries et al., 2016; Vanneste et al., 2015), bryophyte sequences, other than M. polymorpha (Phytozome) and P. patens (Genbank), from available transcriptomes (Wickett et al., 2014; Dong et al., 2019; Radhakrishnan et al., 2020), Klebsormidium and Caulerpa sequences derived from genome sequences (Hori et al., 2014; Arimoto, 2019), and other charophyte and chlorophyte sequences from Genbank and additional transcriptomes (Cooper and Delwiche, 2016; Ju et al., 2015).Complete or partial coding nucleotide sequences were manually aligned as amino acid translations using Se-Al v2.0a11 for Macintosh (http://tree.bio.ed.ac.uk/software/seal/). We excluded ambiguously aligned sequences to produce an alignment of 231 nucleotides (77 amino acids) for 132 BELL sequences. Alignments of KNOX genes included the homeodomain, MEINOX (KNOX) and ELK domains (Joo et al., 2018), comprising 630 nucleotides (210 amino acids) for 124 KNOX sequences. Alignments of nucleotides and amino acids were employed in the subsequent Bayesian analysis. Bayesian phylogenetic analysis was performed using Mr Bayes 3.2.1 (Huelsenbeck and Ronquist, 2001; Huelsenbeck et al., 2001). The Bayesian analyses for the nucleotide data sets were run for 15,000,000 (BELL) or 5,000,000 (KNOX) generations, which was sufficient for convergence of the two simultaneous runs (BELL, 0.0357; KNOX, 0.0423). In both cases, to allow for the burn-in phase, 50% of the total number of saved trees were discarded. The graphic representation of the trees was generated using the FigTree (version 1.4.0) software (http://tree.bio.ed.ac.uk/software/figtree/). Sequence alignments and command files used to run the Bayesian phylogenetic analyses can be provided upon request.In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.Acceptance summary:This work is an exciting and critical study that speaks to the field of evolutionary developmental biology beyond plant biology. It focuses on life cycle evolution in green plants, where the haploid-dominant life cycle in green algae evolved into the diploid-dominant one in land plants. The thorough analysis of molecular phylogeny, and a detailed investigation of developmental expression provides a comprehensive account of entire TALE class proteins, and is more than expected in a single study.Decision letter after peer review:Thank you for submitting your article "Gamete-specific expression of TALE class HD genes activates the diploid sporophyte program in Marchantia polymorpha" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Sheila McCormick as Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Christian Hardtke as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Jae-Hyeok Lee (Reviewer #2).The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.Summary:The manuscript of Dierschke et al., represents an exciting and critical study that speaks to the field of evolutionary developmental biology beyond plant biology. It focuses on life cycle evolution in green plants, where the haploid-dominant life cycle in green algae evolved into the diploid-dominant one in land plants. TALE class HD TF KNOX and BELL subfamilies have been implicated in the haploid-to-diploid transition in the green lineage, in the unicellular Chlamydomonas and in Physcomitrella. The authors therefore investigated the function of TALE class homeodomain transcription factors in the liverwort M. polymorpha. The authors use a range of molecular genetics and genomic approaches to demonstrate that MpBELL234 are expressed in male gametophyte cells (antheridia), while MpKNOX1 is expressed in the female egg cell before fertilization and during sporophyte development. They provide evidence that maternally supplied MpKNOX1 is required for later sporophyte development and that paternally supplied MpBELL234 is required for zygotic development. They found two distinct heterodimeric complexes of TALE class proteins, MpKNOX1 plus MpBELL2/3/4 and MpKNOX2 plus MpBELL1. They used BiFC to demonstrate the interaction between BELL and KNOX. Finally, they suggest that Polycomb, a repressive chromatin modifying complex, is required to repress the MpBELL1 and MpKNOX2 genetic program during gametophyte development. The thorough analysis of molecular phylogeny, and a detailed investigation of developmental expression by combining RNA-seq and ChIP-seq, provides a comprehensive account of entire TALE class proteins, and is more than expected in a single study.Essential revisions:We have major concerns about the ChIP-seq approach that was used to demonstrate the absence of H3K27me3 on the target loci. It appears that the ChIP-seq was not replicated, and the presence/absent of histone PTMs on the loci under study has not been verified by other means. Normally the claims made in the manuscript would need to be backed up by either another replication of the ChIP-seq experiment or with validation using ChIP-qPCR. More details are presented below. In the event that this is not possible now, the phrasing should be changed to emphasize that this needed data is missing (i.e. qualify the claims).Furthermore, the first part of the paper (p. 5-19) is about gene structure and phylogeny and we don't get to the knockdown experiment until p. 20. As these aspects are not represented in the abstract, we suggest that this section of the paper is minimized. In the discussion the phylogeny is not discussed until the end, again giving the impression that it is not that important. The lengthy molecular phylogeny section would need to be further elaborated into a significant finding or provided as a brief supplement for improving the narrative's focus. The discussion is about 8 pages and probably could be condensed to the main points.1. This manuscript reports that KNOX2 and BELL1 form heterodimers, whose functional study may await an analysis of their knockout plants. While the authors consider KNOX2/BELL1 unlikely to be involved in the haploid-to-diploid transition, KNOX2 in the moss, Physcomitrella, plays a critical role in preventing gametophyte development from resuming in diploid sporophytes. Thereby, a de-repressed KNOX2/BELL1 heterodimer may produce aberrant phenotypes during gametophyte development in Marchantia.2. The authors investigated if the repression of TALE-HD (MpBELL1 and MpKNOX2) in the vegetative gametophyte occurs via PRC2. They use an inducible system with miR to knockdown MpE(z). The expectation is that E(z) knockdown will decrease PRC2 activity and hence reduce overall H3K27me3 marks within the genome. It's a rather 'brutal' approach because it is likely that many loci genome-wide will be de-repressed. But the manuscript gives no information about the overall genomic landscape after the E(z) miR knockdown. If available, describe phenotypes of the amiR-MpE(z)1 plants following 17-b-estradiol treatment in addition to the TALE gene expression. For the two experiments using transgenic plants harboring conditional amiR-MpE(z) or ectopically expressed MpKNOX1 or MpBELL3, the effects of the transgenes may differ between male and female gametophytes. Please indicate whether both or single-sex plants were analyzed.3. Detailed analysis of MpBELL2/3/4 knockouts in the male gametophytes indicates their essential role in post-fertilization sporophyte development. However, the two observed phenotypes, zygotic, or early embryonic arrest, may be due to the expression of BELL2/3/4 from the female nucleus, as indicated by the authors (P44 L6). Since both Mpbell2/3/4 male and female plants were reported, it is reasonable to ask for the phenotypes of the sporophytes produced by Mpbell2/3/4 knockout male and female.4. The authors state on page 21 "Compared to the overall H3K27me3 landscape in wild-type gametophytic tissue, the number and amplitude of methylation peaks are consistently reduced after 48h of down-regulation of MpE(z)1." What is meant by "consistently"? This statement needs to be backed up with data and proper statistical analysis. This requires at least one replicate experiment, and certainly no claims can be made about "amplitude" of the peaks (to compare ChIP-seq peaks quantitatively 'spike in' experiments are needed, etc., etc, – not a trivial task). One option could be to take a few dozen loci and carry out ChIP-qPCR experiments.It is stated that MpKNOX2 and MpBELL1 show altered H3K27me3 patterns (which provides the basis for the conclusion that polycomb is required to repress BELL1 and KNOX2 but not KNOX1 and BELL234). The validation of ChIP-seq data using another replicate of ChIP-seq or ChIP-qPCR is also crucial to back up the following sentence "The three primarily gametophyte-expressed KNOX1 genes were neither marked with H3K27me3 nor induced in the amiR-MpE(z)1 background".5. Details of the ChIP-seq data (sequencing stats, coverage, %uniquely mapped reads, n. reads, was more than one peak caller used, FRip values, n. peaks, etc). Is the GEO reference for the ChIP-seq data missing? In the methods section the authors mention MACS as a peak caller. A table should be provided with loci/peaks that have been called before/after induction.Revisions expected in follow-up work:Another replication of the ChIP-seq experiment or validation using ChIP-qPCR.Essential revisions:We have major concerns about the ChIP-seq approach that was used to demonstrate the absence of H3K27me3 on the target loci. It appears that the ChIP-seq was not replicated, and the presence/absent of histone PTMs on the loci under study has not been verified by other means. Normally the claims made in the manuscript would need to be backed up by either another replication of the ChIP-seq experiment or with validation using ChIP-qPCR. More details are presented below. In the event that this is not possible now, the phrasing should be changed to emphasize that this needed data is missing (i.e. qualify the claims).Similar H3K27me3 data for Marchantia polymorpha have been published in Montgomery et al., (2020 Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin. Current Biology 30 573-588), including peaks at MpKNOX2 and MpBELL1 (This publication actually used our unpublished results for corroboration of their results). Thus, given we are (1) unlikely to repeat these experiments and (2) the data were only included into this story to provide an explanation of why we did not include MpKNOX2 and MpBELL1 in further analyses, we have removed all chromatin state data from the present paper. The omission is further supported by the rather minimal advancement in scientific understanding in the experiments and the monetary limitations or repeating the experiments.Furthermore, the first part of the paper (p. 5-19) is about gene structure and phylogeny and we don't get to the knockdown experiment until p. 20. As these aspects are not represented in the abstract, we suggest that this section of the paper is minimized. In the discussion the phylogeny is not discussed until the end, again giving the impression that it is not that important. The lengthy molecular phylogeny section would need to be further elaborated into a significant finding or provided as a brief supplement for improving the narrative's focus. The discussion is about 8 pages and probably could be condensed to the main points.We believe the phylogenetic aspect, a good fraction of which is supplemental, is important to provide the context of which genes we are examining — while MpKNOX1 is a typical land plant KNOX1 gene (updated Figure 1— figure Supplement 2), the MpBELL genes we describe are phylogenetically distinct from any land plant BELL gene that has been described previously. During the long lockdown of the pandemic, we have revisited this analysis by including more sequences from other liverworts and mosses. We can confidently state that the ancestral land plant had three BELL paralogs, only one of which was retained by flowering plants. However, all three paralogs were retained by liverworts and at least two were retained by at least some mosses (updated Figure 2). Further, sequence comparisons suggest that each of the three BELL paralogs might recognize different target sequences (new Figure 2— figure Supplement 1). Thus, we think it necessary to include such data describing which genes we are analysing prior to discussing phenotypes of loss-of-function alleles. This is now alluded to in the abstract.1. This manuscript reports that KNOX2 and BELL1 form heterodimers, whose functional study may await an analysis of their knockout plants. While the authors consider KNOX2/BELL1 unlikely to be involved in the haploid-to-diploid transition, KNOX2 in the moss, Physcomitrella, plays a critical role in preventing gametophyte development from resuming in diploid sporophytes. Thereby, a de-repressed KNOX2/BELL1 heterodimer may produce aberrant phenotypes during gametophyte development in Marchantia.Since MpKNOX2 and MpBELL1 are not expressed in either female or male reproductive tissues and only commence their expression after sporophyte development has progressed we did not consider these genes as candidates for the haploid to diploid transition. However, we included data on their protein-protein interactions (updated Figure5 — figure Supplement 1) to provide further evidence that they are not involved in the haploid to diploid transition as they only interact with one another, and not with MpKNOX1 or MpBELL3/4.2. The authors investigated if the repression of TALE-HD (MpBELL1 and MpKNOX2) in the vegetative gametophyte occurs via PRC2. They use an inducible system with miR to knockdown MpE(z). The expectation is that E(z) knockdown will decrease PRC2 activity and hence reduce overall H3K27me3 marks within the genome. It's a rather 'brutal' approach because it is likely that many loci genome-wide will be de-repressed. But the manuscript gives no information about the overall genomic landscape after the E(z) miR knockdown. If available, describe phenotypes of the amiR-MpE(z)1 plants following 17-b-estradiol treatment in addition to the TALE gene expression. For the two experiments using transgenic plants harboring conditional amiR-MpE(z) or ectopically expressed MpKNOX1 or MpBELL3, the effects of the transgenes may differ between male and female gametophytes. Please indicate whether both or single-sex plants were analyzed.These data have been removed from the new version of the manuscript, but to satisfy the reviewers' curiosity, constitutive expression of amiR-MpE(z)1 is lethal.3. Detailed analysis of MpBELL2/3/4 knockouts in the male gametophytes indicates their essential role in post-fertilization sporophyte development. However, the two observed phenotypes, zygotic, or early embryonic arrest, may be due to the expression of BELL2/3/4 from the female nucleus, as indicated by the authors (P44 L6). Since both Mpbell2/3/4 male and female plants were reported, it is reasonable to ask for the phenotypes of the sporophytes produced by Mpbell2/3/4 knockout male and female.In the limited time available for experimental work due to pandemic-related shutdowns, we have focussed our attention on the question posed in this comment. Furthermore, in our communications with the authors of the manuscript cosubmitted with our manuscript (Hisanaga et al.,), we have sought to reconcile differences in our results. There were two key differences between our analyses. First, in Hisanaga et al., Mpknox1 sporophytes were obtained at a moderate frequency. Second, in Hisanaga et al., Mpbell3/4 mutations affected both males and females, whereas we had only observed an effect in males. These two discrepancies were explored and are discussed below with respect to the new data that has been added to the present manuscript.First, we explored whether some of the discrepancies may be attributed to differences in our observations due to having two very different techniques to generate sporophytes. In Hisanaga et al., the fertilization event is in vitro, allowing a detailed analysis of early events, but precluding observations beyond two weeks in these plants. In contrast, we have analysed the results of crosses in vivo. In our new analyses (new Figure 5A,B, leftmost panel) we made detailed observation of the results of crosses between wild-type males and wild-type females, which uncovered hitherto undescribed embryo abortion, mostly at about the 1-week developmental stage. We speculate that this is due to signalling within the female bearing multiple sporophytes, as such systems are common in angiosperms, where fruits are aborted if there is an 'over'-abundance of fruit set. While this is an interesting phenomenon, we did not explore it further except to acknowledge that this provides a baseline to which other crosses need to be compared.MpKNOX1: We performed additional crosses to quantify the maternal effects Mpknox1 mutants. In no case were we able to observe any mature sporophytes, and only in rare cases was any development beyond the zygote stage observed (new Figure 3—figure Supplement 1). Thus, we believe the difference between our results and those of Hisanaga et al., are that our allele is null (large deletion), while the allele of Hisanaga et al., is a hypomorph in which translational read-through may produce functional protein at a low level. However, we made an additional observation that implies a function of MpKNOX1 maintenance of female reproductive tissues in the absence of fertilization (new Figure 3M-N).MpBELL234: As described above for wild type and Mpknox1, we also performed a large number of reciprocal crosses between wild type and Mpbell34 mutants. We confirmed the result of Hisanaga et al., that Mpbell34 mutants have defects when in either male or female, and when Mpbell34 is mutant in both, most sporophytes abort at the zygote stage (new Figure 5AB). Thus, our results are now concordant. However, this brought up additional questions as our reporter genes for MpBELL2, MpBELL3 and MpBELL4 were expressed only in the male (Figure 4A-F). These reporters were constructed using sequences 5' to the longest transcripts produced from the loci (new Figure 4G). Thus, we performed additional experiments to detect MpBELL3 and MpBELL4 expression by in situ hybridization (new Figure 4H and new Figure 4— figure Supplement 2F). We then explored why our reporter genes exhibited male specific expression. Strand-specific analyses of available transcriptome data (Figure 4— figure Supplement 1B) suggested female and sporophyte expressed transcripts at both the MpBELL3 and MpBELL loci; however, these transcripts were much shorter than the transcripts expressed in males, but at least in the case of MpBELL3 would encode the carboxyl-most homeodomain (Figure 4— figure Supplement 1B). We further investigated these alternative transcripts by searching for potential alternative transcription start sites and using a combination of 5’- and 3’-RACE and rtPCR, we defined the transcripts present in the male, female, and sporophyte (new Figure 4— figure Supplement 2 and new Figure 4— figure Supplement 3). Using this data, we constructed a reporter gene using putative regulatory sequences internal to MpBELL4 (new Figure 4G). This reporter line drove expression in the female, specifically in the archegonia and egg cells (new Figure 4I), but was not expressed in sporophytes (new Figure 4J), perhaps due to some regulatory sequences lacking in our construct. These observations have been incorporated into the discussion of how TALE-HD genes in Marchantia activate the diploid genetic program.4. The authors state on page 21 "Compared to the overall H3K27me3 landscape in wild-type gametophytic tissue, the number and amplitude of methylation peaks are consistently reduced after 48h of down-regulation of MpE(z)1." What is meant by "consistently"? This statement needs to be backed up with data and proper statistical analysis. This requires at least one replicate experiment, and certainly no claims can be made about "amplitude" of the peaks (to compare ChIP-seq peaks quantitatively 'spike in' experiments are needed, etc., etc, – not a trivial task). One option could be to take a few dozen loci and carry out ChIP-qPCR experiments.It is stated that MpKNOX2 and MpBELL1 show altered H3K27me3 patterns (which provides the basis for the conclusion that polycomb is required to repress BELL1 and KNOX2 but not KNOX1 and BELL234). The validation of ChIP-seq data using another replicate of ChIP-seq or ChIP-qPCR is also crucial to back up the following sentence "The three primarily gametophyte-expressed KNOX1 genes were neither marked with H3K27me3 nor induced in the amiR-MpE(z)1 background".These data have been removed from the new version of the manuscript.5. Details of the ChIP-seq data (sequencing stats, coverage, %uniquely mapped reads, n. reads, was more than one peak caller used, FRip values, n. peaks, etc). Is the GEO reference for the ChIP-seq data missing? In the methods section the authors mention MACS as a peak caller. A table should be provided with loci/peaks that have been called before/after induction.These data have been removed from the new version of the manuscript.Revisions expected in follow-up work:Another replication of the ChIP-seq experiment or validation using ChIP-qPCR.This data has been removed from the new version of the manuscript.
Authors: Marcus B Noyes; Ryan G Christensen; Atsuya Wakabayashi; Gary D Stormo; Michael H Brodsky; Scot A Wolfe Journal: Cell Date: 2008-06-27 Impact factor: 41.582
Authors: Jan de Vries; Angela Melanie Fischer; Mayo Roettger; Sophie Rommel; Henriette Schluepmann; Andrea Bräutigam; Annelie Carlsbecker; Sven Bernhard Gould Journal: New Phytol Date: 2015-09-11 Impact factor: 10.151
Figure 1.
Figure 1—figure supplement 1.
Figure 3—figure supplement 3.
Figure 2.
Figure 2—figure supplement 1.
Figure 3.
Figure 3—figure supplement 1.
Figure 3—figure supplement 2.
Figure 4.
Figure 4—figure supplement 1.
Figure 4—figure supplement 2.
Figure 5.
Figure 5—figure supplement 1.
Figure 5—figure supplement 2.
Figure 6.
Figure 6—figure supplement 1.