A simplified explanation for the frameshift mutation that created a novel C-terminal motif in the APETALA3 gene lineage

Structural and Functional Analysis of the MADS-Box Genes Reveals Their Functions in Cold Stress Responses and Flower Development in Tea Plant (Camellia sinensis)

MADS-box genes comprise a large family of transcription factors that play crucial roles in all aspects of plant growth and development. However, no detailed information on the evolutionary relationship and functional characterization of MADS-box genes is currently available for some representative lineages, such as the Camellia plant. In this study, 136 MADS-box genes were detected from a reference genome of the tea plant (Camellia sinensis) by employing a 569 bp HMM (Hidden Markov Model) developed using nucleotide sequencing including 73 type I and 63 type II genes. An additional twenty-seven genes were identified, with five MIKC-type genes. Truncated and/or inaccurate gene models were manually verified and curated to improve their functional characterization. Subsequently, phylogenetic relationships, chromosome locations, conserved motifs, gene structures, and gene expression profiles were systematically investigated. Tea plant MIKC genes were divided into all 14 major eudicot subfamilies, and no gene was found in Mβ. The expansion of MADS-box genes in the tea plant was mainly contributed by WGD/fragment and tandem duplications. The expression profiles of tea plant MADS-box genes in different tissues and seasons were analyzed, revealing widespread evolutionary conservation and genetic redundancy. The expression profiles linked to cold stress treatments suggested the wide involvement of MADS-box genes from the tea plant in response to low temperatures. Moreover, a floral 'ABCE' model was proposed in the tea plant and proved to be both conserved and ancient. Our analyses offer a detailed overview of MADS-box genes in the tea plant, allowing us to hypothesize the potential functions of unknown genes and providing a foundation for further functional characterizations.

AT-HOOK MOTIF NUCLEAR LOCALIZED (AHL) proteins of ancient origin radiate new functions

AHL (AT-HOOK MOTIF NUCLEAR LOCALIZED) protein is an important transcription factor in plants that regulates a wide range of biological process. It is considered to have evolved from an independent PPC domain in prokaryotes to a complete protein in modern plants. AT-hook motif and PPC conserved domains are the main functional domains of AHL. Since the discovery of AHL, their evolution and function have been continuously studied. The AHL gene family has been identified in multiple species and the functions of several members of the gene family have been studied. Here, we summarize the evolution and structural characteristics of AHL genes, and emphasize their biological functions. This review will provide a basis for further functional study and crop breeding.

Whole-transcriptome analysis of differentially expressed genes between ray and disc florets and identification of flowering regulatory genes in Chrysanthemum morifolium

Chrysanthemum morifolium has ornamental and economic values. However, there has been minimal research on the morphology of the chrysanthemum florets and related genes. In this study, we used the leaves as a control to screen for differentially expressed genes between ray and disc florets in chrysanthemum flowers. A total of 8,359 genes were differentially expressed between the ray and disc florets, of which 3,005 were upregulated and 5,354 were downregulated in the disc florets. Important regulatory genes that control flower development and flowering determination were identified. Among them, we identified a TM6 gene (CmTM6-mu) that belongs to the Class B floral homeotic MADS-box transcription factor family, which was specifically expressed in disc florets. We isolated this gene and found it was highly similar to other typical TM6 lineage genes, but a single-base deletion at the 3' end of the open reading frame caused a frame shift that generated a protein in which the TM6-specific paleoAP3 motif was missing at the C terminus. The CmTM6-mu gene was ectopically expressed in Arabidopsis thaliana. Petal and stamen developmental processes were unaffected in transgenic A. thaliana lines; however, the flowering time was earlier than in the wild-type control. Thus, the C-terminal of paleoAP3 appears to be necessary for the functional performance in regulating the development of petals or stamens and CmTM6-mu may be involved in the regulation of flowering time in chrysanthemum. The results of this study will be useful for future research on flowering molecular mechanisms and for the breeding of novel flower types.

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.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

Insights Into the Molecular Evolution of AT-Hook Motif Nuclear Localization Genes in Brassica napus

AT-hook motif nuclear localization (AHL) proteins belong to a family of transcription factors, and play important roles in plant growth and development and response to various stresses through protein-DNA and protein-protein interactions. To better understand the Brassica napus AHL gene family, AHL genes in B. napus and related species were analyzed. Using Arabidopsis as a reference, 122 AHL gene family members were first identified in B. napus. According to the phylogenetic tree and gene organization, the BnaAHLs were classified into two clades (Clade-A and Clade-B) and three types (Type-I, Type-II, and Type-III). Gene organization and motif distribution analysis suggested that the AHL gene family is relatively conserved during evolution. These BnaAHLs are unevenly distributed on 38 chromosomes and expanded by whole-genome duplication (WGD) or segmental duplication. And large-scale loss events have also occurred in evolution. All types of BnaAHLs are subject to purification or neutral selection, while some positive selection sites are also identified in Type-II and Type-III groups. At the same time, the purification effect of Type-I members are stronger than that of the others. In addition, RNA-seq data and cis-acting element analysis also suggested that the BnaAHLs play important roles in B. napus growth and development, as well as in response to some abiotic and biotic stresses. Protein-protein interaction analysis identified some important BnaAHL-binding proteins, which also play key roles in plant growth and development. This study is helpful to fully understand the origin and evolution of the AHL gene in B. napus, and lays the foundation for their functional studies.

Resurrected Protein Interaction Networks Reveal the Innovation Potential of Ancient Whole-Genome Duplication

The evolution of plants is characterized by whole-genome duplications, sometimes closely associated with the origin of large groups of species. The gamma (γ) genome triplication occurred at the origin of the core eudicots, which comprise ∼75% of flowering plants. To better understand the impact of whole-genome duplication, we studied the protein interaction network of MADS domain transcription factors, which are key regulators of reproductive development. We reconstructed, synthesized, and tested the interactions of ancestral proteins immediately before and closely after the triplication and directly compared these ancestral networks to the extant networks of Arabidopsis thaliana and tomato (Solanum lycopersicum). We found that gamma expanded the MADS domain interaction network more strongly than subsequent genomic events. This event strongly rewired MADS domain interactions and allowed for the evolution of new functions and installed robustness through new redundancy. Despite extensive rewiring, the organization of the network was maintained through gamma. New interactions and protein retention compensated for its potentially destructive impact on network organization. Post gamma, the network evolved from an organization around the single hub SEP3 to a network organized around multiple hubs and well-connected proteins lost, rather than gained, interactions. The data provide a resource for comparative developmental biology in flowering plants.

Evolution by duplication: paleopolyploidy events in plants reconstructed by deciphering the evolutionary history of VOZ transcription factors

BACKGROUND

Facilitated by the rapid progress of sequencing technology, comparative genomic studies in plants have unveiled recurrent whole genome duplication (i.e. polyploidization) events throughout plant evolution. The evolutionary past of plant genes should be analyzed in a background of recurrent polyploidy events in distinctive plant lineages. The Vascular Plant One Zinc-finger (VOZ) gene family encode transcription factors associated with a number of important traits including control of flowering time and photoperiodic pathways, but the evolutionary trajectory of this gene family remains uncharacterized.

RESULTS

In this study, we deciphered the evolutionary history of the VOZ gene family by analyses of 107 VOZ genes in 46 plant genomes using integrated methods: phylogenic reconstruction, Ks-based age estimation and genomic synteny comparisons. By scrutinizing the VOZ gene family phylogeny the core eudicot γ event was well circumscribed, and relics of the precommelinid τ duplication event were detected by incorporating genes from oil palm and banana. The more recent T and ρ polyploidy events, closely coincident with the species diversification in Solanaceae and Poaceae, respectively, were also identified. Other important polyploidy events captured included the "salicoid" event in poplar and willow, the "early legume" and "soybean specific" events in soybean, as well as the recent polyploidy event in Physcomitrella patens. Although a small transcription factor gene family, the evolutionary history of VOZ genes provided an outstanding record of polyploidy events in plants. The evolutionary past of VOZ gene family demonstrated a close correlation with critical plant polyploidy events which generated species diversification and provided answer to Darwin's "abominable mystery".

CONCLUSIONS

We deciphered the evolutionary history of VOZ transcription factor family in plants and ancestral polyploidy events in plants were recapitulated simultaneously. This analysis allowed for the generation of an idealized plant gene tree demonstrating distinctive retention and fractionation patterns following polyploidy events.

Genome-wide survey of potato MADS-box genes reveals that StMADS1 and StMADS13 are putative downstream targets of tuberigen StSP6A

Background

MADS-box genes encode transcription factors that are known to be involved in several aspects of plant growth and development, especially in floral organ specification. To date, the comprehensive analysis of potato MADS-box gene family is still lacking after the completion of potato genome sequencing. A genome-wide characterization, classification, and expression analysis of MADS-box transcription factor gene family was performed in this study.

Results

A total of 153 MADS-box genes were identified and categorized into MIKC subfamily (MIKC^C and MIKC^*) and M-type subfamily (Mα, Mβ, and Mγ) based on their phylogenetic relationships to the Arabidopsis and rice MADS-box genes. The potato M-type subfamily had 114 members, which is almost three times of the MIKC members (39), indicating that M-type MADS-box genes have a higher duplication rate and/or a lower loss rate during potato genome evolution. Potato MADS-box genes were present on all 12 potato chromosomes with substantial clustering that mainly contributed by the M-type members. Chromosomal localization of potato MADS-box genes revealed that MADS-box genes, mostly MIKC, were located on the duplicated segments of the potato genome whereas tandem duplications mainly contributed to the M-type gene expansion. The potato MIKC subfamily could be further classified into 11 subgroups and the TT16-like, AGL17-like, and FLC-like subgroups found in Arabidopsis were absent in potato. Moreover, the expressions of potato MADS-box genes in various tissues were analyzed by using RNA-seq data and verified by quantitative real-time PCR, revealing that the MIKC^C genes were mainly expressed in flower organs and several of them were highly expressed in stolon and tubers. StMADS1 and StMADS13 were up-regulated in the StSP6A-overexpression plants and down-regulated in the StSP6A-RNAi plant, and their expression in leaves and/or young tubers were associated with high level expression of StSP6A.

Conclusion

Our study identifies the family members of potato MADS-box genes and investigate the evolution history and functional divergence of MADS-box gene family. Moreover, we analyze the MIKC^C expression patterns and screen for genes involved in tuberization. Finally, the StMADS1 and StMADS13 are most likely to be downstream target of StSP6A and involved in tuber development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5113-z) contains supplementary material, which is available to authorized users.

Identification and Characterization of the MADS-Box Genes and Their Contribution to Flower Organ in Carnation (Dianthus caryophyllus L.)

Dianthus is a large genus containing many species with high ornamental economic value. Extensive breeding strategies permitted an exploration of an improvement in the quality of cultivated carnation, particularly in flowers. However, little is known on the molecular mechanisms of flower development in carnation. Here, we report the identification and description of MADS-box genes in carnation (DcaMADS) with a focus on those involved in flower development and organ identity determination. In this study, 39 MADS-box genes were identified from the carnation genome and transcriptome by the phylogenetic analysis. These genes were categorized into four subgroups (30 MIKC^c, two MIKC*, two Mα, and five Mγ). The MADS-box domain, gene structure, and conserved motif compositions of the carnation MADS genes were analysed. Meanwhile, the expression of DcaMADS genes were significantly different in stems, leaves, and flower buds. Further studies were carried out for exploring the expression of DcaMADS genes in individual flower organs, and some crucial DcaMADS genes correlated with their putative function were validated. Finally, a new expression pattern of DcaMADS genes in flower organs of carnation was provided: sepal (three class E genes and two class A genes), petal (two class B genes, two class E genes, and one SHORT VEGETATIVE PHASE (SVP)), stamen (two class B genes, two class E genes, and two class C), styles (two class E genes and two class C), and ovary (two class E genes, two class C, one AGAMOUS-LIKE 6 (AGL6), one SEEDSTICK (STK), one B sister, one SVP, and one Mα). This result proposes a model in floral organ identity of carnation and it may be helpful to further explore the molecular mechanism of flower organ identity in carnation.

Evolution of Gene Duplication in Plants

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

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.

Prevalent Exon-Intron Structural Changes in the APETALA1/FRUITFULL, SEPALLATA, AGAMOUS-LIKE6, and FLOWERING LOCUS C MADS-Box Gene Subfamilies Provide New Insights into Their Evolution

AP1/FUL, SEP, AGL6, and FLC subfamily genes play important roles in flower development. The phylogenetic relationships among them, however, have been controversial, which impedes our understanding of the origin and functional divergence of these genes. One possible reason for the controversy may be the problems caused by changes in the exon-intron structure of genes, which, according to recent studies, may generate non-homologous sites and hamper the homology-based sequence alignment. In this study, we first performed exon-by-exon alignments of these and three outgroup subfamilies (SOC1, AG, and STK). Phylogenetic trees reconstructed based on these matrices show improved resolution and better congruence with species phylogeny. In the context of these phylogenies, we traced evolutionary changes of exon-intron structures in each subfamily. We found that structural changes have occurred frequently following gene duplication and speciation events. Notably, exons 7 and 8 (if present) suffered more structural changes than others. With the knowledge of exon-intron structural changes, we generated more reasonable alignments containing all the focal subfamilies. The resulting trees showed that the SEP subfamily is sister to the monophyletic group formed by AP1/FUL and FLC subfamily genes and that the AGL6 subfamily forms a sister group to the three abovementioned subfamilies. Based on this topology, we inferred the evolutionary history of exon-intron structural changes among different subfamilies. Particularly, we found that the eighth exon originated before the divergence of AP1/FUL, FLC, SEP, and AGL6 subfamilies and degenerated in the ancestral FLC-like gene. These results provide new insights into the origin and evolution of the AP1/FUL, FLC, SEP, and AGL6 subfamilies.

Ectopic expression of FaesAP3, a Fagopyrum esculentum (Polygonaceae) AP3 orthologous gene rescues stamen development in an Arabidopsis ap3 mutant

Arabidopsis thaliana APETALA3 (AP3) and Antirrhinum majus DEFICIENS (DEF) MADS box genes are required to specify petal and stamen identity. AP3 and DEF are members of the euAP3 lineage, which arose by gene duplication coincident with radiation of the core eudicots. In order to investigate the molecular mechanisms underlying organ development in early diverging clades of core eudicots, we isolated and identified an AP3 homolog, FaesAP3, from Fagopyrum esculentum (buckwheat, Polygonaceae), a multi-food-use pseudocereal with healing benefits. Protein sequence alignment and phylogenetic analyses revealed that FaesAP3 grouped into the euAP3 lineage. Expression analysis showed that FaesAP3 was transcribed only in developing stamens, and differed from AP3 and DEF, which expressed in developing petals and stamens. Moreover, ectopic expression of FaesAP3 rescued stamen development without complementation of petal development in an Arabidopsis ap3 mutant. Our results suggest that FaesAP3 is involved in the development of stamens in buckwheat. These results also suggest that FaesAP3 holds some potential for biotechnical engineering to create a male sterile line of F. esculentum.

Analysis of the APETALA3- and PISTILLATA-like genes in Hedyosmum orientale (Chloranthaceae) provides insight into the evolution of the floral homeotic B-function in angiosperms

BACKGROUND AND AIMS

According to the floral ABC model, B-function genes appear to play a key role in the origin and diversification of the perianth during the evolution of angiosperms. The basal angiosperm Hedyosmum orientale (Chloranthaceae) has unisexual inflorescences associated with a seemingly primitive reproductive morphology and a reduced perianth structure in female flowers. The aim of this study was to investigate the nature of the perianth and the evolutionary state of the B-function programme in this species.

METHODS

A series of experiments were conducted to characterize B-gene homologues isolated from H. orientale, including scanning electron microscopy to observe the development of floral organs, phylogenetic analysis to reconstruct gene evolutionary history, reverse transcription-PCR, quantitative real-time PCR and in situ hybridization to identify gene expression patterns, the yeast two-hybrid assay to explore protein dimerization affinities, and transgenic analyses in Arabidopsis thaliana to determine activities of the encoded proteins.

KEY RESULTS

The expression of HoAP3 genes was restricted to stamens, whereas HoPI genes were broadly expressed in all floral organs. HoAP3 was able to partially restore the stamen but not petal identity in Arabidopsis ap3-3 mutants. In contrast, HoPI could rescue aspects of both stamen and petal development in Arabidopsis pi-1 mutants. When the complete C-terminal sequence of HoPI was deleted, however, no or weak transgenic phenotypes were observed and homodimerization capability was completely abolished.

CONCLUSIONS

The results suggest that Hedyosmum AP3-like genes have an ancestral function in specifying male reproductive organs, and that the activity of the encoded PI-like proteins is highly conserved between Hedyosmum and Arabidopsis. Moreover, there is evidence that the C-terminal region is important for the function of HoPI. Our findings indicate that the development of the proposed perianth in Hedyosmum does not rely on the B homeotic function.

Protein change in plant evolution: tracing one thread connecting molecular and phenotypic diversity

Proteins change over the course of evolutionary time. New protein-coding genes and gene families emerge and diversify, ultimately affecting an organism's phenotype and interactions with its environment. Here we survey the range of structural protein change observed in plants and review the role these changes have had in the evolution of plant form and function. Verified examples tying evolutionary change in protein structure to phenotypic change remain scarce. We will review the existing examples, as well as draw from investigations into domestication, and quantitative trait locus (QTL) cloning studies searching for the molecular underpinnings of natural variation. The evolutionary significance of many cloned QTL has not been assessed, but all the examples identified so far have begun to reveal the extent of protein structural diversity tolerated in natural systems. This molecular (and phenotypic) diversity could come to represent part of natural selection's source material in the adaptive evolution of novel traits. Protein structure and function can change in many distinct ways, but the changes we identified in studies of natural diversity and protein evolution were predicted to fall primarily into one of six categories: altered active and binding sites; altered protein-protein interactions; altered domain content; altered activity as an activator or repressor; altered protein stability; and hypomorphic and hypermorphic alleles. There was also variability in the evolutionary scale at which particular changes were observed. Some changes were detected at both micro- and macroevolutionary timescales, while others were observed primarily at deep or shallow phylogenetic levels. This variation might be used to determine the trajectory of future investigations in structural molecular evolution.

Functional specialization of duplicated AP3-like genes in Medicago truncatula

The B-class of MADS box genes has been studied in a wide range of plant species, but has remained largely uncharacterized in legumes. Here we investigate the evolutionary fate of the duplicated AP3-like genes of a legume species. To obtain insight into the extent to which B-class MADS box gene functions are conserved or have diversified in legumes, we isolated and characterized the two members of the AP3 lineage in Medicago truncatula: MtNMH7 and MtTM6 (euAP3 and paleoAP3 genes, respectively). A non-overlapping and complementary expression pattern of both genes was observed in petals and stamens. MtTM6 was expressed predominantly in the outer cell layers of both floral organs, and MtNMH7 in the inner cell layers of petals and stamens. Functional analyses by reverse genetics approaches (RNAi and Tnt1 mutagenesis) showed that the contribution of MtNMH7 to petal identity is more important than that of MtTM6, whereas MtTM6 plays a more important role in stamen identity than its paralog MtNMH7. Our results suggest that the M. truncatula AP3-like genes have undergone a functional specialization process associated with complete partitioning of gene expression patterns of the ancestral gene lineage. We provide information regarding the similarities and differences in petal and stamen development among core eudicots.

Rate heterogeneity in six protein-coding genes from the holoparasite Balanophora (Balanophoraceae) and other taxa of Santalales

BACKGROUND AND AIMS

The holoparasitic flowering plant Balanophora displays extreme floral reduction and was previously found to have enormous rate acceleration in the nuclear 18S rDNA region. So far, it remains unclear whether non-ribosomal, protein-coding genes of Balanophora also evolve in an accelerated fashion and whether the genes with high substitution rates retain their functionality. To tackle these issues, six different genes were sequenced from two Balanophora species and their rate variation and expression patterns were examined.

METHODS

Sequences including nuclear PI, euAP3, TM6, LFY and RPB2 and mitochondrial matR were determined from two Balanophora spp. and compared with selected hemiparasitic species of Santalales and autotrophic core eudicots. Gene expression was detected for the six protein-coding genes and the expression patterns of the three B-class genes (PI, AP3 and TM6) were further examined across different organs of B. laxiflora using RT-PCR.

KEY RESULTS

Balanophora mitochondrial matR is highly accelerated in both nonsynonymous (d(N)) and synonymous (d(S)) substitution rates, whereas the rate variation of nuclear genes LFY, PI, euAP3, TM6 and RPB2 are less dramatic. Significant d(S) increases were detected in Balanophora PI, TM6, RPB2 and d(N) accelerations in euAP3. All of the protein-coding genes are expressed in inflorescences, indicative of their functionality. PI is restrictively expressed in tepals, synandria and floral bracts, whereas AP3 and TM6 are widely expressed in both male and female inflorescences.

CONCLUSIONS

Despite the observation that rates of sequence evolution are generally higher in Balanophora than in hemiparasitic species of Santalales and autotrophic core eudicots, the five nuclear protein-coding genes are functional and are evolving at a much slower rate than 18S rDNA. The mechanism or mechanisms responsible for rapid sequence evolution and concomitant rate acceleration for 18S rDNA and matR are currently not well understood and require further study in Balanophora and other holoparasites.

Divergences of MPF2-like MADS-domain proteins have an association with the evolution of the inflated calyx syndrome within Solanaceae

The inflated calyx syndrome (ICS) is a post-floral novelty within Solanaceae. Previous work has shown that MPF2-like MADS-box genes have been recruited for the development and evolution of ICS through heterotopic expression from vegetative to floral organs. ICS seems to be a plesiomorphic trait in Physaleae, but it has been secondarily lost in some lineages during evolution. We hypothesized that molecular and functional divergences of MPF2-like proteins might play a role in the loss of ICS. In this study we analyzed the phylogeny, selection and various functions of MPF2-like proteins with respect to the evolution of ICS. Directional selection of MPF2-like orthologs toward evolution of ICS was detected. While auto-activation capacity between proteins varies in yeast, MPF2-like interaction with floral MADS-domain proteins is robustly detected, hence substantiating their integration into the floral developmental programs. Dimerization with A- (MPF3) and E-function (PFSEP1/3) proteins seems to be essential for ICS development within Solanaceae. Moreover, the occurrence of the enlarged sepals, reminiscent of ICS, and MPF2-like interactions with these specific partners were observed in transgenic Arabidopsis. The interaction spectrum relevant to ICS seems to be plesiomorphic, reinforcing the plesiomorphy of this trait. The inability of some MPF2-like to interact with either the A-function or any of the E-function partners characterized is correlated with the loss of ICS in the lineages that showed a MPF2-like expression in the calyx. Our findings suggest that, after recruitment of MPF2-like genes for floral development, diversification in their coding region due to directional selection leads to a modification of the MADS-domain protein interacting spectrum, which might serve as a constraint for the evolution of ICS within Solanaceae.

The MADS and the Beauty: Genes Involved in the Development of Orchid Flowers

Since the time of Darwin, biologists have studied the origin and evolution of the Orchidaceae, one of the largest families of flowering plants. In the last two decades, the extreme diversity and specialization of floral morphology and the uncoupled rate of morphological and molecular evolution that have been observed in some orchid species have spurred interest in the study of the genes involved in flower development in this plant family. As part of the complex network of regulatory genes driving the formation of flower organs, the MADS-box represents the most studied gene family, both from functional and evolutionary perspectives. Despite the absence of a published genome for orchids, comparative genetic analyses are clarifying the functional role and the evolutionary pattern of the MADS-box genes in orchids. Various evolutionary forces act on the MADS-box genes in orchids, such as diffuse purifying selection and the relaxation of selective constraints, which sometimes reveals a heterogeneous selective pattern of the coding and non-coding regions. The emerging theory regarding the evolution of floral diversity in orchids proposes that the diversification of the orchid perianth was a consequence of duplication events and changes in the regulatory regions of the MADS-box genes, followed by sub- and neo-functionalization. This specific developmental-genetic code is termed the "orchid code."

Gamma paleohexaploidy in the stem lineage of core eudicots: significance for MADS-box gene and species diversification

Comparative genome biology has unveiled the polyploid origin of all angiosperms and the role of recurrent polyploidization in the amplification of gene families and the structuring of genomes. Which species share certain ancient polyploidy events, and which do not, is ill defined because of the limited number of sequenced genomes and transcriptomes and their uneven phylogenetic distribution. Previously, it has been suggested that most, but probably not all, of the eudicots have shared an ancient hexaploidy event, referred to as the gamma triplication. In this study, detailed phylogenies of subfamilies of MADS-box genes suggest that the gamma triplication has occurred before the divergence of Gunnerales but after the divergence of Buxales and Trochodendrales. Large-scale phylogenetic and K(S)-based approaches on the inflorescence transcriptomes of Gunnera manicata (Gunnerales) and Pachysandra terminalis (Buxales) provide further support for this placement, enabling us to position the gamma triplication in the stem lineage of the core eudicots. This triplication likely initiated the functional diversification of key regulators of reproductive development in the core eudicots, comprising 75% of flowering plants. Although it is possible that the gamma event triggered early core eudicot diversification, our dating estimates suggest that the event occurred early in the stem lineage, well before the rapid speciation of the earliest core eudicot lineages. The evolutionary significance of this paleopolyploidy event may thus rather lie in establishing a species lineage that was resilient to extinction, but with the genomic potential for later diversification. We consider that the traits generated from this potential characterize extant core eudicots both chemically and morphologically.

Expression of floral MADS-box genes in Sinofranchetia chinensis (Lardizabalaceae): implications for the nature of the nectar leaves

BACKGROUND AND AIMS

The perianths of the Lardizabalaceae are diverse. The second-whorl floral organs of Sinofranchetia chinensis (Lardizabalaceae) are nectar leaves. The aim of this study was to explore the nature of this type of floral organ, and to determine its relationship to nectar leaves in other Ranunculales species, and to other floral organs in Sinofranchetia chinensis.

METHODS

Approaches of evolutionary developmental biology were used, including 3' RACE (rapid amplification of cDNA ends) for isolating floral MADS-box genes, phylogenetic analysis for reconstructing gene evolutionary history, in situ hybridization and tissue-specific RT-PCR for identifying gene expression patterns and SEM (scanning electron microscopy) for observing the epidermal cell morphology of floral organs.

KEY RESULTS

Fourteen new floral MADS-box genes were isolated from Sinofranchetia chinensis and from two other species of Lardizabalaceae, Holboellia grandiflora and Decaisnea insignis. The phylogenetic analysis of AP3-like genes in Ranunculales showed that three AP3 paralogues from Sinofranchetia chinensis belong to the AP3-I, -II and -III lineages. In situ hybridization results showed that SIchAP3-3 is significantly expressed only in nectar leaves at the late stages of floral development, and SIchAG, a C-class MADS-box gene, is expressed not only in stamens and carpels, but also in nectar leaves. SEM observation revealed that the adaxial surface of nectar leaves is covered with conical epidermal cells, a hallmark of petaloidy.

CONCLUSIONS

The gene expression data imply that the nectar leaves in S. chinensis might share a similar genetic regulatory code with other nectar leaves in Ranunculales species. Based on gene expression and morphological evidence, it is considered that the nectar leaves in S. chinensis could be referred to as petals. Furthermore, the study supports the hypothesis that the nectar leaves in some Ranunculales species might be derived from stamens.

Phylogenetic and evolutionary analysis of A-, B-, C- and E-class MADS-box genes in the basal eudicot Platanus acerifolia

London plane tree (Platanus acerifolia Willd.) is an important member of the Platanaceae family, being popular as an urban landscaping tree. Here, we report the isolation of five MADS-box genes from the basal angiosperm, Platanus acerifolia. Sequence and phylogenetic analyses identified FRUITFUL-like, APETELA3-like, AGAMOUS-like, SEPALLATA1-like and SEPALLATA3-like sequences and, hence, we term the respective Platanus acerifolia genes as PlacFUL, PlacAP3, PlacAG, PlacSEP1 and PlacSEP3. From these identities we infer that they represent candidate A-, B-, C-class and two E-class genes, respectively. The conserved MIK or MIKC domains from the nucleotide and protein sequences of PlacFUL, PlacAP3, PlacAG, PlacSEP1 and PlacSEP3 were analyzed using the maximum-likelihood, MrBayes and neighbor-joining methods. The results confirmed P. acerifolia as a basal eudicot. Expression pattern was determined by reverse transcriptase PCR, which showed all paralogous genes have distinct expression patterns, suggesting that they had undergone functional divergence.

An approach for searching insertions in bacterial genes leading to the phase shift of triplet periodicity

The concept of the phase shift of triplet periodicity (TP) was used for searching potential DNA insertions in genes from 17 bacterial genomes. A mathematical algorithm for detection of these insertions has been developed. This approach can detect potential insertions and deletions with lengths that are not multiples of three bases, especially insertions of relatively large DNA fragments (>100 bases). New similarity measure between triplet matrixes was employed to improve the sensitivity for detecting the TP phase shift. Sequences of 17,220 bacterial genes with each consisting of more than 1,200 bases were analyzed, and the presence of a TP phase shift has been shown in ~16% of analysed genes (2,809 genes), which is about 4 times more than that detected in our previous work. We propose that shifts of the TP phase may indicate the shifts of reading frame in genes after insertions of the DNA fragments with lengths that are not multiples of three bases. A relationship between the phase shifts of TP and the frame shifts in genes is discussed.

Analysis of the formation of flower shapes in wild species and cultivars of tree peony using the MADS-box subfamily gene

Tree peony (Paeonia suffricotisa) cultivars have a unique character compared with wild species; the stamen petalody results in increased whorls of petals and generates different flower forms, which are one of the most important traits for cultivar classification. In order to investigate how petaloid stamens are formed, we obtained the coding sequence (666 bp) and genomic DNA sequence of the PsTM6 genes (belongs to B subfamily of MADS-box gene family) from 23 tree peony samples, Five introns and six exons consisted of the genomic DNA sequence. The analysis of cis-acting regulatory elements in the third and fourth intron indicated that they were highly conserved in all samples. Partial putative amino acids were analyzed and the results suggested that functional differentiation of PsTM6 paralogs apparently affected stamen petalody and flower shape formation due to due to amino acid substitution caused by differences in polarity and electronic charge. Sliding window analysis indicated that the different regions of PsTM6 were subjected to different selection forces, especially in the K domain. This is the first attempt to investigate genetic control of the stamen petalody based on the PsTM6 sequence. This will provide a basis for understanding the evolution of PsTM6 and its the function of in determining stamen morphology of tree peony.

Molecular aspects of flower development in grasses

The grass family (Poaceae) of the monocotyledons includes about 10,000 species and represents one of the most important taxa among angiosperms. Their flower morphology is remarkably different from those of other monocotyledons and higher eudicots. The peculiar floral structure of grasses is the floret, which contains carpels and stamens, like eudicots, but lacks petals and sepals. The reproductive organs are surrounded by two lodicules, which correspond to eudicot petals, and by a palea and lemma, whose correspondence to eudicot organs remains controversial. The molecular and genetic analysis of floral morphogenesis and organ specification, primarily performed in eudicot model species, led to the ABCDE model of flower development. Several genes required for floral development in grasses correspond to class A, B, C, D, and E genes of eudicots, but others appear to have unique and diversified functions. In this paper, we outline the present knowledge on the evolution and diversification of grass genes encoding MIKC-type MADS-box transcription factors, based on information derived from studies in rice, maize, and wheat. Moreover, we review recent advances in studying the genes involved in the control of flower development and the extent of structural and functional conservation of these genes between grasses and eudicots.

Functional analysis of the two Brassica AP3 genes involved in apetalous and stamen carpelloid phenotypes

The Arabidopsis homeotic genes APETALA3 (AP3) and PISTILLATA (PI) are B genes which encode MADS-box transcription factors and specify petal and stamen identities. In the current study, the stamen carpelloid (SC) mutants, HGMS and AMS, of B. rapa and B. napus were investigated and two types of AP3 genes, B.AP3.a and B.AP3.b, were functional characterized. B.AP3.a and B.AP3.b share high similarity in amino acid sequences except for 8 residues difference located at the C-terminus. Loss of this 8 residues in B.AP3.b led to the change of PI-derived motifs. Meanwhile, B.AP3.a specified petal and stamen development, whereas B.AP3.b only specified stamen development. In B. rapa, the mutations of both genes generated the SC mutant HGMS. In B. napus that contained two B.AP3.a and two B.AP3.b, loss of the two B.AP3.a functions was the key reason for the apetalous mutation, however, the loss-of-function in all four AP3 was related to the SC mutant AMS. We inferred that the 8 residues or the PI-derived motif in AP3 gene probably relates to petal formation.

Robustness and evolvability in the B-system of flower development

BACKGROUND

Gene duplication has often been invoked as a key mechanism responsible for evolution of new morphologies. The floral homeotic B-group gene family has undergone a number of gene duplication events, and yet the functions of these genes appear to be largely conserved. However, detailed comparative analysis has indicated that such duplicate genes have considerable cryptic variability in their functions. In the Solanaceae, two duplicate B-group gene lineages have been retained in three subfamilies. Comparisons of orthologous genes across members of the Solanaceae have demonstrated that the combined function of all four B-gene members is to establish petal and stamen identity, but that this function was partitioned differently in each species. These observations emphasize both the robustness and the evolvability of the B-system.

SCOPE

We provide an overview of how the B-function genes can robustly specify petal and stamen identity and at the same time evolve through changes in protein-protein interaction, gene expression patterns, copy number variation or alterations in the downstream genes they control. By using mathematical models we explore regulatory differences between species and how these impose constraints on downstream gene regulation.

CONCLUSIONS

Evolvability of the B-genes can be understood through the multiple ways in which the B-system can be robust. Quantitative approaches should allow for the incorporation of more biological realism in the representations of these regulatory systems and this should contribute to understanding the constraints under which different B-systems can function and evolve. This, in turn, can provide a better understanding of the ways in which B-genes have contributed to flower diversity.

Molecular, genetic and transcriptional evidence for a role of VvAGL11 in stenospermocarpic seedlessness in grapevine

BACKGROUND

Stenospermocarpy is a mechanism through which certain genotypes of Vitis vinifera L. such as Sultanina produce berries with seeds reduced in size. Stenospermocarpy has not yet been characterized at the molecular level.

RESULTS

Genetic and physical maps were integrated with the public genomic sequence of Vitis vinifera L. to improve QTL analysis for seedlessness and berry size in experimental progeny derived from a cross of two seedless genotypes. Major QTLs co-positioning for both traits on chromosome 18 defined a 92-kb confidence interval. Functional information from model species including Vitis suggested that VvAGL11, included in this confidence interval, might be the main positional candidate gene responsible for seed and berry development.Characterization of VvAGL11 at the sequence level in the experimental progeny identified several SNPs and INDELs in both regulatory and coding regions. In association analyses performed over three seasons, these SNPs and INDELs explained up to 78% and 44% of the phenotypic variation in seed and berry weight, respectively. Moreover, genetic experiments indicated that the regulatory region has a larger effect on the phenotype than the coding region. Transcriptional analysis lent additional support to the putative role of VvAGL11's regulatory region, as its expression is abolished in seedless genotypes at key stages of seed development. These results transform VvAGL11 into a functional candidate gene for further analyses based on genetic transformation.For breeding purposes, intragenic markers were tested individually for marker assisted selection, and the best markers were those closest to the transcription start site.

CONCLUSION

We propose that VvAGL11 is the major functional candidate gene for seedlessness, and we provide experimental evidence suggesting that the seedless phenotype might be caused by variations in its promoter region. Current knowledge of the function of its orthologous genes, its expression profile in Vitis varieties and the strong association between its sequence variation and the degree of seedlessness together indicate that the D-lineage MADS-box gene VvAGL11 corresponds to the Seed Development Inhibitor locus described earlier as a major locus for seedlessness. These results provide new hypotheses for further investigations of the molecular mechanisms involved in seed and berry development.

Evolution of prokaryotic genes by shift of stop codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene's reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene's end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3'UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3'UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.

Genomic tools development for Aquilegia: construction of a BAC-based physical map

BACKGROUND

The genus Aquilegia, consisting of approximately 70 taxa, is a member of the basal eudicot lineage, Ranuculales, which is evolutionarily intermediate between monocots and core eudicots, and represents a relatively unstudied clade in the angiosperm phylogenetic tree that bridges the gap between these two major plant groups. Aquilegia species are closely related and their distribution covers highly diverse habitats. These provide rich resources to better understand the genetic basis of adaptation to different pollinators and habitats that in turn leads to rapid speciation. To gain insights into the genome structure and facilitate gene identification, comparative genomics and whole-genome shotgun sequencing assembly, BAC-based genomics resources are of crucial importance.

RESULTS

BAC-based genomic resources, including two BAC libraries, a physical map with anchored markers and BAC end sequences, were established from A. formosa. The physical map was composed of a total of 50,155 BAC clones in 832 contigs and 3939 singletons, covering 21X genome equivalents. These contigs spanned a physical length of 689.8 Mb (~2.3X of the genome) suggesting the complex heterozygosity of the genome. A set of 197 markers was developed from ESTs induced by drought-stress, or involved in anthocyanin biosynthesis or floral development, and was integrated into the physical map. Among these were 87 genetically mapped markers that anchored 54 contigs, spanning 76.4 Mb (25.5%) across the genome. Analysis of a selection of 12,086 BAC end sequences (BESs) from the minimal tiling path (MTP) allowed a preview of the Aquilegia genome organization, including identification of transposable elements, simple sequence repeats and gene content. Common repetitive elements previously reported in both monocots and core eudicots were identified in Aquilegia suggesting the value of this genome in connecting the two major plant clades. Comparison with sequenced plant genomes indicated a higher similarity to grapevine (Vitis vinifera) than to rice and Arabidopsis in the transcriptomes.

CONCLUSIONS

The A. formosa BAC-based genomic resources provide valuable tools to study Aquilegia genome. Further integration of other existing genomics resources, such as ESTs, into the physical map should enable better understanding of the molecular mechanisms underlying adaptive radiation and elaboration of floral morphology.

Hidden variability of floral homeotic B genes in Solanaceae provides a molecular basis for the evolution of novel functions

B-class MADS box genes specify petal and stamen identities in several core eudicot species. Members of the Solanaceae possess duplicate copies of these genes, allowing for diversification of function. To examine the changing roles of such duplicate orthologs, we assessed the functions of B-class genes in Nicotiana benthamiana and tomato (Solanum lycopersicum) using virus-induced gene silencing and RNA interference approaches. Loss of function of individual duplicates can have distinct phenotypes, yet complete loss of B-class gene function results in extreme homeotic transformations of petal and stamen identities. We also show that these duplicate gene products have qualitatively different protein-protein interaction capabilities and different regulatory roles. Thus, compensatory changes in B-class MADS box gene duplicate function have occurred in the Solanaceae, in that individual gene roles are distinct, but their combined functions are equivalent. Furthermore, we show that species-specific differences in the stamen regulatory network are associated with differences in the expression of the microRNA miR169. Whereas there is considerable plasticity in individual B-class MADS box transcription factor function, there is overall conservation in the roles of the multimeric MADS box B-class protein complexes, providing robustness in the specification of petal and stamen identities. Such hidden variability in gene function as we observe for individual B-class genes can provide a molecular basis for the evolution of regulatory functions that result in novel morphologies.

Expression divergence of the AGL6 MADS domain transcription factor lineage after a core eudicot duplication suggests functional diversification

BACKGROUND

Because of their known role as transcriptional regulators of key plant developmental processes, the diversification of MADS-box gene function is thought to be a major driving force in the developmental evolution of plants. Yet the function of some MADS-box gene subfamilies has remained elusive thus far. One such lineage, AGL6, has now been functionally characterized in three angiosperm species, but a phylogenetic framework for comparison of AGL6 gene function is currently missing.

RESULTS

Based on phylogenetic analyses of newly isolated and EST-based sequences, we describe the duplication history of the AGL6 subfamily in angiosperms. Our analyses provide support for four ancient duplications in the evolution of the AGL6 lineage: one at the base of core eudicots resulting in euAGL6 and AGL6-like gene clades, one during basal angiosperm diversification and two in monocot evolution. To investigate whether the spatial domains in which AGL6 genes function have diverged after duplication, we use quantitative Real Time PCR. We show that the core eudicot AGL6-like clade acquired expression in vegetative tissues, while its paralog euAGL6 remains predominantly confined to reproductive tissues.

CONCLUSIONS

These and previous data lead us to propose that the AGL6 lineage in core eudicots, in addition to functions related to the expression in reproductive structures, may have acquired a function in developmental transitions of vegetative shoots.

Aquilegia as a model system for the evolution and ecology of petals

The ranunculid genus Aquilegia holds extraordinary promise as a model system for investigating a wide range of questions relating to the evolution and ecology of petals. New genetic and genomic resources, including an extensive EST database, BAC libraries and physical maps, as well as virus-induced gene silencing are facilitating this research on multiple fronts. At the developmental genetic level, Aquilegia has been important for elucidating the developmental programme for specifying petals and petaloid characteristics. Data suggest that duplication events among the petal and stamen identity genes have resulted in sub- and neofunctionalization. This expansion of gene function does not include the petaloidy of Aquilegia sepals, however, which does not depend on the same loci that control identity of the second whorl petals. Of special interest is the elaboration of the petal into a nectar spur, a major innovation for the genus. Intra- and interspecific variation in the shape and colour of petals, especially the spurs, has been shown to be adaptative for different pollinators. Thus, understanding the genetic basis of these traits will help us connect the ecological interactions driving speciation with the genetic changes responsible for remodelling morphology. Progress in this area has focused on the multiple, parallel transitions in flower colour and nectar spur length across the genus. For flower colour, upstream transcription factors appear to be primarily targets of natural selection. Thus research in Aquilegia spans the initial evolution of petals and petaloidy to the diversification of petal morphology to the ecological basis of petal form, thereby providing a comprehensive picture of the evolutionary biology of this critical angiosperm feature.

Functional characterization of B class MADS-box transcription factors in Gerbera hybrida

According to the classical ABC model, B-function genes are involved in determining petal and stamen development. Most core eudicot species have B class genes belonging to three different lineages: the PI, euAP3, and TM6 lineages, although both Arabidopsis and Antirrhinum appear to have lost their TM6-like gene. Functional studies were performed for three gerbera (Gerbera hybrida) B class MADS-box genes--PI/GLO-like GGLO1, euAP3 class GDEF2, and TM6-like GDEF1--and data are shown for a second euAP3-like gene, GDEF3. In phylogenetic analysis, GDEF3 is a closely related paralogue of GDEF2, and apparently stems from a duplication common to all Asteraceae. Expression analysis and transgenic phenotypes confirm that GGLO1 and GDEF2 mediate the classical B-function since they determine petal and stamen identities. However, based on assays in yeast, three B class heterodimer combinations are possible in gerbera. In addition to the interaction of GGLO1 and GDEF2 proteins, GGLO1 also pairs with GDEF1 and GDEF3. This analysis of GDEF1 represents the first functional characterization of a TM6-like gene in a core eudicot species outside Solanaceae. Similarly to its relatives in petunia and tomato, the expression pattern and transgenic phenotypes indicate that GDEF1 is not involved in determination of petal identity, but has a redundant role in regulating stamen development.

Analysis of B function in legumes: PISTILLATA proteins do not require the PI motif for floral organ development in Medicago truncatula

The B-class gene PISTILLATA (PI) codes for a MADS-box transcription factor required for floral organ identity in angiosperms. Unlike Arabidopsis, it has been suggested that legume PI genes contribute to a variety of processes, such as the development of floral organs, floral common petal-stamen primordia, complex leaves and N-fixing root nodules. Another interesting feature of legume PI homologues is that some of them lack the highly conserved C-terminal PI motif suggested to be crucial for function. Therefore, legume PI genes are useful for addressing controversial questions on the evolution of B-class gene function, including how they may have diverged in both function and structure to affect different developmental processes. However, functional analysis of legume PI genes has been hampered because no mutation in any B-class gene has been identified in legumes. Here we fill this gap by studying the PI function in the model legume species Medicago truncatula using mutant and RNAi approaches. Like other legume species, M. truncatula has two PI homologues. The expression of the two genes, MtPI and MtNGL9, has strongly diverged, suggesting differences in function. Our analyses show that these genes are required for petal and stamen identity, where MtPI appears to play a predominant role. However, they appear not to be required for development of the nodule, the common primordia or the complex leaf. Moreover, both M. truncatula PI homologues lack the PI motif, which indicates that the C-terminal motif is not essential for PI activity.

Pistillata--duplications as a mode for floral diversification in (Basal) asterids

Basal asterid families, and to a lesser extent the asterids as a whole, are characterized by a high variation in petal and stamen morphology. Moreover, the stamen number, the adnation of stamens to petals, and the degree of sympetaly vary considerably among basal asterid taxa. The B group genes, members of the APETALA3 (AP3) and PISTILLATA (PI) gene lineages, have been shown to specify petal and stamen identities in several core eudicot species. Duplicate genes in these lineages have been shown in some cases to have diversified in their function; for instance in Petunia, a PI paralog is required for the fusion of stamens to the corolla tube, illustrating that such genes belonging to this lineage are not just involved in specifying the identity of the stamens and petals but can also specify novel floral morphologies. This motivated us to study the duplication history of class B genes throughout asterid lineages, which comprise approximately one-third of all flowering plants. The evolutionary history of the PI gene subfamily indicates that the two genes in Petunia result from an ancient duplication event, coinciding with the origin of core asterids. A second duplication event occurred before the speciation of basal asterid Ericales families. These and other duplications in the PI lineage are not correlated with duplications in the AP3 lineage. To understand the molecular evolution of the Ericales PI genes after duplication, we have described their expression patterns using reverse transcription polymerase chain reaction and in situ hybridization, reconstructed how selection shaped their protein sequences and tested their protein interaction specificity with other class B proteins. We find that after duplication, PI paralogs have acquired multiple different expression patterns and negative selective pressure on their codons is relaxed, whereas substitutions in sites putatively involved in protein-protein interactions show positive selection, allowing for a change in the interaction behavior of the PI paralogs after duplication. Together, these observations suggest that the asterids have preferentially recruited PI duplicate genes to diverse and potentially novel roles in asterid flower development.

Evolution of plant MADS box transcription factors: evidence for shifts in selection associated with early angiosperm diversification and concerted gene duplications

Phylogenomic analyses show that gene and genome duplication events have led to the diversification of transcription factor gene families throughout the evolutionary history of land plants and that gene duplications have played an important role in shaping regulatory networks influencing key phenotypic characters including floral development and flowering time. A molecular evolutionary investigation of the mode and tempo of selection acting on the angiosperm MADS box AP1/SQUA, AP3/PI, AG/AGL11, and SEP gene subfamilies revealed site-specific patterns of shifting evolutionary constraint throughout angiosperm history. Specific positions in the four canonical MADS box gene regions, especially K domains and C-terminal regions of all four of these MADS box gene subfamilies exhibited clade-specific shifts in selective constraint following concerted duplication events. Moreover, the frequency of site-specific shifts in constraint was correlated with gene duplications and early angiosperm diversification. We hypothesize that coevolution among interacting MADS box proteins may be responsible for simultaneous increases in the ratio of nonsynonymous to synonymous substitutions (d(N)/d(S) = omega) early in angiosperm history and following concerted duplication events.

Positive selection and ancient duplications in the evolution of class B floral homeotic genes of orchids and grasses

BACKGROUND

Positive selection is recognized as the prevalence of nonsynonymous over synonymous substitutions in a gene. Models of the functional evolution of duplicated genes consider neofunctionalization as key to the retention of paralogues. For instance, duplicate transcription factors are specifically retained in plant and animal genomes and both positive selection and transcriptional divergence appear to have played a role in their diversification. However, the relative impact of these two factors has not been systematically evaluated. Class B MADS-box genes, comprising DEF-like and GLO-like genes, encode developmental transcription factors essential for establishment of perianth and male organ identity in the flowers of angiosperms. Here, we contrast the role of positive selection and the known divergence in expression patterns of genes encoding class B-like MADS-box transcription factors from monocots, with emphasis on the family Orchidaceae and the order Poales. Although in the monocots these two groups are highly diverse and have a strongly canalized floral morphology, there is no information on the role of positive selection in the evolution of their distinctive flower morphologies. Published research shows that in Poales, class B-like genes are expressed in stamens and in lodicules, the perianth organs whose identity might also be specified by class B-like genes, like the identity of the inner tepals of their lily-like relatives. In orchids, however, the number and pattern of expression of class B-like genes have greatly diverged.

RESULTS

The DEF-like genes from Orchidaceae form four well-supported, ancient clades of orthologues. In contrast, orchid GLO-like genes form a single clade of ancient orthologues and recent paralogues. DEF-like genes from orchid clade 2 (OMADS3-like genes) are under less stringent purifying selection than the other orchid DEF-like and GLO-like genes. In comparison with orchids, purifying selection was less stringent in DEF-like and GLO-like genes from Poales. Most importantly, positive selection took place before the major organ reduction and losses in the floral axis that eventually yielded the zygomorphic grass floret.

CONCLUSION

In DEF-like genes of Poales, positive selection on the region mediating interactions with other proteins or DNA could have triggered the evolution of the regulatory mechanisms behind the development of grass-specific reproductive structures. Orchidaceae show a different trend, where gene duplication and transcriptional divergence appear to have played a major role in the canalization and modularization of perianth development.

Using triplet periodicity of nucleotide sequences for finding potential reading frame shifts in genes

We introduce a novel approach for the detection of possible mutations leading to a reading frame (RF) shift in a gene. Deletions and insertions of DNA coding regions are considerable events for genes because an RF shift results in modifications of the extensive region of amino acid sequence coded by a gene. The suggested method is based on the phenomenon of triplet periodicity (TP) in coding regions of genes and its relative resistance to substitutions in DNA sequence. We attempted to extend 326 933 regions of continuous TP found in genes from the KEGG databank by considering possible insertions and deletions. We revealed totally 824 genes where such extension was possible and statistically significant. Then we generated amino acid sequences according to active (KEGG's) and hypothetically ancient RFs in order to find confirmation of a shift at a protein level. Consequently, 64 sequences have protein similarities only for ancient RF, 176 only for active RF, 3 for both and 581 have no protein similarity at all. We aimed to have revealed lower bound for the number of genes in which a shift between RF and TP is possible. Further ways to increase the number of revealed RF shifts are discussed.

Aquilegia: a new model for plant development, ecology, and evolution

The lower eudicot genus Aquilegia holds enormous potential for investigating aspects of development, ecology, and evolution that are otherwise unrepresented among existing model systems. Its evolutionary history is of particular interest because it represents a phylogenetic midpoint between models such as Arabidopsis and Oryza but, at the same time, has experienced a recent adaptive radiation within the genus. To take advantage of these features, a collaborative group has developed a number of genetic and genomic resources for Aquilegia that have facilitated the study of its distinct morphology. This work has demonstrated that although the petaloid sepals of Aquilegia do not depend on B-class genes for their identity, these loci do control development of the petals, stamens, and novel staminodium. Overall, Aquilegia stands as a key example of the potential utility and speed of developing new genetic model systems.

Chapter 4. New model systems for the study of developmental evolution in plants

The number of genetically tractable plant model systems is rapidly increasing, thanks to the decreasing cost of sequencing and the wide amenability of plants to stable transformation and other functional approaches. In this chapter, I discuss emerging model systems from throughout the land plant phylogeny and consider how their unique attributes are contributing to our understanding of development, evolution, and ecology. These new models are being developed using two distinct strategies: in some cases, they are selected because of their close relationship to the established models, while in others, they are chosen with the explicit intention of exploring distantly related plant lineages. Such complementary approaches are yielding exciting new results that shed light on both micro- and macroevolutionary processes in the context of developmental evolution.

Genome-wide analyses of alternative splicing in plants: opportunities and challenges

Alternative splicing (AS) creates multiple mRNA transcripts from a single gene. While AS is known to contribute to gene regulation and proteome diversity in animals, the study of its importance in plants is in its early stages. However, recently available plant genome and transcript sequence data sets are enabling a global analysis of AS in many plant species. Results of genome analysis have revealed differences between animals and plants in the frequency of alternative splicing. The proportion of plant genes that have one or more alternative transcript isoforms is approximately 20%, indicating that AS in plants is not rare, although this rate is approximately one-third of that observed in human. The majority of plant AS events have not been functionally characterized, but evidence suggests that AS participates in important plant functions, including stress response, and may impact domestication and trait selection. The increasing availability of plant genome sequence data will enable larger comparative analyses that will identify functionally important plant AS events based on their evolutionary conservation, determine the influence of genome duplication on the evolution of AS, and discover plant-specific cis-elements that regulate AS. This review summarizes recent analyses of AS in plants, discusses the importance of further analysis, and suggests directions for future efforts.

Extraction of high-quality tissue-specific RNA from London plane trees (Platanus acerifolia), permitting the construction of a female inflorescence cDNA library

The London plane tree (Platanus acerifolia Willd.) has global importance as an urban landscaping tree and is the subject of genetic-improvement programs for productive sterility, disease and/or insect resistance. Molecular analysis techniques are crucial to such programs, but may be impeded by specific difficulties encountered during nucleic acid isolation. A detailed RNA isolation and purification protocol, based on established cetyltrimethyl-ammonium bromide (CTAB) extraction techniques combined with additional purification steps using butanol and the ionic detergent CTAB, which overcomes these problems in the woody species P. acerifolia, was conducted. In short, phenolic compounds are bound to soluble polyvinylpyrrolidone and then separated out through LiCl precipitation of the RNA. Subsequently, protein- and carbohydrate-contaminants are removed by chloroform partitioning followed by LiCl-mediated precipitation. The resulting isolates of RNA were found to be of sufficient quality for successful use in reverse transcription PCR analysis. Furthermore, RNA isolates from female inflorescences were used for the construction of a cDNA library. This library was found to contain several full-length cDNA clones of MADS-box genes, consistent with the library being representative of inflorescence expression profiles.

B-class MADS-box genes in trioecious papaya: two paleoAP3 paralogs, CpTM6-1 and CpTM6-2, and a PI ortholog CpPI

In the ABC model of flower development, B function organ-identity genes act in the second and third whorls of the flower to control petal and stamen identity. The trioecious papaya has male, female, and hermaphrodite flowers and is an ideal system for testing the B-class gene expression patterns in trioecious plants. We cloned papaya B-class genes, CpTM6-1, CpTM6-2, and CpPI, using MADS box gene specific degenerate primers followed by cDNA library screening and sequencing of positive clones. While phylogenetic analyses show that CpPI is the ortholog of the Arabidopsis gene PI, the CpTM6-1 and CpTM6-2 loci are representatives of the paralogous TM6 lineage that contain paleoAP3 motifs unlike the euAP3 gene observed in Arabidopsis. These two paralogs appeared to have originated from a tandem duplication occurred approximately 13.4 million year ago (mya) (bootstrap range 13.36 +/- 2.42). In-situ hybridization and RT-PCR showed that the papaya B-class genes were highly expressed in young flowers across all floral organ primordia. As the flower organs developed, all three B-class genes were highly expressed in petals of all three-sex types and in stamens of hermaphrodite and male flowers. CpTM6-1 expressed at low levels in sepals and carpels, whereas CpTM6-2 expressed at a low level in sepals and at a high level in leaves. Our results showed that B-class gene homologs could function as predicted by the ABC model in trioecous flowers but differential expressions of CpTM6-1, and CpTM6-2, and CpPI suggested the diversification of their functions after the duplication events.

The MIK region rather than the C-terminal domain of AP3-like class B floral homeotic proteins determines functional specificity in the development and evolution of petals

In core eudicots, euAP3-type MADS-box genes encode a PISTILLATA (PI)-derived motif, as well as a C-terminal euAP3 motif that originated from a paleoAP3 motif of an ancestral APETALA3 (AP3)-like protein through a translational frameshift mutation. To determine the functional and evolutionary relevance of these motifs, a series of point mutation and domain-swap constructs were generated, involving CsAP3, a paleoAP3-type gene from the basal angiosperm Chloranthus spicatus encoding a truncated paleoAP3 motif, and AtAP3, a euAP3-type gene from the core eudicot Arabidopsis thaliana. The chimeric constructs were expressed in A. thaliana under the control of the AP3 promoter or the CaMV 35S promoter in an ap3 mutant or wild-type background, respectively. Significant recovery of AP3 function was obtained in both complementation and ectopic expression experiments whenever the region upstream of the C-terminal motifs (MIK region) from A. thaliana was taken, even when the PI-derived motif and the truncated paleoAP3 motif of CsAP3 substituted for the corresponding sequences from AtAP3. However, no or very weak complementation or gain-of-function was seen when the MIK region was from CsAP3. Our data suggest that changes in the MIK region rather than mutations in the C-terminal domain were of crucial importance for the evolution of the functional specificity of euAP3-type proteins in stamen and petal development.

Conserved C-terminal motifs of the Arabidopsis proteins APETALA3 and PISTILLATA are dispensable for floral organ identity function

The B-class genes APETALA3 (AP3) and PISTILLATA (PI) in Arabidopsis (Arabidopsis thaliana) and their orthologs in other species have been the focus of studies to elucidate the development of petals and stamens in angiosperm flowers. Evolutionary analysis indicates that B-class genes have undergone multiple gene duplication events in angiosperms. The resultant B-class lineages are characterized by short, conserved amino acid sequences at the extreme C-terminal end of the B-class proteins. AP3 is a member of the euAP3 lineage that contains both the euAP3 and PI-derived motifs at the C terminus. PI is a member of the PI lineage that contains the C-terminal PI motif at the C terminus. Despite conservation over a wide evolutionary distance, the function of C-terminal motifs is not well understood. In this study, we demonstrate that truncated forms of AP3 and PI, which lack the conserved C-terminal motifs, function to direct floral organ identity specification in Arabidopsis plants. By contrast, larger truncations, which remove the third putative amphipathic alpha-helix in the K domain of AP3 or PI, are nonfunctional. We conclude that the euAP3 and PI-derived motifs of AP3 and the PI motif of PI are not essential for floral organ identity function of AP3 and PI in Arabidopsis.

Are petals sterile stamens or bracts? The origin and evolution of petals in the core eudicots

BACKGROUND

The aim of this paper is to discuss the controversial origins of petals from tepals or stamens and the links between the morphological expression of petals and floral organ identity genes in the core eudicots.

SCOPE

I challenge the widely held classical view that petals are morphologically derived from stamens in the core eudicots, and sepals from tepals or bracts. Morphological data suggest that tepal-derived petals have evolved independently in the major lineages of the core eudicots (i.e. asterids, Santalales and rosids) from Berberidopsis-like prototypes, and that staminodial petals have arisen only in few isolated cases where petals had been previously lost (Caryophyllales, Rosales). The clear correlation between continuous changes in petal morphology, and a scenario that indicates numerous duplications to have taken place in genes controlling floral organ development, can only be fully understood within a phylogenetic context. B-gene expression plays a fundamental role in the evolution of the petals by controlling petaloidy, but it does not clarify petal homology.

CONCLUSIONS

An increased synorganization of the flower in the core eudicots linked with the establishment of floral whorls restricts the petaloid gene expression to the second whorl, reducing the similarities of petals with tepals from which they were originally derived. An increased flower size linked with secondary polyandry or polycarpelly may lead to a breakdown of the restricted gene expression and a reversal to ancestral characteristics of perianth development. An altered 'sliding boundary' hypothesis is proposed for the core eudicots to explain shifts in petaloidy of the perianth and the event of staminodial petals. The repetitive changes of function in the perianth of the core eudicots are linked with shifts in petaloidy to the outer perianth whorl, or losses of petal or sepal whorls that can be secondarily compensated for by the inclusion of bracts in the flower. The origin and evolution of petals appears to be as complex on a molecular basis as it is from a morphological point of view.