Molecular control of normal and acrocona mutant seed cone development in Norway spruce (Picea abies) and the evolution of conifer ovule-bearing organs

Reflections on the ABC model of flower development

The formulation of the ABC model by a handful of pioneer plant developmental geneticists was a seminal event in the quest to answer a seemingly simple question: how are flowers formed? Fast forward 30 years and this elegant model has generated a vibrant and diverse community, capturing the imagination of developmental and evolutionary biologists, structuralists, biochemists and molecular biologists alike. Together they have managed to solve many floral mysteries, uncovering the regulatory processes that generate the characteristic spatio-temporal expression patterns of floral homeotic genes, elucidating some of the mechanisms allowing ABC genes to specify distinct organ identities, revealing how evolution tinkers with the ABC to generate morphological diversity, and even shining a light on the origins of the floral gene regulatory network itself. Here we retrace the history of the ABC model, from its genesis to its current form, highlighting specific milestones along the way before drawing attention to some of the unsolved riddles still hidden in the floral alphabet.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

Thinking outside the F-box: how UFO controls angiosperm development

The formation of inflorescences and flowers is essential for the successful reproduction of angiosperms. In the past few decades, genetic studies have identified the LEAFY transcription factor and the UNUSUAL FLORAL ORGANS (UFO) F-box protein as two major regulators of flower development in a broad range of angiosperm species. Recent research has revealed that UFO acts as a transcriptional cofactor, redirecting the LEAFY floral regulator to novel cis-elements. In this review, we summarize the various roles of UFO across species, analyze past results in light of new discoveries and highlight the key questions that remain to be solved.

Phytohormone biosynthesis and transcriptional analyses provide insight into the main growth stage of male and female cones Pinus koraiensis

The cone is a crucial component of the whole life cycle of gymnosperm and an organ for sexual reproduction of gymnosperms. In Pinus koraiensis, the quantity and development process of male and female cones directly influence seed production, which in turn influences the tree's economic value. There are, however, due to the lack of genetic information and genomic data, the morphological development and molecular mechanism of female and male cones of P. koraiensis have not been analyzed. Long-term phenological observations were used in this study to document the main process of the growth of both male and female cones. Transcriptome sequencing and endogenous hormone levels at three critical developmental stages were then analyzed to identify the regulatory networks that control these stages of cones development. The most significant plant hormones influencing male and female cones growth were discovered to be gibberellin and brassinosteroids, according to measurements of endogenous hormone content. Additionally, transcriptome sequencing allowed the identification of 71,097 and 31,195 DEGs in male and female cones. The synthesis and control of plant hormones during cones growth were discovered via enrichment analysis of key enrichment pathways. FT and other flowering-related genes were discovered in the coexpression network of flower growth development, which contributed to the growth development of male and female cones of P. koraiensis. The findings of this work offer a cutting-edge foundation for understanding reproductive biology and the molecular mechanisms that control the growth development of male and female cones in P. koraiensis.

Cone-setting in spruce is regulated by conserved elements of the age-dependent flowering pathway

Reproductive phase change is well characterized in angiosperm model species, but less studied in gymnosperms. We utilize the early cone-setting acrocona mutant to study reproductive phase change in the conifer Picea abies (Norway spruce), a gymnosperm. The acrocona mutant frequently initiates cone-like structures, called transition shoots, in positions where wild-type P. abies always produces vegetative shoots. We collect acrocona and wild-type samples, and RNA-sequence their messenger RNA (mRNA) and microRNA (miRNA) fractions. We establish gene expression patterns and then use allele-specific transcript assembly to identify mutations in acrocona. We genotype a segregating population of inbred acrocona trees. A member of the SQUAMOSA BINDING PROTEIN-LIKE (SPL) gene family, PaSPL1, is active in reproductive meristems, whereas two putative negative regulators of PaSPL1, miRNA156 and the conifer specific miRNA529, are upregulated in vegetative and transition shoot meristems. We identify a mutation in a putative miRNA156/529 binding site of the acrocona PaSPL1 allele and show that the mutation renders the acrocona allele tolerant to these miRNAs. We show co-segregation between the early cone-setting phenotype and trees homozygous for the acrocona mutation. In conclusion, we demonstrate evolutionary conservation of the age-dependent flowering pathway and involvement of this pathway in regulating reproductive phase change in the conifer P. abies.

Transcriptional Regulation of Pine Male and Female Cone Initiation and Development: Key Players Identified Through Comparative Transcriptomics

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly been identified through homology and phylogenetic reconstruction with their angiosperm counterparts. We used RNA-sequencing to generate transcriptomes of the early morphological stages of cone development in the conifer species Pinus densiflora and used these to gain a deeper insight into the transcriptional changes during male and female cone development. Paired-end Illumina sequencing was used to generate transcriptomes from non-reproductive tissue and male and female cones at four time points with a total of 382.82 Gbp of data generated. After assembly and stringent filtering, a total of 37,164 transcripts were retrieved, of which a third were functionally annotated using the Mercator plant pipeline. Differentially expressed gene (DEG) analysis resulted in the identification of 172,092 DEGs in the nine tissue types. This, alongside GO gene enrichment analyses, pinpointed transcripts putatively involved in conifer reproductive structure development, including co-orthologs of several angiosperm flowering genes and several that have not been previously reported in conifers. This study provides a comprehensive transcriptome resource for male and early female cone development in the gymnosperm species Pinus densiflora. Characterisation of this resource has allowed the identification of potential key players and thus provides valuable insights into the molecular regulation of reproductive structure development in monoecious conifers.

Phytohormone profiles and related gene expressions after endodormancy release in developing Pinus tabuliformis male strobili

Development after endo-dormancy release ensures perennial plants, such as forest trees, proper response to environmental changes and enhances their adaptability. In northern hemisphere, megasporophore and microsporophore of conifers undergo dormancy to complete their development. Here combined with transcriptome data, we used high-performance liquid chromatography/electrospray ionization tandem mass spectrometry (ESI-HPLC-MS/MS) to quantitatively analyse the various hormones (Abscisic Acid (ABA), 3-Indoleacetic acid (IAA), Gibberellins (GAs), Cytokinin (CTK), Jasmonic acid (JA) and Salicylic acid (SA)) of Chinese pine (Pinus tabuliformis Carr.) male strobili after endo-dormancy release. More specifically, we analysed endogenous hormones and their related-genes and verified the important role of ABA in plants growth and development. We observed rapid decrease in ABA content after dormancy release, resulting in reducing the inhibitory effect on male strobili growth. Similarly, rapid drop in ABA/GA ratio was observed and was associated with the start of male strobili growth and development. Combined with transcriptome data, we found that HAB2-SnRK2.10 played a central role in the ABA pathway in the entire network of hormones regulating male strobili development. Due to external environment warming, the differentially expressed HAB2-SnRK gene led to ABA content rapid decline, thus initiating male strobili growth. We constructed a network of hormone-regulated development to understand the interactions between hormones after male strobili dormancy release of male strobili. This study provided essential foundations for studying megasporophore and microsporophore growth mechanism after endo-dormancy and offered new ideas for flower development in gymnosperms and angiosperms.

Reconstructing development of the earliest seed integuments raises a new hypothesis for the evolution of ancestral seed-bearing structures

How plant seeds originated remains unresolved, in part due to disconnects between fossil intermediates and developmental genetics in extant species. The Carboniferous fossil Genomosperma is considered among the most primitive known seeds, with highly lobed integument and exposed nucellus. We have used this key fossil taxon to investigate the evolutionary origins of seed development. We examined sectioned Genomosperma specimens using modern digital 3D reconstruction techniques and established population-level measurements of Genomosperma ovules for quantitative analysis. Genomosperma ovules show significant variation in integumentary lobe fusion and curvature. Our analysis suggests that this variation represents a single species with significant variations in lobe number and fusion, reminiscent of floral development in extant species. We conclude that changes in lobe flexure occurred late in development, consistent with a previously hypothesized function in pollen guidance/retention. We also identify seeds of Genomosperma within cupules for the first time. The presence of a cupule adds evidence towards the plesiomorphy of cupules within seed plants. Together with the similarities identified between the Genomosperma lobed integument and floral organs, we propose that the cupule, integument and nucellus together developed in a shoot-like fashion, potentially ancestral to extant seed plant reproductive shoots.

The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development

BACKGROUND

Cunninghamia lanceolata (Chinese fir), a member of the conifer family Cupressaceae, is one of the most popular cultivated trees for wood production in China. Continuous research is being performed to improve C. lanceolata breeding values. Given the high rate of seed abortion (one of the reasons being the failure of ovule and pollen development) in C. lanceolata, the proper formation of female/male cones could theoretically increase the number of offspring in future generations. MIKC MADS-box genes are well-known for their roles in the flower/cone development and comprise the typical/atypical floral development model for both angiosperms and gymnosperms.

RESULTS

We performed a transcriptomic analysis to find genes differentially expressed between female and male cones at a single, carefully determined developmental stage, focusing on the MIKC MADS-box genes. We finally obtained 47 unique MIKC MADS-box genes from C. lanceolata and divided these genes into separate branches. 27 out of the 47 MIKC MADS-box genes showed differential expression between female and male cones, and most of them were not expressed in leaves. Out of these 27 genes, most B-class genes (AP3/PI) were up-regulated in the male cone, while TM8 genes were up-regulated in the female cone. Then, with no obvious overall preference for AG (class C + D) genes in female/male cones, it seems likely that these genes are involved in the development of both cones. Finally, a small number of genes such as GGM7, SVP, AGL15, that were specifically expressed in female/male cones, making them candidate genes for sex-specific cone development.

CONCLUSIONS

Our study identified a number of MIKC MADS-box genes showing differential expression between female and male cones in C. lanceolata, illustrating a potential link of these genes with C. lanceolata cone development. On the basis of this, we postulated a possible cone development model for C. lanceolata. The gene expression library showing differential expression between female and male cones shown here, can be used to discover unknown regulatory networks related to sex-specific cone development in the future.

Colourful cones: how did flower colour first evolve?

Angiosperms that are biotically pollinated typically produce flowers with bright and contrasting colours that help to attract pollinators and hence contribute to the reproductive success of the species. This colourful array contrasts with the much less multicoloured reproductive structures of the four living gymnosperm lineages, which are mostly wind pollinated, though cycads and Gnetales are predominantly pollinated by insects that feed on surface fluids from the pollination drops. This review examines the possible evolutionary pathways and cryptic clues for flower colour in both living and fossil seed plants. It investigates how the ancestral flowering plants could have overcome the inevitable trade-off that exists between attracting pollinators and minimizing herbivory, and explores the possible evolutionary and biological inferences from the colours that occur in some living gymnosperms. The red colours present in the seed-cone bracts of some living conifers result from accumulation of anthocyanin pigments; their likely primary function is to help protect the growing plant tissues under particular environmental conditions. Thus, the visual cue provided by colour in flower petals could have first evolved as a secondary effect, probably post-dating the evolution of bee colour vision but occurring before the subsequent functional accumulation of a range of different flower pigments.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

A Dead Gene Walking: Convergent Degeneration of a Clade of MADS-Box Genes in Crucifers

Genes are "born," and eventually they "die." These processes shape the phenotypic evolution of organisms and are hence of great biological interest. If genes die in plants, they generally do so quite rapidly. Here, we describe the fate of GOA-like genes that evolve in a dramatically different manner. GOA-like genes belong to the subfamily of Bsister genes of MIKC-type MADS-box genes. Typical MIKC-type genes encode conserved transcription factors controlling plant development. We show that ABS-like genes, a clade of Bsister genes, are indeed highly conserved in crucifers (Brassicaceae) maintaining the ancestral function of Bsister genes in ovule and seed development. In contrast, their closest paralogs, the GOA-like genes, have been undergoing convergent gene death in Brassicaceae. Intriguingly, erosion of GOA-like genes occurred after millions of years of coexistence with ABS-like genes. We thus describe Delayed Convergent Asymmetric Degeneration, a so far neglected but possibly frequent pattern of duplicate gene evolution that does not fit classical scenarios. Delayed Convergent Asymmetric Degeneration of GOA-like genes may have been initiated by a reduction in the expression of an ancestral GOA-like gene in the stem group of Brassicaceae and driven by dosage subfunctionalization. Our findings have profound implications for gene annotations in genomics, interpreting patterns of gene evolution and using genes in phylogeny reconstructions of species.

Gene duplication led to divergence of expression patterns, protein-protein interaction patterns and floral development functions of AGL6-like genes in the basal angiosperm Magnolia wufengensis (Magnoliaceae)

The MADS-box family genes play critical roles in the regulation of growth and development of flowering plants. AGAMOUS-LIKE 6 (AGL6)-like genes are one of the most enigmatic subfamilies of the MADS-box family because of highly variable expression patterns and ambiguous functions, which have long puzzled researchers. A lot of AGL6 homologs have been identified from gymnosperms and angiosperms. However, only a few have been characterized, especially for basal angiosperm taxa. Magnolia wufengensis is a woody basal angiosperm from the family Magnoliaceae. In the current study, the phylogenesis, expression and protein-protein interaction (PPI) patterns, and functions of two AGL6 homologs from M. wufengensis, MawuAGL6-1 and MawuAGL6-2, were analyzed. Phylogenetic analysis indicated that the two AGL6 duplicates may have arisen by gene duplication before the divergence of Magnoliaceae and Lauraceae, with the diversification of their expression and PPI patterns after gene duplication. Functional analysis revealed that, in addition to common functions in accelerating flowering, MawuAGL6-1 might be responsible for flower meristem determinacy, while MawuAGL6-2 is preferentially recruited to regulate tepal morphogenesis. These findings further advance our understanding of the evolution of phylogenesis, expression, interaction and functions of AGL6 lineage genes from basal angiosperms, as well as the entire AGL6 lineage genes, and the significance of AGL6 lineage genes in the evolution and biological diversity.

Origination and selection of ABCDE and AGL6 subfamily MADS-box genes in gymnosperms and angiosperms

BACKGROUND

The morphological diversity of flower organs is closely related to functional divergence within the MADS-box gene family. Bryophytes and seedless vascular plants have MADS-box genes but do not have ABCDE or AGAMOUS-LIKE6 (AGL6) genes. ABCDE and AGL6 genes belong to the subgroup of MADS-box genes. Previous works suggest that the B gene was the first ABCDE and AGL6 genes to emerge in plant but there are no mentions about the probable origin time of ACDE and AGL6 genes. Here, we collected ABCDE and AGL6 gene 381 protein sequences and 361 coding sequences from gymnosperms and angiosperms and reconstructed a complete Bayesian phylogeny of these genes. In this study, we want to clarify the probable origin time of ABCDE and AGL6 genes is a great help for understanding the role of the formation of the flower, which can decipher the forming order of MADS-box genes in the future.

RESULTS

These genes appeared to have been under purifying selection and their evolutionary rates are not significantly different from each other. Using the Bayesian evolutionary analysis by sampling trees (BEAST) tool, we estimated that: the mutation rate of the ABCDE and AGL6 genes was 2.617 × 10-3 substitutions/site/million years, and that B genes originated 339 million years ago (MYA), CD genes originated 322 MYA, and A genes shared the most recent common ancestor with E/AGL6 296 MYA, respectively.

CONCLUSIONS

The phylogeny of ABCDE and AGL6 genes subfamilies differed. The APETALA1 (AP1 or A gene) subfamily clustered into one group. The APETALA3/PISTILLATA (AP3/PI or B genes) subfamily clustered into two groups: the AP3 and PI clades. The AGAMOUS/SHATTERPROOF/SEEDSTICK (AG/SHP/STK or CD genes) subfamily clustered into a single group. The SEPALLATA (SEP or E gene) subfamily in angiosperms clustered into two groups: the SEP1/2/4 and SEP3 clades. The AGL6 subfamily clustered into a single group. Moreover, ABCDE and AGL6 genes appeared in the following order: AP3/PI → AG/SHP/STK → AGL6/SEP/AP1. In this study, we collected candidate sequences from gymnosperms and angiosperms. This study highlights important events in the evolutionary history of the ABCDE and AGL6 gene families and clarifies their evolutionary path.

Pinus massoniana Introgression Hybrids Display Differential Expression of Reproductive Genes

Pinus massoniana and P. hwangshanensis are two conifer species located in southern China, which are of both economic and ornamental value. Around the middle and lower reaches of the Yangtze River, P. massoniana occurs mainly at altitudes below 700 m, while P. hwangshanensis can be found above 900 m. At altitudes where the distribution of both pines overlaps, a natural introgression hybrid exists, which we will further refer to as the Z pine. This pine has a morphological character that shares attributes of both P. massoniana and P. hwangshanensis. However, compared to the other two pines, its reproductive structure, the pinecone, has an ultra-low ripening rate with seeds that germinate poorly. In this study, we aimed to find the reason for the impaired cone maturation by comparing transcriptome libraries of P. massoniana and Z pine cones at seven successive growth stages. After sequencing and assembly, we obtained unigenes and then annotated them against NCBI’s non-redundant nucleotide and protein sequences, Swiss-Prot, Clusters of Orthologous Groups, Gene Ontology and KEGG Orthology databases. Gene expression levels were estimated and differentially expressed genes (DEGs) of the two pines were mined and analyzed. We found that several of them indeed relate to reproductive process. At every growth stage, these genes are expressed at a higher level in P. massoniana than in the Z pine. These data provide insight into understanding which molecular mechanisms are altered between P. massoniana and the Z pine that might cause changes in the reproductive process.

The Evolution of Sexual Fluids in Gymnosperms From Pollination Drops to Nectar

A current synthesis of data from modern and fossil plants paints a new picture of sexual fluids, including nectar, as a foundational component of gymnosperm reproductive evolution. We review the morpho-anatomical adaptations, their accompanying secretions, and the functional compounds involved. We discuss two types of secretions: (1) those involved in fertilization fluids produced by gametophytes and archegonia of zooidogamous gymnosperms, i.e., Ginkgo and cycads, and (2) those involved in pollen capture mechanisms (PCMs), i.e., pollination drops. Fertilization fluids provide both liquid in which sperm swim, as well as chemotactic signals that direct sperm to the egg. Such fertilization fluids were probably found among many extinct plants such as ancient cycads and others with swimming sperm, but were subsequently lost upon the evolution of siphonogamy (direct delivery of sperm to the egg by pollen tubes), as found in modern gnetophytes, conifers, and Pinaceae. Pollination drops are discussed in terms of three major types of PCMs and the unique combinations of morphological and biochemical adaptations that define each. These include their amino acids, sugars, calcium, phosphate and proteins. The evolution of PCMs is also discussed with reference to fossil taxa. The plesiomorphic state of extant gymnosperms is a sugar-containing pollination drop functioning as a pollen capture surface, and an in ovulo pollen germination medium. Additionally, these drops are involved in ovule defense, and provide nectar for pollinators. Pollination drops in anemophilous groups have low sugar concentrations that are too low to provide insects with a reward. Instead, they appear to be optimized for defense and microgametophyte development. In insect-pollinated modern Gnetales a variety of tissues produce sexual fluids that bear the biochemical signature of nectar. Complete absence of fluid secretions is restricted to a few, poorly studied modern conifers, and is presumably derived. Aspects of pollination drop dynamics, e.g., regulation of secretion and retraction, are reviewed. Lastly, we discuss pollination drops' control of pollen germination. Large gaps in our current knowledge include the composition of fertilization fluids, the pollination drops of Podocarpaceae, and the overall hydrodynamics of sexual fluids in general.

LEAFY maintains apical stem cell activity during shoot development in the fern Ceratopteris richardii

During land plant evolution, determinate spore-bearing axes (retained in extant bryophytes such as mosses) were progressively transformed into indeterminate branching shoots with specialized reproductive axes that form flowers. The LEAFY transcription factor, which is required for the first zygotic cell division in mosses and primarily for floral meristem identity in flowering plants, may have facilitated developmental innovations during these transitions. Mapping the LEAFY evolutionary trajectory has been challenging, however, because there is no functional overlap between mosses and flowering plants, and no functional data from intervening lineages. Here, we report a transgenic analysis in the fern Ceratopteris richardii that reveals a role for LEAFY in maintaining cell divisions in the apical stem cells of both haploid and diploid phases of the lifecycle. These results support an evolutionary trajectory in which an ancestral LEAFY module that promotes cell proliferation was progressively co-opted, adapted and specialized as novel shoot developmental contexts emerged.

The first plants colonized land around 500 million years ago. These plants had simple shoots with no branches, similar to the mosses that live today. Later on, some plants evolved more complex structures including branched shoots and flowers (collectively known as the “flowering plants”). Ferns are a group of plants that evolved midway between the mosses and flowering plants and have branched shoots but no flowers.

The gradual transition from simple to more complex plant structures required changes to the way in which cells divide and grow within plant shoots. Whereas animals produce new cells throughout their body, most plant cells divide in areas known as meristems. All plants grow from embryos, which contain meristems that will form the roots and shoots of the mature plant. A gene called LEAFY is required for cells in moss embryos to divide. However, in flowering plants LEAFY does not carry out this role, instead it is only required to make the meristems that produce flowers.

How did LEAFY transition from a general role in embryos to a more specialized role in making flowers? To address this question, Plackett, Conway et al. studied the two LEAFY genes in a fern called Ceratopteris richardii. The experiments showed that at least one of these LEAFY genes was active in the meristems of fern shoots throughout the lifespan of the plant. The shoots of ferns with less active LEAFY genes could not form the leaves seen in normal C. richardii plants. This suggests that as land plants evolved, the role of LEAFY changed from forming embryos to forming complex shoot structures.

Most of our major crops are flowering plants. By understanding how the role of LEAFY has changed over the evolution of land plants, it might be possible to manipulate LEAFY genes in crop plants to alter shoot structures to better suit specific environments.

Incorporating fossils into the Pinaceae tree of life

PREMISE OF THE STUDY

Pinaceae have a rich but enigmatic early fossil record, much of which is represented by permineralized seed cones. Our incomplete knowledge of morphology and anatomy in living and extinct species poses an important barrier to understanding their phylogenetic relationships and timing of diversification.

METHODS

We expanded a morphology matrix to 46 fossil and 31 extant Pinaceae species, mainly adding characters from stem and leaf anatomy and seed cones. Using parsimony and Bayesian inference, we compared phylogenetic relationships for extant taxa with and without fossils from the morphology matrix combined with an alignment of plastid gene sequences.

KEY RESULTS

Combined analysis of morphological and molecular characters resulted in a phylogeny of extant Pinaceae that was robust at all nodes except those relating to the interrelationships of Pinus, Picea, and Cathaya and the position of Cedrus. Simultaneous analysis of all fossil and extant species did not result in changes in the relationships among the extant species but did greatly reduce branch support. We found that the placement of most fossils was sensitive to the method of phylogenetic reconstruction when analyzing them singly with the extant species.

CONCLUSIONS

A robust phylogenetic hypothesis for the main lineages of Pinaceae is emerging. Most Early Cretaceous fossils are stem or crown lineages of Pinus, but close relationships also were found between fossils and several other extant genera. The phylogenetic position of fossils broadly supports the existence of extant genera in the Lower Cretaceous.

Understanding the cone scale in Cupressaceae: insights from seed-cone teratology in Glyptostrobus pensilis

Both wild-type and teratological seed cones are described in the monoecious conifer Glyptostrobus pensilis and compared with those of other Cupressaceae sensu lato and other conifers. Some Cupressaceae apparently possess a proliferation of axillary structures in their cone scales. In our interpretation, in Glyptostrobus each bract of both typical and atypical seed cones bears two descending accessory shoots, interpreted here as seed scales (ovuliferous scales). The primary seed scale is fertile and forms the ovules, the second is sterile and forms characteristic tooth-like structures. The bract and the two axillary seed scales are each supplied with a single distinct vascular bundle that enters the cone axis as a separate strand; this vasculature also characterises the descending accessory short shoots in the vegetative parts of the crown. In wild-type seed cones, the fertile seed scale is reduced to its ovules, and the ovules are always axillary. In contrast, the ovules of some of the teratological seed cones examined were located at the centre of the cone scale. An additional tissue found on the upper surface of the sterile lower seed scale is here interpreted as the axis of the fertile seed scale. Thus, the central position of the ovules can be explained by recaulescent fusion of the upper fertile and lower sterile seed scales. In several teratological cone scales, the ovules were enveloped by an additional sterile tissue that is uniseriate and represents an epidermal outgrowth of the fertile seed scale. Close to the ovules, the epidermis was detached from lower tissue and surrounded the ovule completely, except at the micropyle. These teratological features are potentially significant in understanding seed-cone homologies among extant conifers.

Integrative Analysis of Three RNA Sequencing Methods Identifies Mutually Exclusive Exons of MADS-Box Isoforms During Early Bud Development in Picea abies

Recent efforts to sequence the genomes and transcriptomes of several gymnosperm species have revealed an increased complexity in certain gene families in gymnosperms as compared to angiosperms. One example of this is the gymnosperm sister clade to angiosperm TM3-like MADS-box genes, which at least in the conifer lineage has expanded in number of genes. We have previously identified a member of this sub-clade, the conifer gene DEFICIENS AGAMOUS LIKE 19 (DAL19), as being specifically upregulated in cone-setting shoots. Here, we show through Sanger sequencing of mRNA-derived cDNA and mapping to assembled conifer genomic sequences that DAL19 produces six mature mRNA splice variants in Picea abies. These splice variants use alternate first and last exons, while their four central exons constitute a core region present in all six transcripts. Thus, they are likely to be transcript isoforms. Quantitative Real-Time PCR revealed that two mutually exclusive first DAL19 exons are differentially expressed across meristems that will form either male or female cones, or vegetative shoots. Furthermore, mRNA in situ hybridization revealed that two mutually exclusive last DAL19 exons were expressed in a cell-specific pattern within bud meristems. Based on these findings in DAL19, we developed a sensitive approach to transcript isoform assembly from short-read sequencing of mRNA. We applied this method to 42 putative MADS-box core regions in P. abies, from which we assembled 1084 putative transcripts. We manually curated these transcripts to arrive at 933 assembled transcript isoforms of 38 putative MADS-box genes. 152 of these isoforms, which we assign to 28 putative MADS-box genes, were differentially expressed across eight female, male, and vegetative buds. We further provide evidence of the expression of 16 out of the 38 putative MADS-box genes by mapping PacBio Iso-Seq circular consensus reads derived from pooled sample sequencing to assembled transcripts. In summary, our analyses reveal the use of mutually exclusive exons of MADS-box gene isoforms during early bud development in P. abies, and we find that the large number of identified MADS-box transcripts in P. abies results not only from expansion of the gene family through gene duplication events but also from the generation of numerous splice variants.

Chicken Toes-Like Leaf and Petalody Flower (CTP) is a novel regulator that controls leaf and flower development in soybean

A soybean mutant displaying chicken toes-like leaves and petalody flowers was identified as being caused by a single recessive gene, termed ctp. Using heterozygous-inbred recombinants, this gene was fine-mapped to a 37-kb region harbouring three predicted genes on chromosome 5. The gene Glyma05g022400.1 was detected to have a 33-bp deletion in its first exon that was responsible for ctp. Validation for this gene was provided by the fact that the 33-bp deletion-derived marker I2 completely co-segregated with the phenotypes of 96 F10-derived residual heterozygous lines and 2200 fine-mapping individuals, and by the fact that the orthologue NbCTP in Nicotiana benthamiana also influenced leaf and flower development under virus-induced gene silencing. In terms of characterization, the CTPs shared highly conserved domains specifically in higher plants, GmCTP was alternatively spliced, and it was expressed in multiple organs, especially in lateral meristems. GmCTP was localized to the nucleus and other regions and performed transcriptional activity. In ctp, the expression levels and splicing patterns were dramatically disrupted, and many key regulators in leaf and/or floral developmental pathways were interrupted. Thus, CTP is a novel and critical pleiotropic regulator of leaf and flower development that participates in multiple regulation pathways, and may play key roles in lateral organ differentiation as a putative novel transcription regulator.

A link between LEAFY and B-gene homologues in Welwitschia mirabilis sheds light on ancestral mechanisms prefiguring floral development

Flowering plants evolved from an unidentified gymnosperm ancestor. Comparison of the mechanisms controlling development in angiosperm flowers and gymnosperm cones may help to elucidate the mysterious origin of the flower. We combined gene expression studies with protein behaviour characterization in Welwitschia mirabilis to test whether the known regulatory links between LEAFY and its MADS-box gene targets, central to flower development, might also contribute to gymnosperm reproductive development. We found that WelLFY, one of two LEAFY-like genes in Welwitschia, could be an upstream regulator of the MADS-box genes APETALA3/PISTILLATA-like (B-genes). We demonstrated that, even though their DNA-binding domains are extremely similar, WelLFY and its paralogue WelNDLY exhibit distinct DNA-binding specificities, and that, unlike WelNDLY, WelLFY shares with its angiosperm orthologue the capacity to bind promoters of Welwitschia B-genes. Finally, we identified several cis-elements mediating these interactions in Welwitschia and obtained evidence that the link between LFY homologues and B-genes is also conserved in two other gymnosperms, Pinus and Picea. Although functional approaches to investigate cone development in gymnosperms are limited, our state-of-the-art biophysical techniques, coupled with expression studies, provide evidence that crucial links, central to the control of floral development, may already have existed before the appearance of flowers.

An Evolutionary Framework for Carpel Developmental Control Genes

Carpels are the female reproductive organs of flowering plants (angiosperms), enclose the ovules, and develop into fruits. The presence of carpels unites angiosperms, and they are suggested to be the most important autapomorphy of the angiosperms, e.g., they prevent inbreeding and allow efficient seed dispersal. Many transcriptional regulators and coregulators essential for carpel development are encoded by diverse gene families and well characterized in Arabidopsis thaliana. Among these regulators are AGAMOUS (AG), ETTIN (ETT), LEUNIG (LUG), SEUSS (SEU), SHORT INTERNODE/STYLISH (SHI/STY), and SEPALLATA1, 2, 3, 4 (SEP1, 2, 3, 4). However, the timing of the origin and their subsequent molecular evolution of these carpel developmental regulators are largely unknown. Here, we have sampled homologs of these carpel developmental regulators from the sequenced genomes of a wide taxonomic sampling of the land plants, such as Physcomitrella patens, Selaginella moellendorfii, Picea abies, and several angiosperms. Careful phylogenetic analyses were carried out that provide a phylogenetic background for the different gene families and provide minimal estimates for the ages of these developmental regulators. Our analyses and published work show that LUG-, SEU-, and SHI/STY-like genes were already present in the Most Recent Common Ancestor (MRCA) of all land plants, AG- and SEP-like genes were present in the MRCA of seed plants and their origin may coincide with the ξ Whole Genome Duplication. Our work shows that the carpel development regulatory network was, in part, recruited from preexisting network components that were present in the MRCA of angiosperms and modified to regulate gynoecium development.

Evolutionary Developmental Biology (Evo-Devo) Research in Latin America

Famous for its blind cavefish and Darwin's finches, Latin America is home to some of the richest biodiversity hotspots of our planet. The Latin American fauna and flora inspired and captivated naturalists from the nineteenth and twentieth centuries, including such notable pioneers such as Fritz Müller, Florentino Ameghino, and Léon Croizat who made a significant contribution to the study of embryology and evolutionary thinking. But, what are the historical and present contributions of the Latin American scientific community to Evo-Devo? Here, we provide the first comprehensive overview of the Evo-Devo laboratories based in Latin America and describe current lines of research based on endemic species, focusing on body plans and patterning, systematics, physiology, computational modeling approaches, ecology, and domestication. Literature searches reveal that Evo-Devo in Latin America is still in its early days; while showing encouraging indicators of productivity, it has not stabilized yet, because it relies on few and sparsely distributed laboratories. Coping with the rapid changes in national scientific policies and contributing to solve social and health issues specific to each region are among the main challenges faced by Latin American researchers. The 2015 inaugural meeting of the Pan-American Society for Evolutionary Developmental Biology played a pivotal role in bringing together Latin American researchers eager to initiate and consolidate regional and worldwide collaborative networks. Such networks will undoubtedly advance research on the extremely high genetic and phenotypic biodiversity of Latin America, bound to be an almost infinite source of amazement and fascinating findings for the Evo-Devo community.

Insights from the pollination drop proteome and the ovule transcriptome of Cephalotaxus at the time of pollination drop production

BACKGROUND AND AIMS

Many gymnosperms produce an ovular secretion, the pollination drop, during reproduction. The drops serve as a landing site for pollen, but also contain a suite of ions and organic compounds, including proteins, that suggests diverse roles for the drop during pollination. Proteins in the drops of species of Chamaecyparis, Juniperus, Taxus, Pseudotsuga, Ephedra and Welwitschia are thought to function in the conversion of sugars, defence against pathogens, and pollen growth and development. To better understand gymnosperm pollination biology, the pollination drop proteomes of pollination drops from two species of Cephalotaxus have been characterized and an ovular transcriptome for C. sinensis has been assembled.

METHODS

Mass spectrometry was used to identify proteins in the pollination drops of Cephalotaxus sinensis and C. koreana RNA-sequencing (RNA-Seq) was employed to assemble a transcriptome and identify transcripts present in the ovules of C. sinensis at the time of pollination drop production.

KEY RESULTS

About 30 proteins were detected in the pollination drops of both species. Many of these have been detected in the drops of other gymnosperms and probably function in defence, polysaccharide metabolism and pollen tube growth. Other proteins appear to be unique to Cephalotaxus, and their putative functions include starch and callose degradation, among others. Together, the proteins appear either to have been secreted into the drop or to occur there due to breakdown of ovular cells during drop production. Ovular transcripts represent a wide range of gene ontology categories, and some may be involved in drop formation, ovule development and pollen-ovule interactions.

CONCLUSIONS

The proteome of Cephalotaxus pollination drops shares a number of components with those of other conifers and gnetophytes, including proteins for defence such as chitinases and for carbohydrate modification such as β-galactosidase. Proteins likely to be of intracellular origin, however, form a larger component of drops from Cephalotaxus than expected from studies of other conifers. This is consistent with the observation of nucellar breakdown during drop formation in Cephalotaxus The transcriptome data provide a framework for understanding multiple metabolic processes that occur within the ovule and the pollination drop just before fertilization. They reveal the deep conservation of WUSCHEL expression in ovules and raise questions about whether any of the S-locus transcripts in Cephalotaxus ovules might be involved in pollen-ovule recognition.

A transcriptomics investigation into pine reproductive organ development

The development of reproductive structures in gymnosperms is still poorly studied because of a lack of genomic information and useful genetic tools. The hermaphroditic reproductive structure derived from unisexual gymnosperms is an even less studied aspect of seed plant evolution. To extend our understanding of the molecular mechanism of hermaphroditism and the determination of sexual identity of conifer reproductive structures in general, unisexual and bisexual cones from Pinus tabuliformis were profiled for gene expression using 60K microarrays. Expression patterns of genes during progression of sexual cone development were analysed using RNA-seq. The results showed that, overall, the transcriptomes of male structures in bisexual cones were more similar to those of female cones. However, the expression of several MADS-box genes in the bisexual cones was similar to that of male cones at the more juvenile developmental stage, while despite these expression shifts, male structures of bisexual cones and normal male cones were histologically indistinguishable and cone development was continuous. This study represents a starting point for in-depth analysis of the molecular regulation of cone development and also the origin of hermaphroditism in pine.

Sequenced genomes and rapidly emerging technologies pave the way for conifer evolutionary developmental biology

Conifers, Ginkgo, cycads and gnetophytes comprise the four groups of extant gymnosperms holding a unique position of sharing common ancestry with the angiosperms. Comparative studies of gymnosperms and angiosperms are the key to a better understanding of ancient seed plant morphologies, how they have shifted over evolution to shape modern day species, and how the genes governing these morphologies have evolved. However, conifers and other gymnosperms have been notoriously difficult to study due to their long generation times, inaccessibility to genetic experimentation and unavailable genome sequences. Now, with three draft genomes from spruces and pines, rapid advances in next generation sequencing methods for genome wide expression analyses, and enhanced methods for genetic transformation, we are much better equipped to address a number of key evolutionary questions relating to seed plant evolution. In this mini-review we highlight recent progress in conifer developmental biology relevant to evo-devo questions. We discuss how genome sequence data and novel techniques might allow us to explore genetic variation and naturally occurring conifer mutants, approaches to reduce long generation times to allow for genetic studies in conifers, and other potential upcoming research avenues utilizing current and emergent techniques. Results from developmental studies of conifers and other gymnosperms in comparison to those in angiosperms will provide information to trace core molecular developmental control tool kits of ancestral seed plants, but foremost they will greatly improve our understanding of the biology of conifers and other gymnosperms in their own right.

The Genome of Dendrobium officinale Illuminates the Biology of the Important Traditional Chinese Orchid Herb

Dendrobium officinale Kimura et Migo is a traditional Chinese orchid herb that has both ornamental value and a broad range of therapeutic effects. Here, we report the first de novo assembled 1.35 Gb genome sequences for D. officinale by combining the second-generation Illumina Hiseq 2000 and third-generation PacBio sequencing technologies. We found that orchids have a complete inflorescence gene set and have some specific inflorescence genes. We observed gene expansion in gene families related to fungus symbiosis and drought resistance. We analyzed biosynthesis pathways of medicinal components of D. officinale and found extensive duplication of SPS and SuSy genes, which are related to polysaccharide generation, and that the pathway of D. officinale alkaloid synthesis could be extended to generate 16-epivellosimine. The D. officinale genome assembly demonstrates a new approach to deciphering large complex genomes and, as an important orchid species and a traditional Chinese medicine, the D. officinale genome will facilitate future research on the evolution of orchid plants, as well as the study of medicinal components and potential genetic breeding of the dendrobe.

Evolution of the Plant Reproduction Master Regulators LFY and the MADS Transcription Factors: The Role of Protein Structure in the Evolutionary Development of the Flower

Understanding the evolutionary leap from non-flowering (gymnosperms) to flowering (angiosperms) plants and the origin and vast diversification of the floral form has been one of the focuses of plant evolutionary developmental biology. The evolving diversity and increasing complexity of organisms is often due to relatively small changes in genes that direct development. These "developmental control genes" and the transcription factors (TFs) they encode, are at the origin of most morphological changes. TFs such as LEAFY (LFY) and the MADS-domain TFs act as central regulators in key developmental processes of plant reproduction including the floral transition in angiosperms and the specification of the male and female organs in both gymnosperms and angiosperms. In addition to advances in genome wide profiling and forward and reverse genetic screening, structural techniques are becoming important tools in unraveling TF function by providing atomic and molecular level information that was lacking in purely genetic approaches. Here, we summarize previous structural work and present additional biophysical and biochemical studies of the key master regulators of plant reproduction - LEAFY and the MADS-domain TFs SEPALLATA3 and AGAMOUS. We discuss the impact of structural biology on our understanding of the complex evolutionary process leading to the development of the bisexual flower.

Flower Development and Perianth Identity Candidate Genes in the Basal Angiosperm Aristolochia fimbriata (Piperales: Aristolochiaceae)

Aristolochia fimbriata (Aristolochiaceae: Piperales) exhibits highly synorganized flowers with a single convoluted structure forming a petaloid perianth that surrounds the gynostemium, putatively formed by the congenital fusion between stamens and the upper portion of the carpels. Here we present the flower development and morphology of A. fimbriata, together with the expression of the key regulatory genes that participate in flower development, particularly those likely controlling perianth identity. A. fimbriata is a member of the magnoliids, and thus gene expression detected for all ABCE MADS-box genes in this taxon, can also help to elucidate patterns of gene expression prior the independent duplications of these genes in eudicots and monocots. Using both floral development and anatomy in combination with the isolation of MADS-box gene homologs, gene phylogenetic analyses and expression studies (both by reverse transcription PCR and in situ hybridization), we present hypotheses on floral organ identity genes involved in the formation of this bizarre flower. We found that most MADS-box genes were expressed in vegetative and reproductive tissues with the exception of AfimSEP2, AfimAGL6, and AfimSTK transcripts that are only found in flowers and capsules but are not detected in leaves. Two genes show ubiquitous expression; AfimFUL that is found in all floral organs at all developmental stages as well as in leaves and capsules, and AfimAG that has low expression in leaves and is found in all floral organs at all stages with a considerable reduction of expression in the limb of anthetic flowers. Our results indicate that expression of AfimFUL is indicative of pleiotropic roles and not of a perianth identity specific function. On the other hand, expression of B-class genes, AfimAP3 and AfimPI, suggests their conserved role in stamen identity and corroborates that the perianth is sepal and not petal-derived. Our data also postulates an AGL6 ortholog as a candidate gene for sepal identity in the Aristolochiaceae and provides testable hypothesis for a modified ABCE model in synorganized magnoliid flowers.

Hughmillerites vancouverensis sp. nov. and the Cretaceous diversification of Cupressaceae

•

PREMISE OF THE STUDY

Two ovulate conifer cones, one of which is attached terminally to a short leafy shoot, reveal the presence of a new species of Hughmillerites in the Early Cretaceous Apple Bay flora of Vancouver Island, British Columbia, Canada. This ancient conifer expands the diversity of Cupressaceae in the Mesozoic and reveals details about the evolution of Subfamily: Cunninghamioideae.•

METHODS

Specimens were studied from anatomical sections prepared using the cellulose acetate peel technique.•

KEY RESULTS

Vegetative shoots have helically arranged leaves that are Cunninghamia-like. Seed cones have many helically arranged bract/scale complexes in which the bract is larger than the ovuliferous scale. Each ovuliferous scale has three free tips that separate from the bract immediately distal to an inverted seed. Several ovuliferous scales show interseminal ridges between seeds.•

CONCLUSIONS

This study documents a new extinct species of cunninghamioid conifers, Hughmillerites vancouverensis, expanding the record of the genus from the Late Jurassic to the Early Cretaceous. This new extinct species emphasizes the important role that conifers from subfamily Cunninghamioideae played in the initial evolutionary radiation of Cupressaceae. In light of recent findings in conifer regulatory genetics, we use H. vancouverensis to hypothesize that variations of expression in certain gene homologues played an important role in the evolution of the cupressaceous ovuliferous scale.

From leaf and branch into a flower: Magnolia tells the story

BACKGROUND

In the classical doctrines, Magnolia was frequently considered the archetype among flowering plants, and its conduplicate carpel with marginal placentation was assumed to be derived from a leaf-like organ bearing ovules along its margins. Although the robustness of this concept has been seriously questioned by advances in botanical research, especially the emergence of Magnolia deeper in the angiosperm tree of life in molecular systematics, it remains the most-taught interpretation for the origin of carpels.

RESULTS

To test the validity of this classical doctrine, we performed comparative anatomical analyses of the vascular bundles in the flowers of Magnolia using fine (8-μm) paraffin -sections. We document the presence of two independent vascular systems in the carpels: the collateral bundles of the dorsal and ventral veins arising from the stelar bundle, and the amphicribral ovular bundles arising from the cortical bundles. This observation in conjunction with data from other fields concurrently suggests that the ovary wall is equivalent to a foliar organ whereas the placenta represents an ovule-bearing shoot.

CONCLUSIONS

Our observation on the former model plant, Magnolia, nullifies the classical doctrine of carpel evolution and supports the Unifying Theory. This conclusion prompts a reconsideration of the concept of angiosperm flower evolution.

MADS goes genomic in conifers: towards determining the ancestral set of MADS-box genes in seed plants

BACKGROUND AND AIMS

MADS-box genes comprise a gene family coding for transcription factors. This gene family expanded greatly during land plant evolution such that the number of MADS-box genes ranges from one or two in green algae to around 100 in angiosperms. Given the crucial functions of MADS-box genes for nearly all aspects of plant development, the expansion of this gene family probably contributed to the increasing complexity of plants. However, the expansion of MADS-box genes during one important step of land plant evolution, namely the origin of seed plants, remains poorly understood due to the previous lack of whole-genome data for gymnosperms.

METHODS

The newly available genome sequences of Picea abies, Picea glauca and Pinus taeda were used to identify the complete set of MADS-box genes in these conifers. In addition, MADS-box genes were identified in the growing number of transcriptomes available for gymnosperms. With these datasets, phylogenies were constructed to determine the ancestral set of MADS-box genes of seed plants and to infer the ancestral functions of these genes.

KEY RESULTS

Type I MADS-box genes are under-represented in gymnosperms and only a minimum of two Type I MADS-box genes have been present in the most recent common ancestor (MRCA) of seed plants. In contrast, a large number of Type II MADS-box genes were found in gymnosperms. The MRCA of extant seed plants probably possessed at least 11-14 Type II MADS-box genes. In gymnosperms two duplications of Type II MADS-box genes were found, such that the MRCA of extant gymnosperms had at least 14-16 Type II MADS-box genes.

CONCLUSIONS

The implied ancestral set of MADS-box genes for seed plants shows simplicity for Type I MADS-box genes and remarkable complexity for Type II MADS-box genes in terms of phylogeny and putative functions. The analysis of transcriptome data reveals that gymnosperm MADS-box genes are expressed in a great variety of tissues, indicating diverse roles of MADS-box genes for the development of gymnosperms. This study is the first that provides a comprehensive overview of MADS-box genes in conifers and thus will provide a framework for future work on MADS-box genes in seed plants.

Evolution of fruit development genes in flowering plants

The genetic mechanisms regulating dry fruit development and opercular dehiscence have been identified in Arabidopsis thaliana. In the bicarpellate silique, valve elongation and differentiation is controlled by FRUITFULL (FUL) that antagonizes SHATTERPROOF1-2 (SHP1/SHP2) and INDEHISCENT (IND) at the dehiscence zone where they control normal lignification. SHP1/2 are also repressed by REPLUMLESS (RPL), responsible for replum formation. Similarly, FUL indirectly controls two other factors ALCATRAZ (ALC) and SPATULA (SPT) that function in the proper formation of the separation layer. FUL and SHP1/2 belong to the MADS-box family, IND and ALC belong to the bHLH family and RPL belongs to the homeodomain family, all of which are large transcription factor families. These families have undergone numerous duplications and losses in plants, likely accompanied by functional changes. Functional analyses of homologous genes suggest that this network is fairly conserved in Brassicaceae and less conserved in other core eudicots. Only the MADS box genes have been functionally characterized in basal eudicots and suggest partial conservation of the functions recorded for Brassicaceae. Here we do a comprehensive search of SHP, IND, ALC, SPT, and RPL homologs across core-eudicots, basal eudicots, monocots and basal angiosperms. Based on gene-tree analyses we hypothesize what parts of the network for fruit development in Brassicaceae, in particular regarding direct and indirect targets of FUL, might be conserved across angiosperms.

Heterochronic genes in plant evolution and development

Evolution of morphology includes evolutionary shifts of developmental processes in space or in time. Heterochronic evolution is defined as a temporal shift. The concept of heterochrony has been very rewarding to investigators of both animal and plant developmental evolution, because it has strong explanatory power when trying to understand morphological diversity. While for animals, extensive literature on heterochrony developed along with the field of evolution of development, in plants the concept has been applied less often and is less elaborately developed. Yet novel genetic findings highlight heterochrony as a developmental and evolutionary process in plants. Similar to what has been found for the worm Caenorhabditis, a heterochronic gene pathway controlling developmental timing has been elucidated in flowering plants. Two antagonistic microRNA's miR156 and miR172 target two gene families of transcription factors, SQUAMOSA PROMOTOR BINDING PROTEIN-LIKE and APETALA2-like, respectively. Here, we propose that this finding now allows the molecular investigation of cases of heterochronic evolution in plants. We illustrate this point by examining microRNA expression patterns in the Antirrhinum majus incomposita and choripetala heterochronic mutants. Some of the more beautiful putative cases of heterochronic evolution can be found outside flowering plants, but little is known about the extent of conservation of this flowering plant pathway in other land plants. We show that the expression of an APETALA2-like gene decreases with age in a fern species. This contributes to the idea that ferns share some heterochronic gene functions with flowering plants.