Identification of the sex-determining factor in the liverwort Marchantia polymorpha reveals unique evolution of sex chromosomes in a haploid system.

Sex determination is a central process for sexual reproduction and is often regulated by a sex determinant encoded on a sex chromosome. Rules that govern the evolution of sex chromosomes via specialization and degeneration following the evolution of a sex determinant have been well studied in diploid organisms. However, distinct predictions apply to sex chromosomes in organisms where sex is determined in the haploid phase of the life cycle: both sex chromosomes, female U and male V, are expected to maintain their gene functions, even though both are non-recombining. This is in contrast to the X-Y (or Z-W) asymmetry and Y (W) chromosome degeneration in XY (ZW) systems of diploids. Here, we provide evidence that sex chromosomes diverged early during the evolution of haploid liverworts and identify the sex determinant on the Marchantia polymorpha U chromosome. This gene, Feminizer, encodes a member of the plant-specific BASIC PENTACYSTEINE transcription factor family. It triggers female differentiation via regulation of the autosomal sex-determining locus of FEMALE GAMETOPHYTE MYB and SUPPRESSOR OF FEMINIZATION. Phylogenetic analyses of Feminizer and other sex chromosome genes indicate dimorphic sex chromosomes had already been established 430 mya in the ancestral liverwort. Feminizer also plays a role in reproductive induction that is shared with its gametolog on the V chromosome, suggesting an ancestral function, distinct from sex determination, was retained by the gametologs. This implies ancestral functions can be preserved after the acquisition of a sex determination mechanism during the evolution of a dominant haploid sex chromosome system.

Introduction

Sex-specific gene expression leading to sexual differentiation is generally triggered by a sex-determining gene, often carried by a sex chromosome. Since the early 20th century, different sex determination systems by sex chromosomes, including XX/XY, ZW/ZZ, XX/XO, and ZO/ZZ, have been described in animals. Similarly, some flowering plants possess XX/XY (or ZW/ZZ) systems directing sex-specific differentiation, leading to flower development during the diploid life phase.2, 3, 4 Sex chromosomes in some flowering plants, including persimmon and poplar, have a single master gene that directs sex determination, as in some animal XX/XY systems. In other flowering plants, such as kiwifruit and asparagus, however, sex determination involves allelic differences in two linked genes on the heterogametic sex chromosome, which act to promote male development and suppress female development. It is also notable that sex chromosomes in a number of dioecious angiosperm species are homomorphic or show minor structural divergence. In diploid systems, suppression of recombination between the sex chromosomes in heterozygotes restricts structural and functional differentiation to the sex chromosome specific to either males (Y) or females (W). The other sex chromosome (X or Z) remains similar to the ancestral state, carrying functional copies of ancestrally present genes, whose presence allows genetic degeneration of the Y- or W-linked regions, sometimes including loss of most ancestral genes. Dosage compensation may evolve to compensate for males’ low X-linked gene expression in XY systems (or females’ low Z expression in ZW systems). Understanding of sex determination mechanisms and sex chromosome evolution has largely come from studying diploid organisms.

In contrast to animals and flowering plants, sexual differentiation takes place in the dominant haploid life phase of bryophytes (e.g., liverworts and mosses) and most algae,, with species with separate male and female individuals termed dioicous. In these groups, males and females are defined by their gamete size, with females producing large immotile egg cells and males producing small motile sperm., Sex chromosomes in haploid systems show fundamentally distinct genetic and evolutionary behavior from those in diploid systems in which one sex is hemizygous (XY or ZW) and one homozygous (XX or ZZ)., In haploid systems, each individual carries only a single sex chromosome. The female and male sex chromosomes are termed U and V, respectively, to emphasize the differences from X and Y (or Z and W) chromosomes. Fertilization results in diploid UV sporophytes, a brief non-sexually differentiated phase of the life cycle, another difference from diploid systems. Due to their similar inheritance, the U and V chromosomes are predicted to evolve similarly. Specifically, both chromosomes are predicted to undergo little or no degeneration, and any degeneration is predicted to affect both the U and V chromosomes similarly, in contrast to the diploid systems outlined above.

Many liverwort species are dioicous and often have heteromorphic sex chromosomes, and dioicy is the likely ancestral condition for liverworts. In Marchantia polymorpha, plants possess either a U or a V sex chromosome19, 20, 21 and develop into females or males, respectively. Genetic analyses of phylogenetically diverse liverworts indicate that the U chromosome carries a feminizing locus and that the V chromosome carries multiple loci required for sperm motility. Further, haploid plants spontaneously containing both U and V chromosomes develop as females, supporting the presence of a dominant U-linked feminizing locus, in whose absence plants develop as males. Genomic analyses of M. polymorpha have shown that both U and V chromosomes have lower gene densities, and higher repetitive content, than the autosomes,19, 20, 21 similar to the symmetrical changes in other non-recombining U/V systems, in the mosses Ceratodon purpureus and Syntrichia caninervis and the brown alga Ectocarpus. However, the sex-determining genes have yet to be identified in liverworts or any other complex multicellular haploid organism, and only recently has it become possible to study the evolution of sex chromosomes in haploid systems and test their predicted differences from diploid systems.

Here, we identify a U-linked transcription factor gene, BASIC PENTACYSTEINE ON THE U CHROMOSOME (BPCU), as the master sex determinant, Feminizer, in M. polymorpha. BPCU has a homolog (gametolog) on the V chromosome, BASIC PENTACYSTEINE ON THE V CHROMOSOME (BPCV), and both share a conserved function in reproductive induction, suggesting retention of an ancestral function despite the BPCU gametolog evolving an additional function in female sex determination. Phylogenetic analyses of BPCU and BPCV, and their orthologs in other liverworts suggest that these gametologs, and thus the U and V chromosomes, had been already established 430 mya in the ancestral liverwort.

Results

Chromosome-level assembly of the U chromosome

As a basis for the search for Feminizer, we constructed a chromosome-level assembly of the U chromosome using PacBio long-read sequencing and Hi-C technologies, as was done previously for the V chromosome. The U chromosome carries rDNA repeats at each of its ends, and, between these, includes a 4.5-Mb region containing genes. Previous analyses identified 74 and 154 genes on the U and V chromosomes, respectively,, but the chromosome-level assembly of the U chromosome and manually refined annotation of both chromosomes (excluding transposable element genes and merging some split genes) result in smaller numbers of 47 and 93, respectively (Data S1A). As was noted previously, most U chromosome genes that do not have a V gametolog appear to be fragmentary or have highly similar autosomal paralogs, suggesting that they may be recent immigrants to the U chromosome (Data S1A). The eight M. polymorpha autosomes carry much larger numbers of genes, between 1,900 and 3,000.

BPCU on the U chromosome is necessary and sufficient for female sex determination

We confirmed the aforementioned predictions of a dominant U-linked sex-determining Feminizer locus in M. polymorpha, because plants carrying both U and V sex chromosomes that arose spontaneously via meiotic non-disjunction (wild type [WT] [UV], the sex chromosome background is given in brackets for clarity) developed female sex organs (Figures 1A–1C). To identify the U-linked Feminizer gene, we screened, via genome editing of U-chromosome candidate genes, for mutations that impact sex organ morphology (Data S1A). Two independently obtained mutants of BPCU formed male sex organs (Figures 1D–1F), suggesting that BPCU is necessary for development of female sex organs. Surprisingly, both of these bpcU mutants possessed not only a U chromosome harboring the mutated bpcU gene but also a V chromosome (Figure 1G). This observation and analysis by flow cytometry (Figure S1) suggest that the two bpcU mutants were derived from spores that had inherited both sex chromosomes, presumably via non-disjunction. To further test BPCU’s dominant feminizing activity, we introduced a genomic fragment of BPCU spanning 2 kb upstream to 3 kb downstream of the BPCU transcribed region into the genome of WT male plants that carry only the V chromosome. The genetically male plants transformed with BPCU (gBPCU [V]) displayed a male-to-female sex conversion phenotype, although no egg cells developed within their sex organs (Figure 1H), a phenotype previously reported for an autosomal sex conversion mutant. The loss of egg cell differentiation in these genotypes is likely due to the lack of U-linked egg-cell differentiation gene(s). We concluded that BPCU is necessary and sufficient to induce the development program of female sex organs.

Figure 1

BPCU on the U chromosome is necessary and sufficient for female sex determination

(A–E) Sexual development of wild-type (WT) and mutant plants. Sex organs of female and male WT plants are shown in (A) and (B), respectively. Photographs were taken from female WT [U] (Tak-2; A), male WT [V] (Tak-1; B), UV plant (C), genome-editing mutants of bpcU (bpcU-101 [UV], D; bpcU-102 [UV], E). Left and right panels: receptacles and gamete-containing organs are shown, respectively. Arrowheads indicate egg cells. Scale bars, 1 mm (left panels); 50 μm (right panels).

(F) Mutations in bpcU. The nucleotide sequences of the WT and the mutants (bpcU [UV]) are aligned. The sequence for guide RNA (gRNA) for genome editing is underlined. The protospacer adjacent motif (PAM) sequence for CRISPR/Cas9 is shown in bold. Deletions are shown in red.

(G) Diagnosis of genetic sex using U-chromosome- and V-chromosome-linked markers.

(H) Feminization of the genetic male lines transformed with BPCU (gBPCU [V]). Sex organs 1 month after reproductive induction are shown for gBPCU#4 [V]. Left and right panels; receptacle and gamete-containing organ are shown, respectively. Scale bars, 1 mm (left panel); 50 μm (right panel).

See also Figure S1.

BPCU acts as Feminizer by regulating the autosomal sex-determining locus, FGMYB/SUF

We next analyzed the expression profiles of sex-specific genes in the sex-converted plants and detected gene expression patterns reflecting the morphological sex phenotype (Figure 2A). Critically, expression of the autosomal locus FEMALE GAMETOPHYTE MYB (FGMYB), which promotes female sex differentiation, was correlated with the morphological sex observed in different genotypes. In males, the expression of FGMYB is suppressed by expression of an antisense long non-coding RNA called SUPPRESSOR OF FEMINIZATION (SUF), which is transcribed from a promoter located 3′ of the FGMYB locus. Consistent with the sex conversion phenotype, the mutants (bpcU-101 and bpcU-102) had decreased FGMYB and increased SUF expression compared with WT females (Figure 2B, left). Furthermore, the BPCU-carrying genetically male plants (gBPCU [V]) had decreased SUF and increased FGMYB expression compared with WT males (Figure 2B, right).

Figure 2

BPCU acts as Feminizer by regulating the autosomal sex-determining locus, FGMYB/SUF

(A) RT-PCR analysis of sex-specific gene expression in bpcU [UV] and gBPCU [V]. Total RNAs from WT [U] (Tak-2), WT [V] (Tak-1), UV plant, bpcU (bpcU-101 [UV] and bpcU-102 [UV]), and gBPCU [V] (lines no. 3 and no. 4) were used for detecting expression of the genes indicated above. Mp3g14390 and Mp6g01560 are genes expressed specifically in males and Mp5g13880 and Mp7g02820 in females. Mp3g23400 (EF1) is a ubiquitously expressed gene.

(B) Expression of FGMYB and SUF. qRT-PCR was performed for RNAs from sex organs. Experiments shown in the left and right panels were performed separately for the loss-of-function lines (bpcU-101 [UV]; bpcU-102 [UV]) and the gain-of-function lines (gBPCU [V] no. 3 and no. 4), respectively, with WT [U] (Tak-2) and WT [V] (Tak-1) as controls. Bars represent mean ± SE. Symbols above the bars indicate grouping by p < 0.05 in a Tukey-Kramer test (n = 3).

(C) Genetic suppression of BPCU-dependent feminization by fgmyb mutation. Sex organs 1 month after reproductive induction are shown for two independently obtained fgmyb mutants. Left and right panels: receptacles and gamete-containing organs are shown, respectively. Scale bars, 1 mm (left panels); 50 μm (right panels).

(D) H3K27me3 status at the FGMYB/SUF locus in the vegetative stage in WT [U], WT [V], and bpcU-1 [U]. See Figure 3 for the genotype and the phenotype of bpcU-1 [U]. IGV browser screenshot is shown. H3K27me3 and RNA-seq tracks show bigwig files normalized to 1 × genomic coverage (reads per genome coverage [RPGCs]) and bins per million mapped reads (BPMs), respectively (window size = 10 bp). Peak calling was done using SICER2 with the default parameters (window size = 200 bp; false discovery rate [FDR] = 0.01). H3 was used as a control for peak calling.

See also Figures S2 and S3.

Figure 3

BPCU and its gametolog BPCV are required for induction of sexual reproduction, but only BPCU confers the feminization

(A) Gene organizations of BPCU and BPCV loci and schematic illustrations of bpcU and bpcV mutations. The coding regions are colored. Gray boxes indicate 5′ and 3′ untranslated regions. Arrowheads indicate the positions of gRNAs. The deleted region in BPCU is shown by a bent broken line.

(B) No sex organ formation in bpcU [U] and bpcV [V]. Female WT [U] (Tak-2; top left), male WT [V] (Tak-1; bottom left), bpcU-1 [U] (top right), and bpcV-1 [V] (bottom right) plants are shown. Images 1 month after reproductive induction are shown. Arrowheads indicate sex organs. Scale bar, 10 mm.

(C) Genetic complementation of the reproductive induction phenotype of bpcU with BPCU or BPCV. Genomic BPCU-transformed bpcU [U] (gBPCU#2bpcU-1 [U]; top) and genomic BPCV-transformed bpcU [U] (gBPCV#7bpcU-1 [U]; bottom) plants are shown. Arrowheads indicate sex organs. Scale bar, 10 mm.

(D) Magnified images of sex organs in gBPCU#2bpcU-1 [U] (top) and gBPCV#7bpcU-1 [U] (bottom). Sex organs 1 month after reproductive induction are shown. Left and right panels: receptacles and gamete-containing organs are shown, respectively. Scale bars, 1 mm (left panels); 50 μm (right panels).

(E) Loss of BPCU or BPCV impairs reproductive induction. Gemmae were grown under white light for 10 days and transferred to the reproductive induction condition (STAR Methods). Days at which a visible receptacle was first formed after reproductive induction are shown. No sex organs were formed in bpcU (bpcU-1, bpcU-2, and bpcU-3) and bpcV (bpcV-1, bpcV-2, and bpcV-3) over 60 days. Bars represent mean ± SD (n = 3).

See also Figures S4 and S5.

BPCU belongs to an alanine-zipper-containing group of the BARLEY B RECOMBINANT/BASIC PENTACYSTEINE (BBR/BPC) proteins that bind to GAGA elements, and repress gene expression via epigenetic regulation. We demonstrated that, in vitro, BPCU binds to GAGA motifs enriched in the promoter of SUF (Figure S2). Thus, we hypothesized that, in females, BPCU directly suppresses SUF expression, consequently relieving FGMYB from suppression by SUF. In support of this hypothesis, fgmyb mutations in the gBPCU [V] background suppressed the feminization induced by gBPCU (Figures 2C and S3). We also observed differences in chromatin status with the repressive H3K27me3 peaks expanded into the gene region of SUF in WT [U], but not in WT [V] and bpcU [U] (Figure 2D), which is consistent with a model that BPCU directly represses SUF through chromatin modification. These observations indicate that BPCU acts as the U-linked Feminizer and that it directly regulates the autosomal sex-determining locus complex, FGMYB/SUF, in M. polymorpha.

BPCU and its gametolog, BPCV, are required for induction of sexual reproduction, but only BPCU confers the feminization

In contrast with other sex-determining genes identified in diploid organisms and also in the haploid alga Volvox carteri, the Feminizer, BPCU, has a paired homolog, BPCV, on the male sex chromosome (Figure S4A). Such gene pairs are called “gametologs” and originate from a common ancestral autosomal gene. Gametologs generally conserve essential functions that are maintained during sex chromosome evolution. Indeed, the highly conserved C-terminal domain of BPCV, which is known as a DNA-binding domain of BPCs, also bound the same target sequences as BPCU in vitro (Figure S2A). BPCU and BPCV are both expressed primarily in the female and male reproductive organs, respectively, although they are expressed in all tissues examined (Figure S4B).

Because the bpcU mutants (bpcU-101 and bpcU-102) obtained by the initial screening were aneuploid with the V chromosomes, we generated bpcU as well as bpcV mutants in a haploid background (Figures 3A and S5). Both bpcU and bpcV haploid mutants failed to induce sex organ formation under growth conditions that induce reproductive development in WT (Figures 3B and 3E). The introduction of a genomic fragment of BPCU or BPCV complemented the reproductive induction phenotype of bpcU (Figures 3C–3E). However, bpcU mutants complemented with BPCV developed male sex organs instead of female sex organs (Figures 3C and 3D). These observations indicate that BPCU and BPCV share a function required to induce sex organ development.

The sex chromosomes in liverworts were established around 430 mya, which is the oldest among the known sex chromosomes

We first examined the antiquity of the BPCU and BPCV gametologs within liverworts and noted that both gametologs are present in species representing the three extant lineages of liverworts, indicating that the divergence of these gametologs began in the ancestral liverwort (Figure 4). Further, the ancestral land plant possessed a single BPC transcription factor gene,, a condition likely retained in common ancestors of each of the bryophyte lineages. The identification of BPCU as Feminizer in M. polymorpha sheds light on liverwort U/V sex chromosome evolution. Sex chromosomes have their evolutionary origin when an ancestral gene on an autosomal pair mutates or when a gene duplicates onto one member of this chromosome pair and functions as a sex determiner. During the evolution of a U/V sex chromosome system, the sex-determining gene might have arisen in a non-recombining region or a non-recombining region might have evolved subsequently. Thereafter, successive chromosomal rearrangements could create evolutionary strata of U-V sequence divergence and potential genetic degeneration, which, as outlined above, should affect both members of the pair.,, However, essential genes required in the haploid generation could be retained on both the U and V and over time diverge into pairs of gametologs that differ in sequence but retain a conserved function. In contrast, genes necessary in only one sex should be retained on the respective sex chromosome and could be lost from the other.

Figure 4

Phylogenetic tree of BPCU and BPCV

The phylogram was reconstructed using nucleotide alignments as described in the STAR Methods section. The liverwort sequences are shown in purple, and the clades containing BPCU and BPCV are highlighted in yellow and turquoise, respectively. The phylogenetic positions of the BPCU and BPCV of M. polymorpha are demarcated by arrows with the pS value between the two M. polymorpha gametologs indicated (0.77). Taxa are color coded as follows: liverworts, purple; Haplomitriopsida, dark purple; Marchantiopsida, bright purple; Jungermanniopsida, violet; hornworts, dark green; mosses, light green; lycophytes, orange; and ferns, brown. Numbers at branches indicate posterior probability values. See also Data S1B and S2.

Because the sex-determining gene or genes should be located in the oldest stratum, the M. polymorpha BPCU probably identifies the oldest non-recombining region of the U chromosome. In this plant, divergence of U-V genes could not be estimated using synonymous site substitution values (pS), because these values are at or near saturation for almost half of the 19 gametolog pairs (Data S1B). Instead, we therefore used phylogenetic analyses based on all alignable coding sequences of each of the gametologs (Figure 4; Data S2). The timing of the incorporation of each pair of M. polymorpha gametologs into non-recombining regions of the M. polymorpha sex chromosomes can be inferred by the tree topology. If the M. polymorpha U and V gametologs are less diverged from Haplomitriopsida sequences than from each other, the divergence of the M. polymorpha gametologs was inferred to have occurred before the divergence of the Haplomitriopsida and the Marchantiopsida. Conversely, if the M. polymorpha U and V gametologs are less diverged from one another than they are to some other liverwort sequences, then their divergence time was inferred based on the relative positions in the tree topology of the Jungermanniopsida and other Marchantiopsida sequences. These analyses suggest that the M. polymorpha U and V chromosomes consist of a mosaic of at least five evolutionary strata, with 7 pairs of “ancient” gametologs, including BPCU and BPCV, pre-dating the divergence of extant liverworts (purple horizontal regions in Figure 5). Other possible strata can be defined based on the other patterns of M. polymorpha U and V gametolog divergence relative to orthologs from other liverwort species (Figure 4; Data S2; see the Figure 5 legend). No pseudoautosomal regions were identified.

Figure 5

Feminizer and sex chromosome evolution in M. polymorpha

A schematic of the genomic architecture of the M. polymorpha ruderalis U and V chromosomes is displayed (vertical bars on the left) with gametologs positioned between the chromosomes and U- and V-specific genes represented by short lines (for expression patterns, see Data S1A). The U chromosome sequence was obtained from a chromosome-scale assembly of the female accession, Tak-2; an updated list of its annotated genes is in Data S1A along with V-linked gametologs and previous (v3.1) annotations where applicable. Results of phylogenetic analyses of Marchantiophyta orthologs are presented as a matrix on the right. Open circles indicate outgroup sequences that diverged before the M. polymorpha U and V divergence. These were inferred when sequence divergence between the M. polymorpha U and V gametologs was smaller than that between either gametolog and the ortholog in the other species indicated. Ingroup sequences (solid circles) were defined as orthologs whose sequence divergence from one of the M. polymorpha gametologs was smaller than the divergence between the two M. polymorpha gametologs, indicating divergence after the M. polymorpha U-V gametolog split. Specific designations of each ortholog as either U or V gametolog related are presented in Data S1C. See Data S2 and Figure 4 for phylogenetic trees of each of the gametologs. For the Marchantiopsida and Jungermanniopsida, we included only dioicous species where the transcriptome (∗) or genome (°) was apparently derived from a single sex (i.e., where only single orthologs were detected for all the genes that we infer to be in the oldest evolutionary stratum [shaded purple]). Predicted nodes of the origins (color coded) of the non-recombining regions of the M. polymorpha U and V chromosomes are indicated on the chromosome schematic and the presently accepted liverwort phylogenetic tree (above the matrix)., Geologic ages of numbered nodes as previously estimated, are listed at left of the tree. See also Data S1A, S1C, and S2.

Discussion

The chromosome level assemblies of the M. polymorpha U and V chromosomes allow a test of the theoretical predictions that UV sex chromosomes should evolve symmetrically and that both should show minimal degeneration. While the M. polymorpha U and V chromosomes have evolved symmetrically, the low gene numbers (<5% of each autosome) suggest that they have undergone substantial degeneration.19, 20, 21 The low gene density is consistent with this, as repetitive sequence accumulation can occur if there are few genes.19, 20, 21 No obvious synteny is detectable between the U and V chromosomes, and the evolutionary strata inferred are not contiguous, implying that each sex chromosome has undergone independent chromosomal rearrangements.

Our results demonstrate that a U-linked transcription factor gene, BPCU, is the sex determinant Feminizer in M. polymorpha. Unexpectedly, however, BPCU is also essential for reproductive induction in females, a function which, in males, is carried out by its gametolog on the V chromosome, BPCV. The conserved function in the reproductive transition of both BPCU and BPCV suggests that this presumably ancestral function was retained despite only BPCU evolved a sex-determining (feminizing) function promoting female organ development (Figure 6). That BPCU and BPCV have apparently equivalent, and required, roles in the reproductive transition explains why the originally isolated bpcU mutants were aneuploid, with BPCV enabling the transition to reproductive development in the bpcU [UV] background.

Figure 6

Genetic model of feminization by BPCU and reproductive induction by BPCU and BPCV

The female and male plants carry BPCU on the U chromosome and BPCV on the V chromosome, respectively. BPCU and BPCV share a function involved in inducing sexual reproduction. Following the transition to reproductive development, the female thallus initiates expression of the transcription factor FGMYB that promotes development of female reproductive organs. In gametophytes that do not carry BPCU, expression of the antisense long non-coding (lnc) RNA SUF prevents FGMYB expression, leading to male reproductive development. The U-linked Feminizer, BPCU functions to repress antisense SUF expression in females, allowing FGMYB expression and leading to female reproductive development.

Including BPCU and BPCV, seven ancient gametolog pairs are present on the M. polymorpha sex chromosomes, implying that large non-recombining regions were already established in the ancestral liverwort, which is estimated to have existed in the Mid-Silurian, ca. 430 mya. This age precedes any other known sex chromosome system. The oldest XY and ZW angiosperm sex chromosomes evolved in the Cenozoic (30–40 mya), the U/V sex chromosome systems in green and brown algae arose in the Cretaceous (75 and 110 mya, respectively), and the mammalian X and Y chromosomes date to the Jurassic (180 mya). The oldest previously described UV sex chromosomes of the moss C. purpureus evolved in the Carboniferous (300 mya), and even the oldest mating-type system known thus far, pheromone receptor alleles in Microbotryum fungal species, has been dated to the Devonian (370 mya).

Relative to BPCU and BPCV, each of the other six ancient pairs of gametologs has experienced recombination either between the sex chromosomes, or between a sex chromosome and an autosome, in various liverwort lineages (Figure 5). These data are consistent with the possibility that BPCU acted as the Feminizer in the ancestral liverwort. Given that FGMYB orthologs promote female gametophyte development in both M. polymorpha and Arabidopsis thaliana, it probably promoted female gametophyte development in the common ancestor of land plants. In extant M. polymorpha, FGMYB acts as a binary switch, resulting in female differentiation when it is active, with differentiation as males as a “default” when FGMYB is repressed. The genetic pathway(s) promoting maleness are as yet unknown. It is thus of interest whether the BPCU and FGMYB/SUF regulatory module is conserved in other liverworts and other bryophyte lineages with sex chromosomes (Figure 6). While there is no solid evidence of liverwort gametologs diverging prior to extant liverwort divergence, this does not preclude an earlier origin of dioicy, a condition possibly present in the ancestral land plant.

One plausible scenario is that the ancestral BPC gene functioned in the transition from vegetative to reproductive growth, and in the ancestral liverwort, an allele of BPC acquired a role in specifying female sex organ development, perhaps via its regulation of FGMYB. This allele evolved into BPCU, while the alternative allele (now the gametolog), BPCV, retained the ancestral function in the reproductive transition but did not acquire a sex-determining function. This can potentially explain the fact that, in addition to its role in sex determination, BPCU functions in the initiation of reproductive development, like its gametolog BPCV. This scenario is distinct from that described for other sex-determining systems. For example, in therian mammals, the male sex determinant SRY does not share the non-sex determination function of its X-linked gametolog, SOX3, which was presumably retained from the common ancestral gene. Likewise, the sex determinant MID in the chlorophyte alga Volvox carteri evolved directly from an already male-specific gene within the ancestral mating-type locus. Our findings provide a distinct trajectory for the evolution of a sex determinant (the Feminizer) in another haploid system and further implies that a pair of gametologs that initially appear functionally equivalent may have been differentiated to serve females and males unequally.

STAR★Methods

Key resources table

Resource availability

Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Takayuki Kohchi (tkohchi@lif.kyoto-u.ac.jp).

Materials availability

Plasmids and transgenic lines generated during this study are available from the Lead Contact without restriction.

Experimental model and subject details

Plant materials and culture

The female accession, Takaragaike-2 (Tak-2), and the male accession, Takaragaike-1 (Tak-1), and their backcrossed female line BC3-38 of Marchantia polymorpha ssp. ruderalis,, were used in this study. The UV plant (Figure 1C) was isolated among progeny sporelings derived from a cross between BC3-38 and Tak-1. These plants were grown primarily on half-strength Gamborg’s B5 medium containing 1% agar, or on vermiculite supplemented by 2,000-fold-diluted HYPONeX (Hyponex Japan, Osaka, Japan) under continuous white light (50-60 μmol photons m-2 s-1; CCFL OPT-40C-N-L, Optorom) at 22°C. To induce reproductive growth, 10-day-old gemmalings or older thalli were transferred to white light conditions supplemented with far-red light (30 μmol photons m-2 s-1, VBL-TFL600-IR730∗, Valore, Kyoto, Japan).,

Method details

Plant transformation

Agrobacterium-mediated transformation of M. polymorpha was performed using standard protocols for sporelings for the initial screening or thalli for all the other lines. Transformants were selected on the media containing 100 μg ml-1 cefotaxime and either 10 μg ml-1 hygromycin or 0.5 μM chlorsulfuron according to the selection marker genes in the binary vectors. After co-cultivation with Agrobacterium, plants were washed with sterile water and plated on the selective media. Regenerating antibiotic-resistant transformants (T1) were transferred to new selective media and grown until gemmae (G1) were produced.

U chromosome assembly

Draft genome assembly was performed from PacBio reads using Redbean ver. 2.5 with default settings. The draft contigs were polished by NextPolish ver v1.3.1 using PacBio long reads and Illumina short reads, followed by a further round of polishing by Pilon ver. 1.23 with Illumina short reads. Separately, a draft genome was also reconstructed using NextDenovo v2.2-beta.0 (https://github.com/Nextomics/NextDenovo), NextPolish, and Pilon, which was then merged into the draft genome obtained from Redbean to compensate missing regions based on comparisons with the reference genomes obtained from MarpolBase (MpTak1v5.1 and JGI3.1, https://marchantia.info). The resultant genome sequences were subjected to Hi-C scaffolding with the 3D-DNA pipeline. The in situ Hi-C library preparation of 2-week-old Tak-2 thalli, sequencing, and data processing were performed as described previously. Finally, the order and orientation of the contigs within scaffolds were manually corrected by comparing them to the v3.1 genome sex-chromosomes. Gene annotation was lifted-over from the reference genomes by aligning transcript sequences using GMAP ver. 2019.09.12 s. Full-length transcriptome sequences (Iso-seq) and RNA-seq reads were obtained from INSDC (BioProject IDs: PRJDB8530 and PRJNA251267), mapped to the newly obtained Tak-2 genome, and used as evidence for gene annotations. All U chromosome-linked genes were manually curated and spurious gene models that did not show any expression level or share homology to known functional domains were eliminated. New gene models were also identified based on evidence from Iso-seq and RNA-seq. The updated gene IDs for the sex chromosome-linked genes are listed in Data S1A.

Genetic nomenclature

BPCU and BPCV correspond to the locus codes, MpUg00370 (this study) and MpVg00350, respectively. Genetic nomenclature in M. polymorpha is as outlined previously. Prefix “Mp” in gene symbols was omitted for simplicity since all genes used in this study were derived from M. polymorpha. The gene symbols, BPCU and BPCV, were renamed from BPC1 and BPC2, respectively, from the reported genome analysis to clearly designate their location on specific sex chromosomes. Locus codes and gene symbols are also available through MarpolBase (https://marchantia.info/).

Microscopy

The stalk of the developing sex organs (archegoniophores or antheridiophores) was cut off, placed on agar medium, and the entire receptacles were photographed with a stereomicroscope SZX16 (OLYMPUS, Tokyo, Japan). The inner gamete-containing organs were removed using a needle on a glass slide and photographed with a microscope Axiophot (ZEISS, Oberkochen, Germany).

RNA extraction and reverse-transcription

Total RNA was extracted from stages 1-2 sex organs or 10-day-old thalli with TRIzol reagent (Thermo Fisher Scientific) according to the manufacture’s instruction. DNA in the solution was digested using RQ1 RNase-Free DNase (Promega, Madison, USA) by incubation at 37°C for 15 min. The quantity and quality of extracted RNA were assessed using a NanoDrop 2000 spectrophotometer (Thermo Scientific). 500 ng of total RNA was reverse-transcribed to synthesize the first-strand cDNA using ReverTra Ace (Toyobo Life Science) with a 20-mer oligo dT primer. The synthesized cDNA was diluted 10 times with H2O and used for the experiments.

Plasmid construction for expressing BPCU and BPCV

Plasmids used in this study were constructed using Gateway cloning system (Thermo Fisher Scientific, Waltham, MA) for M. polymorpha,, as described below. Primers used for PCR amplification are listed in Table S1. The constructs were verified by Sanger sequencing.

To create the vectors for expressing BPCU and BPCV in plants, the gene sequences of BPCU (2-kb upstream to 3-kb downstream of the transcribed region) and BPCV (1.5-kb upstream to 2-kb downstream of the transcribed region) were amplified from Tak-2 and Tak-1 genomic DNAs, respectively, using KOD FX Neo DNA polymerase (Toyobo Life Science, Osaka, Japan), and cloned into pENTR/D-TOPO (Thermo Fisher Scientific). The resulting entry vectors (BPCU-gene-pENTRdTOPO and BPCV-gene-pENTRdTOPO) were recombined with pMpGWB101 and pMpGWB301 using Gateway LR clonase II Enzyme mix (Thermo Fisher Scientific) to create BPCU-gene-pMpGWB101, BPCU-gene-pMpGWB301, and BPCV-gene-pMpGWB301.

To create the vectors for the expression of BPCU and BPCV proteins in Escherichia coli, the coding sequences (CDSs) of BPCU and BPCV were amplified by RT-PCR from Tak-2 and Tak-1 thallus RNA using KOD FX Neo DNA polymerase and KOD One PCR Master Mix (Toyobo Life Science), respectively. These PCR products were cloned into pENTR D-TOPO (Thermo Fisher Scientific) to generate BPCU-cds-pENTRdTOPO and BPCV-cds-pENTRdTOPO. These entry vectors were recombined with pET-DEST42 using Gateway LR clonase II Enzyme mix to create BPCU-cds-pET-DEST42 and BPCV-cds-pET-DEST42. These vectors were introduced into E. coli BL21-RIL cells for the protein expression described below.

Genome editing

To generate mutants for BPCU (bpcU-101, bpcU-102, and bpcU-3), BPCV (bpcV-1, bpcV-2 and bpcV-3) and FGMYB (fgmyb-31 and fgmyb-41), oligonucleotides listed in Table S1 were used as guide RNA (gRNA) sequences for CRISPR/Cas9-mediated genome editing with pMpGE_En03 (Addgene 71353), and pMpGE010 (Addgene 71536) or pMpGE011 (Addgene 71537). To generate large deletion mutants for BPCU (bpcU-1 and bpcU-2), four pairs of oligonucleotides (Table S1) were cloned as gRNA sequences into pMpGE_En04, pBC-GE12, pBC-GE23, and pBC-GE34, and further combined to generate an entry clone having four gRNAs. The resultant gRNA sequences were subcloned into pMpGE017 by the Gateway LR reaction as described previously. For the initial Feminizer screening, the gene-edited mutants were generated via transformation of F1 spores from a cross between Tak-2 and Tak-1 with the above constructs. bpcU mutants (bpcU-1 and bpcU-2) and bpcV mutants (bpcV-1, bpcV-2, and bpcV-3) were generated in Tak-2 and Tak-1, respectively, and bpcU-3 in the BC3-38 background, via thallus transformation.

Genotyping and sex diagnosis

Genomic DNA was extracted from a small piece of thallus (3 mm x 3 mm) crushed with micropestle in 100 μL of genotyping buffer containing 100 mM Tris-HCl, pH 9.5, 1M KCl, and 100 mM EDTA as described previously, and used for PCR using KOD FX Neo DNA polymerase. For gene-edited mutant identification, gRNA-targeted regions were amplified from the genomic DNA of transformants, and sequenced using the primer pairs listed in Table S1. Sex diagnosis was performed using the primer sets for rhf73 and rbm27 that are specific to sex chromosomes as described previously.

Generation of transgenic lines expressing genomic BPCU and genomic BPCV

To create the transgenic plants carrying the BPCU and BPCV transgenes, thallus transformation was used to introduce BPCU gene pMpGWB101 and BPCU gene pMpGWB301 into Tak-1 and bpcU-1, respectively, and also BPCV gene pMpGWB301 into bpcU-1. To confirm genomic integration of the transgenes, G1 gemmae were used for PCR amplification using KOD One PCR Master Mix and primers listed in Table S1, and progeny of G1 plants were used for the experiments. More than 10 lines were independently isolated for each construct, and the data from multiple representative lines were shown.

RT-PCR

RT-PCR was performed in 10 μL reaction solution, containing KOD One PCR Master Mix, 1 μL diluted cDNA and 0.3 μL each of 10 μM forward and reverse primers listed in Table S1 PRM (Mp3g14390) and LC7 (Mp6g01560) were amplified under the following cycle conditions: initial incubation at 94°C for 2 min, followed by 37 cycles of 98°C for 10 s, 58°C for 5 s and then 68°C for 1 s. Mp5g13880 was amplified under the following cycle conditions: initial incubation at 94°C for 2 min, followed by 35 cycles of 98°C for 10 s, 60°C for 5 s and then 68°C for 1 s. Mp7g02820 was amplified under the following cycle conditions: initial incubation at 94°C for 2 min, followed by 35 cycles of 98°C for 10 s, 58°C for 5 s and then 68°C for 1 s. EF1 (Mp3g23400) was amplified under the following cycle conditions: initial incubation at 94°C for 2 min, followed by 31 cycles of 98°C for 10 s, 58°C for 5 s and then 68°C for 1 s.

Quantitative real-time RT-PCR

Quantitative real-time RT-PCR (qRT-PCR) was performed in 25 μL reaction solution, containing homemade Taq DNA polymerase, SYBR Green I Nucleic Acid Gel Stain (Lonza, Basel, Switzerland), 5 μL diluted cDNA, and 1 μL each of 10 μM forward and reverse primers listed in Table S1. Reactions were carried out with a CFX96 real-time PCR detection system (Bio-Rad Laboratories, Hercules, CA), under the following cycle conditions: initial incubation at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, and 60°C for 30 s. Three technical and three biological replicates were measured for each reaction. EF1 (Mp3g23400) was used as internal control. Relative expression levels were calculated by the 2−ΔΔCT method.

DPI-ELISA with BPCU and selected DNA-probes

E. coli BL21-RIL cells harboring BPCU-cds-pET-DEST42 or BPCV-cds-pET-DEST42 were used for the protein expression. Crude protein extraction was performed with HEPES buffer (4 mM HEPES pH 7.5, 100 mM KCl, 8% (v/v) glycerol, 0.2% biotin-free BSA, 5 mM dithiothreithol (DTT), supplemented with 1x cOmplete Protease Inhibitor without EDTA (Roche, Basel, Switzerland). Protein extracts from non-transformed E. coli BL21-RIL cells served as controls. Semiquantitative DPI-ELISA experiments were performed with 5 μg of crude protein extracts (in 30 μl buffered solution) and biotinylated double-stranded (ds) DNA-probes. Promoter probes were PCR amplified by using biotinylated (sense) and non-biotinylated (antisense) primer combinations. Residual biotinylated primers were removed with a NucleoSpin Gel and PCR Clean-up Kit (Macherey & Nagel, Duren, Germany). DPI-ELISA with purified promoter probes was performed according to the general protocol. The use of streptavidin-coated ELISA plates (5 pmol/well) allows for the semiquantitative and comparable readout within the same protein extract., Pre-blocked clear flat-bottom Streptavidin Coated Plates (96-wells) (Pierce Thermo Fischer Scientific) were used with Safire II multiplate reader (Tecan, Mannedorf, Switzerland) for photometric measurements. Detection of 6xHis-epitope tagged BPCU was carried out with either Penta·His HRP Conjugate Kit (QIAGEN, Hilden, Germany) or HRP Anti-6X His tag antibody (Abcam, Cambridge, UK) and ortho-phenylenediamine (OPD) as a substrate. The HRP reaction was stopped in all wells by addition of sulfuric acid (2 N H2SO4) as soon as a visible coloration appeared. Average signal intensities and standard deviation were computed from two technical replicates on the same plate and from two independent protein extracts. Relative signal intensities were calculated relative to the highest photometric absorbance of a reference protein-DNA combination.

Flow cytometry analysis

Flow cytometry was performed according as previously described with slight modifications. One hundred mg of 9-day-old M. polymorpha thalli and 37-day-old A. thaliana (Col-0) rosette leaves were chopped in 1.5 mL of Galbraith Buffer (45 mM MgCl2, 30 mM Trisodium citrate dihydrate, 20 mM MOPS, pH 7.0) supplemented with 0.1% v/v Triton X-100 using a single edge blade (AccuTec Blades, Verona, USA) for 5 min on ice, and consecutively passed through 40 μm (Corning, NY, USA), 20 μm and 10 μm cell strainers (SYSMEX, Kobe, Japan). RNaseA (Sigma-Aldrich, St. Louis, USA) was added to a final concentration of 50 μg ml-1 and incubated on ice for 15 min, then propidium iodide (Nacalai Tesque, Kyoto, Japan) was added to the same final concentration and incubated on ice for 30 min. The prepared samples were analyzed by a BD Accuri C6 Flow Cytometer (BD Biosciences, Franklin Lakes, USA) for 300,000 events each.

Chromatin profiling analysis

The CUT&RUN (Cleavage Under Targets & Release Using Nuclease) assays were performed as previously described with slight modifications. DAPI stained 40,000 nuclei were sorted by BD FACSARIA III (BD Biosciences). 10 μL of BioMagPlus Concanavalin A bead slurry (Polysciences, #86057) was washed twice in 700 μL of binding buffer (20 mM HEPES pH 7.9, 10 mM KCl, 1 mM CaCl2 and 1 mM MnCl2) and resuspended in 10 μl of binding buffer, and added to the nuclear suspension. The mixture was thoroughly mixed on a rotator for 10 min at room temperature to allow binding of nuclei to the beads. The nuclei were pelleted with a magnet stand and resuspended with 50 μl of antibody buffer (20 mM HEPES pH 7.5, 150 mM NaCl, 0.5 mM spermidine, 1 × protease inhibitor cocktail EDTA free (Roche) and 2 mM EDTA pH8.0). 0.5 μg of antibody was thoroughly mixed and incubated overnight at 4°C. Nuclei were washed twice on a magnet stand with 1 mL of wash buffer and incubated with 700 ng μl-1 of Protein A–micrococcal nuclease fusion protein (pA-MN) in 50 μl of wash buffer at room temperature. After 10 minutes, the mixture was washed twice with 1 mL wash buffer to remove unbound pA-MN. Nuclei were resuspended in 150 μl of wash buffer and chilled on ice at 0°C. 3 μl of 100 mM CaCl2 was added to activate pA-MN with incubation at 0°C for 120 min. The reaction was stopped by the addition of 100 μl of 2 × STOP buffer (340 mM NaCl, 20 mM EDTA, 4 mM EGTA, 50 μg ml–1 RNase A, 50 μg ml–1 glycogen and MNase digested spike-in DNA (HEK293 cells)) and incubated in thermo-mixier at 37°C for 10 minutes at 500 rpm. The protein–DNA complex was released by centrifugation and then digested by 0.05 mg proteinase K and 0.1% SDS at 70 °C for 10 min. DNA was extracted by Phenol: Chloroform: Isoamyl-alcohol (25:24:1) with phase lock gel tubes. Library preparation of CUT&RUN was performed using a NEB Ultra II library preparation kit. 150-500 bp of DNA was selected by SPRI beads. DNA amount and fragment size were analyzed by Nanodrop 3300 and Fragment analyzer respectively. Libraries were pooled at similar molar amounts and sequenced using the Nextseq 500 platform. Paired-end sequencing was performed (read length, 42 bp × 2; index, 6 bp).

Analysis of sex chromosome evolution

Predicted orthologs of M. polymorpha gametolog were collected from available liverwort genome and transcriptome sequence data.82, 83, 84, 85, 86 Since the available liverwort transcriptomes were not generated to investigate sex chromosome evolution, we initially applied two filters. First, sequences derived from monooicous species were removed as monoicy is associated with a loss of sex chromosomes. Second, for analyses of gametologs other than BPC, we included only dioicous liverwort species for which only a single sex appeared to be sampled, because the transcriptome data appeared to be derived from a single sex (based on the presence of a single gametolog for each gene that we inferred to be in the oldest evolutionary stratum). An exception was the Haplomitriopsida species, where we included all identified sequences since these sequences were critical for identifying the oldest evolutionary stratum. For the phylogenetic analysis of BPC orthologs, all identified sequences from dioicous liverwort species were included. Moss and hornwort orthologs were also included to ascertain whether any non-recombining regions pre-date the origin of liverworts, and vascular plant (lycophyte and fern) orthologs were used to root the trees.,

For phylogenetic analyses of the gametologs (Figure 4; Data S2), all site types in coding regions (including non-synonymous as well as synonymous sites) were used. 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/). Poorly alignable sequences were removed, and alignments of nucleotides were employed in subsequent Bayesian analysis using MrBayes 3.2.1., The analysis was run for 1,000,000 generations (except for BPCU/MpUg00370 which was run for 10,000,000 generations), which was sufficient for convergence of the two simultaneous runs (split frequencies < 0.05). To allow a burn-in phase, the first 50% of the trees generated 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.

The proportions of synonymous substitutions per synonymous site (Ks values) were calculated with SNAP v2.1.1 (https://hcv.lanl.gov/content/sequence/SNAP/SNAP.html) as described previously., Ks values for the 19 pairs of M. polymorpha gametologs were previously reported. To relate these values to the liverwort phylogeny, we aligned sequences of eleven autosomal orthologs of liverwort genes across the same phylogenetic groups as used for studying gametologs, and estimated Ks values separately for each gene, with the standard deviation calculated as a measure of variation in values of different genes at each phylogenetic node (Data S1B).

Data analysis of CUT&RUN

CUT&RUN reads were preprocessed to filter out low-quality reads with Fastp v0.20.1, and mapped to the MpTakv6 genome presented in this paper using Bowtie2 v2.4.1 with the “--very --sensitive” flag. Reads with MapQ less than 10 were removed with Samtools v1.9.0, and duplicates or inserts less than 140 bp reads were removed with alignmentSieve in Deeptools v3.5.0. Deduplicated reads from biological replicates were merged. We called peaks for chromatin marks by SICER2 with default parameters. Bigwig files were generated using “bamcoverage” in Deeptools with “--extendReads --exactScaling --binSize 10 --normalizeUsing RPGC” flags. The effective genome size was calculated from MpTakv6 with faCount.

Transcriptome Mapping

RNA-seq reads were preprocessed to filter out low-quality reads with Fastp v0.20.1, and mapped to the MpTakv6 genome presented in this paper using STAR v2.7.3a with default parameters. Bigwig files were generated using “bamcoverage” in Deeptools with “--binSize 10 --normalizeUsing BPM” flags.

Quantification and statistical analysis

General statistical analyses

All statistical analyses were conducted in R (version 3.6.0). The specific tests, sample, size, and p value are displayed in the corresponding legends.

Gene expression analysis of BPCU and BPCV

Gene expression data from female sex organs: DRR050351, DRR050352, DRR050353, male sex organs: DRR050346, DRR050347, DRR050348, female thalli: DRR118943, DRR118944, DRR118945, male thalli: DRR118949, DRR118950, DRR118951, spores: SRR4450260, SRR4450261, SRR4450262, and sporelings: SRR4450254, SRR4450255, SRR4450256 were downloaded from Sequence Read Archive. Transcript abundances were quantified by Salmon package v1.40 with–validateMappings option. Transcripts Per Million values of transcripts were integrated into values corresponding to genome annotations by tximport package.

REAGENT or RESOURCE	SOURCE	IDENTIFIER
Antibodies
Rabbit polyclonal H3	Abcam	Cat# ab1791; RRID: AB_302613
Rabbit polyclonal H3K27me3	Millipore	Cat# 07-449; RRID: AB_310624
HRP Anti-6X His tag antibody	Abcam	Cat# ab1187, RRID: AB_298652
Bacterial and virus strains
Escherichia coli DH5a	Widely distributed	N/A
Agrobacterium tumefaciens GV2260	Deblaere et al.45	N/A
Escherichia coli BL21-CodonPlus-RIL	Agilent	Cat# 230240
Biological samples
HEK293 DNA	Danhua Jiang, Beijing, China	N/A
Chemicals, peptides, and recombinant proteins
Gamborg’s B5 salts	Gamborg et al.46	N/A
HYPONeX	Hyponex Japan	Cat# 4977517180036
hygromycin B	Nacalai Tesque	Cat# 07296-24
cefotaxime (CLAFORAN)	Sanofi	Cat# 6132409D1050
TRIzol reagent	ThermoFisher Scientific	Cat# 15596018
ReverTra Ace	Toyobo Life Science	Cat# TRT-101
RNase A	ThermoFisher Scientific	Cat# EN0531
RNase A	Sigma-Aldrich	Cat# R5125
RQ1 RNase-Free DNase	Promega	Cat# M610A
Proteinase K	Thermo Fisher Scientific	Cat# EO0491
Taq DNA polymerase	Pluthero47	N/A
Bio-Mag Plus Concanavalin A coated beads	polysciences	Cat# 86057
cOmplete Protease Inhibitor Cocktail	Roche	Cat# 11697498001
pA-MNase	Skene and Henikoff48	Henikoff lab batch #6 purified 11.01.2017
propidium iodide	Nacalai Tesque	Cat# 29037-76
ortho-phenylenediamine (OPD)	Merck-Sigma/Aldrich	Cat# P5412
sulfuric acid (2 N H2SO4)	Roth	Cat# 2609.1
SYBR Green Nucleic Acid Gel Stain	Lonza	Cat# 50513
Critical commercial assays
KOD FX Neo DNA polymerase	Toyobo Life Science	Cat# KFX-201
KOD One PCR Master Mix	Toyobo Life Science	Cat# KMM-101
pENTR/D-TOPO Cloning kit	Thermo Fisher Scientific	Cat# K240020
Gateway LR clonase II Enzyme mix	Thermo Fisher Scientific	Cat# 11791020
NucleoSpin Gel and PCR Clean-up Kit	Macherey & Nagel	Cat# 740609.50
Pre-blocked clear flat-bottom Streptavidin Coated Plates (96-wells)	Pierce Thermo Fischer Scientific	Cat# 15124
Penta·His HRP Conjugate Kit	QIAGEN	Cat# 34460
Deposited data
Marchantia polymorpha genome v3.1	Bowman et al.19	https://marchantia.info; SRA: SRR1800537
Marchantia polymorpha genome v5.1	Montgomery et al.20	https://marchantia.info; SRA: PRJNA553138
Marchantia polymorpha genome v6.1	This paper	https://marchantia.info; SRA: PRJDB11173
Iso-seq and RNA-seq for gene annotations	Montgomery et al.20	SRA: PRJDB8530 and PRJNA251267
Tak-1 thallus CUT&RUN	Montgomery et al.20	SRA: PRJNA553138
Tak-2 and bpcU thallus CUT&RUN	This paper	SRA: PRJNA757234
Tak-1, Tak-2 and bpcU thallus RNA-seq	This paper	SRA: PRJNA757569
RNA-seq for expression analyses in various tissues	Bowman et al.,19 Higo et al.,49 and Yamaoka et al.50	SRA: DRR050351-53; DRR050346-48; DRR118943-45; DRR118949-51; SRR4450260-62; SRR4450254-56
Experimental models: Organisms/strains
Marchantia polymorpha Tak-1	Bowman et al.19	N/A
Marchantia polymorpha Tak-2	Bowman et al.19	N/A
Marchantia polymorpha BC3-38	Yamaoka et al.50	N/A
Marchantia polymorpha UV plant	This paper	N/A
Marchantia polymorpha bpcU-101	This paper	N/A
Marchantia polymorpha bpcU-102	This paper	N/A
Marchantia polymorpha gBPCU#3	This paper	N/A
Marchantia polymorpha gBPCU#4	This paper	N/A
Marchantia polymorpha gBPCU#4fgmyb-31	This paper	N/A
Marchantia polymorpha gBPCU#4fgmyb-41	This paper	N/A
Marchantia polymorpha bpcU-1	This paper	N/A
Marchantia polymorpha bpcU-2	This paper	N/A
Marchantia polymorpha bpcU-3	This paper	N/A
Marchantia polymorpha bpcV-1	This paper	N/A
Marchantia polymorpha bpcV-2	This paper	N/A
Marchantia polymorpha bpcV-3	This paper	N/A
Marchantia polymorpha gBPCU#2bpcU-1	This paper	N/A
Marchantia polymorpha gBPCU#4bpcU-1	This paper	N/A
Marchantia polymorpha gBPCU#6bpcU-1	This paper	N/A
Marchantia polymorpha gBPCV#7bpcU-1	This paper	N/A
Marchantia polymorpha gBPCV#8bpcU-1	This paper	N/A
Marchantia polymorpha gBPCV#9bpcU-1	This paper	N/A
Oligonucleotides
See Table S1	This paper	N/A
Recombinant DNA
pENTR D-TOPO	Thermo Fisher Scientific	Cat# 45-0218
BPCU gene pENTRdTOPO	This paper	N/A
BPCV gene pENTRdTOPO	This paper	N/A
pMpGWB101	Ishizaki et al.51	GenBank: LC057443
pMpGWB301	Ishizaki et al.51	GenBank: LC057517
BPCU gene pMpGWB101	This paper	N/A
BPCU gene pMpGWB301	This paper	N/A
BPCV gene pMpGWB301	This paper	N/A
pMpGE_En03	Sugano et al.52	GenBank: LC090755
pMpGE_En04	Koide et al.53	N/A
pBC-GE12	Koide et al.53	N/A
pBC-GE23	Koide et al.53	N/A
pBC-GE34	Koide et al.53	N/A
BPCU.1 gRNA pMpGE_En03	This paper	N/A
BPCU gRNALDs pMpGE_En04	This paper	N/A
BPCV.1 gRNA pMpGE_En03	This paper	N/A
BPCV.2 gRNA pMpGE_En03	This paper	N/A
FGMYB gRNA3 pMpGE_En03	This paper	N/A
FGMYB gRNA4 pMpGE_En03	This paper	N/A
pMpGE010	Sugano et al.52	GenBank: LC090756
pMpGE011	Sugano et al.52	GenBank: LC090757
pMpGE017	Koide et al.53	N/A
BPCU.1 gRNA pMpGE010	This paper	N/A
BPCU gRNALDs pMpGE017	This paper	N/A
BPCV.1 gRNA pMpGE010	This paper	N/A
BPCV.2 gRNA pMpGE010	This paper	N/A
FGMYB gRNA3 pMpGE011	This paper	N/A
FGMYB gRNA4 pMpGE011	This paper	N/A
BPCU cds pENTRdTOPO	This paper	N/A
BPCV cds pENTRdTOPO	This paper	N/A
pET-DEST42	Thermo Fisher Scientific	Cat# 12276010
BPCU cds pET-DEST42	This paper	N/A
BPCV cds pET-DEST42	This paper	N/A
Software and algorithms
fastp v0.20.1	Chen et al.54	https://github.com/OpenGene/fastp
Salmon v0.14.1	Patro et al.55	https://github.com/COMBINE-lab/salmon
bowtie2 v2.4.1	Langmead and Salzberg56	https://github.com/BenLangmead/bowtie2
STAR v2.7.3a	Dobin et al.57	https://github.com/alexdobin/STAR
samtools v1.9.0	Danecek et al.58	https://github.com/samtools/
deeptools v3.5.0	Ramírez et al.59	https://github.com/deeptools/deepTools
faCount	UCSC Genome Browser	https://hgdownload.cse.ucsc.edu/admin/exe/
SICER2	Zang et al.30	https://zanglab.github.io/SICER2/
tximport	Soneson et al.60	https://bioconductor.org/packages/release/bioc/html/tximport.html
IGV v2.8.0	Thorvaldsdóttir et al.29	https://software.broadinstitute.org/software/igv/
Redbean v2.5	Ruan and Li61	https://github.com/ruanjue/wtdbg2
NextPolish v1.3.1	Hu et al.62	https://github.com/Nextomics/NextPolish
Pilon v1.23	Walker et al.63	https://github.com/broadinstitute/pilon
NextDenovo v2.2-beta.0	GrandOmics	https://github.com/Nextomics/NextDenovo
3D-DNA pipeline	Dudchenko et al.64	https://github.com/aidenlab/3d-dna
GMAP ver. 2019.09.12	Wu and Watanabe65	http://research-pub.gene.com/gmap/
Se-Al v2.0a11	Institute of Evolutionary Biology University of Edinburgh	http://tree.bio.ed.ac.uk/software/seal/
MrBayes 3.2.1	Huelsenbeck and Ronquist66 and Huelsenbeck et al.67	https://nbisweden.github.io/MrBayes/index.html
FigTree v1.4.0	Institute of Evolutionary Biology University of Edinburgh	http://tree.bio.ed.ac.uk/software/figtree/
SNAP v2.1.1	Korber68	https://hcv.lanl.gov/content/sequence/SNAP/SNAP.html
Other
BD FACSARIA III	BD Biosciences	https://www.bdbiosciences.com/en-us
BD Accuri C6 Flow Cytometer	BD Biosciences	https://www.bdbiosciences.com/en-us
CFX96 real-time PCR detection system	Bio-Rad Laboratories	https://www.bio-rad.com
NanoDrop 2000 spectrophotometer	Thermo Scientific	https://www.thermofisher.com/us/en/home.html