Analysis of the floral transcriptome uncovers new regulators of organ determination and gene families related to flower organ differentiation in Gerbera hybrida (Asteraceae)

Polyketide reductases in defense-related parasorboside biosynthesis in Gerbera hybrida share processing strategies with microbial polyketide synthase systems

Plant polyketides are well-known for their crucial functions in plants and their importance in the context of human health. They are synthesized by type III polyketide synthases (PKSs) and their final functional diversity is determined by post-PKS tailoring enzymes. Gerbera hybrida is rich in two defense-related polyketides: gerberin and parasorboside. Their synthesis is known to be initiated by GERBERA 2-PYRONE SYNTHASE 1 (G2PS1), but the polyketide reductases (PKRs) that determine their final structure have not yet been identified. We identified two PKR candidates in the pathway, GERBERA REDUCTASE 1 (GRED1) and GRED2. Gene expression and metabolite analysis of different gerbera tissues, cultivars, and transgenic gerbera plants, and in vitro enzyme assays, were performed for functional characterization of the enzymes. GRED1 and GRED2 catalyze the second reduction step in parasorboside biosynthesis. They reduce the proximal keto domain of the linear CoA bound intermediate before lactonization. We identified a crucial tailoring step in an important gerbera PKS pathway and show that plant polyketide biosynthesis shares processing strategies with fungi and bacteria. The two tailoring enzymes are recruited from the ancient sporopollenin biosynthetic pathway to a defense-related PKS pathway in gerbera. Our data provide an example of how plants recruit conserved genes to new functions in secondary metabolism that are important for environmental adaptation.

Integrative Analysis of miRNAs and Their Targets Involved in Ray Floret Growth in Gerbera hybrida

MicroRNAs (miRNAs) are involved in regulating many aspects of plant growth and development at the post-transcriptional level. Gerbera (Gerbera hybrida) is an important ornamental crop. However, the role of miRNAs in the growth and development of gerbera is still unclear. In this study, we used high-throughput sequencing to analyze the expression profiles of miRNAs in ray floret during inflorescence opening. A total of 164 miRNAs were obtained, comprising 24 conserved miRNAs and 140 novel miRNAs. Ten conserved and 15 novel miRNAs were differentially expressed during ray floret growth, and 607 differentially expressed target genes of these differentially expressed miRNAs were identified using psRNATarget. We performed a comprehensive analysis of the expression profiles of the miRNAs and their targets. The changes in expression of five miRNAs (ghy-miR156, ghy-miR164, ghy-miRn24, ghy-miRn75 and ghy-miRn133) were inversely correlated with the changes in expression of their eight target genes. The miRNA cleavage sites in candidate target gene mRNAs were determined using 5′-RLM-RACE. Several miRNA-mRNA pairs were predicted to regulate ray floret growth and anthocyanin biosynthesis. In conclusion, the results of small RNA sequencing provide valuable information to reveal the mechanisms of miRNA-mediated ray floret growth and anthocyanin accumulation in gerbera.

Sequencing and analysis of gerbera daisy leaf transcriptomes reveal disease resistance and susceptibility genes differentially expressed and associated with powdery mildew resistance

BACKGROUND

RNA sequencing has been widely used to profile genome-wide gene expression and identify candidate genes controlling disease resistance and other important traits in plants. Gerbera daisy is one of the most important flowers in the global floricultural trade, and powdery mildew (PM) is the most important disease of gerbera. Genetic improvement of gerbera PM resistance has become a crucial goal in gerbera breeding. A better understanding of the genetic control of gerbera resistance to PM can expedite the development of PM-resistant cultivars.

RESULTS

The objectives of this study were to identify gerbera genotypes with contrasting phenotypes in PM resistance and sequence and analyze their leaf transcriptomes to identify disease resistance and susceptibility genes differentially expressed and associated with PM resistance. An additional objective was to identify SNPs and SSRs for use in future genetic studies. We identified two gerbera genotypes, UFGE 4033 and 06-245-03, that were resistant and susceptible to PM, respectively. De novo assembly of their leaf transcriptomes using four complementary pipelines resulted in 145,348 transcripts with a N50 of 1124 bp, of which 67,312 transcripts contained open reading frames and 48,268 were expressed in both genotypes. A total of 494 transcripts were likely involved in disease resistance, and 17 and 24 transcripts were up- and down-regulated, respectively, in UFGE 4033 compared to 06-245-03. These gerbera disease resistance transcripts were most similar to the NBS-LRR class of plant resistance genes conferring resistance to various pathogens in plants. Four disease susceptibility transcripts (MLO-like) were expressed only or highly expressed in 06-245-03, offering excellent candidate targets for gene editing for PM resistance in gerbera. A total of 449,897 SNPs and 19,393 SSRs were revealed in the gerbera transcriptomes, which can be a valuable resource for developing new molecular markers.

CONCLUSION

This study represents the first transcriptomic analysis of gerbera PM resistance, a highly important yet complex trait in a globally important floral crop. The differentially expressed disease resistance and susceptibility transcripts identified provide excellent targets for development of molecular markers and genetic maps, cloning of disease resistance genes, or targeted mutagenesis of disease susceptibility genes for PM resistance in gerbera.

Cis-regulatory differences in isoform expression associate with life history strategy variation in Atlantic salmon

A major goal in biology is to understand how evolution shapes variation in individual life histories. Genome-wide association studies have been successful in uncovering genome regions linked with traits underlying life history variation in a range of species. However, lack of functional studies of the discovered genotype-phenotype associations severely restrains our understanding how alternative life history traits evolved and are mediated at the molecular level. Here, we report a cis-regulatory mechanism whereby expression of alternative isoforms of the transcription co-factor vestigial-like 3 (vgll3) associate with variation in a key life history trait, age at maturity, in Atlantic salmon (Salmo salar). Using a common-garden experiment, we first show that vgll3 genotype associates with puberty timing in one-year-old salmon males. By way of temporal sampling of vgll3 expression in ten tissues across the first year of salmon development, we identify a pubertal transition in vgll3 expression where maturation coincided with a 66% reduction in testicular vgll3 expression. The late maturation allele was not only associated with a tendency to delay puberty, but also with expression of a rare transcript isoform of vgll3 pre-puberty. By comparing absolute vgll3 mRNA copies in heterozygotes we show that the expression difference between the early and late maturity alleles is largely cis-regulatory. We propose a model whereby expression of a rare isoform from the late allele shifts the liability of its carriers towards delaying puberty. These results exemplify the potential importance of regulatory differences as a mechanism for the evolution of life history traits.

Alternative life history strategies are an important source of diversity within populations and promote the maintenance of adaptive capacity and population resilience. However, in many cases the molecular basis of different life history strategies remains elusive. Age at maturity is a key adaptive life history trait in Atlantic salmon and has a relatively simple genetic basis. Using salmon age at maturity as a model, we report a mechanism whereby different transcript isoforms of the key age at maturity gene, vestigial-like 3 (vgll3), associate with variation in the timing of male puberty. Our results show how gene regulatory differences in conjunction with variation in gene transcript structure can encode for complex alternative life histories.

Analyzing Differentially Expressed Genes and Pathways Associated with Pistil Abortion in Japanese Apricot via RNA-Seq

Reproduction is a critical stage in the flower development process, and its failure causes serious problems affecting fruit quality and yield. Pistil abortion is one of the main factors in unsuccessful reproduction and occurs in many fruit plants. In Japanese apricot, the problem of pistil abortion is very common and affects fruit quality and plant yield; however, its molecular mechanism is not clearly understood. Therefore, in the current study, we used RNA-Seq to identify the differentially expressed genes (DEGs) and pathways actively involved in pistil abortion. A total of 3882 differentially expressed genes were found after cutoff and pairwise comparison analysis. According to KEGG pathway analysis, plant hormone signaling transduction and metabolic pathways were found most significantly enriched in this study. A total of 60 transcription factor families such as MADS-box, NAC and TCP showed their role in this process. RT-qPCR assays confirmed that the expression levels were consistent with RNA-Seq results. This study provides an alternative to be considered for further studies and understanding of pistil abortion processes in Japanese apricot, and it provides a reference related to this issue for other deciduous fruit crops.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

Genome-Wide Association Mapping of Floral Traits in Cultivated Sunflower (Helianthus annuus)

Floral morphology and pigmentation are both charismatic and economically relevant traits associated with cultivated sunflower (Helianthus annuus L.). Recent work has linked floral morphology and pigmentation to pollinator efficiency and seed yield. Understanding the genetic architecture of such traits is essential for crop improvement, and gives insight into the role of genetic constraints in shaping floral diversity. A diversity panel of 288 sunflower genotypes was phenotyped for a variety of morphological, phenological, and color traits in both a greenhouse and a field setting. Association mapping was performed using 5788 SNP markers using a mixed linear model approach. Several dozen markers across 10 linkage groups were significantly associated with variation in morphological and color trait variation. Substantial trait plasticity was observed between greenhouse and field phenotyping, and associations differed between environments. Color traits mapped more strongly than morphology in both settings, with markers together explaining 16% of petal carotenoid content in the greenhouse, and 17% and 24% of variation in disc anthocyanin presence in the field and greenhouse, respectively. Morphological traits like disc size mapped more strongly in the field, with markers together explaining up to 19% of disc size variation. Loci identified here through association mapping within cultivated germplasm differ from those identified through biparental crosses between modern cultivated sunflower and either its wild progenitor or domesticated landraces. Several loci lie within genomic regions involved in domestication. Differences between phenotype expression under greenhouse and field conditions highlight the importance of plasticity in determining floral morphology and pigmentation.

Non-coding RNAs having strong positive interaction with mRNAs reveal their regulatory nature during flowering in a wild relative of pigeonpea (Cajanus scarabaeoides)

In higher plants, flower development is a result of crosstalk between many factors like photoperiod, vernalization, hormone concentration, epigenetic modification etc. and is also regulated by non-coding RNAs (ncRNAs). In the present study, we are reporting the involvement of long non-coding RNAs (lncRNAs) and miRNAs during the process of flower development in Cajanus scarabaeoides, an important wild relative of pigeonpea. The transcriptome of floral and leaf tissues revealed a total of 1672 lncRNAs and 57 miRNAs being expressed during flower development. Prediction analysis of identified lncRNAs showed that 1593 lncRNAs were targeting 3420 mRNAs and among these, 98 were transcription factors (TFs) belonging to 48 groups. All the identified 57 miRNAs were novel, suggesting their genera specificity. Prediction of the secondary structure of lncRNAs and miRNAs followed by interaction analysis revealed that 199 lncRNAs could interact with 47 miRNAs where miRNAs were acting in the root of interaction. Gene Ontology of the ncRNAs and their targets showed the potential role of lncRNAs and miRNAs in the flower development of C. scarabaeoides. Among the identified interactions, 17 lncRNAs were endogenous target mimics (eTMs) for miRNAs that target flowering-related transcription factors. Expression analysis of identified transcripts revealed that higher expression of Csa-lncRNA_1231 in the bud sequesters Csa-miRNA-156b by indirectly mimicking the miRNA and leading to increased expression of flower-specific SQUAMOSA promoter-binding protein-like (SPL-12) TF indicating their potential role in flower development. The present study will help in understanding the molecular regulatory mechanism governing the induction of flowering in C. scarabaeoides.

Transcriptomic response to cold of thermophilous medicinal plant Marsdenia tenacissima

Extracts from Marsdenia tenacissima, involving tenacissoside H, I and G, have been used as remedies of cancer, inflammation and asthma. Low temperature serves as one of the main factors constrain the planting expansion and quality of M. tenacissima, but its functional mechanism has been known scarcely for the lack of genomic information and transcriptional profile. Here we investigated the transcriptomic responses of M. tenacissima under cold stress to gain insight into the molecular mechanism of low temperature sensitivity. Total RNAs were collected from samples obtained at 4-time points (after 0, 3, 6 and 48 h cold treatments with 4 °C, respectively), then used for library construction and sequenced on the Illumina Hiseq™ 4000 platform. Passing quality assessments, 500794 transcripts, and 206137 unigenes were de novo assembly out in Trinity v2.4.0, holding contig N50 of 2566 bp and unigene mean length of 754 bp. 44.20% of assembled unigenes were annotated to the well-known public protein database on a basis of sequence similarity. Using statistical comparison of the fragments per kilo base of transcript per million reads mapped (FPKM) values between conditions, 6082 group-specific differentially expressed genes (DEGs) were identified and considered as cold-responsive genes, which contained copious transcription factors and active secondary metabolism. Among them, 43 unigenes were constantly up-regulated expression along with cold time, which mainly implicated in the biosynthesis of secondary metabolites, carbon metabolism, RNA and DNA metabolism. Conversely, 21 unigenes involved in photosynthesis, cell wall, protein degradation, and transporters were downregulated continually with cold timescale. Experimentally, MtEF1α was chosen as the best housekeeping gene. Functional enrichments found that damaging of cold stress on M. tenacissima may be ascribed to inability of photosynthesis, ribsome processing, flavonoid biosynthesis and terpenoids degradation. Correlation analysis between cold induced transcription factors and tenacissoside biosynthesis-related genes indicated that 3β-HSD significant positively correlated with bHLH51, and 4-MSO with NF-YB, GRAS3, Trihelix, FAR1, MYB60, MYBS1, bZIP43. Further promoter clone found MYB-binding site in the promoter of 4-MSO. In view of the reported cold tolerance of MYB60, it is recommended as a potential candidate suitable for future molecular design of exaptation cultivation with high bioactive constituents.

The clock gene Gigantea 1 from Petunia hybrida coordinates vegetative growth and inflorescence architecture

The gene GIGANTEA (GI) appeared early in land plants. It is a single copy gene in most plants and is found in two to three copies in Solanaceae. We analyzed the silencing of one GI copy, Petunia hybrida GI1 (PhGI1), by hairpin RNAs in Petunia in order to gain knowledge about its range of functions. Decreased transcript levels of PhGI1 were accompanied also by a reduction of PhGI2. They were further associated with increased time period between two consecutive peaks for PhGI1 and CHANEL (PhCHL), the orthologue of the blue light receptor gene ZEITLUPE (ZTL), confirming its role in maintaining circadian rhythmicity. Silenced plants were bigger with modified internode length and increased leaf size while flowering time was not altered. We uncovered a new function for PhGI1 as silenced plants showed reduction of flower bud number and the appearance of two flower buds in the bifurcation point, were normally one flower bud and the inflorescence meristem separate. Furthermore, one of the flower buds consistently showed premature flower abortion. Flowers that developed fully were significantly smaller as a result of decreased cell size. Even so the circadian pattern of volatile emission was unchanged in the silenced lines, flowers emitted 20% less volatiles on fresh weight basis over 24 hours and showed changes in the scent profile. Our results indicate a novel role of PhGI1 in the development of reproductive organs in Petunia. PhGI1 therefore represses growth in vegetative plant parts, maintains the typical cymose inflorescence structure, and inhibits premature flower abortion.

RNA-Seq analysis reveals an essential role of tyrosine metabolism pathway in response to root-rot infection in Gerbera hybrida

Gerbera hybrida is one of the top five cut flowers across the world, it is host for the root rot causing parasite called Phytophthora cryptogea. In this study, plantlets of healthy and root-rot pathogen-infected G. hybrida were used as plant materials for transcriptome analyis using high-throughput Illumina sequencing technique. A total 108,135 unigenes were generated with an average length of 727 nt and N50 equal to 1274 nt out of which 611 genes were identified as DEGs by DESeq analyses. Among DEGs, 228 genes were up-regulated and 383 were down-regulated. Through this annotated data and Kyoto encyclopedia of genes and genomes (KEGG), molecular interaction network, transcripts accompanying with tyrosine metabolism, phenylalanine, tyrosine, and tryptophan biosynthesis, phenylpropanoid and flavonoid biosynthesis, and plant hormone signal transduction pathways were thoroughly observed considering expression pattern. The involvement of DEGs in tyrosine metabolism pathway was validated by real-time qPCR. We found that genes related with tyrosine metabolism were activated and up-regulated against stress response. The expression of GhTAT, GhAAT, GhHPD, GhHGD and GhFAH genes was significantly increased in the leaves and petioles at four and six dpi (days post inoculation) as compared with control. The study predicts the gene sequences responsible for the tyrosine metabolism pathway and its responses against root-rot resistance in gerbera plant. In future, identification of such genes is necessary for the better understanding of rot resistance mechanism and to develop a root rot resistance strategy for ornamental plants.

Differential expression of anthocyanin structural genes and transcription factors determines coloration patterns in gerbera flowers

We investigated the expression of anthocyanin structural genes and transcription factors (TFs) associated with varying anthocyanin content during different developmental stages (S1-S4) of the gerbera cultivars 'Nathasha' and 'Rosalin'. Accumulation of anthocyanin started at S1 and reached a maximum at S3 in both cultivars. Enhancement of anthocyanin content in 'Nathasha' was associated with upregulation of ANS and MYB10, whereas in 'Rosalin', upregulation was associated with CHS1, MYB10, and MYC1. Low-temperature exposure (6 °C) enhanced anthocyanin content to a greater extent than that at 22 °C via stronger upregulation of CHS1 and MYB10 in 'Nathasha' and CHS1 in 'Rosalin', irrespective of flower developmental stage. However, differences in anthocyanin content between the two cultivars were found to be influenced by the expression levels of all structural genes and TFs, irrespective of flower developmental stage and temperature conditions. We suggest that differences in the regulation mechanisms of anthocyanin biosynthesis and coloration pattern between 'Nathasha' and 'Rosalin' are related to differences in the expression patterns of structural genes and TFs; however, further functional studies of the key genes in anthocyanin biosynthesis are needed.

Transcriptional control of anthocyanin biosynthesis genes and transcription factors associated with flower coloration patterns in Gerbera hybrida

We analyzed the expression of anthocyanin biosynthesis genes and transcription factors (TFs) in the Gerbera hybrida cultivars 'Bintang' and 'Alliance' that exhibit different coloration patterns. Differential expression of biosynthesis genes and TFs was associated with variable anthocyanin content at different flower developmental stages (S1-S3) in both cultivars; higher anthocyanin content was correlated with higher levels of gene expression. Exposure to different temperatures (6 and 22 °C) also resulted in different anthocyanin content levels: the lower temperature (6 °C) enhanced anthocyanin content compared to the higher temperature (22 °C). However, the increased anthocyanin content of 'Bintang' compared to 'Alliance' was the result of higher levels of expression of all detected genes, regardless of flower stage and temperature conditions. Therefore, we conclude that transcriptional control of the detected genes is associated with the mechanisms of anthocyanin biosynthesis and coloration patterns in gerberas; however, further studies of the key genes are needed.

A Mini Zinc-Finger Protein (MIF) from Gerbera hybrida Activates the GASA Protein Family Gene, GEG, to Inhibit Ray Petal Elongation

Petal appearance is an important horticultural trail that is generally used to evaluate the ornamental value of plants. However, knowledge of the molecular regulation of petal growth is mostly derived from analyses of Arabidopsis thaliana, and relatively little is known about this process in ornamental plants. Previously, GEG (Gerbera hybrida homolog of the gibberellin [GA]-stimulated transcript 1 [GAST1] from tomato), a gene from the GA stimulated Arabidopsis (GASA) family, was reported to be an inhibitor of ray petal growth in the ornamental species, G. hybrida. To explore the molecular regulatory mechanism of GEG in petal growth inhibition, a mini zinc-finger protein (MIF) was identified using yeast one-hybrid (Y1H) screen. The direct binding of GhMIF to the GEG promoter was verified by using an electrophoretic mobility shift assay and a dual-luciferase assay. A yeast two-hybrid (Y2H) revealed that GhMIF acts as a transcriptional activator. Transient transformation assay indicated that GhMIF is involved in inhibiting ray petal elongation by activating the expression of GEG. Spatiotemporal expression analyses and hormone treatment assay showed that the expression of GhMIF and GEG is coordinated during petal development. Taken together, these results suggest that GhMIF acts as a direct transcriptional activator of GEG, a gene from the GASA protein family to regulate the petal elongation.

Two polyketide synthases are necessary for 4-hydroxy-5-methylcoumarin biosynthesis in Gerbera hybrida

Gerbera (Gerbera hybrida) is an economically important ornamental species and a model plant of the Asteraceae family for flower development and secondary metabolism. Gerberin and parasorboside, two bitter tasting glucosidic lactones, are produced in high amounts in nearly all gerbera tissues. Gerbera and its close relatives also produce a rare coumarin, 4-hydroxy-5-methylcoumarin (HMC). Unlike most coumarins, 5-methylcoumarins have been suggested to be derived through the acetate-malonate pathway. All of these polyketide-derived glucosylated molecules are considered to have a role in defense against herbivores and phytopathogens in gerbera. Gerbera expresses three genes encoding 2-pyrone synthases (G2PS1-3). The enzymes are chalcone synthase-like polyketide synthases with altered starter substrate specificity. We have shown previously that G2PS1 is responsible for the synthesis of 4-hydroxy-6-methyl-2-pyrone (triacetolactone), a putative precursor of gerberin and parasorboside. Here we show that polyketide synthases G2PS2 and G2PS3 are necessary for the biosynthesis of HMC in gerbera, and that a reductase enzyme is likely required to complete the pathway to HMC. G2PS2 is expressed in the leaf blade and inflorescences of gerbera, while G2PS3 is strictly root specific. Heterologous expression of G2PS2 or G2PS3 in tobacco leads to the formation of 4,7-dihydroxy-5-methylcoumarin, apparently an unreduced precursor of HMC, while ectopic expression in gerbera leads to HMC formation in tissues where nontransgenic tissue does not express the genes and does not accumulate the compound. Using protein modelling and site-directed mutagenesis we identified the residues I203 and T344 in G2PS2 and G2PS3 to be critical for pentaketide synthase activity.

Proteomic insights into floral biology

The flower is the most important biological structure for ensuring angiosperms reproductive success. Not only does the flower contain critical reproductive organs, but the wide variation in morphology, color, and scent has evolved to entice specialized pollinators, and arguably mankind in many cases, to ensure the successful propagation of its species. Recent proteomic approaches have identified protein candidates related to these flower traits, which has shed light on a number of previously unknown mechanisms underlying these traits. This review article provides a comprehensive overview of the latest advances in proteomic research in floral biology according to the order of flower structure, from corolla to male and female reproductive organs. It summarizes mainstream proteomic methods for plant research and recent improvements on two dimensional gel electrophoresis and gel-free workflows for both peptide level and protein level analysis. The recent advances in sequencing technologies provide a new paradigm for the ever-increasing genome and transcriptome information on many organisms. It is now possible to integrate genomic and transcriptomic data with proteomic results for large-scale protein characterization, so that a global understanding of the complex molecular networks in flower biology can be readily achieved. This article is part of a Special Issue entitled: Plant Proteomics--a bridge between fundamental processes and crop production, edited by Dr. Hans-Peter Mock.

Transcriptomic Analysis of Differentially Expressed Genes during Flower Organ Development in Genetic Male Sterile and Male Fertile Tagetes erecta by Digital Gene-Expression Profiling

Tagetes erecta is an important commercial plant of Asteraceae family. The male sterile (MS) and male fertile (MF) two-type lines of T. erecta have been utilized in F1 hybrid production for many years, but no report has been made to identify the genes that specify its male sterility that is caused by homeotic conversion of floral organs. In this study, transcriptome assembly and digital gene expression profiling were performed to generate expression profiles of MS and MF plants. A cDNA library was generated from an equal mixture of RNA isolated from MS and MF flower buds (1 mm and 4 mm in diameter). Totally, 87,473,431 clean tags were obtained and assembled into 128,937 transcripts among which 65,857 unigenes were identified with an average length of 1,188 bp. About 52% of unigenes (34,176) were annotated in Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KO (KEGG Ortholog database) and/or GO. Taking the above transcriptome as reference, 125 differentially expressed genes were detected in both developmental stages of MS and MF flower buds. MADS-box genes were presumed to be highly related to male sterility in T. erecta based on histological and cytological observations. Twelve MADS-box genes showed significantly different expression levels in flower buds 4 mm in diameter, whereas only one gene expressed significantly different in flower buds 1 mm in diameter between MS and MF plants. This is the first transcriptome analysis in T. erecta and will provide a valuable resource for future genomic studies, especially in flower organ development and/or differentiation.

Transcriptome Analysis of Gerbera hybrida Including in silico Confirmation of Defense Genes Found

For the ornamental crop Gerbera hybrida, breeding at the moment is done using conventional methods. As this has drawbacks in breeding speed and efficiency, especially for complex traits like disease resistance, we set out to develop genomic resources. The leaf and flower bud transcriptomes of four parents, used to generate two gerbera populations, were sequenced using Illumina paired-end sequencing. In total, 36,770 contigs with an average length of 1397 bp were generated and these have been the starting point for SNP identification and annotation. The consensus contig sequences were used to map reads of individual parents, to identify genotype specific SNPs, and to assess the presence of common SNPs between genotypes. Comparison with the non-redundant protein database (nr) showed that 29,146 contigs gave BLAST hits. Of sequences with blast results, 73.3% obtained a clear gene ontology (GO) annotation. EST contigs coding for enzymes were found in Kyoto Encyclopedia of Genes and Genomes maps (KEGG). Through, these annotated data and KEGG molecular interaction network, transcripts associated with the phenylpropanoid metabolism, other secondary metabolite biosynthesis pathways, phytohormone biosynthesis and signal transduction were analyzed in more detail. Identifying genes involved in these processes could provide genetic and genomic resources for studying the mechanism of disease resistance in gerbera.

Functional diversification of duplicated chalcone synthase genes in anthocyanin biosynthesis of Gerbera hybrida

• Chalcone synthase (CHS) is the key enzyme in the first committed step of the flavonoid biosynthetic pathway and catalyzes the stepwise condensation of 4-coumaroyl-CoA and malonyl-CoA to naringenin chalcone. In plants, CHS is often encoded by a small family of genes that are temporally and spatially regulated. Our earlier studies have shown that GCHS4 is highly activated by ectopic expression of an MYB-type regulator GMYB10 in gerbera (Gerbera hybrida). • The tissue- and development-specific expression patterns of three gerbera CHS genes were examined. Virus-induced gene silencing (VIGS) was used to knock down GCHS1 and GCHS4 separately in gerbera inflorescences. • Our data show that GCHS4 is the only CHS encoding gene that is expressed in the cyanidin-pigmented vegetative tissues of gerbera cv Terraregina. GCHS3 expression is pronounced in the pappus bristles of the flowers. Expression of both GCHS1 and GCHS4 is high in the epidermal cells of gerbera petals, but only GCHS1 is contributing to flavonoid biosynthesis. • Gerbera contains a family of three CHS encoding genes showing different spatial and temporal regulation. GCHS4 expression in gerbera petals is regulated post-transcriptionally, at the level of either translation elongation or protein stability.

Establishing genomic tools and resources for Guizotia abyssinica (L.f.) Cass.-the development of a library of expressed sequence tags, microsatellite loci, and the sequencing of its chloroplast genome

We present an EST library, chloroplast genome sequence, and nuclear microsatellite markers that were developed for the semi-domesticated oilseed crop noug (Guizotia abyssinica) from Ethiopia. The EST library consists of 25 711 Sanger reads, assembled into 17 538 contigs and singletons, of which 4781 were functionally annotated using the Arabidopsis Information Resource (TAIR). The age distribution of duplicated genes in the EST library shows evidence of two paleopolyploidizations-a pattern that noug shares with several other species in the Heliantheae tribe (Compositae family). From the EST library, we selected 43 microsatellites and then designed and tested primers for their amplification. The number of microsatellite alleles varied between 2 and 10 (average 4.67), and the average observed and expected heterozygosities were 0.49 and 0.54, respectively. The chloroplast genome was sequenced de novo using Illumina's sequencing technology and completed with traditional Sanger sequencing. No large re-arrangements were found between the noug and sunflower chloroplast genomes, but 1.4% of sites have indels and 1.8% show sequence divergence between the two species. We identified 34 tRNAs, 4 rRNA sequences, and 80 coding sequences, including one region (trnH-psbA) with 15% sequence divergence between noug and sunflower that may be particularly useful for phylogeographic studies in noug and its wild relatives.

A global view of transcriptome dynamics during flower development in chickpea by deep sequencing

Measurement of gene expression can provide important clues about gene function and molecular basis of developmental processes. Here, we have analysed the chickpea transcriptome in vegetative and flower tissues by exploiting the potential of high-throughput sequencing to measure gene expression. We mapped more than 295 million reads to quantify the transcript abundance during flower development. We detected the expression of more than 90% genes in at least one tissue analysed. We found quite a large number of genes were differentially expressed during flower development as compared to vegetative tissues. Further, we identified several genes expressed in a stage-specific manner. Various transcription factor families and metabolic pathways involved in flower development were elucidated. The members of MADS-box family were most represented among the transcription factor genes up-regulated during various stages of flower development. The abundant expression of several well-known genes implicated in flower development in chickpea flower development stages confirmed our results. In addition, we detected the expression specificities of lineage-specific genes during flower development. The expression data presented in this study is the most comprehensive dataset available for chickpea as of now and will serve as resource for unraveling the functions of many specific genes involved in flower development in chickpea and other legumes.

Next-generation sequencing of the Chrysanthemum nankingense (Asteraceae) transcriptome permits large-scale unigene assembly and SSR marker discovery

Background

Simple sequence repeats (SSRs) are ubiquitous in eukaryotic genomes. Chrysanthemum is one of the largest genera in the Asteraceae family. Only few Chrysanthemum expressed sequence tag (EST) sequences have been acquired to date, so the number of available EST-SSR markers is very low.

Methodology/Principal Findings

Illumina paired-end sequencing technology produced over 53 million sequencing reads from C. nankingense mRNA. The subsequent de novo assembly yielded 70,895 unigenes, of which 45,789 (64.59%) unigenes showed similarity to the sequences in NCBI database. Out of 45,789 sequences, 107 have hits to the Chrysanthemum Nr protein database; 679 and 277 sequences have hits to the database of Helianthus and Lactuca species, respectively. MISA software identified a large number of putative EST-SSRs, allowing 1,788 primer pairs to be designed from the de novo transcriptome sequence and a further 363 from archival EST sequence. Among 100 primer pairs randomly chosen, 81 markers have amplicons and 20 are polymorphic for genotypes analysis in Chrysanthemum. The results showed that most (but not all) of the assays were transferable across species and that they exposed a significant amount of allelic diversity.

Conclusions/Significance

SSR markers acquired by transcriptome sequencing are potentially useful for marker-assisted breeding and genetic analysis in the genus Chrysanthemum and its related genera.

Non-cell-autonomous postmortem lignification of tracheary elements in Zinnia elegans

Postmortem lignification of xylem tracheary elements (TEs) has been debated for decades. Here, we provide evidence in Zinnia elegans TE cell cultures, using pharmacological inhibitors and in intact Z. elegans plants using Fourier transform infrared microspectroscopy, that TE lignification occurs postmortem (i.e., after TE programmed cell death). In situ RT-PCR verified expression of the lignin monomer biosynthetic cinnamoyl CoA reductase and cinnamyl alcohol dehydrogenase in not only the lignifying TEs but also in the unlignified non-TE cells of Z. elegans TE cell cultures and in living, parenchymatic xylem cells that surround TEs in stems. These cells were also shown to have the capacity to synthesize and transport lignin monomers and reactive oxygen species to the cell walls of dead TEs. Differential gene expression analysis in Z. elegans TE cell cultures and concomitant functional analysis in Arabidopsis thaliana resulted in identification of several genes that were expressed in the non-TE cells and that affected lignin chemistry on the basis of pyrolysis-gas chromatography/mass spectrometry analysis. These data suggest that living, parenchymatic xylem cells contribute to TE lignification in a non-cell-autonomous manner, thus enabling the postmortem lignification of TEs.

Transcriptome analysis of Gerbera hybrida ray florets: putative genes associated with gibberellin metabolism and signal transduction

In this study, the transcriptome of the Gerbera hybrida ray floret was constructed using a high-throughput Illumina sequencing platform. All 47,104 UniGenes with an average length of 845 nt and an N50 equaling 1321 nt were generated from 72,688,546 total primary reads after filtering and assembly. A total of 36,693 transcripts were annotated by comparison with non-redundant National Center for Biotechnology Information (NCBI) protein (Nr), non-redundant NCBI nucleotide (Nt), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases after removing exogenous contaminated sequences. The majority of the genes that are associated with gibberellin metabolism and signal transduction were identified. The targets for signal transduction of other plant hormones were also enumerated. Our study provides a systematic overview of the hormone signal transduction genes that are involved in ray floral development in Asteraceae and should facilitate further understanding of the crucial roles of phytohormones in plant growth.

Virus-induced gene silencing for Asteraceae--a reverse genetics approach for functional genomics in Gerbera hybrida

Virus-induced gene silencing (VIGS) is a natural defence mechanism in plants which leads to sequence-specific degradation of viral RNA. For identifying gene functions, Tobacco rattle virus (TRV)-based VIGS has been applied for silencing of endogenous genes in many plant species. Gerbera hybrida (Asteraceae) has emerged as a novel model for studies in flower development and secondary metabolism. For this highly heterozygous species, functional studies have been conducted through reverse genetic methods by producing stable transgenic lines, which, however, is labour-intensive and time-consuming. For the development of TRV-based VIGS system for gerbera, and for the first time for an Asteraceaeous species, we screened several gerbera cultivars and optimized the agroinfiltration methods for efficient silencing. Gene fragments for gerbera phytoene desaturase (GPDS) and Mg-chelatase subunits (GChl-H and GChl-I), expressed from a TRV vector, induced silencing phenotypes in leaves, scapes, and involucral bracts indicating their feasibility as markers for green tissues. In addition, robust silencing symptoms were achieved in gerbera floral tissues by silencing the anthocyanin pathway gene for chalcone synthase (GCHS1) and a gerbera B-type MADS-box gene globosa (GGLO1), confirming the phenotypes previously observed in stable transgenic lines. Unexpectedly, photobleaching induced by GPDS and GChl-H or GChl-I silencing, or by the herbicide norflurazon, resulted in silencing of the polyketide synthase gene G2PS1, which has no apparent connections to carotenoid or chlorophyll biosynthesis. We have shown feasibility of VIGS for functional studies in gerbera, but our results also show that selection of the marker gene for silencing must be critically evaluated.

Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Carefully analyzed expression profiles can serve as a valuable reference for deciphering gene functions. We exploited the potential of whole genome microarrays to measure the spatial and temporal expression profiles of rice genes in 19 stages of vegetative and reproductive development. We could verify expression of 22,980 genes in at least one of the tissues. Differential expression analysis with respect to five vegetative tissues and preceding stages of development revealed reproductive stage-preferential/-specific genes. By using subtractive logic, we identified 354 and 456 genes expressing specifically during panicle and seed development, respectively. The metabolic/hormonal pathways and transcription factor families playing key role in reproductive development were elucidated after overlaying the expression data on the public databases and manually curated list of transcription factors, respectively. During floral meristem differentiation (P1) and male meiosis (P3), the genes involved in jasmonic acid and phenylpropanoid biosynthesis were significantly upregulated. P6 stage of panicle, containing mature gametophytes, exhibited enrichment of transcripts involved in homogalacturonon degradation. Genes regulating auxin biosynthesis were induced during early seed development. We validated the stage-specificity of regulatory regions of three panicle-specific genes, OsAGO3, OsSub42, and RTS, and an early seed-specific gene, XYH, in transgenic rice. The data generated here provides a snapshot of the underlying complexity of the gene networks regulating rice reproductive development.

Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Carefully analyzed expression profiles can serve as a valuable reference for deciphering gene functions. We exploited the potential of whole genome microarrays to measure the spatial and temporal expression profiles of rice genes in 19 stages of vegetative and reproductive development. We could verify expression of 22,980 genes in at least one of the tissues. Differential expression analysis with respect to five vegetative tissues and preceding stages of development revealed reproductive stage-preferential/-specific genes. By using subtractive logic, we identified 354 and 456 genes expressing specifically during panicle and seed development, respectively. The metabolic/hormonal pathways and transcription factor families playing key role in reproductive development were elucidated after overlaying the expression data on the public databases and manually curated list of transcription factors, respectively. During floral meristem differentiation (P1) and male meiosis (P3), the genes involved in jasmonic acid and phenylpropanoid biosynthesis were significantly upregulated. P6 stage of panicle, containing mature gametophytes, exhibited enrichment of transcripts involved in homogalacturonon degradation. Genes regulating auxin biosynthesis were induced during early seed development. We validated the stage-specificity of regulatory regions of three panicle-specific genes, OsAGO3, OsSub42, and RTS, and an early seed-specific gene, XYH, in transgenic rice. The data generated here provides a snapshot of the underlying complexity of the gene networks regulating rice reproductive development.

PASSIOMA: Exploring Expressed Sequence Tags during Flower Development in Passiflora spp

The genus Passiflora provides a remarkable example of floral complexity and diversity. The extreme variation of Passiflora flower morphologies allowed a wide range of interactions with pollinators to evolve. We used the analysis of expressed sequence tags (ESTs) as an approach for the characterization of genes expressed during Passiflora reproductive development. Analyzing the Passiflora floral EST database (named PASSIOMA), we found sequences showing significant sequence similarity to genes known to be involved in reproductive development such as MADS-box genes. Some of these sequences were studied using RT-PCR and in situ hybridization confirming their expression during Passiflora flower development. The detection of these novel sequences can contribute to the development of EST-based markers for important agronomic traits as well as to the establishment of genomic tools to study the naturally occurring floral diversity among Passiflora species.

Os11Gsk gene from a wild rice, Oryza rufipogon improves yield in rice

Chromosomal segments from wild rice species Oryza rufipogon, introgressed into an elite indica rice restorer line (KMR3) using molecular markers, resulted in significant increase in yield. Here we report the transcriptome analysis of flag leaves and fully emerged young panicles of one of the high yielding introgression lines IL50-7 in comparison to KMR3. A 66-fold upregulated gene Os11Gsk, which showed no transcript in KMR3 was highly expressed in O. rufipogon and IL50-7. A 5-kb genomic region including Os11Gsk and its flanking regions could be PCR amplified only from IL50-7, O. rufipogon, japonica varieties of rice-Nipponbare and Kitaake but not from the indica varieties, KMR3 and Taichung Native-1. Three sister lines of IL50-7 yielding higher than KMR3 showed presence of Os11Gsk, whereas the gene was absent in three other ILs from the same cross having lower yield than KMR3, indicating an association of the presence of Os11Gsk with high yield. Southern analysis showed additional bands in the genomic DNA of O. rufipogon and IL50-7 with Os11Gsk probe. Genomic sequence analysis of ten highly co-expressed differentially regulated genes revealed that two upregulated genes in IL50-7 were derived from O. rufipogon and most of the downregulated genes were either from KMR3 or common to KMR3, IL50-7, and O. rufipogon. Thus, we show that Os11Gsk is a wild rice-derived gene introduced in KMR3 background and increases yield either by regulating expression of functional genes sharing homology with it or by causing epigenetic modifications in the introgression line.

Over-expression of the Gerbera hybrida At-SOC1-like1 gene Gh-SOC1 leads to floral organ identity deterioration

BACKGROUND AND AIMS

The family of MADS box genes is involved in a number of processes besides controlling floral development. In addition to supplying homeotic functions defined by the ABC model, they influence flowering time and transformation of vegetative meristem into inflorescence meristem, and have functions in roots and leaves. Three Gerbera hybrida At-SOC1-like genes (Gh-SOC1-Gh-SOC3) were identified among gerbera expressed sequence tags.

METHODS

Evolutionary relationships between SOC1-like genes from gerbera and other plants were studied by phylogenetic analysis. The function of the gerbera gene Gh-SOC1 in gerbera floral development was studied using expression analysis, protein-protein interaction assays and reverse genetics. Transgenic gerbera lines over-expressing or downregulated for Gh-SOC1 were obtained using Agrobacterium transformation and investigated for their floral phenotype.

KEY RESULTS

Phylogenetic analysis revealed that the closest paralogues of At-SOC1 are Gh-SOC2 and Gh-SOC3. Gh-SOC1 is a more distantly related paralogue, grouping together with a number of other At-SOC1 paralogues from arabidopsis and other plant species. Gh-SOC1 is inflorescence abundant and no expression was seen in vegetative parts of the plant. Ectopic expression of Gh-SOC1 did not promote flowering, but disturbed the development of floral organs. The epidermal cells of ray flower petals appeared shorter and their shape was altered. The colour of ray flower petals differed from that of the wild-type petals by being darker red on the adaxial side and greenish on the abaxial surface. Several protein-protein interactions with other gerbera MADS domain proteins were identified.

CONCLUSIONS

The At-SOC1 paralogue in gerbera shows a floral abundant expression pattern. A late petal expression might indicate a role in the final stages of flower development. Over-expression of Gh-SOC1 led to partial loss of floral identity, but did not affect flowering time. Lines where Gh-SOC1 was downregulated did not show a phenotype. Several gerbera MADS domain proteins interacted with Gh-SOC1.

Large scale interaction analysis predicts that the Gerbera hybrida floral E function is provided both by general and specialized proteins

BACKGROUND

The ornamental plant Gerbera hybrida bears complex inflorescences with morphologically distinct floral morphs that are specific to the sunflower family Asteraceae. We have previously characterized several MADS box genes that regulate floral development in Gerbera. To study further their behavior in higher order complex formation according to the quartet model, we performed yeast two- and three-hybrid analysis with fourteen Gerbera MADS domain proteins to analyze their protein-protein interaction potential.

RESULTS

The exhaustive pairwise interaction analysis showed significant differences in the interaction capacity of different Gerbera MADS domain proteins compared to other model plants. Of particular interest in these assays was the behavior of SEP-like proteins, known as GRCDs in Gerbera. The previously described GRCD1 and GRCD2 proteins, which are specific regulators involved in stamen and carpel development, respectively, showed very limited pairwise interactions, whereas the related GRCD4 and GRCD5 factors displayed hub-like positions in the interaction map. We propose GRCD4 and GRCD5 to provide a redundant and general E function in Gerbera, comparable to the SEP proteins in Arabidopsis. Based on the pairwise interaction data, combinations of MADS domain proteins were further subjected to yeast three-hybrid assays. Gerbera B function proteins showed active behavior in ternary complexes. All Gerbera SEP-like proteins with the exception of GRCD1 were excellent partners for B function proteins, further implicating the unique role of GRCD1 as a whorl- and flower-type specific C function partner.

CONCLUSIONS

Gerbera MADS domain proteins exhibit both conserved and derived behavior in higher order protein complex formation. This protein-protein interaction data can be used to classify and compare Gerbera MADS domain proteins to those of Arabidopsis and Petunia. Combined with our reverse genetic studies of Gerbera, these results reinforce the roles of different genes in the floral development of Gerbera. Building up the elaborate capitulum of Gerbera calls for modifications and added complexity in MADS domain protein behavior compared to the more simple flowers of, e.g., Arabidopsis.

Characterization of SQUAMOSA-like genes in Gerbera hybrida, including one involved in reproductive transition

BACKGROUND

The flowering process in plants proceeds through the induction of an inflorescence meristem triggered by several pathways. Many of the genes associated with both the flowering process and floral architecture encode transcription factors of the MADS domain family. Gerbera, a member of the sunflower family, Asteraceae, bears compressed inflorescence heads (capitula) with three different flower types characterized by differences in both sexuality and floral symmetry. To understand how such a complex inflorescence structure is achieved at the molecular level, we have characterized the array of Gerbera MADS box genes. The high number of SQUAMOSA-like genes in Gerbera compared to other model species raised the question as to whether they may relate to Gerbera's complex inflorescence structure and whether or not a homeotic A function is present.

RESULTS

In this paper we describe six Gerbera genes related to the SQUAMOSA/APETALA1/FRUITFULL genes of snapdragon and Arabidopsis. Based on phylogenetic analysis of the entire gene lineage, our data indicates that GSQUA1 and GSQUA3 are members of the SQUA/AP1 clade, while GSQUA2, GSQUA4, GSQUA5 and GSQUA6 are co-orthologs of the Arabidopsis FUL gene. GSQUA1/GSQUA3 and GSQUA4/GSQUA5/GSQUA6, respectively, represent several gene duplication events unknown in the model systems that may be specific to either Gerbera or Asteraceae. GSQUA genes showed specific expression profiles. GSQUA1, GSQUA2, and GSQUA5 were inflorescence abundant, while GSQUA3, GSQUA4, and GSQUA6 expression was also detected in vegetative organs. Overexpression of GSQUA2 in Gerbera led to accelerated flowering, dwarfism and vegetative abnormalities, all new and specific phenomena observed in transgenic Gerbera plants with modified MADS box gene expression.

CONCLUSIONS

Based on expression patterns, none of the Gerbera SQUA-like genes are likely to control flower organ identity in the sense of the floral A function. However, our data shows that the FUL-like gene GSQUA2 plays a vital role in meristem transition. The roles of other GSQUA-genes in Gerbera floral development are intriguing, but require still further study.

Analysis of floral transcription factors from Lycoris longituba

Transcription factors (TFs) are proteins that bind to specific promoter regions of their target genes and regulate gene transcription. Many of these factors have been found to influence flowering. Lycoris longituba exhibits a great deal of diversity in flower color and flower form, making it a suitable model for the study of floral development. We have identified 338 putative TFs from more than thirty thousand ESTs sequenced from the floral tissue of L. longituba, and validated them using real-time RT-PCR. Fifty-one of the TFs were recognized as being potentially flower-specific, and the expression patterns of some of them during six flowering phases have been elucidated. Homolog annotation and phylogenetic analysis revealed that some TFs that belong to several TF families, such as MADS, MYB-related, NAC, and ABI3-VP1, were suggested to play important roles in the flowering process. Our dataset may be used to identify priority target TF genes for further study.

Within and between whorls: comparative transcriptional profiling of Aquilegia and Arabidopsis

Background

The genus Aquilegia is an emerging model system in plant evolutionary biology predominantly because of its wide variation in floral traits and associated floral ecology. The anatomy of the Aquilegia flower is also very distinct. There are two whorls of petaloid organs, the outer whorl of sepals and the second whorl of petals that form nectar spurs, as well as a recently evolved fifth whorl of staminodia inserted between stamens and carpels.

Methodology/Principal Findings

We designed an oligonucleotide microarray based on EST sequences from a mixed tissue, normalized cDNA library of an A. formosa x A. pubescens F2 population representing 17,246 unigenes. We then used this array to analyze floral gene expression in late pre-anthesis stage floral organs from a natural A. formosa population. In particular, we tested for gene expression patterns specific to each floral whorl and to combinations of whorls that correspond to traditional and modified ABC model groupings. Similar analyses were performed on gene expression data of Arabidopsis thaliana whorls previously obtained using the Ath1 gene chips (data available through The Arabidopsis Information Resource).

Conclusions/Significance

Our comparative gene expression analyses suggest that 1) petaloid sepals and petals of A. formosa share gene expression patterns more than either have organ-specific patterns, 2) petals of A. formosa and A. thaliana may be independently derived, 3) staminodia express B and C genes similar to stamens but the staminodium genetic program has also converged on aspects of the carpel program and 4) staminodia have unique up-regulation of regulatory genes and genes that have been implicated with defense against microbial infection and herbivory. Our study also highlights the value of comparative gene expression profiling and the Aquilegia microarray in particular for the study of floral evolution and ecology.

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.

Comparative analyses of the proteomes of leaves and flowers at various stages of development reveal organ-specific functional differentiation of proteins in soybean

The functional differentiation of protein networks in individual organs and tissues of soybean at various developmental stages was investigated by proteomic approach. Protein extraction by Mg/NP-40 buffer followed by alkaline phenol-based method was optimized for proteomic analysis. Proteome analyses of leaves at various developmental stages showed 26 differentially expressed proteins, wherein proteins in translocon at the outer/inner envelope membrane of chloroplast protein-transport machineries increased significantly at the first trifoliate. Immunoblot analysis showed chaperonin-60 expressed abundantly in young leaves, whereas HSP 70 and ATP-synthase beta were constitutively expressed in all tissues. The net photosynthesis rate and chlorophyll content showed an age-dependent correlation in leaves. These results suggest that proteins involved in carbon assimilation, folding and assembly, and energy may work synchronously and show a linear correlation to photosynthesis at developmental stages of leaves. Comparison of flower bud and flower proteome reveals 29 differentially expressed proteins, wherein proteins involved in mitochondrial protein transport and assembly, secondary metabolism, and pollen-tube growth were up-regulated during flower development. Together, these results suggest that during developmental stages, each type of tissue is associated with a specific group of proteins; wherein proteins involved in energy, sugar metabolism, and folding, assembly, and destination may play pivotal roles in the maturation process of each organ or tissue.

Identification of flowering genes in strawberry, a perennial SD plant

BACKGROUND

We are studying the regulation of flowering in perennial plants by using diploid wild strawberry (Fragaria vesca L.) as a model. Wild strawberry is a facultative short-day plant with an obligatory short-day requirement at temperatures above 15 degrees C. At lower temperatures, however, flowering induction occurs irrespective of photoperiod. In addition to short-day genotypes, everbearing forms of wild strawberry are known. In 'Baron Solemacher' recessive alleles of an unknown repressor, SEASONAL FLOWERING LOCUS (SFL), are responsible for continuous flowering habit. Although flower induction has a central effect on the cropping potential, the molecular control of flowering in strawberries has not been studied and the genetic flowering pathways are still poorly understood. The comparison of everbearing and short-day genotypes of wild strawberry could facilitate our understanding of fundamental molecular mechanisms regulating perennial growth cycle in plants.

RESULTS

We have searched homologs for 118 Arabidopsis flowering time genes from Fragaria by EST sequencing and bioinformatics analysis and identified 66 gene homologs that by sequence similarity, putatively correspond to genes of all known genetic flowering pathways. The expression analysis of 25 selected genes representing various flowering pathways did not reveal large differences between the everbearing and the short-day genotypes. However, putative floral identity and floral integrator genes AP1 and LFY were co-regulated during early floral development. AP1 mRNA was specifically accumulating in the shoot apices of the everbearing genotype, indicating its usability as a marker for floral initiation. Moreover, we showed that flowering induction in everbearing 'Baron Solemacher' and 'Hawaii-4' was inhibited by short-day and low temperature, in contrast to short-day genotypes.

CONCLUSION

We have shown that many central genetic components of the flowering pathways in Arabidopsis can be identified from strawberry. However, novel regulatory mechanisms exist, like SFL that functions as a switch between short-day/low temperature and long-day/high temperature flowering responses between the short-day genotype and the everbearing 'Baron Solemacher'. The identification of putative flowering gene homologs and AP1 as potential marker gene for floral initiation will strongly facilitate the exploration of strawberry flowering pathways.

Expressed sequence tags from the flower pathogen Claviceps purpurea

SUMMARY The ascomycete Claviceps purpurea (ergot) is a biotrophic flower pathogen of rye and other grasses. The deleterious toxic effects of infected rye seeds on humans and grazing animals have been known since the Middle Ages. To gain further insight into the molecular basis of this disease, we generated about 10 000 expressed sequence tags (ESTs)-about 25% originating from axenic fungal culture and about 75% from tissues collected 6-20 days after infection of rye spikes. The pattern of axenic vs. in planta gene expression was compared. About 200 putative plant genes were identified within the in planta library. A high percentage of these were predicted to function in plant defence against the ergot fungus and other pathogens, for example pathogenesis-related proteins. Potential fungal pathogenicity and virulence genes were found via comparison with the pathogen-host interaction database (PHI-base; http://www.phi-base.org) and with genes known to be highly expressed in the haustoria of the bean rust fungus. Comparative analysis of Claviceps and two other fungal flower pathogens (necrotrophic Fusarium graminearum and biotrophic Ustilago maydis) highlighted similarities and differences in their lifestyles, for example all three fungi have signalling components and cell wall-degrading enzymes in their arsenal. In summary, the analysis of axenic and in planta ESTs yielded a collection of candidate genes to be evaluated for functional roles in this plant-microbe interaction.

The MYB305 transcription factor regulates expression of nectarin genes in the ornamental tobacco floral nectary

We have isolated and characterized the cDNA encoding the ornamental tobacco (Nicotiana langsdorffii X N. sanderae) homolog of the antirrhinum (Antirrhinum majus) MYB305. This transcription factor was robustly expressed at Stage 12 of nectary development but was only weakly expressed in the earlier Stage 6 nectaries. The ornamental tobacco MYB305 contains a conserved R2R3 MYB DNA binding domain with 76 amino acids in the activation domain. A green fluorescent protein-MYB305 fusion localized to nucleus of tobacco protoplasts and yeast one-hybrid assays demonstrated that it functions as a transcription activator. A conserved 23-amino acid C-terminal domain is required to activate gene expression. The coding region of the myb305 cDNA was expressed in Escherichia coli as a glutathione S-transferase fusion protein and was purified to homogeneity. This protein shows binding to two consensus MYB binding sites on the ornamental tobacco Nectarin I (nec1) promoter as well as to the single site located on the Nectarin V (nec5) promoter. Deletions of either of the binding sites from the nec1 promoter significantly reduced expression in nectary tissues. Temporally, MYB305 expression precedes that of nec1 and nec5, as would be expected if the MYB305 factor regulates expression of the nec1 and nec5 genes. Ectopic expression of MYB305 in foliage was able to induce expression of both nec1 and nec5, as well as two flavonoid biosynthetic genes in the foliage. Finally, RNA interference knockdown of MYB305 resulted in reduced expression of both nectarins and flavonoid biosynthetic genes. We conclude that expression of MYB305 regulates expression of the major nectarin genes in the floral nectary.

Identification of transcription factors predominantly expressed in soybean flowers and characterization of GmSEP1 encoding a SEPALLATA1-like protein

By microarray analysis and real-time RT-PCR, we identified 28 soybean flower-enriched transcription factors such as MADS-box proteins, zinc finger proteins, and MYB proteins. Among them, one MADS-box protein GmSEP1 was chosen for further analysis. GmSEP1 contains 8 exons and 7 introns, showing similar exon-intron structure with Arabidopsis SEP genes. Phylogenetic analysis also suggested that GmSEP1 belonged to AGL2/SEP subfamily and likely a soybean orthologue of Arabidopsis SEP1. Subcellular localization assay suggested that GmSEP1 was localized in nucleus. Semi-quantitative RT-PCR showed that GmSEP1 was predominantly expressed in reproductive organs including petal, sepal, stamen, carpel and seed. By real-time RT-PCR, we found that GmSEP1 was expressed with a low level in stamens of soybean mutant NJS-10Hfs in which some stamens are converted to petal-like structure. Further we found that expression of GmSEP1 was altered during the course of seed development and more accumulated in outer and inner integuments, epidermis and endothelium at the globular stage, heart stage and cotyledon stage and in seed coat parenchyma at the early maturation stage of seed development. Taken together, GmSEP1 might play important role in soybean reproductive development, in particular, in petal and seed coat development.

Genome-wide analysis of gene expression profiles during ear development of maize

In order to unravel the molecular mechanism of maize ear development, a microarray containing approximately 56,000 probes was used to monitor the gene expression profiles of ears at four developmental stages. The results showed that 2,794 genes, accounting for 5.0% of the total probes, changed significantly during ear development. Among the 2,794 genes, 1,844 genes differentially expressed during the spikelet differentiation phase, 836 genes during the floret primordium differentiation phase and 645 genes during the floret organ differentiation phase. Hierarchical clustering revealed that the differentially expressed genes had 9 major expression patterns. Based on Mips Functional Catalogue, 684 differentially expressed genes were grouped into at least one functional category, including metabolism (30.4%), protein related function (29.2%), biogenesis of cellular components (15.4%) and transcription (13.7%). The analysis revealed that the auxin signaling pathway play an important role in ear development. Moreover, regulation of some transcription factors may play a key role during ear development. RT-PCR and in situ hybridization for some selected genes validated our microarray data and supplied additional information on ear developmental processes.

Non-targeted analysis of spatial metabolite composition in strawberry (Fragariaxananassa) flowers

Formation of flower organs and the subsequent pollination process require a coordinated spatial and temporal regulation of particular metabolic pathways. In this study a comparison has been made between the metabolite composition of individual flower organs of strawberry (Fragariaxananassa) including the petal, sepal, stamen, pistil and the receptacle that gives rise to the strawberry fruit. Non-targeted metabolomics analysis of the semi-polar secondary metabolites by the use of UPLC-qTOF-MS was utilized in order to localize metabolites belonging to various chemical classes (e.g. ellagitannins, proanthocyanidins, flavonols, terpenoids, and spermidine derivatives) to the different flower organs. The vast majority of the tentatively identified metabolites were ellagitannins that accumulated in all five parts of the flower. Several metabolite classes were detected predominantly in certain flower organs, as for example spermidine derivatives were present uniquely in the stamen and pistil, and the proanthocyanidins were almost exclusively detected in the receptacle and sepals. The latter organ was also rich in terpenoids (i.e. triterpenoid and sesquiterpenoid derivatives) whereas phenolic acids and flavonols were the predominant classes of compounds detected in the petals. Furthermore, we observed extensive variation in the accumulation of metabolites from the same class in a single organ, particularly in the case of ellagitannins, and the flavonols quercetin, kaempferol and isorhamnetin. These results allude to spatially-restricted production of secondary metabolite classes and specialized derivatives in flowers that take part in implementing the unique program of individual organs in the floral life cycle.

Multiple paleopolyploidizations during the evolution of the Compositae reveal parallel patterns of duplicate gene retention after millions of years

Of the approximately 250,000 species of flowering plants, nearly one in ten are members of the Compositae (Asteraceae), a diverse family found in almost every habitat on all continents except Antarctica. With an origin in the mid Eocene, the Compositae is also a relatively young family with remarkable diversifications during the last 40 My. Previous cytologic and systematic investigations suggested that paleopolyploidy may have occurred in at least one Compositae lineage, but a recent analysis of genomic data was equivocal. We tested for evidence of paleopolyploidy in the evolutionary history of the family using recently available expressed sequence tag (EST) data from the Compositae Genome Project. Combined with data available on GenBank, we analyzed nearly 1 million ESTs from 18 species representing seven genera and four tribes. Our analyses revealed at least three ancient whole-genome duplications in the Compositae-a paleopolyploidization shared by all analyzed taxa and placed near the origin of the family just prior to the rapid radiation of its tribes and independent genome duplications near the base of the tribes Mutisieae and Heliantheae. These results are consistent with previous research implicating paleopolyploidy in the evolution and diversification of the Heliantheae. Further, we observed parallel retention of duplicate genes from the basal Compositae genome duplication across all tribes, despite divergence times of 33-38 My among these lineages. This pattern of retention was also repeated for the paleologs from the Heliantheae duplication. Intriguingly, the categories of genes retained in duplicate were substantially different from those in Arabidopsis. In particular, we found that genes annotated to structural components or cellular organization Gene Ontology categories were significantly enriched among paleologs, whereas genes associated with transcription and other regulatory functions were significantly underrepresented. Our results suggest that paleopolyploidy can yield strikingly consistent signatures of gene retention in plant genomes despite extensive lineage radiations and recurrent genome duplications but that these patterns vary substantially among higher taxonomic categories.

Persea americana (avocado): bringing ancient flowers to fruit in the genomics era

The avocado (Persea americana) is a major crop commodity worldwide. Moreover, avocado, a paleopolyploid, is an evolutionary "outpost" among flowering plants, representing a basal lineage (the magnoliid clade) near the origin of the flowering plants themselves. Following centuries of selective breeding, avocado germplasm has been characterized at the level of microsatellite and RFLP markers. Nonetheless, little is known beyond these general diversity estimates, and much work remains to be done to develop avocado as a major subtropical-zone crop. Among the goals of avocado improvement are to develop varieties with fruit that will "store" better on the tree, show uniform ripening and have better post-harvest storage. Avocado transcriptome sequencing, genome mapping and partial genomic sequencing will represent a major step toward the goal of sequencing the entire avocado genome, which is expected to aid in improving avocado varieties and production, as well as understanding the evolution of flowers from non-flowering seed plants (gymnosperms). Additionally, continued evolutionary and other comparative studies of flower and fruit development in different avocado strains can be accomplished at the gene expression level, including in comparison with avocado relatives, and these should provide important insights into the genetic regulation of fruit development in basal angiosperms.

Identification of target genes for a MYB-type anthocyanin regulator in Gerbera hybrida

Genetic modification of the flavonoid pathway has been used to produce novel colours and colour patterns in ornamental plants as well as to modify the nutritional and pharmaceutical properties of food crops. It has been suggested that co-ordinate control of multiple steps of the pathway with the help of regulatory genes would lead to a more predictable control of metabolic flux. Regulation of anthocyanin biosynthesis has been studied in a common ornamental plant, Gerbera hybrida (Asteraceae). An R2R3-type MYB factor, GMYB10, shares high sequence similarity and is phylogenetically grouped together with previously characterized regulators of anthocyanin pigmentation. Ectopic expression of GMYB10 leads to strongly enhanced accumulation of anthocyanin pigments as well as to an altered pigmentation pattern in transgenic gerbera plants. Anthocyanin analysis indicates that GMYB10 specifically induces cyanidin biosynthesis in undifferentiated callus and in vegetative tissues. Furthermore, in floral tissues enhanced pelargonidin production is detected. Microarray analysis using the gerbera 9K cDNA array revealed a highly predicted set of putative target genes for GMYB10 including new gene family members of both early and late biosynthetic genes of the flavonoid pathway. However, completely new candidate targets, such as a serine carboxypeptidase-like gene as well, as two new MYB domain factors, GMYB11 and GMYB12, whose exact function in phenylpropanoid biosynthesis is not clear yet, were also identified.

Transcriptional analysis of petal organogenesis in Gerbera hybrida

Understanding of the molecular interplay, which determines early steps of flower formation has grown considerably during last years. In contrast, genetic actions responsible for how flower organs acquire their size and shape at later phases of organogenesis are still poorly understood. We have exploited the large and anatomically simple Gerbera (Gerbera hybrida var. Terra regina) ray flower petals to describe transcriptional changes during organogenesis. Gerbera 9 K cDNA microarray was utilized to profile gene expression at six different developmental stages of petal organogenesis, at the earliest stage expansion of petals is starting and at the latest stage petals have reached their final size and shape. Genes potentially participating in petal opening were identified based on the similarity in expression with a known marker gene. Our results showed characteristic sets of genes expressed during the cell division and cell expansion phases of petal development. Interestingly, there was a transition stage during which neither cell division nor cell expansion marker genes were abundantly expressed. Moreover, constitutive expression of late petal specific genes indicates that they participate in petal organogenesis throughout the development and they are not involved in stage specific switch points.