Developmental Genetics of Corolla Tube Formation: Role of the tasiRNA-ARF Pathway and a Conceptual Model

Differential Gene Expression and Unbalanced Homeolog Expression Bias in 4 Million-Year-Old Allopolyploids of Nicotiana Section Repandae

Allopolyploidy, a phenomenon prevalent in angiosperms involving hybridization and whole-genome duplication, results in species with multiple subgenomes, altering genome structure and gene expression, leading to novel phenotypes. Allopolyploids often experience unbalanced homeolog expression bias, the preferential expression of homeologs from one of the two progenitor genomes. To explore the consequences of allopolyploidy and unbalanced homeolog expression bias, we investigate global gene expression and the fate of homeologs in Nicotiana (Solanaceae). We focus on Nicotiana section Repandae, including three allotetraploid species, Nicotiana nudicaulis, N. repanda, and N. stocktonii, derived from diploid progenitors N. sylvestris and N. obtusifolia ∼4.3 Ma. We identify genes with differential expression and investigate expression of candidate genes for flower size variation. Our results show expression differences with the allopolyploids intermediate between the two progenitor species, with a slight bias toward N. obtusifolia. Moreover, we demonstrate unbalanced homeolog expression bias toward the N. obtusifolia subgenome across developmental stages in the allopolyploids, with a stronger bias in N. nudicaulis. In contrast, unbalanced homeolog expression bias shifts toward N. sylvestris for flower size genes in N. nudicaulis, showing that genes involved in particular phenotypes can display different patterns of unbalanced homeolog expression than the overall transcriptome. We also see differential expression of several known flower size genes across corolla developmental stages. Our results highlight the role of unbalanced homeolog expression bias in shaping the evolutionary trajectory of Nicotiana species and provide a foundation for future research into the ecological and evolutionary implications of allopolyploidy in flowering plants.

Building beauty: Understanding how hormone signaling regulates petal patterning and morphogenesis

The corolla of flowering plants provides pivotal functions for the reproduction of angiosperms, directly impacting the fitness of individuals. Different petal shapes and patterns contribute to these functions and, thus, participate in the production of morphological diversity and the emergence of new species. During petal morphogenesis, the coordination of cell fate specification, cell division, and cell expansion is coherent and robust across the petal blade and is set according to proximo-distal, medio-lateral, and abaxial-adaxial axes. However, the mechanisms specifying petal polarity and controlling cell behavior in a position-dependent manner as petals develop remain poorly understood. In this review, we draw parallels with other evolutionarily related plant lateral organs such as leaves to argue that hormones likely play central, yet largely unexplored, roles in such coordination. By examining petal development in Arabidopsis and other angiosperms, we frame what are the knowns and the unknowns of hormones contributions to petal morphogenesis and patterning. Finally, we argue that using emerging model organisms can provide invaluable information to tackle questions that have long remained unanswered, broadening our understanding by allowing us to investigate petal morphogenesis and the tinkering of phytohormone signaling through an evolutionary lens.

BLADE-ON-PETIOLE interacts with CYCLOIDEA to fine-tune CYCLOIDEA-mediated flower symmetry in monkeyflowers (Mimulus)

Morphological novelties, or key innovations, are instrumental to the diversification of the organisms. In plants, one such innovation is the evolution of zygomorphic flowers, which is thought to promote outcrossing and increase flower morphological diversity. We isolated three allelic mutants from two Mimulus species displaying altered floral symmetry and identified the causal gene as the ortholog of Arabidopsis BLADE-ON-PETIOLE. We found that MlBOP and MlCYC2A physically interact and this BOP-CYC interaction module is highly conserved across the angiosperms. Furthermore, MlBOP self-ubiquitinates and suppresses MlCYC2A self-activation. MlCYC2A, in turn, impedes MlBOP ubiquitination. Thus, this molecular tug-of-war between MlBOP and MlCYC2A fine-tunes the expression of MlCYC2A, contributing to the formation of bilateral symmetry in flowers, a key trait in angiosperm evolution.

CYCLOIDEA paralogs function partially redundantly to specify dorsal flower development in Mimulus lewisii

PREMISE

Duplicated genes (paralogs) are abundant in plant genomes, and their retention may influence the function of genetic programs and contribute to evolutionary novelty. How gene duplication affects genetic modules and what forces contribute to paralog retention are outstanding questions. The CYCLOIDEA(CYC)-dependent flower symmetry program is a model for understanding the evolution of gene duplication, providing multiple examples of paralog partitioning and novelty. However, a novel CYC gene lineage duplication event near the origin of higher core Lamiales (HCL) has received little attention.

METHODS

To understand the evolutionary fate of duplicated HCL CYC2 genes, we determined the effects on flower symmetry by suppressing MlCYC2A and MlCYC2B expression using RNA interference (RNAi). We determined the phenotypic effects on flower symmetry in single- and double-silenced backgrounds and coupled our functional analyses with expression surveys of MlCYC2A, MlCYC2B, and a putative downstream RADIALIS (MlRAD5) ortholog.

RESULTS

MlCYC2A and MlCYC2B jointly contribute to bilateral flower symmetry. MlCYC2B exhibits a clear dorsal flower identity function and may additionally function in carpel development. MlCYC2A functions in establishing dorsal petal shape. Further, our results suggest an MlCYC2A-MlCYC2B regulatory interaction, which may affect pathway homeostasis.

CONCLUSIONS

Our results suggest that CYC paralogs specific to higher core Lamiales may be selectively retained for their joint contribution to bilateral flower symmetry, similar to the independently derived CYC paralogs in the Lamiales model for bilateral flower symmetry research, Antirrhinum majus (snapdragon).

A BLADE-ON-PETIOLE orthologue regulates corolla differentiation in the proximal region in Torenia fournieri

The three-dimensional shape of a flower is integrated by morphogenesis along different axes. Differentiation along the petal proximodistal axis is tightly linked to the specification of pollinators; however, it is still unclear how a petal patterns this axis. The corolla of Torenia fournieri exhibits strong differentiation along the proximodistal axis, and we previously found a proximal regulator, TfALOG3, controlling corolla neck differentiation. Here, we report another gene, TfBOP2, which is predominantly expressed in the proximal region of the corolla. TfBOP2 mutants have shorter proximal corolla tubes and longer distal lobe, demonstrating its function as a proximal regulator. Arabidopsis BOPs mutant shows similar defects, favouring a shared role of BOPs homologues. Genetic analysis demonstrates the interaction between TfBOP2 and TfALOG3, and we further found that TfALOG3 physically interacts with TfBOP2 and can recruit TfBOP2 to the nuclear region. Our study favours a hypothetical shared BOP-ALOG complex that is recruited to regulate corolla differentiation in the proximal region axis of T. fournieri.

An insight into the roles of regulatory ncRNAs in plants: An abiotic stress and developmental perspective

Different environmental cues lead to changes in physiology, biochemistry and molecular status of plant's growth. Till date, various genes have been accounted for their role in regulating plant development and response to abiotic stress. Excluding genes that code for a functional protein in a cell, a large chunk of the eukaryotic transcriptome consists of non-coding RNAs (ncRNAs) which lack protein coding capacity but are still functional. Recent advancements in Next Generation Sequencing (NGS) technology have led to the unearthing of different types of small and large non-coding RNAs in plants. Non-coding RNAs are broadly categorised into housekeeping ncRNAs and regulatory ncRNAs which work at transcriptional, post-transcriptional and epigenetic levels. Diverse ncRNAs play different regulatory roles in nearly all biological processes including growth, development and response to changing environments. This response can be perceived and counteracted by plants using diverse evolutionarily conserved ncRNAs like miRNAs, siRNAs and lncRNAs to participate in complex molecular regimes by activating gene-ncRNA-mRNA regulatory modules to perform the downstream function. Here, we review the current understanding with a focus on recent advancements in the functional studies of the regulatory ncRNAs at the nexus of abiotic stresses and development. Also, the potential roles of ncRNAs in imparting abiotic stress tolerance and yield improvement in crop plants are also discussed with their future prospects.

Floral Development Stage-Specific Transcriptomic Analysis Reveals the Formation Mechanism of Different Shapes of Ray Florets in Chrysanthemum

The formation mechanism of different ray floret shapes of chrysanthemum (Chrysanthemum × morifolium) remains elusive due to its complex genetic background. C. vestitum, with the basic ray floret shapes of the flat, spoon, and tubular types, is considered a model material for studying ray floret morphogenesis. In this study, the flat and tubular type lines of C. vestitum at specific stages were used to investigate the key genes that regulate morphological differences in ray florets. We found that the expression levels of genes related to auxin synthesis, transport, and response were generally higher in the tubular type than in the flat type. CvARF3 was highly expressed in the flat type, while CvARF5 and CvARF6 were highly expressed in the tubular type. Additionally, the transcription levels of Class B and E genes closely related to petal development, including CvPI, CvAP3, Cvdefh21, CvSEP3, and CvCDM77, were expressed at higher levels in the tubular type than the flat type. Based on the results, it is proposed that auxin plays a key role in the development of ray florets, and auxin-related genes, especially CvARFs, may be key genes to control the morphological difference of ray florets. Simultaneously, MADS-box genes are involved in the co-regulation of ray floret morphogenesis. The results provide novel insights into the molecular mechanism of different petal type formation and lay a theoretical foundation for the directional breeding of petal type in chrysanthemums.

Evolutionarily inspired solutions to the crop pollination crisis

The global decline in insect diversity threatens pollination services, potentially impacting crop production and food security. Here, we argue that this looming pollination crisis is generally approached from an ecological standpoint, and that consideration of evolutionary principles offers a novel perspective. First, we outline that wild plant species have overcome 'pollination crises' throughout evolutionary history, and show how associated principles can be applied to crop pollination. We then highlight technological advances that can be used to adapt crop flowers for optimal pollination by local wild pollinators, especially by increasing generalization in pollination systems. Thus, synergies among fundamental evolutionary research, genetic engineering, and agro-ecological science provide a promising template for addressing a potential pollination crisis, complementing much-needed strategies focused on pollinator conservation.

Testing candidate genes linked to corolla shape variation of a pollinator shift in Rhytidophyllum (Gesneriaceae)

Floral adaptations to specific pollinators like corolla shape variation often result in reproductive isolation and thus speciation. But despite their ecological importance, the genetic bases of corolla shape transitions are still poorly understood, especially outside model species. Hence, our goal was to identify candidate genes potentially involved in corolla shape variation between two closely related species of the Rhytidophyllum genus (Gesneriaceae family) from the Antilles with contrasting pollination strategies. Rhytidophyllum rupincola has a tubular corolla and is strictly pollinated by hummingbirds, whereas R. auriculatum has more open flowers and is pollinated by hummingbirds, bats, and insects. We surveyed the literature and used a comparative transcriptome sequence analysis of synonymous and non-synonymous nucleotide substitutions to obtain a list of genes that could explain floral variation between R. auriculatum and R. rupincola. We then tested their association with corolla shape variation using QTL mapping in a F2 hybrid population. Out of 28 genes tested, three were found to be good candidates because of a strong association with corolla shape: RADIALIS, GLOBOSA, and JAGGED. Although the role of these genes in Rhytidophyllum corolla shape variation remains to be confirmed, these findings are a first step towards identifying the genes that have been under selection by pollinators and thus involved in reproductive isolation and speciation in this genus.

Genome-Wide Characterization, Evolutionary Analysis of ARF Gene Family, and the Role of SaARF4 in Cd Accumulation of Sedum alfredii Hance

Auxin response factors (ARFs) play important roles in plant development and environmental adaption. However, the function of ARFs in cadmium (Cd) accumulation are still unknown. Here, 23 SaARFs were detected in the genome of hyperaccumulating ecotype of Sedum alfredii Hance (HE), and they were not evenly distributed on the chromosomes. Their protein domains remained highly conservative. SaARFs in the phylogenetic tree can be divided into three groups. Genes in the group Ⅰ contained three introns at most. However, over ten introns were found in other two groups. Collinearity relationships were exhibited among ten SaARFs. The reasons for generating SaARFs may be segmental duplication and rearrangements. Collinearity analysis among different species revealed that more collinear genes of SaARFs can be found in the species with close relationships of HE. A total of eight elements in SaARFs promoters were related with abiotic stress. The qRT-PCR results indicated that four SaARFs can respond to Cd stress. Moreover, that there may be functional redundancy among six SaARFs. The adaptive selection and functional divergence analysis indicated that SaARF4 may undergo positive selection pressure and an adaptive-evolution process. Overexpressing SaARF4 effectively declined Cd accumulation. Eleven single nucleotide polymorphism (SNP) sites relevant to Cd accumulation can be detected in SaARF4. Among them, only one SNP site can alter the sequence of the SaARF4 protein, but the SaARF4 mutant of this site did not cause a significant difference in cadmium content, compared with wild-type plants. SaARFs may be involved in Cd-stress responses, and SaARF4 may be applied for decreasing Cd accumulation of plants.

Plant and animal small RNA communications between cells and organisms

Since the discovery of eukaryotic small RNAs as the main effectors of RNA interference in the late 1990s, diverse types of endogenous small RNAs have been characterized, most notably microRNAs, small interfering RNAs (siRNAs) and PIWI-interacting RNAs (piRNAs). These small RNAs associate with Argonaute proteins and, through sequence-specific gene regulation, affect almost every major biological process. Intriguing features of small RNAs, such as their mechanisms of amplification, rapid evolution and non-cell-autonomous function, bestow upon them the capacity to function as agents of intercellular communications in development, reproduction and immunity, and even in transgenerational inheritance. Although there are many types of extracellular small RNAs, and despite decades of research, the capacity of these molecules to transmit signals between cells and between organisms is still highly controversial. In this Review, we discuss evidence from different plants and animals that small RNAs can act in a non-cell-autonomous manner and even exchange information between species. We also discuss mechanistic insights into small RNA communications, such as the nature of the mobile agents, small RNA signal amplification during transit, signal perception and small RNA activity at the destination.

Transcriptomic Analysis Suggests Auxin Regulation in Dorsal-Ventral Petal Asymmetry of Wild Progenitor Sinningia speciosa

The establishment of dorsal-ventral (DV) petal asymmetry is accompanied by differential growth of DV petal size, shape, and color differences, which enhance ornamental values. Genes involved in flower symmetry in Sinningia speciosa have been identified as CYCLOIDEA (SsCYC), but which gene regulatory network (GRN) is associated with SsCYC to establish DV petal asymmetry is still unknown. To uncover the GRN of DV petal asymmetry, we identified 630 DV differentially expressed genes (DV-DEGs) from the RNA-Seq of dorsal and ventral petals in the wild progenitor, S. speciosa 'ES'. Validated by qRT-PCR, genes in the auxin signaling transduction pathway, SsCYC, and a major regulator of anthocyanin biosynthesis were upregulated in dorsal petals. These genes correlated with a higher endogenous auxin level in dorsal petals, with longer tube length growth through cell expansion and a purple dorsal color. Over-expression of SsCYC in Nicotiana reduced petal size by regulating cell growth, suggesting that SsCYC also controls cell expansion. This suggests that auxin and SsCYC both regulate DV petal asymmetry. Transiently over-expressed SsCYC, however, could not activate most major auxin signaling genes, suggesting that SsCYC may not trigger auxin regulation. Whether auxin can activate SsCYC or whether they act independently to regulate DV petal asymmetry remains to be explored in the future.

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA-target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA-target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA-target interactions.

Nuclear and cytoplasmic RNA exosomes and PELOTA1 prevent miRNA-induced secondary siRNA production in Arabidopsis

Amplification of short interfering RNA (siRNAs) via RNA-dependent RNA polymerases (RdRPs) is of fundamental importance in RNA silencing. Plant microRNA (miRNA) action generally does not involve engagement of RdRPs, in part thanks to a poorly understood activity of the cytoplasmic exosome adaptor SKI2. Here, we show that inactivation of the exosome subunit RRP45B and SKI2 results in similar patterns of miRNA-induced siRNA production. Furthermore, loss of the nuclear exosome adaptor HEN2 leads to secondary siRNA production from miRNA targets largely distinct from those producing siRNAs in ski2. Importantly, mutation of the Release Factor paralogue PELOTA1 required for subunit dissociation of stalled ribosomes causes siRNA production from miRNA targets overlapping with, but distinct from, those affected in ski2 and rrp45b mutants. We also show that in exosome mutants, miRNA targets can be sorted into producers and non-producers of illicit secondary siRNAs based on trigger miRNA levels and miRNA:target affinity rather than on presence of 5′-cleavage fragments. We propose that stalled RNA-Induced Silencing Complex (RISC) and ribosomes, but not mRNA cleavage fragments released from RISC, trigger siRNA production, and that the exosome limits siRNA amplification by reducing RISC dwell time on miRNA target mRNAs while PELOTA1 does so by reducing ribosome stalling.

MicroRNAs Are Involved in Regulating Plant Development and Stress Response through Fine-Tuning of TIR1/AFB-Dependent Auxin Signaling

Auxin, primarily indole-3-acetic acid (IAA), is a versatile signal molecule that regulates many aspects of plant growth, development, and stress response. Recently, microRNAs (miRNAs), a type of short non-coding RNA, have emerged as master regulators of the auxin response pathways by affecting auxin homeostasis and perception in plants. The combination of these miRNAs and the autoregulation of the auxin signaling pathways, as well as the interaction with other hormones, creates a regulatory network that controls the level of auxin perception and signal transduction to maintain signaling homeostasis. In this review, we will detail the miRNAs involved in auxin signaling to illustrate its in planta complex regulation.

The leaf polarity factors SGS3 and YABBYs regulate style elongation through auxin signaling in Mimulus lewisii

Style length is a major determinant of breeding strategies in flowering plants and can vary dramatically between and within species. However, little is known about the genetic and developmental control of style elongation. We characterized the role of two classes of leaf adaxial-abaxial polarity factors, SUPPRESSOR OF GENE SILENCING3 (SGS3) and the YABBY family transcription factors, in the regulation of style elongation in Mimulus lewisii. We also examined the spatiotemporal patterns of auxin response during style development. Loss of SGS3 function led to reduced style length via limiting cell division, and downregulation of YABBY genes by RNA interference resulted in shorter styles by decreasing both cell division and cell elongation. We discovered an auxin response minimum between the stigma and ovary during the early stages of pistil development that marks style differentiation. Subsequent redistribution of auxin response to this region was correlated with style elongation. Auxin response was substantially altered when both SGS3 and YABBY functions were disrupted. We suggest that auxin signaling plays a central role in style elongation and that the way in which auxin signaling controls the different cell division and elongation patterns underpinning natural style length variation is a major question for future research.

Comparative analysis of corolla tube development across three closely related Mimulus species with different pollination syndromes

Fusion of petals to form a corolla tube is considered a key innovation contributing to the diversification of many flowering plant lineages. Corolla tube length often varies dramatically among species and is a major determinant of pollinator preference. However, our understanding of the developmental dynamics underlying corolla tube length variation is very limited. Here we examined corolla tube growth in the Mimulus lewisii species complex, an emerging model system for studying the developmental genetics and evo-devo of pollinator-associated floral traits. We compared developmental and cellular processes associated with corolla tube length variation among the bee-pollinated M. lewisii, the hummingbird-pollinated Mimulus verbenaceus, and the self-pollinated Mimulus parishii. We found that in all three species, cell size is non-uniformly distributed along the mature tube, with the longest cells just distal to the stamen insertion site. Differences in corolla tube length among the three species are not associated with processes of organogenesis or early development but are associated with variation in multiple processes occurring later in development, including the location and duration of cell division and cell elongation. The tube growth curves of the small-flowered M. parishii and large-flowered M. lewisii are essentially indistinguishable, except that M. parishii tubes stop growing earlier at a smaller size, suggesting a critical role of heterochrony in the shift from outcrossing to selfing. These results not only highlight the developmental process associated with corolla tube variation among species but also provide a baseline reference for future developmental genetic analyses of mutants or transgenic plants with altered corolla tube morphology in this emerging model system.