A eudicot MIXTA family ancestor likely functioned in both conical cells and trichomes

The MIXTA family of MYB transcription factors modulate the development of diverse epidermal features in land plants. This study investigates the evolutionary history and function of the MIXTA gene family in the early-diverging eudicot model lineage Thalictrum (Ranunculaceae), with R2R3 SBG9-A MYB transcription factors representative of the pre-core eudicot duplication and thus hereby referred to as “paleoMIXTA” (PMX). Cloning and phylogenetic analysis of Thalictrum paleoMIXTA (ThPMX) orthologs across 23 species reveal a genus-wide duplication coincident with a whole-genome duplication. Expression analysis by qPCR confirmed that the highest expression is found in carpels, while newly revealing high expression in leaves and nuanced differences between paralogs in representative polyploid species. The single-copy ortholog from the diploid species T. thalictroides (TthPMX, previously TtMYBML2), which has petaloid sepals with conical–papillate cells and trichomes on leaves, was functionally characterized by virus-induced gene silencing (VIGS), and its role in leaves was also assessed from heterologous overexpression in tobacco. Another ortholog from a species with conical–papillate cells on stamen filaments, TclPMX, was also targeted for silencing. Overexpression assays in tobacco provide further evidence that the paleoMIXTA lineage has the potential for leaf trichome function in a core eudicot. Transcriptome analysis by RNA-Seq on leaves of VIGS-treated plants suggests that TthPMX modulates leaf trichome development and morphogenesis through microtubule-associated mechanisms and that this may be a conserved pathway for eudicots. These experiments provide evidence for a combined role for paleoMIXTA orthologs in (leaf) trichomes and (floral) conical–papillate cells that, together with data from other systems, makes the functional reconstruction of a eudicot ancestor most likely as also having a combined function

Virus-Induced Gene Silencing as a Tool for Comparative Functional Studies in Thalictrum

Perennial woodland herbs in the genus Thalictrum exhibit high diversity of floral morphology, including four breeding and two pollination systems. Their phylogenetic position, in the early-diverging eudicots, makes them especially suitable for exploring the evolution of floral traits and the fate of gene paralogs that may have shaped the radiation of the eudicots. A current limitation in evolution of plant development studies is the lack of genetic tools for conducting functional assays in key taxa spanning the angiosperm phylogeny. We first show that virus-induced gene silencing (VIGS) of a PHYTOENE DESATURASE ortholog (TdPDS) can be achieved in Thalictrum dioicum with an efficiency of 42% and a survival rate of 97%, using tobacco rattle virus (TRV) vectors. The photobleached leaf phenotype of silenced plants significantly correlates with the down-regulation of endogenous TdPDS (P<0.05), as compared to controls. Floral silencing of PDS was achieved in the faster flowering spring ephemeral T. thalictroides. In its close relative, T. clavatum, silencing of the floral MADS box gene AGAMOUS (AG) resulted in strong homeotic conversions of floral organs. In conclusion, we set forth our optimized protocol for VIGS by vacuum-infiltration of Thalictrum seedlings or dormant tubers as a reference for the research community. The three species reported here span the range of floral morphologies and pollination syndromes present in Thalictrum. The evidence presented on floral silencing of orthologs of the marker gene PDS and the floral homeotic gene AG will enable a comparative approach to the study of the evolution of flower development in this group

Gene Duplication and Transference of Function in the paleoAP3 Lineage of Floral Organ Identity Genes

The floral organ identity gene APETALA3 (AP3) is a MADS-box transcription factor involved in stamen and petal identity that belongs to the B-class of the ABC model of flower development. Thalictrum (Ranunculaceae), an emerging model in the non-core eudicots, has AP3 homologs derived from both ancient and recent gene duplications. Prior work has shown that petals have been lost repeatedly and independently in Ranunculaceae in correlation with the loss of a specific AP3 paralog, and Thalictrum represents one of these instances. The main goal of this study was to conduct a functional analysis of the three AP3 orthologs present in Thalictrum thalictroides, representing the paleoAP3 gene lineage, to determine the degree of redundancy versus divergence after gene duplication. Because Thalictrum lacks petals, and has lost the petal-specific AP3, we also asked whether heterotopic expression of the remaining AP3 genes contributes to the partial transference of petal function to the first whorl found in insect-pollinated species. To address these questions, we undertook functional characterization by virus-induced gene silencing (VIGS), protein–protein interaction and binding site analyses. Our results illustrate partial redundancy among Thalictrum AP3s, with deep conservation of B-class function in stamen identity and a novel role in ectopic petaloidy of sepals. Certain aspects of petal function of the lost AP3 locus have apparently been transferred to the other paralogs. A novel result is that the protein products interact not only with each other, but also as homodimers. Evidence presented here also suggests that expression of the different ThtAP3 paralogs is tightly integrated, with an apparent disruption of B function homeostasis upon silencing of one of the paralogs that codes for a truncated protein. To explain this result, we propose two testable alternative scenarios: that the truncated protein is a dominant negative mutant or that there is a compensational response as part of a back-up circuit. The evidence for promiscuous protein–protein interactions via yeast two-hybrid combined with the detection of AP3 specific binding motifs in all B-class gene promoters provide partial support for these hypotheses

Ferns:the missing link in shoot evolution and development

Shoot development in land plants is a remarkably complex process that gives rise to an extreme diversity of forms. Our current understanding of shoot developmental mechanisms comes almost entirely from studies of angiosperms (flowering plants), the most recently diverged plant lineage. Shoot development in angiosperms is based around a layered multicellular apical meristem that produces lateral organs and/or secondary meristems from populations of founder cells at its periphery. In contrast, non-seed plant shoots develop from either single apical initials or from a small population of morphologically distinct apical cells. Although developmental and molecular information is becoming available for non-flowering plants, such as the model moss Physcomitrella patens, making valid comparisons between highly divergent lineages is extremely challenging. As sister group to the seed plants, the monilophytes (ferns and relatives) represent an excellent phylogenetic midpoint of comparison for unlocking the evolution of shoot developmental mechanisms, and recent technical advances have finally made transgenic analysis possible in the emerging model fern Ceratopteris richardii. This review compares and contrasts our current understanding of shoot development in different land plant lineages with the aim of highlighting the potential role that the fern C. richardii could play in shedding light on the evolution of underlying genetic regulatory mechanisms

Arabidopsis and Tobacco SUPERMAN regulate hormone signalling and mediate cell proliferation and differentiation

Arabidopsis thaliana SUPERMAN (SUP) plays an important role during flower development by maintaining the boundary between stamens and carpels in the inner two whorls. It was proposed that SUP maintains this boundary by regulating cell proliferation in both whorls, as loss-of-function superman mutants produce more stamens at the expense of carpels. However, the cellular mechanism that underlies SUP function remains unknown. Here Arabidopsis or tobacco (Nicotiana tabacum) SUP was overexpressed in tobacco plants to substantiate SUP's role as a regulator of cell proliferation and boundary definition and provide evidence that its biological role may be mediated via hormonal changes. It was found that moderate levels of SUP stimulated cell growth and proliferation, whereas high levels were inhibitory. SUP stimulated auxin- and cytokinin-regulated processes, and cells overexpressing SUP displayed reduced hormone dependency for proliferation and regeneration into plants. SUP also induced proliferation of female traits in the second and third flower whorls and promoted differentiation of petaloid properties in sepals, further supporting a role for SUP as a boundary regulator. Moreover, cytokinin suppressed stamen development and promoted differentiation of carpeloid tissues, suggesting that SUP may regulate male and female development via its effect on cytokinin signalling. Taken together, these observations suggest a model whereby the effect of SUP on cell growth and proliferation involves the modulation of auxin- and cytokinin-regulated processes. Furthermore, differential SUP expression or different sensitivities of different cell types to SUP may determine whether SUP stimulates or suppresses their proliferation

Gene Duplication and Differential Expression of Flower Symmetry Genes in Rhododendron (Ericaceae)

Bilaterally symmetric flowers have evolved over a hundred times in angiosperms, yet orthologs of the transcription factors CYCLOIDEA (CYC), RADIALIS (RAD), and DIVARICATA (DIV) are repeatedly implicated in floral symmetry changes. We examined these candidate genes to elucidate the genetic underpinnings of floral symmetry changes in florally diverse Rhododendron, reconstructing gene trees and comparing gene expression across floral organs in representative species with radial and bilateral flower symmetries. Radially symmetric R. taxifolium Merr. and bilaterally symmetric R. beyerinckianum Koord. had four and five CYC orthologs, respectively, from shared tandem duplications. CYC orthologs were expressed in the longer dorsal petals and stamens and highly expressed in R. beyerinckianum pistils, whereas they were either ubiquitously expressed, lost from the genome, or weakly expressed in R. taxifolium. Both species had two RAD and DIV orthologs uniformly expressed across all floral organs. Differences in gene structure and expression of Rhododendron RAD compared to other asterids suggest that these genes may not be regulated by CYC orthologs. Our evidence supports CYC orthologs as the primary regulators of differential organ growth in Rhododendron flowers, while also suggesting certain deviations from the typical asterid gene regulatory network for flower symmetry

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

Data from: Fluidigm2PURC: automated processing and haplotype inference for double-barcoded PCR amplicons

Premise of the Study: Targeted enrichment strategies for phylogenomic inference are a time‐ and cost‐efficient way to collect DNA sequence data for large numbers of individuals at multiple, independent loci. Automated and reproducible processing of these data is a crucial step for researchers conducting phylogenetic studies. Methods and Results: We present Fluidigm2PURC, an open source Python utility for processing paired‐end Illumina data from double‐barcoded PCR amplicons. In combination with the program PURC (Pipeline for Untangling Reticulate Complexes), our scripts process raw FASTQ files for analysis with PURC and use its output to infer haplotypes for diploids, polyploids, and samples with unknown ploidy. We demonstrate the use of the pipeline with an example data set from the genus Thalictrum (Ranunculaceae). Conclusions: Fluidigm2PURC is freely available for Unix‐like operating systems on GitHub (https://github.com/pblischak/fluidigm2purc) and for all operating systems through Docker (https://hub.docker.com/r/pblischak/fluidigm2purc)

Downregulation of <i>TdPDS</i> and detection of TRV transcripts in VIGS photobleached leaves of <i>Thalictrum dioicum</i>.

<p><b>A: Expression of </b><b><i>TdACTIN</i></b><b> control, native </b><b><i>TdPDS</i></b><b> and viral transcripts in leaves by Reverse Transcriptase (RT)-PCR.</b> Untreated and mock-treated (empty TRV2) <i>T. dioicum</i> plants are compared to TRV2-<i>TdPDS</i> treated plants showing photobleached (white), variegated (green/white) and green leaf tissue. RT-PCR was performed with locus-specific primers to the housekeeping gene <i>ACTIN</i> (loading control), to endogenous <i>TdPDS</i> and to the viral transcripts TRV1/TRV2. Approximate band size indicated for TRV2: larger band results from the presence of the <i>TdPDS</i> insert, smaller band from an empty TRV2 (mock control). <b>B: Comparative expression of </b><b><i>TdPDS</i></b><b> normalized with </b><b><i>TdACTIN</i></b><b> among treatments and resulting phenotypes of </b><b><i>Thalictrum dioicum</i></b><b>.</b> Values based on quantification of RT-PCR gel bands in part A using ImageJ (see text for details). Different letters indicate statistical significance in a one-way ANOVA followed by Tuckey test (p<0.05), same letters indicate no statistical difference. Average and standard error bars are shown. Sample sizes are n = 4 for untreated and mock-treated, n = 6 for treated bleached and n = 3 for treated variegated or green.</p