Deep evolutionary origin of gamete-directed zygote activation by KNOX/BELL transcription factors in green plants

  1. Tetsuya Hisanaga
  2. Shota Fujimoto
  3. Yihui Cui
  4. Katsutoshi Sato
  5. Ryosuke Sano
  6. Shohei Yamaoka
  7. Takayuki Kohchi
  8. Frédéric Berger
  9. Keiji Nakajima  Is a corresponding author
  1. Graduate School of Science and Technology, Nara Institute of Science and Technology, Japan
  2. Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter, Austria
  3. Graduate School of Biostudies, Kyoto University, Japan
7 figures, 1 table and 3 additional files

Figures

Figure 1 with 1 supplement
Comparative transcriptome analysis of wild-type and Mprkd archegonia and identification of MpKNOX1 as an egg-specific gene.

(A) Schematic illustration of RNA-seq analysis comparing the archegonia transcriptomes from wild-type females and egg-deficient Mprkd mutant females. About 4000 and 2000 archegonia were collected from wild-type and each of the two Mprkd mutant lines, and randomly allocated into four and two replicates, respectively. (B) Comparison of the domain arrangements of MpKNOX1 vs. representative class I KNOX proteins from Arabidopsis thaliana (AtSTM; GenBank accession number AEE33958.1), Zea mays (ZmKN1; AAP21616.1), Ceratopteris richardii (CerKNOX1, BAB18582.1), Physcomitrium patens (PpMKN2; AAK61308.2), and Chlamydomonas reinhardtii (ChrGSM1; ABJ15867.1).

Figure 1—figure supplement 1
Comparative transcriptome analysis of wild-type and Mprkd archegonia and identification of MpKNOX1 as an egg-specific gene in M. polymorpha.

(A) M-A plot of differential expression analysis between wild-type and Mprkd archegonia. log2 differences between wild-type and Mprkd (Y-axis) are plotted against log2 average between wild-type and Mprkd (X-axis). Differentially expressed genes (DEGs) are plotted in magenta. (B) Captured Genome Browser image showing the expression levels along the MpKNOX1 locus in the indicated tissue types (archegoniophore: SRX301555; antheridiophore: SRX301553; sporophyte: SRX301556; sporeling: SRX301559; thallus: SRX301557). MpKNOX1 is specifically expressed in archegoniophores. The unit for Y axis is a read count normalized against 1× genome coverage.

MpKNOX1 is specifically expressed in egg cells.

(A) RT-PCR analysis of MpKNOX1. Lanes are labeled as follows: M: size markers; Th ♂: male thalli; Th ♀: female thalli; Ant: antheridiophores; Arc: archegoniophores; Spo: sporophytes of 3-week-old plants. Constitutively expressed MpEF1α was used as a control. Red arrow indicates the expected size of PCR products from spliced MpKNOX1 mRNA. Bands at the top of the gel likely correspond to unspliced MpKNOX1 transcripts. Shown is a representative result from the experiments using three independently collected plant samples each with two technical replicates (two PCRs from each cDNA pool). See Figure 4—figure supplement 1A for the primer position. (B, C) Expression of the MpKNOX1 transcriptional reporter. Magenta: chlorophyll autofluorescence; green: GFP fluorescence. Lower panels are merged photographs of fluorescence and bright-field images. Bars, 10 μm.

Figure 3 with 2 supplements
Time-course observation of subcellular dynamics during M. polymorpha fertilization.

(A) Schematic representation of sexual reproduction in M. polymorpha. Female plants develop umbrella-shaped sexual branches termed archegoniophores that form egg-containing archegonia. Male plants develop disc-shaped sexual branches termed antheridiophores that form antheridia, which produce numerous motile sperm cells. Upon soaking antheridiophores in water, the sperm cells are released from the antheridia and swim to egg cells in the archegonia. After fertilization, each zygote undergoes embryogenesis by dividing and differentiating into a sporophyte body consisting of a capsule containing haploid spores and a short supportive stalk called the seta. Black bars, 1 mm. White bars, 10 µm. (B) Illustration of the in vitro fertilization method used in this study. Excised archegoniophores and antheridiophores were co-cultured in water for 1 hr to allow fertilization to take place. The archegoniophores were transferred to a fresh tube containing water for further culturing. The tube lids were left open to allow gas exchange to occur. Archegoniophores containing sporophytes were cultured for up to 2 weeks. (C) A DAPI-stained zygote after 1 hr of co-culture. Most zygotes contained sperm nuclei at this time (green arrowhead). Bars, 10 µm. (D–I) DAPI-stained zygotes and sporophytes at the indicated days after fertilization (DAF). Male pronuclei (green arrowheads) were visible at 1–3 DAF (D, E) in wild-type fertilized eggs. In most zygotes, karyogamy was completed, and cells were cleaved at 4 DAF (F). Sporophyte cells continued to divide at 5–14 DAF (G–I), as visualized by the presence of multiple nuclei (blue arrowheads; not all nuclei are labeled in H and I). Ca: calyptra. Bars, 10 µm.

Figure 3—figure supplement 1
Cell wall regeneration during zygote development.

(A–F) Confocal images of archegonia containing an egg cell, zygotes, or embryos. Cell walls stained with SCRI Renaissance 2200 were detected only after embryogenesis (arrowheads), but not in mature eggs or zygotes at 2 days after fertilization or earlier. Bars, 10 µm.

Figure 3—figure supplement 2
Cellular dynamics of zygotes and embryos generated by in planta crossing.

(A–H) DAPI-stained zygotes and embryos derived from in planta crosses at the indicated stages, showing that the timing of subcellular dynamics was equivalent to that generated by the in vitro fertilization method shown in Figure 3. At 4 days after fertilization, most zygotes completed karyogamy (D), whereas some showed an M-phase image (E) or completed the first division (F). Green arrowheads: male pronuclei; blue arrowheads: sporophytic nuclei. Bars, 10 µm.

Figure 4 with 6 supplements
Maternally inherited MpKNOX1 is required for nuclear fusion.

(A, E) Bright-field images of wild-type (WT; A) and Mpknox1-1ge (E) archegonia. (B–D and F–H). 1 day after fertilization (DAF) (B, F), 5 DAF (C, G), and 7 DAF (D, H) zygotes from a cross between WT female and male plants (B–D), and a cross between Mpknox1-1ge female and WT male plants (F–H), indicating that maternal MpKNOX1 is dispensable for fertilization (B, F) but is required for embryogenesis (C, D, G, H). (I–K) Bar graphs showing the ratios of zygotes containing male pronuclei (mp+) vs. those not containing male pronuclei (mp-) at 1 DAF (I), and developed vs. arrested zygotes at 5 DAF (J) and 7 DAF (K) derived from the crosses of either WT or Mpknox1 females with WT males. Number of observed zygotes is shown above each bar. (L–Q) Egg cells of MpSUN-GFP marker lines in the WT (L) or Mpknox1-2ge (O) female background were crossed with WT males. At 3 DAF, male and female pronuclei were in contact with each other in both WT (M) and Mpknox1 (P) eggs. At 5 DAF, zygotes derived from a WT egg started to divide (N), while those from an Mpknox1 female (Q) were arrested without nuclear membrane fusion. Green arrowhead: male pronucleus; orange arrowhead: female pronucleus; blue arrowhead: embryo nucleus. Bars, 10 µm.

Figure 4—figure supplement 1
Generation of loss-of-function mutant lines by CRISPR/Cas9.

Gene organization and the locations of CRISPR/Cas9-introduced mutations in MpKNOX1 (A), MpBELL3 (B), and MpBELL4 (C) are indicated by the following symbols: gray line: 5′- and 3′-flanking sequences; gray box: 5′- and 3′-UTR; black box: coding regions; arrowhead: mutation positions; black arrow: transcriptional direction; dotted line: splice pattern; gray arrow: primers used in RT-PCR shown in Figures 2A and 5A. Sequence alignments of wild-type and mutant alleles are shown below each gene model. Mismatched nucleotides and gaps are shown in red.

Figure 4—figure supplement 2
MpKNOX1 and MpBELL3/4 are dispensable for gametophyte development.

Three samples of vegetative thalli and gametangiophores from wild-type (WT), Mpknox1, and Mpbell3/4 mutant lines are shown. Bars, 5 mm.

Figure 4—figure supplement 3
Expression of MpKNOX1-GFP by the MpKNOX1 promoter complements the karyogamy defects of Mpknox1 mutants.

(A–C) 5 days after fertilization (DAF) zygotes obtained by crossing the indicated female lines with wild-type males. Green and blue arrowheads indicate a male pronucleus and embryo nuclei, respectively. Bars, 10 µm.

Figure 4—figure supplement 4
Sporangium and spore formation in wild-type and Mpknox1 plants.

(A, B) Archegoniophores of wild-type (A) and Mpknox1-1ge (B) female plants 4 weeks after crossing with wild-type males. Blue arrowheads indicate sporangia harboring expanded capsules. (C–F) Capsule (C, D) and spores (E, F) in wild-type (C, E) and Mpknox1-1ge (D, F) female plants crossed with wild-type males. (G) A box plot showing the number of sporangia per archegoniophore from indicated crosses. Five gametangiophores were analyzed for each genotype. *p=0.355* and **0.128, Student’s t-test. Bars, 5 mm (A, B), 0.5 mm (C, D), and 50 µm (E, F).

Figure 4—figure supplement 5
MpKNOX1 is dispensable for sperm differentiation and embryogenesis.

(A–C) DAPI-stained sperm cells from wild-type (WT; A) and Mpknox1 male plants (B, C) show indistinguishable morphology. (D–F) 5 days after fertilization (DAF) sporophytes in WT female plants crossed with WT (D) or Mpknox1males of distinct alleles (E, F). Blue arrowheads indicate nuclei of embryo cells. (G) Proportion of developed and arrested zygotes in WT female plants crossed with WT or Mpknox1 males of distinct alleles. Numbers of observed zygotes are indicated above each bar. Bars, 5 µm (A–C) and 10 µm (D–F).

Figure 4—figure supplement 6
Paternally inherited Mpknox1 does not enhance the zygote arrest phenotype caused by maternally inherited Mpknox1.

(A–D) 7 days after fertilization (DAF) sporophytes in Mpknox1-2ge female plants crossed with wild-type (A, B) or Mpknox1-2ge male plants (C, D). Shown are representative images of zygotes arrested at karyogamy (A, C), developing sporophytes (B), and developing sporophytes under mitosis (D). Green and blue arrowheads indicate a male pronucleus and embryo nuclei, respectively. Bars, 10 µm. (E) A bar graph showing the ratios of developed vs. arrested zygotes derived from the crosses of Mpknox1 eggs with either wild-type or Mpknox1 sperm. Numbers of observed zygotes are indicated above each bar. *p=0.681 and **1, Fisher’s exact test.

Figure 5 with 3 supplements
Both paternally and maternally inherited MpBELL genes are required for karyogamy.

(A) RT-PCR analysis indicating that MpBELL3 and MpBELL4 are specifically expressed in antheridiophores. The lanes are labeled as in Figure 2A. Shown is a representative result from the experiments using three independently collected samples each with two technical replicates. See Figure 4—figure supplement 1A for the primer position. (B–K) Zygotes at 1 day after fertilization (DAF; B, D, G, J) and 5 DAF (C, E, F, H, I, K) from crosses between a wild-type female and wild-type male (B, C) or Mpbell3/4-2ge male (D–F) and a Mpbell3/4-2ge female and a wild-type male (G–I) or Mpbell3/4-2ge male (J, K). The presence of male pronuclei (green arrowheads) in zygotes of all genotypes at 1 DAF (B, D, G, J) indicates that MpBELL3 and MpBELL4 are dispensable for plasmogamy. Note that zygotes produced from both or one Mpbell3/4 parent exhibit a variable degree of karyogamy arrest as visualized by the retention of male pronuclei (green arrowheads) among those starting embryonic division (nuclei labeled with blue arrowheads). Bars, 10 µm. (L) A bar graph showing the ratios of zygotes containing male pronuclei (mp+) vs. those not containing male pronuclei (mp-) from indicated crosses at 1 DAF. Numbers of observed zygotes are shown above each bar. (M) A bar graph showing the ratios of developed vs. arrested zygotes at 5 DAF from the indicated crosses. Numbers of observed zygotes are shown above each bar.

Figure 5—figure supplement 1
MpBELL3 and MpBELL4 are preferentially expressed in antheridiophores.

(A, B) Bar graphs showing the expression levels of MpBELL genes in the indicated organs, constructed from publicly available transcriptome data (Bowman et al., 2017) (A) and the RNA-seq data obtained in this study (B). (C) Snapshots of Genome Browser views of the MpBELL3 and MpBELL4 loci displaying their preferential expression in antheridiophores. The unit for Y axis is a read count normalized against 1× genome coverage.

Figure 5—figure supplement 2
MpBELL3 and MpBELL4 are dispensable for gamete differentiation.

( A–C) DAPI-stained sperm cells from wild-type (A), Mpbell3/4-1ge (B), and Mpbell3/4-2ge (C) male plants. (D–F) Differential Interference contrast (DIC) images of archegonia in wild-type (D), Mpbell3/4-1ge (E), and Mpbell3/4-2ge (F) female plants. Bars, 5 µm (A–C) and 10 µm (D–F).

Figure 5—figure supplement 3
MpBELL3 and MpBELL4 are dispensable for gamete differentiation.

(A–F) 7 days after fertilization (DAF) sporophytes from the indicated crosses. Note that zygotes inheriting Mpbell3/4 from either or both parent(s) were arrested at karyogamy with variable ratios. Arrested zygotes were recognized with retention of male pronuclei (green arrowheads), whereas developing zygotes with divided nuclei (blue arrowheads). Bars, 10 µm. (G) A bar graph showing the ratios of developed vs. arrested zygotes at 7 DAF from the indicated crosses. Numbers of observed zygotes are shown above each bar. (H) A box plot showing the number of sporangia per archegoniophore from indicated crosses. Five gametangiophores were analyzed for each genotype.

MpKNOX1 transiently localizes to female and male pronuclei prior to karyogamy.

(A–F) GFP (upper panels) and DAPI (lower panels) signals from gMpKNOX1-GFP/Mpknox1-1ge eggs (A) and zygotes obtained by crossing a gMpKNOX1-GFP/Mpknox1-1ge female with a wild-type (B, C) or Mpbell3/4-2ge (D–F) male. Note that before fertilization MpKNOX1-GFP was exclusively localized to the cytosol (A). At 12 HAF, MpKNOX1-GFP signals were enriched in female (orange arrowhead) and male (green arrowhead) pronuclei in the wild-type background (B). In the absence of paternally inherited MpBELL3 and MpBELL4, about a half of zygotes did not exhibit nuclear-enriched MpKNOX1-GFP signals at 12 HAF (D, E). At 24 HAF, weak and exclusively cytosolic GFP signals were detected in both genotypes (C, F). Bars, 10 µm. (G) A bar graph showing the ratios of zygotes with and without pronuclei-enriched MpKNOX1-GFP signals from the indicated crosses. Numbers of observed zygotes are shown above each bar.

Functions and expression patterns of KNOX/BELL transcription factors in green plants.

(A) Expression patterns of KNOX and BELL proteins and their predicted role in zygote activation in M. polymorpha. ‘K’ and ‘B’ represent MpKNOX1 and MpBELL protein subunits, respectively. Orange, green, and blue circles represent female pronuclei, male pronuclei, and nuclei of embryo cells, respectively. Purple curved arrows indicate auto-amplification of BELL levels by KNOX/BELL-mediated transcriptional control. X indicates unknown karyogamy-promoting factor(s) whose expression and/or functions are activated by KNOX/BELL-mediated transcription. (B) Predicted evolutionary trajectory of KNOX/BELL expression patterns along the green plant lineages. Orange/green circles and blue rectangles represent gametes and sporophyte bodies, respectively. ‘K‘ and ‘B’ indicate the expression of KNOX and BELL proteins, respectively. Purple stars indicate the predicted positions at which the functional transition of KNOX/BELL from zygote activation to sporophyte morphogenesis occurred. Also indicated are the expression patterns of evolutionarily conserved regulators of sexual differentiation; a male-determinant factor MID of volvocine algae, and female- and male-differentiation factors, FGMYB and DUO1, respectively, of land plants. Note that the KNOX/BELL expression patterns in ancestral plants at the bottom of the tree are an inference.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (Marchantia polymorpha)MpKNOX1MalpolBaseMp5g01600
Gene (Marchantia polymorpha)MpBELL1MalpolBaseMp8g18310
Gene (Marchantia polymorpha)MpBELL2MalpolBaseMp4g09650
Gene (Marchantia polymorpha)MpBELL3MalpolBaseMp8g02970
Gene (Marchantia polymorpha)MpBELL4MalpolBaseMp8g07680
Gene (Marchantia polymorpha)MpBELL5MalpolBaseMp5g11060
Gene (Marchantia polymorpha)MpSUNMalpolBaseMp5g02400
Gene (Marchantia polymorpha)MpECproMalpolBaseMp5g18000
Strain, strain background (Marchantia polymorpha, male)Tak-1DOI:10.1016/j.cell.2017.09.030Male wild-type strain of Marchantia
Strain, strain background (Marchantia polymorpha, female)Tak-2DOI:10.1016/j.cell.2017.09.030Female wild-type strain of Marchantia
Recombinant DNA reagentpMpGE_En03AddgeneRRID:Addgene_71535Entry plasmid for gRNA
Recombinant DNA reagentpMpGE010AddgeneRRID:Addgene_71536Destination plasmid for CRISPR/Cas9
Recombinant DNA reagentpMpGE_En04DOI:10.15252/embj.2018100240Entry plasmid for gRNA
Recombinant DNA reagentpBC-GE12DOI:10.15252/embj.2018100240Entry plasmid for gRNA
Recombinant DNA reagentpBC-GE23DOI:10.15252/embj.2018100240Entry plasmid for gRNA
Recombinant DNA reagentpBC-GE34DOI:10.15252/embj.2018100240Entry plasmid for gRNA
Recombinant DNA reagentpMpGE017DOI:10.15252/embj.2018100240Destination plasmid for Cas9 nickase
Recombinant DNA reagentpMpGE010_MpKNOX1geThis paperPlasmid to create Mpknox1 mutants
Recombinant DNA reagentpMpGE017_MpBELL4ge-MpBELL3geThis paperPlasmid to create Mpbell3/4 mutants
Recombinant DNA reagentpMpSL30DOI:10.15252/embj.2018100240Destination plasmid
Recombinant DNA reagentpMpSL30_MpKNOX1pro-H2B-GFP-3’MpKNOX1This paperMpKNOX1 promoter reporter
Recombinant DNA reagentpMpSL30_gMpKNOX1-GFPThis paperMpKNOX1 complementation
Recombinant DNA reagentpMpSL30_ ECpro-MpSUN-GFPThis paperEgg cell-specific nuclear envelop marker
Genetic reagent (Marchantia polymorpha)Mpknox1-1geThis paperMpknox1 CRISPR mutant obtained from sporeling transformation
Genetic reagent (Marchantia polymorpha)Mpknox1-2geThis paperMpknox1 CRISPR mutant obtained from sporeling transformation
Genetic reagent (Marchantia polymorpha)Mpknox1-3geThis paperMpknox1 CRISPR mutant obtained from sporeling transformation
Genetic reagent (Marchantia polymorpha)Mpbell3/4-1geThis paperMpbell3/4 CRISPR mutant obtained from sporeling transformation
Genetic reagent (Marchantia polymorpha)Mpbell3/4-2geThis paperMpbell3/4 CRISPR mutant obtained from sporeling transformation
Genetic reagent (Marchantia polymorpha)MpKNOX1pro:H2B-GFP/WT ♀This paperpMpSL30_MpKNOX1pro-H2B-GFP-3′MpKNOX1 was transformed into Tak-2
Genetic reagent (Marchantia polymorpha)gMpKNOX1-GFP/Mpknox1-1geThis paperpMpSL30_gMpKNOX1-GFP was transformed into Mpknox1-1ge
Genetic reagent (Marchantia polymorpha)ECpro:MpSUN-GFPThis paperpMpSL30_ ECpro-MpSUN-GFP was transformed into Tak-2 or Mpknox1-2ge
Genetic reagent (Marchantia polymorpha)Mprkd-1DOI:10.1016/j.cub.2016.05.013
Genetic reagent (Marchantia polymorpha)Mprkd-3DOI:10.1016/j.cub.2016.05.013
Commercial assay or kitRNeasy Plant Mini KitQiagen74904
Commercial assay or kitTruSeq RNA Sample Prep Kit v2IlluminaRS-122-2001
Commercial assay or kitGateway LR Clonase II Enzyme mixThermo Fisher Scientific11791020
Software, algorithmTopHat ver. 2.0.14doi:10.1093/bioinformatics/btp120RRID:SCR_013035
Software, algorithmDESeq2doi:10.1186/s13059-014-0550-8RRID:SCR_015687
Software, algorithmTCCDOI:10.1186/1471-2105-14-219RRID:SCR_001779

Additional files

Supplementary file 1

List of genes with changed mRNA levels in Mprkd archegonia.

Genes with reduced and increased expression levels in Mprkd archegonia are listed in sheet 1 and sheet 2, respectively.

https://cdn.elifesciences.org/articles/57090/elife-57090-supp1-v2.xlsx
Supplementary file 2

Primers used in this study.

Lowercase and uppercase letters indicate synthetic adaptor and target DNA sequences, respectively.

https://cdn.elifesciences.org/articles/57090/elife-57090-supp2-v2.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/57090/elife-57090-transrepform1-v2.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Tetsuya Hisanaga
  2. Shota Fujimoto
  3. Yihui Cui
  4. Katsutoshi Sato
  5. Ryosuke Sano
  6. Shohei Yamaoka
  7. Takayuki Kohchi
  8. Frédéric Berger
  9. Keiji Nakajima
(2021)
Deep evolutionary origin of gamete-directed zygote activation by KNOX/BELL transcription factors in green plants
eLife 10:e57090.
https://doi.org/10.7554/eLife.57090