Identifying the genetic mechanisms that underpin the evolution of new organ and tissue systems is an aim of evolutionary developmental biology. Comparative functional genetic studies between angiosperms and bryophytes can define those genetic changes that were responsible for developmental innovations. Here, we report the generation of a transcriptome atlas covering most phases in the life cycle of the model bryophyte Physcomitrella patens, including detailed sporophyte developmental progression. We identified a comprehensive set of sporophyte-specific transcription factors, and found that many of these genes have homologs in angiosperms that function in developmental processes such as flowering and shoot branching. Deletion of the PpTCP5 transcription factor results in development of supernumerary sporangia attached to a single seta, suggesting that it negatively regulates branching in the moss sporophyte. Given that TCP genes repress branching in angiosperms, we suggest that this activity is ancient. Finally, comparison of P. patens and Arabidopsis thaliana transcriptomes led us to the identification of a conserved core of transcription factors expressed in tip-growing cells. We identified modifications in the expression patterns of these genes that could account for developmental differences between P. patens tip-growing cells and A. thaliana pollen tubes and root hairs.

Elucidating the genetic basis of morphological changes in evolution remains a major challenge in biology. Repeated independent trait changes are of particular interest because they can indicate adaptation in different lineages or genetic and developmental constraints on generating morphological variation. In animals, changes to "hot spot" genes with minimal pleiotropy and large phenotypic effects underlie many cases of repeated morphological transitions. By contrast, only few such genes have been identified from plants, limiting cross-kingdom comparisons of the principles of morphological evolution. Here, we demonstrate that the REDUCED COMPLEXITY (RCO) locus underlies more than one naturally evolved change in leaf shape in the Brassicaceae. We show that the difference in leaf margin dissection between the sister species Capsella rubella and Capsella grandiflora is caused by cis-regulatory variation in the homeobox gene RCO-A, which alters its activity in the developing lobes of the leaf. Population genetic analyses in the ancestral C. grandiflora indicate that the more-active C. rubella haplotype is derived from a now rare or lost C. grandiflora haplotype via additional mutations. In Arabidopsis thaliana, the deletion of the RCO-A and RCO-B genes has contributed to its evolutionarily derived smooth leaf margin, suggesting the RCO locus as a candidate for an evolutionary hot spot. We also find that temperature-responsive expression of RCO-A can explain the phenotypic plasticity of leaf shape to ambient temperature in Capsella, suggesting a molecular basis for the well-known negative correlation between temperature and leaf margin dissection.

SummaryStreptophytes colonized the land some time before 470 million years ago [1, 2, 3]. The colonization coincided with an increase in morphological and cellular diversity [4, 5, 6, 7]. This increase in diversity is correlated with a proliferation in transcription factors encoded in genomes [8, 9, 10]. This suggests that gene duplication and subsequent diversification of function was instrumental in the generation of land plant diversity. Here, we investigate the diversification of the streptophyte-specific Lotus japonicus ROOTHAIRLESS LIKE (LRL) transcription factor (TF) [11, 12] subfamily of basic loop helix (bHLH) proteins by comparing gene function in early divergent and derived land plant species. We report that the single Marchantia polymorpha LRL gene acts as a general growth regulator required for rhizoid development, a function that has been partially conserved throughout multicellular streptophytes. In contrast, the five relatively derived Arabidopsis thaliana LRL genes comprise two antagonistically acting groups of differentially expressed genes. The diversification of LRL genes accompanied the evolution of an antagonistic regulatory element controlling root hair development.

Lateral inhibition patterns differentiated cell types among equivalent cells during development in bacteria, metazoans, and plants. Tip-growing rhizoid cells develop among flat epidermal cells in the epidermis of the early-diverging land plant Marchantia polymorpha. We show that the majority of rhizoid cells develop individually, but some develop in linear, one-dimensional groups (chains) of between 2 and 7 rhizoid cells in wild-type plants. The distribution of rhizoid cells can be accounted for within a simple cellular automata model of lateral inhibition. The model predicted that in the absence of lateral inhibition, two-dimensional rhizoid cell groups (clusters) form. These can be larger than those formed with lateral inhibition. M. polymorpha rhizoid differentiation is positively regulated by the ROOT HAIR DEFECTIVE SIX-LIKE1 (MpRSL1) basic-helix-loop-helix (bHLH) transcription factor, which is directly repressed by the FEW RHI-ZOIDS1 (MpFRH1) microRNA (miRNA). To test if MpFRH1 miRNA acts during lateral inhibition, we generated loss-of-function (lof) mutants without the MpFRH1 miRNA. Two-dimensional clusters of rhizoids develop in Mpfrh1 lof mutants as predicted by the model for plants that lack lateral inhibition. Furthermore, two-dimensional clusters of up to 9 rhizoid cells developed in the Mpfrh1 lof mutants compared to a maximum number of 7 observed in wild-type groups. The higher steady-state levels of MpRSL1 mRNA in Mpfrh1 lof mutants indicate that MpFRH1-mediated lateral inhibition involves the repression of MpRSL1 activity. Together, the modeling and genetic data indicate that MpFRH1 miRNA mediates lateral inhibition by repressing MpRSL1 during pattern formation in the M. polymorpha epidermis.

Basic helix-loop-helix transcription factors encoded by RSL class I genes control a gene regulatory network that positively regulates the development of filamentous rooting cells – root hairs and rhizoids – in land plants. The GLABRA2 transcription factor negatively regulates these genes in the angiosperm Arabidopsis thaliana. To find negative regulators of RSL class I genes in early diverging land plants we conducted a mutant screen in the liverwort Marchantia polymorpha. This identified FEW RHIZOIDS1 (MpFRH1) microRNA (miRNA) that negatively regulates the RSL class I gene MpRSL1. The miRNA and its mRNA target constitute a feedback mechanism that controls epidermal cell differentiation. MpFRH1 miRNA target sites are conserved among liverwort RSL class I mRNAs but are not present in RSL class I mRNAs of other land plants. These findings indicate that while RSL class I genes are ancient and conserved, independent negative regulatory mechanisms evolved in different lineages during land plant evolution.