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

Phosphate (Pi) is a pivotal nutrient that constraints plant development and productivity in natural ecosystems. Land colonization by plants, more than 470 million years ago, evolved adaptive mechanisms to conquer Pi-scarce environments. However, little is known about the molecular basis underlying such adaptations at early branches of plant phylogeny. To shed light on how early divergent plants respond to Pi limitation, we analyzed the morpho-physiological and transcriptional dynamics of Marchantia polymorpha upon Pi starvation. Our phylogenomic analysis highlights some gene networks present since the Chlorophytes and others established in the Streptophytes (eg. PHR1-SPX1 and STOP1-ALMT1, respectively). At the morpho-physiological level, the response is characterized by the induction of phosphatase activity, media acidification, accumulation of auronidins, reduction of internal Pi concentration and developmental modifications of rhizoids. The transcriptional response involves the induction of MpPHR1, Pi transporters, lipid turnover enzymes and MpMYB14, an essential transcription factor for auronidins biosynthesis. MpSTOP2 up-regulation correlates with expression changes in genes related to organic acid biosynthesis and transport, suggesting preference for citrate exudation. Analysis of MpPHR1 binding sequences (P1BS) shows enrichment of this cis regulatory element in differentially expressed genes. Our study unravels the strategies, at diverse levels of organization, exerted by M. polymorpha to cope with low Pi availability.Significance StatementThis study unravels the transcriptional and morphophysiological mechanisms executed by the non-vascular, and rootless, plant Marchantia polymorpha upon phosphate starvation conditions. The findings in this study shed light on the mechanisms that early land plants may have developed for the conquest of substrates poor in available phosphate, some of which are still conserved by current-day plants. Moreover, our results open several working hypotheses and novel perspectives for the study of Pi-starvation responses along plant evolution.

The most studied DNA methylation pathway in plants is the RNA Directed DNA Methylation (RdDM), a conserved mechanism that involves the role of noncoding RNAs to control the expansion of the noncoding genome. Genome-wide DNA methylation levels have been reported to correlate with genome size. However, little is known about the catalog of noncoding RNAs and the impact on DNA methylation in small plant genomes with reduced noncoding regions. Because of the small length of intergenic regions in the compact genome of the carnivorous plant Utricularia gibba, we investigated its repertoire of noncoding RNA and DNA methylation landscape. Here, we report that, compared to other angiosperms, U. gibba has an unusual distribution of small RNAs and reduced global DNA methylation levels. DNA methylation was determined using a novel strategy based on long-read DNA sequencing with the Pacific Bioscience platform and confirmed by whole-genome bisulfite sequencing. Moreover, some key genes involved in the RdDM pathway may not represented by compensatory paralogs or comprise truncated proteins, for example, U. gibba DICER-LIKE 3 (DCL3), encoding a DICER endonuclease that produces 24-nt small-interfering RNAs, has lost key domains required for complete function. Our results unveil that a truncated DCL3 correlates with a decreased proportion of 24-nt small-interfering RNAs, low DNA methylation levels, and developmental abnormalities during female gametogenesis in U. gibba. Alterations in female gametogenesis are reminiscent of RdDM mutant phenotypes in Arabidopsis thaliana. It would be interesting to further study the biological implications of the DCL3 truncation in U. gibba, as it could represent an initial step in the evolution of RdDM pathway in compact genomes.

No abstract

No abstract

Plant hormones are master regulators of developmental and genetic mechanisms to deal with diverse environmental cues. Upon phosphate (Pi) limitation, vascular plants modify phytohormone metabolism to coordinate diverse mechanisms to overcome such stress. However, the transcriptional program underlying hormonal signaling in response to Pi scarcity in early branches of land plant phylogeny, remains unclear. Therefore, we explored the transcriptional dynamics of key genes involved in auxin, cytokinin, ethylene, jasmonate, gibberellin and abscisic acid metabolism in the early divergent land plant Marchantia polymorpha, upon Pi starvation. Our RNAseq approach revealed major changes in genes associated with auxin and ethylene biosynthesis. Genes involved in cytokinin synthesis are repressed. Interestingly, genes involved in auxin and ethylene signaling such as MpARF1 and MpARF2 are upregulated. In contrast, MpARRb is down-regulated. Moreover, genes involved in the synthesis of jasmonates were highly upregulated, but those related to signaling did not change in expression. Our data suggest that auxin and ethylene act as positive regulators of rhizoid development under Pi-limited conditions, whereas cytokinin may act as a negative regulator. The transcriptional behaviour of some hormonerelated genes in Marchantia is similar to those described in controlling root hair development in Arabidopsis, Maize and Rice, upon Pi scarcity.

The most studied DNA methylation pathway in plants is the RNA Directed DNA Methylation (RdDM), a conserved mechanism that involves the role of noncoding RNAs to control the expansion of the noncoding genome. Genome-wide methylation levels have been reported to correlate with genome size. However, little is known about the catalog of noncoding RNAs and the impact on DNA methylation in compact plant genomes. Because the small genome size of the carnivorous plant Utricularia gibba we investigate the noncoding RNA landscape and global DNA methylation in a compact genome. Here, we report that, compared to other angiosperms, U. gibba has an unusual distribution of noncoding RNAs and reduced global DNA methylation levels, as determined by a novel strategy based on long-read DNA sequencing with the Pacific Bioscience platform and confirmed by whole-genome bisulfite sequencing. Moreover, reduced DNA methylation correlates with lack of a functional RdDM pathway, as U. gibba DICER-LIKE 3 (DCL3), encoding a DICER endonuclease that produces 24-nt small-interfering RNAs lost key domains required for complete function. Our findings unveil that lack of a functional DCL3 in U. gibba correlates with a decreased proportion of 24-nt small-interfering RNAs, low genome methylation levels, and developmental abnormalities during female gametogenesis that are reminiscent of RdDM mutant phenotypes in Arabidopsis thaliana. It would be interesting to further study the biological implications of the DCL3 truncation in U. gibba, as it could represent an initial step in the evolution of apomixis in compact genomes.

SummaryThe most studied DNA methylation pathway in plants is the RNA Directed DNA Methylation (RdDM), which is a conserved mechanism that involves noncoding-RNAs to control the expansion of intergenic regions. However, little is known about relationship between plant genome size reductions and DNA methylation.Because the compact genome size of the carnivorous plant Utricularia gibba, we investigate in this plant the noncoding-RNA landscape and DNA methylation through a combination of cytological, evolutionary, and genome-wide transcriptomic and methylation approaches.We report an unusual distribution of noncoding RNAs in U. gibba in comparison with other characterized angiosperms, which correlated with a lower level of global genome methylation, as determined by a novel strategy based on long-read DNA sequencing and corroborated by whole-genome bisulfite analysis. Moreover, found that genes involved in the RdDM pathway may not be functionally active in U. gibba, including a truncated DICER-LIKE 3 (DCL3), involved in the production of 24-nt small-RNAs.Our findings suggest that selective pressure to conserve a fully functional RdDM pathway might be reduced in compact genomes and a defective DCL3 correlate with a decreased proportion of 24-nt small-RNAs and developmental alterations in U. gibba, which could represent an initial step in the evolution of apomixis.