Distinct subfunctionalization and neofunctionalization of the B-class MADS-box genes in Physalis floridana

Plant Nuclear Factor Y (NF-Y) Transcription Factors: Evolving Insights into Biological Functions and Gene Expansion

Gene expansion is a common phenomenon in plant transcription factor families; however, the underlying molecular mechanisms remain elusive. Examples of gene expansion in transcription factors are found in all eukaryotes. One example is plant nuclear factor Y (NF-Y) transcription factors. NF-Y is ubiquitous to eukaryotes and comprises three independent protein families: NF-YA, NF-YB, and NF-YC. While animals and fungi mostly have one of each NF-Y subunit, NF-Y is greatly expanded in plants. For example, humans have one each of NF-YA, NF-YB, and NF-YC, while the model plant Arabidopsis has ten each of NF-YA, NF-YB, and NF-YC. Our understanding of the plant NF-Y, including its biological roles, molecular mechanisms, and gene expansion, has improved over the past few years. Here we will review its biological roles and focus on studies demonstrating that NF-Y can serve as a model for plant gene expansion. These studies show that NF-Y can be classified into ancestrally related subclasses. Further, the primary structure of each NF-Y contains a conserved core domain flanked by non-conserved N- and C-termini. The non-conserved N- and C-termini, under pressure for diversifying selection, may provide clues to this gene family's retention and functional diversification following gene duplication. In summary, this review demonstrates that NF-Y expansion has the potential to be used as a model to study the gene expansion and retention of transcription factor families.

Genome-wide analysis of the MADS-box gene family of sea buckthorn (Hippophae rhamnoides ssp. sinensis) and their potential role in floral organ development

Sea buckthorn (Hippophae rhamnoides ssp. sinensis) is a deciduous shrub or small tree in the Elaeagnaceae family. It is dioecious, featuring distinct structures in female and male flowers. The MADS-box gene family plays a crucial role in flower development and differentiation of floral organs in plants. However, systematic information on the MADS-box family in sea buckthorn is currently lacking. This study presents a genome-wide survey and expression profile of the MADS-box family of sea buckthorn. We identified 92 MADS-box genes in the H. rhamnoides ssp. Sinensis genome. These genes are distributed across 12 chromosomes and classified into Type I (42 genes) and Type II (50 genes). Based on the FPKM values in the transcriptome data, the expression profiles of HrMADS genes in male and female flowers of sea buckthorn showed that most Type II genes had higher expression levels than Type I genes. This suggesting that Type II HrMADS may play a more significant role in sea buckthorn flower development. Using the phylogenetic relationship between sea buckthorn and Arabidopsis thaliana, the ABCDE model genes of sea buckthorn were identified and some ABCDE model-related genes were selected for qRT-PCR analysis in sea buckthorn flowers and floral organs. Four B-type genes may be involved in the identity determination of floral organs in male flowers, and D-type genes may be involved in pistil development. It is hypothesized that ABCDE model genes may play an important role in the identity of sea buckthorn floral organs. This study analyzed the role of MADS-box gene family in the development of flower organs in sea buckthorn, which provides an important theoretical basis for understanding the regulatory mechanism of sex differentiation in sea buckthorn.

Transcriptomic comparison sheds new light on regulatory networks for dimorphic flower development in response to photoperiod in Viola prionantha

BACKGROUND

Chasmogamous (CH)-cleistogamous (CL) dimorphic flowers are developed in Viola prionantha. However, the environmental and genetic factors necessary for the CH-CL transition are unknown.

RESULTS

In the present work, short-day (SD) conditions induced CH flowers, whereas long days (LDs) triggered CL flowers in V. prionantha. Compared to fully developed CH flowers, CL flowers had less mature stamens, no nectar glands, and immature petals. Comparative transcriptomics revealed differentially expressed genes (DEGs) during CL and CH development. Core genes in the photoperiod pathway, such as V. prionantha orthologs of GIGANTEA (GI), CONSTANS (CO), and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1), which promote floral induction, were highly expressed in CL flowers, whereas UNUSUAL FLORAL ORGANS (UFO) and B-class MADS-box genes for floral organ identity and development showed an opposite alteration. Moreover, genes in the glycolytic process, sucrose metabolic process, and fatty acid biosynthetic process were all highly expressed in CH flowers. Interestingly, V. prionantha orthologs of the B-class MADS-box genes APETALA3 (AP3) and PISTILLATA (PI) might relate to these sugar-fatty acid processes and were co-expressed with GAIP-B-like and YABBY5 (YAB5), which regulate the development of the petal, stamen, and nectary. Compared to CH flowers, DEGs and hub genes in the most significantly correlated modules of the gene co-expression network, which are involved in abiotic and biotic responses, were upregulated in CL flowers.

CONCLUSIONS

We proposed an integrative model for transcription regulation of genes in the photoperiod pathway, floral organ development, stress response, and sugar-fatty acid processes to determine CH-CL flower development in V. prionantha. Particularly, under LDs, activated GI may induce genes involved in the stress-response pathways, and then downregulated AP3 and PI or UFO to inhibit the sugar-fatty acid metabolic processes, together forming CL flowers. In contrast, CH flowers were produced under SDs. This work provides novel insights into the developmental evolution of dimorphic flowers in Viola.

Cloning, Expression, and Tobacco Overexpression Analyses of a PISTILLATA/GLOBOSA-like (OfGLO1) Gene from Osmanthus fragrans

GLOBOSA (GLO), a B-class MADS-box gene, is involved in floral organ determination but has rarely been studied in Osmanthus fragrans, which is a very popular ornamental tree species in China. Here, the full-length cDNA of a homologous GLO1 gene (named OfGLO1) was cloned from a flower bud of O. fragrans using the RACE technique. The OfGLO1 has a 645 bp open reading frame, encoding 214 amino acids. Similar to other PI/GLO proteins, OfGLO1 has two conserved domains, MADS MEF2-like and K-box, and a 16-amino-acid PI motif in the C terminal region. Our phylogeny analysis classified OfGLO1 as a PI-type member of the B-class MADS-box gene family. The qRT-PCR assay showed that the expression of OfGLO1 in O. fragrans was continuously upregulated from the tight bud stage to the full flowering stage but barely expressed in the pistils, sepals, and non-floral organs, such as root, leaf, and stem. The genetic effect of OfGLO1 was assayed by ectopic expression in tobacco plants. Compared with the wild-type, OfGLO1 transformants showed reduced plant size, earlier flowering, shorter stamens, and lower seed setting rates. Furthermore, some stamens were changed into petal-like structures. These findings indicate that OfGLO1 plays an important role in the regulation of flower development. This study improved our understanding of class B gene function in woody plants.

Physalis floridana CRABS CLAW mediates neofunctionalization of GLOBOSA genes in carpel development

Floral B-function MADS-box genes, such as GLOBOSA (GLO), function in corolla and stamen organ identity specification. The functions of these genes outside these floral whorls are rarely reported. DOLL1 is a GLO gene controlling corolla and androecium organ identity. In this study we found that, in Physalis floridana double-layered-lantern 1 (doll1) mutant pollinated with wild-type pollen, fruit set was extremely low, indicating that doll1 females are dysfunctional. Stigma and style structure, stigma receptivity, pollen tube guidance, and embryo sac development were also impaired in doll1. P. floridana CRABS CLAW (PFCRC), predominantly expressed in carpels, was repressed in doll1 native carpels. Loss-of-function of PFCRC altered carpel meristem determinacy, carpel closure, and ovule number, and the resultant 'pistil' consisted of multiple spirally-arranged dorsiventral carpels occasionally with 1-2 naked ovules on the margin and trichomes at each mutated carpel tip, implying an alteration of carpel organ identity. Regulatory and genetic interactions between B-class MADS-box genes and PFCRC were revealed in a context-dependent manner in floral development. Our work reveals a new role for the B-function genes in carpel and ovule development via regulating PFCRC, providing a new understanding of genetic regulatory networks between MADS-domain and CRC transcription factors in mediating carpel organ specification, functionality, and origin.

Isolation and Characterization of APETALA3 Orthologs and Promoters from the Distylous Fagopyrum esculentum

Common buckwheat (Fagopyrum esculentum) produces distylous flowers with undifferentiated petaloid tepals, which makes it obviously different from flowers of model species. In model species Arabidopsis, APETALA3 (AP3) is expressed in petal and stamen and specifies petal and stamen identities during flower development. Combining with our previous studies, we found that small-scale gene duplication (GD) event and alternative splicing (AS) of common buckwheat AP3 orthologs resulted in FaesAP3_1, FaesAP3_2 and FaesAP3_2a. FaesAP3_2 and FaesAP3_2a were mainly expressed in the stamen of thrum and pin flower. Promoters functional analysis suggested that intense GUS staining was observed in the whole stamen in pFaesAP3_2::GUS transgenic Arabidopsis, while intense GUS staining was observed only in the filament of stamen in pFaesAP3_1::GUS transgenic Arabidopsis. These suggested that FaesAP3_1 and FaesAP3_2 had overlapping functions in specifying stamen filament identity and work together to determine normal stamen development. Additionally, FaesAP3_2 and FaesAP3_2a owned the similar ability to rescue stamen development of Arabidopsis ap3-3 mutant, although AS resulted in a frameshift mutation and consequent omission of the complete PI-derived motif and euAP3 motif of FaesAP3_2a. These suggested that the MIK region of AP3-like proteins was crucial for determining stamen identity, while the function of AP3-like proteins in specifying petal identity was gradually obtained after AP3 Orthologs acquiring a novel C-terminal euAP3 motif during the evolution of core eudicots. Our results also provide a clue to understanding the early evolution of the functional specificity of euAP3-type proteins involving in floral organ development in core eudicots, and also suggested that FaesAP3_2 holds the potential application for biotechnical engineering to develop a sterile male line of F. esculentum.

Differential Evolution of α-Glucan Water Dikinase (GWD) in Plants

The alpha-glucan water dikinase (GWD) enzyme catalyzes starch phosphorylation, an integral step in transitory starch degradation. The high phosphate content in stored starch has great industrial value, due to its physio–chemical properties making it more versatile, although the phosphate content of stored starch varies depending on the botanical source. In this study, we used various computational approaches to gain insights into the evolution of the GWD protein in 48 plant species with possible roles in enzyme function and alteration of phosphate content in their stored starch. Our analyses identified deleterious mutations, particularly in the highly conserved 5 aromatic amino acid residues in the dual tandem carbohydrate binding modules (CBM-45) of GWD protein in C. zofingiensis, G. hirsutum, A. protothecoides, P. miliaceum, and C. reinhardtii. These findings will inform experimental designs for simultaneous repression of genes coding for GWD and the predicted interacting proteins to elucidate the role this enzyme plays in starch degradation. Our results reveal significant diversity in the evolution of GWD enzyme across plant species, which may be evolutionarily advantageous according to the varying needs for phosphorylated stored starch between plants and environments.

Sequence analysis and protein interactions of Arabidopsis CIA2 and CIL proteins

BACKGROUND

A previous screening of Arabidopsis thaliana for mutants exhibiting dysfunctional chloroplast protein transport identified the chloroplast import apparatus (cia) gene. The cia2 mutant has a pale green phenotype and reduced rate of protein import into chloroplasts, but leaf shape and size are similar to wild-type plants of the same developmental stage. Microarray analysis showed that nuclear CIA2 protein enhances expression of the Toc75, Toc33, CPN10 and cpRPs genes, thereby up-regulating protein import and synthesis efficiency in chloroplasts. CIA2-like (CIL) shares 65% sequence identity to CIA2, suggesting that CIL and CIA2 are homologous proteins in Arabidopsis. Here, we further assess the protein interactions and sequence features of CIA2 and CIL.

RESULTS

Subcellular localizations of truncated CIA2 protein fragments in our onion transient assay demonstrate that CIA2 contains two nuclear localization signals (NLS) located at amino acids (aa) 62-65 and 291-308, whereas CIL has only one NLS at aa 47-50. We screened a yeast two-hybrid (Y2H) Arabidopsis cDNA library to search for putative CIA2-interacting proteins and identified 12 nuclear proteins, including itself, CIL, and flowering-control proteins (such as CO, NF-YB1, NF-YC1, NF-YC9 and ABI3). Additional Y2H experiments demonstrate that CIA2 and CIL mainly interact with flowering-control proteins via their N-termini, but preferentially form homo- or hetero-dimers through their C-termini. Moreover, sequence alignment showed that the N-terminal sequences of CIA2, CIL and NF-YA are highly conserved. Therefore, NF-YA in the NF-Y complex could be substituted by CIA2 or CIL.

CONCLUSIONS

We show that Arabidopsis CIA2 and CIL can interact with CO and NF-Y complex, so not only may they contribute to regulate chloroplast function but also to modulate flower development.

Rapid functional divergence after small-scale gene duplication in grasses

BACKGROUND

Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events.

RESULTS

We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon.

CONCLUSIONS

Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals.

Transcriptome-wide analysis of MADS-box family genes involved in aluminum and fluoride assimilation in Camellia sinensis

MADS-box transcription factors (TFs) are involved in a variety of processes in flowering plants ranging from root growth to flower and fruit development. However, studies of the tolerance-related functions of MADS-box genes are very limited, and to date no such studies have been conducted on Camellia sinensis. To gain insight into the functions of genes of this family and to elucidate the role they may play in tissue development and Al and F response, we identified 45 MADS-box genes through transcriptomic analysis of C. sinensis. Phylogenetic analysis of these CsMADS-box genes, along with their homologues in Arabidopsis thaliana, enabled us to classify them into distinct groups, including: M-type (Mα), MIKC* and MIKC^c (which contains the SOC1, AGL12, AGL32, SEP, ANR1, SVP, and FLC subgroups). Conserved motif analysis of the CsMADS-box proteins revealed diverse motif compositions indicating a complex evolutionary relationship. Finally, we examined the expression patterns of CsMADS-box genes in various tissues and under different Al and F concentration treatments. Our qPCR results showed that these CsMADS-box genes were involved in Al and F accumulation and root growth in C. sinensis. These findings lay the foundation for future research on the function of CsMADS-box genes and their role in response to Al and F accumulation in root tissues.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity

BACKGROUND

Floral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system.

RESULTS

The transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae.

CONCLUSIONS

We found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

MADS-box genes and crop domestication: the jack of all traits

MADS-box genes are key regulators of virtually every aspect of plant reproductive development. They play especially prominent roles in flowering time control, inflorescence architecture, floral organ identity determination, and seed development. The developmental and evolutionary importance of MADS-box genes is widely acknowledged. However, their role during flowering plant domestication is less well recognized. Here, we provide an overview illustrating that MADS-box genes have been important targets of selection during crop domestication and improvement. Numerous examples from a diversity of crop plants show that various developmental processes have been shaped by allelic variations in MADS-box genes. We propose that new genomic and genome editing resources provide an excellent starting point for further harnessing the potential of MADS-box genes to improve a variety of reproductive traits in crops. We also suggest that the biophysics of MADS-domain protein-protein and protein-DNA interactions, which is becoming increasingly well characterized, makes them especially suited to exploit coding sequence variations for targeted breeding approaches.

Geometric morphometrics reveals shifts in flower shape symmetry and size following gene knockdown of CYCLOIDEA and ANTHOCYANIDIN SYNTHASE

BACKGROUND

While floral symmetry has traditionally been assessed qualitatively, recent advances in geometric morphometrics have opened up new avenues to specifically quantify flower shape and size using robust multivariate statistical methods. In this study, we examine, for the first time, the ability of geometric morphometrics to detect morphological differences in floral dorsoventral asymmetry following virus-induced gene silencing (VIGS). Using Fedia graciliflora Fisch. & Meyer (Valerianaceae) as a model, corolla shape of untreated flowers was compared using canonical variate analysis to knockdown phenotypes of CYCLOIDEA2A (FgCYC2A), ANTHOCYANIDIN SYNTHASE (FgANS), and empty vector controls.

RESULTS

Untreated flowers and all VIGS treatments were morphologically distinct from each other, suggesting that VIGS may cause subtle shifts in floral shape. Knockdowns of FgCYC2A were the most dramatic, affecting the position of dorsal petals in relation to lateral petals, thereby resulting in more actinomorphic-like flowers. Additionally, FgANS knockdowns developed larger flowers with wider corolla tube openings.

CONCLUSIONS

These results provide a method to quantify the role that specific genes play in the developmental pathway affecting the dorsoventral axis of symmetry in zygomorphic flowers. Additionally, they suggest that ANS may have an unintended effect on floral size and shape.

Duplication and Whorl-Specific Down-Regulation of the Obligate AP3-PI Heterodimer Genes Explain the Origin of Paeonia lactiflora Plants with Spontaneous Corolla Mutation

Herbaceous peony (Paeonia lactiflora) is a globally important ornamental plant. Spontaneous floral mutations occur frequently during cultivation, and are selected as a way to release new cultivars, but the underlying evolutionary developmental genetics remain largely elusive. Here, we investigated a collection of spontaneous corolla mutational plants (SCMPs) whose other floral organs were virtually unaffected. Unlike the corolla in normal plants (NPs) that withered soon after fertilization, the transformed corolla (petals) in SCMPs was greenish and persistent similar to the calyx (sepals). Epidermal cellular morphology of the SCMP corolla was also similar to that of calyx cells, further suggesting a sepaloid corolla in SCMPs. Ten floral MADS-box genes from these Paeonia plants were comparatively characterized with respect to sequence and expression. Codogenic sequence variation of these MADS-box genes was not linked to corolla changes in SCMPs. However, we found that both APETALA3 (AP3) and PISTILLATA (PI) lineages of B-class MADS-box genes were duplicated, and subsequent selective expression alterations of these genes were closely associated with the origin of SCMPs. AP3-PI obligate heterodimerization, essential for organ identity of corolla and stamens, was robustly detected. However, selective down-regulation of these duplicated genes might result in a reduction of this obligate heterodimer concentration in a corolla-specific manner, leading to the sepaloid corolla in SCMPs, thus representing a new sepaloid corolla model taking advantage of gene duplication. Our work suggests that modifying floral MADS-box genes could facilitate the breeding of novel cultivars with distinct floral morphology in ornamental plants, and also provides new insights into the functional evolution of the MADS-box genes in plants.

Expression of B-class MADS-box genes in response to variations in photoperiod is associated with chasmogamous and cleistogamous flower development in Viola philippica

BACKGROUND

Some plants develop a breeding system that produces both chasmogamous (CH) and cleistogamous (CL) flowers. However, the underlying molecular mechanism remains elusive.

RESULTS

In the present study, we observed that Viola philippica develops CH flowers with short daylight, whereas an extended photoperiod induces the formation of intermediate CL and CL flowers. In response to long daylight, the respective number and size of petals and stamens was lower and smaller than those of normally developed CH flowers, and a minimum of 14-h light induced complete CL flowers that had no petals but developed two stamens of reduced fertility. The floral ABC model indicates that B-class MADS-box genes largely influence the development of the affected two-whorl floral organs; therefore, we focused on characterizing these genes in V. philippica to understand this particular developmental transition. Three such genes were isolated and respectively designated as VpTM6-1, VpTM6-2, and VpPI. These were differentially expressed during floral development (particularly in petals and stamens) and the highest level of expression was observed in CH flowers; significantly low levels were detected in intermediate CL flowers, and the lowest level in CL flowers. The observed variations in the levels of expression after floral induction and organogenesis apparently occurred in response to variations in photoperiod.

CONCLUSIONS

Therefore, inhibition of the development of petals and stamens might be due to the downregulation of B-class MADS-box gene expression by long daylight, thereby inducing the generation of CL flowers. Our work contributes to the understanding of the adaptive evolutionary formation of dimorphic flowers in plants.

The tomato floral homeotic protein FBP1-like gene, SlGLO1, plays key roles in petal and stamen development

MADS-box transcription factors play important role in plant growth and development, especially floral organ identities. In our study, a MADS-box gene SlGLO1- tomato floral homeotic protein FBP1-like gene was isolated. Its tissue-specific expression profile analysis showed that SlGLO1 was highly expressed in petals and stamens. RNAi (RNA interference) repression of SlGLO1 resulted in floral organ abnormal phenotypes, including green petals with shorter size, and aberrant carpelloid stamens. SlGLO1-silenced lines are male sterile. Total chlorophyll content was increased and chlorophyll biosynthetic genes were significantly up-regulated in SlGLO1-silenced petals and stamens. Furthermore, B-class genes expression analysis indicated that the repressed function of SlGLO1 led to the enhanced expression of TAP3 and the down-regulation of TPI in the petals and stamens, while the expression of TM6 was reduced in petals and increased in stamens and carpels of SlGLO1-RNAi plants. Additionally, pollen grains of transgenic lines were aberrant and failed to germinate and tomato pollen-specific genes were down-regulated by more than 90% in SlGLO1-silenced lines. These results suggest that SlGLO1 plays important role in regulating plant floral organ and pollen development in tomato.

Evolutionary developmental genetics of fruit morphological variation within the Solanaceae

Morphological variations of fruits such as shape and size, and color are a result of adaptive evolution. The evolution of morphological novelties is particularly intriguing. An understanding of these evolutionary processes calls for the elucidation of the developmental and genetic mechanisms that result in particular fruit morphological characteristics, which determine seed dispersal. The genetic and developmental basis for fruit morphological variation was established at a microevolutionary time scale. Here, we summarize the progress on the evolutionary developmental genetics of fruit size, shape and color in the Solanaceae. Studies suggest that the recruitment of a pre-existing gene and subsequent modification of its interaction and regulatory networks are frequently involved in the evolution of morphological diversity. The basic mechanisms underlying changes in plant morphology are alterations in gene expression and/or gene function. We also deliberate on the future direction in evolutionary developmental genetics of fruit morphological variation such as fruit type. These studies will provide insights into plant developmental processes and will help to improve the productivity and fruit quality of crops.