Genome-wide expression profiling and identification of gene activities during early flower development in Arabidopsis

Developmentally dependent reprogramming of the Arabidopsis floral transcriptome under sufficient and limited water availability

BACKGROUND

Environmental stresses negatively impact reproductive development and yield. Drought stress, in particular, has been examined during Arabidopsis reproductive development at morphological and transcriptomic levels. However, drought-responsive transcriptomic changes at different points in reproductive development remain unclear. Additionally, an investigation of the entire transcriptome at various stages during flower development is of great interest.

RESULTS

Here, we treat Arabidopsis plants with well-watered and moderately and severely limiting water amounts when the first flowers reach maturity and generate RNA-seq datasets for early, middle, and late phases during flower development at 5, 6, and 7 days following treatment. Under different drought conditions, flowers in different developmental phases display differential sets of drought-responsive genes (DTGs), including those that are enriched in different GO functional categories, such as transcriptional regulation and response to stresses (early phase), lipid storage (middle phase), and pollen and seed development and metabolic processes (late phase). Some gene families have different members induced at different floral phases, suggesting that similar biochemical functions are carried out by distinct members. Developmentally-regulated genes (DVGs) with differential expression among the three floral phases belong to GO terms that are similar between water conditions, such as development and reproduction, metabolism and transport, and signaling and stress response. However, for different water conditions, such similar GO terms correspond to either distinct gene families or different members of a gene family, suggesting that drought affects the expression of distinct families or family members during reproductive development. A further comparison among transcriptomes of tissues collected on different days after treatment identifies differential gene expression, suggesting age-related genes (ARGs) might reflect the changes in the overall plant physiology in addition to drought response and development.

CONCLUSION

Together, our study provides new insights into global transcriptome reprogramming and candidate genes for drought response, flower development, aging and coordination among these complex biological processes.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Pathways to sex determination in plants: how many roads lead to Rome?

The presence of thousands of independent origins of dioecy in angiosperms provides a unique opportunity to address the parallel evolution of the molecular pathways underlying unisexual flowers. Recent progress towards identifying sex determination genes has identified hormone response pathways, mainly associated with cytokinin and ethylene response pathways, as having been recruited multiple times independently to control unisexuality. Moreover, transcriptomics has begun to identify commonalities among intermediate sections of signal transduction pathways. These recent advances set the stage for development of a comparative evolutionary development research program to identify the shared and unique aspects of the genetic pathways of unisexual flower development in angiosperms.

The Quest for Molecular Regulation Underlying Unisexual Flower Development

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled "Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber" published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.

Next-generation transcriptome analysis in transgenic birch overexpressing and suppressing APETALA1 sheds lights in reproduction development and diterpenoid biosynthesis

Overexpression of BpAP1 could cause early flowering in birch. BpAP1 affected the expression of many flowering-related unigenes and diterpenoid biosynthesis in transgenic birch, and BpPI was a putative target gene of BpAP1. APETALA1 (AP1) is an MADS-box transcription factor that is involved in the flowering process in plants and has been a focus of genetic studies examining flower development. Here, we carried out transcriptome analysis of birch (Betula platyphylla Suk.), including BpAP1 overexpression lines, BpAP1 suppression lines, and non-transgenic line (NT). Compared with NT, we detected 8302 and 7813 differentially expressed unigenes in 35S::BpAP1 and 35S::BpAP1RNAi transgenic lines, respectively. Overexpression and suppression of BpAP1 in birch affected diterpenoid biosynthesis and altered expression of many flowering-related unigenes. Moreover, combining information from the RNA-seq database and the birch genome, we predicted downstream target genes of BpAP1. Among the 166 putative target genes of BpAP1, there was a positive correlation between BpAP1 and BpPI. These results provide references for further examining the relationship between BpAP1 and its target genes, and reveal that BpAP1 functions as a transcription regulator in birch.

Functional characterization of drought-responsive modules and genes in Oryza sativa: a network-based approach

Drought is one of the major environmental stress conditions affecting the yield of rice across the globe. Unraveling the functional roles of the drought-responsive genes and their underlying molecular mechanisms will provide important leads to improve the yield of rice. Co-expression relationships derived from condition-dependent gene expression data is an effective way to identify the functional associations between genes that are part of the same biological process and may be under similar transcriptional control. For this purpose, vast amount of freely available transcriptomic data may be used. In this study, we consider gene expression data for different tissues and developmental stages in response to drought stress. We analyze the network of co-expressed genes to identify drought-responsive genes modules in a tissue and stage-specific manner based on differential expression and gene enrichment analysis. Taking cues from the systems-level behavior of these modules, we propose two approaches to identify clusters of tightly co-expressed/co-regulated genes. Using graph-centrality measures and differential gene expression, we identify biologically informative genes that lack any functional annotation. We show that using orthologous information from other plant species, the conserved co-expression patterns of the uncharacterized genes can be identified. Presence of a conserved neighborhood enables us to extrapolate functional annotation. Alternatively, we show that single 'guide-gene' approach can help in understanding tissue-specific transcriptional regulation of uncharacterized genes. Finally, we confirm the predicted roles of uncharacterized genes by the analysis of conserved cis-elements and explain the possible roles of these genes toward drought tolerance.

Patterns of gene expression during Arabidopsis flower development from the time of initiation to maturation

Background

The formation of flowers is one of the main model systems to elucidate the molecular mechanisms that control developmental processes in plants. Although several studies have explored gene expression during flower development in the model plant Arabidopsis thaliana on a genome-wide scale, a continuous series of expression data from the earliest floral stages until maturation has been lacking. Here, we used a floral induction system to close this information gap and to generate a reference dataset for stage-specific gene expression during flower formation.

Results

Using a floral induction system, we collected floral buds at 14 different stages from the time of initiation until maturation. Using whole-genome microarray analysis, we identified 7,405 genes that exhibit rapid expression changes during flower development. These genes comprise many known floral regulators and we found that the expression profiles for these regulators match their known expression patterns, thus validating the dataset. We analyzed groups of co-expressed genes for over-represented cellular and developmental functions through Gene Ontology analysis and found that they could be assigned specific patterns of activities, which are in agreement with the progression of flower development. Furthermore, by mapping binding sites of floral organ identity factors onto our dataset, we were able to identify gene groups that are likely predominantly under control of these transcriptional regulators. We further found that the distribution of paralogs among groups of co-expressed genes varies considerably, with genes expressed predominantly at early and intermediate stages of flower development showing the highest proportion of such genes.

Conclusions

Our results highlight and describe the dynamic expression changes undergone by a large number of genes during flower development. They further provide a comprehensive reference dataset for temporal gene expression during flower formation and we demonstrate that it can be used to integrate data from other genomics approaches such as genome-wide localization studies of transcription factor binding sites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1699-6) contains supplementary material, which is available to authorized users.

Genetics of dioecy and causal sex chromosomes in plants

Dioecy (separate male and female individuals) ensures outcrossing and is more prevalent in animals than in plants. Although it is common in bryophytes and gymnosperms, only 5% of angiosperms are dioecious. In dioecious higher plants, flowers borne on male and female individuals are, respectively deficient in functional gynoecium and androecium. Dioecy is inherited via three sex chromosome systems: XX/XY, XX/X0 and WZ/ZZ, such that XX or WZ is female and XY, X0 or ZZ are males. The XX/XY system generates the rarer XX/X0 and WZ/ZZ systems. An autosome pair begets XY chromosomes. A recessive loss-of-androecium mutation (ana) creates X chromosome and a dominant gynoecium-suppressing (GYS) mutation creates Y chromosome. The ana/ANA and gys/GYS loci are in the sex-determining region (SDR) of the XY pair. Accumulation of inversions, deleterious mutations and repeat elements, especially transposons, in the SDR of Y suppresses recombination between X and Y in SDR, making Y labile and increasingly degenerate and heteromorphic from X. Continued recombination between X and Y in their pseudoautosomal region located at the ends of chromosomal arms allows survival of the degenerated Y and of the species. Dioecy is presumably a component of the evolutionary cycle for the origin of new species. Inbred hermaphrodite species assume dioecy. Later they suffer degenerate-Y-led population regression. Cross-hybridization between such extinguishing species and heterologous species, followed by genome duplication of segregants from hybrids, give rise to new species.

Identification of tapetum-specific genes by comparing global gene expression of four different male sterile lines in Brassica oleracea

The tapetum plays an important role in anther development by providing necessary enzymes and nutrients for pollen development. However, it is difficult to identify tapetum-specific genes on a large-scale because of the difficulty of separating tapetum cells from other anther tissues. Here, we reported the identification of tapetum-specific genes by comparing the gene expression patterns of four male sterile (MS) lines of Brassica oleracea. The abortive phenotypes of the four MS lines revealed different defects in tapetum and pollen development but normal anther wall development when observed by transmission electron microscopy. These tapetum displayed continuous defective characteristics throughout the anther developmental stages. The transcriptome from flower buds, covering all anther developmental stages, was analyzed and bioinformatics analyses exploring tapetum development-related genes were performed. We identified 1,005 genes differentially expressed in at least one of the MS lines and 104 were non-pollen expressed genes (NPGs). Most of the identified NPGs were tapetum-specific genes considering that anther walls were normally developed in all four MS lines. Among the 104 NPGs, 22 genes were previously reported as being involved in tapetum development. We further separated the expressed NPGs into different developmental stages based on the MS defects. The data obtained in this study are not only informative for research on tapetum development in B. oleracea, but are also useful for genetic pathway research in other related species.

Estimating the proportion of true null hypotheses when the statistics are discrete

MOTIVATION

In high-dimensional testing problems π0, the proportion of null hypotheses that are true is an important parameter. For discrete test statistics, the P values come from a discrete distribution with finite support and the null distribution may depend on an ancillary statistic such as a table margin that varies among the test statistics. Methods for estimating π0 developed for continuous test statistics, which depend on a uniform or identical null distribution of P values, may not perform well when applied to discrete testing problems.

RESULTS

This article introduces a number of π0 estimators, the regression and 'T' methods that perform well with discrete test statistics and also assesses how well methods developed for or adapted from continuous tests perform with discrete tests. We demonstrate the usefulness of these estimators in the analysis of high-throughput biological RNA-seq and single-nucleotide polymorphism data.

AVAILABILITY AND IMPLEMENTATION

implemented in R.

Analysis of the floral transcriptome of Tarenaya hassleriana (Cleomaceae), a member of the sister group to the Brassicaceae: towards understanding the base of morphological diversity in Brassicales

Background

Arabidopsis thaliana, a member of the Brassicaceae family is the dominant genetic model plant. However, while the flowers within the Brassicaceae members are rather uniform, mainly radially symmetrical, mostly white with fixed organ numbers, species within the Cleomaceae, the sister family to the Brassicaceae show a more variable floral morphology. We were interested in understanding the molecular basis for these morphological differences. To this end, the floral transcriptome of a hybrid Tarenaya hassleriana, a Cleomaceae with monosymmetric, bright purple flowers was sequenced, annotated and analyzed in respect to floral regulators.

Results

We obtained a comprehensive floral transcriptome with high depth and coverage close to saturation analyzed using rarefaction analysis a method well known in biodiversity studies. Gene expression was analyzed by calculating reads per kilobase gene model per million reads (RPKM) and for selected genes in silico expression data was corroborated by qRT-PCR analysis. Candidate transcription factors were identified based on differences in expression pattern between A. thaliana and T. hassleriana, which are likely key regulators of the T. hassleriana specific floral characters such as coloration and male sterility in the hybrid plant used. Analysis of lineage specific genes was carried out with members of the fabids and malvids.

Conclusions

The floral transcriptome of T. hassleriana provides insights into key pathways involved in the regulation of late anthocyanin biosynthesis, male fertility, flowering time and organ growth regulation which are unique traits compared the model organism A. thaliana. Analysis of lineage specific genes carried out with members of the fabids and malvids suggests an extensive gene birth rate in the lineage leading to core Brassicales while only few genes were potentially lost during core Brassicales evolution, which possibly reflects the result of the At-β whole genome duplication. Our analysis should facilitate further analyses into the molecular mechanisms of floral morphogenesis and pigmentation and the mechanisms underlying the rather diverse floral morphologies in the Cleomaceae.

On the role of AtDMC1, AtRAD51 and its paralogs during Arabidopsis meiosis

Meiotic recombination plays a critical role in achieving accurate chromosome segregation and increasing genetic diversity. Many studies, mostly in yeast, have provided important insights into the coordination and interplay between the proteins involved in the homologous recombination pathway, especially the recombinase RAD51 and the meiosis-specific DMC1. Here we summarize the current progresses on the function of both recombinases and the CX3 complex encoded by AtRAD51 paralogs, in the plant model species Arabidopsis thaliana. Similarities and differences respect to the function of these proteins in other organisms are also indicated.

Alternative splicing during Arabidopsis flower development results in constitutive and stage-regulated isoforms

Alternative splicing (AS) is a process in eukaryotic gene expression, in which the primary transcript of a multi-exon gene is spliced into two or more different mature transcripts, thereby increasing proteome diversity. AS is often regulated differentially between different tissues or developmental stages. Recent studies suggested that up to 60% of intron-containing genes in Arabidopsis thaliana undergo AS. Yet little is known about this complicated and important process during floral development. To investigate the preferential expression of different isoforms of individual alternatively spliced genes, we used high throughput RNA-Seq technology to explore the transcriptomes of three floral development stages of Arabidopsis thaliana and obtained information of various AS events. We identified approximately 24,000 genes that were expressed at one or more of these stages, and found that nearly 25% of multi-exon genes had two or more spliced variants. This is less frequent than the previously reported 40–60% for multiple organs and stages of A. thaliana, indicating that many genes expressed in floral development function with a single predominant isoform. On the other hand, 1716 isoforms were differentially expressed between the three stages, suggesting that AS might still play important roles in stage transition during floral development. Moreover, 337 novel transcribed regions were identified and most of them have a single exon. Taken together, our analyses provide a comprehensive survey of AS in floral development and facilitate further genomic and genetic studies.

The Arabidopsis RAD51 paralogs RAD51B, RAD51D and XRCC2 play partially redundant roles in somatic DNA repair and gene regulation

The eukaryotic RAD51 gene family has seven ancient paralogs conserved between plants and animals. Among these, RAD51, DMC1, RAD51C and XRCC3 are important for homologous recombination and/or DNA repair, whereas single mutants in RAD51B, RAD51D or XRCC2 show normal meiosis, and the lineages they represent diverged from each other evolutionarily later than the other four paralogs, suggesting possible functional redundancy. The function of Arabidopsis RAD51B, RAD51D and XRCC2 genes in mitotic DNA repair and meiosis was analyzed using molecular genetic, cytological and transcriptomic approaches. The relevant double and triple mutants displayed normal vegetative and reproductive growth. However, the triple mutant showed greater sensitivity than single or double mutants to DNA damage by bleomycin. RNA-Seq transcriptome analysis supported the idea that the triple mutant showed DNA damage similar to that caused by bleomycin. On bleomycin treatment, many genes were altered in the wild-type but not in the triple mutant, suggesting that the RAD51 paralogs have roles in the regulation of gene transcription, providing an explanation for the hypersensitive phenotype of the triple mutant to bleomycin. Our results provide strong evidence that Arabidopsis XRCC2, RAD51B and RAD51D have complex functions in somatic DNA repair and gene regulation, arguing for further studies of these ancient genes that have been maintained in both plants and animals during their long evolutionary history.

Analysis of Arabidopsis floral transcriptome: detection of new florally expressed genes and expansion of Brassicaceae-specific gene families

The flower is essential for sexual reproduction of flowering plants and has been extensively studied. However, it is still not clear how many genes are expressed in the flower. Here, we performed RNA-seq analysis as a highly sensitive approach to investigate the Arabidopsis floral transcriptome at three developmental stages. We provide evidence that at least 23, 961 genes are active in the Arabidopsis flower, including 8512 genes that have not been reported as florally expressed previously. We compared gene expression at different stages and found that many genes encoding transcription factors are preferentially expressed in early flower development. Other genes with expression at distinct developmental stages included DUF577 in meiotic cells and DUF220, DUF1216, and Oleosin in stage 12 flowers. DUF1216 and DUF577 are Brassicaceae specific, and together with other families experienced expansion within the Brassicaceae lineage, suggesting novel/greater roles in Brassicaceae floral development than other plants. The large dataset from this study can serve as a resource for expression analysis of genes involved in flower development in Arabidopsis and for comparison with other species. Together, this work provides clues regarding molecular networks underlying flower development.

Developmental genetics of the perianthless flowers and bracts of a paleoherb species, Saururus chinensis

Saururus chinensis is a core member of Saururaceae, a perianthless (lacking petals or sepals) family. Due to its basal phylogenetic position and unusual floral composition, study of this plant family is important for understanding the origin and evolution of perianthless flowers and petaloid bracts among angiosperm species. To isolate genes involved in S. chinensis flower development, subtracted floral cDNA libraries were constructed by using suppression subtractive hybridization (SSH) on transcripts isolated from developing inflorescences and seedling leaves. The subtracted cDNA libraries contained a total of 1,141 ESTs and were used to create cDNA microarrays to analyze transcript profiles of developing inflorescence tissues. Subsequently, qRT-PCR analyses of eight MADS-box transcription factors and in situ hybridizations of two B-class MADS-box transcription factors were performed to verify and extend the cDNA microarray results. Finally, putative phylogenetic relationships within the B-class MADS-box gene family were determined using the discovered S. chinensis B-class genes to compare K-domain sequences with B genes from other basal angiosperms. Two hundred seventy-seven of the 1,141 genes were found to be expressed differentially between S. chinensis inflorescence tissues and seedling leaves, 176 of which were grouped into at least one functional category, including transcription (14.75%), energy (12.59%), metabolism (9.12%), protein-related function (8.99%), and cellular transport (5.76%). qRT-PCR and in situ hybridization of selected MADS-box genes supported our microarray data. Phylogenetic analysis indicated that a total of six B-class MADS-box genes were isolated from S. chinensis. The differential regulation of S. chinensis B-class MADS-box transcription factors likely plays a role during the development of subtending bracts and perianthless flowers. This study contributes to our understanding of inflorescence development in Saururus, and represents an initial step toward understanding the formation of petaloid bracts in this species.

Regulation of the Arabidopsis anther transcriptome by DYT1 for pollen development

Several genes encoding transcription factors have been shown to be essential for male fertility in plants, suggesting that transcriptional regulation is a major mechanism controlling anther development in Arabidopsis. DYSFUNCTIONAL TAPETUM 1 (DYT1), a putative bHLH transcription factor, plays a critical role in regulating tapetum function and pollen development. Here, we compare the transcriptomes of young anthers of wild-type and the dyt1 mutant, demonstrating that DYT1 is upstream of at least 22 genes encoding transcription factors and regulates the expression of a large number of genes, including genes involved in specific metabolic pathways. We also show that DYT1 can bind to DNA in a sequence-specific manner in vitro, and induction of DYT1 activity in vivo activated the expression of the downstream transcription factor genes MYB35 and MS1. We generated DYT1-SRDX transgenic plants whose fertility was dramatically reduced, implying that DYT1 probably acts as a transcriptional activator. Furthermore, we used yeast two-hybrid assays to show that DYT1 forms homodimers and heterodimers with other bHLH transcription factors. Our results demonstrate the important role of DYT1 in regulating anther transcriptome and function, and supporting normal pollen 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.

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.

AMS-dependent and independent regulation of anther transcriptome and comparison with those affected by other Arabidopsis anther genes

BACKGROUND

In flowering plants, the development of male reproductive organs is controlled precisely to achieve successful fertilization and reproduction. Despite the increasing knowledge of genes that contribute to anther development, the regulatory mechanisms controlling this process are still unclear.

RESULTS

In this study, we analyzed the transcriptome profiles of early anthers of sterile mutants aborted microspores (ams) and found that 1,368 genes were differentially expressed in ams compared to wild type anthers, affecting metabolism, transportation, ubiquitination and stress response. Moreover, the lack of significant enrichment of potential AMS binding sites (E-box) in the promoters of differentially expressed genes suggests both direct and indirect regulation for AMS-dependent regulation of anther transcriptome involving other transcription factors. Combining ams transcriptome profiles with those of two other sterile mutants, spl/nzz and ems1/exs, expression of 3,058 genes were altered in at least one mutant. Our investigation of expression patterns of major transcription factor families, such as bHLH, MYB and MADS, suggested that some closely related homologs of known anther developmental genes might also have similar functions. Additionally, comparison of expression levels of genes in different organs suggested that anther-preferential genes could play important roles in anther development.

CONCLUSION

Analysis of ams anther transcriptome and its comparison with those of spl/nzz and ems1/exs anthers uncovered overlapping and distinct sets of regulated genes, including those encoding transcription factors and other proteins. These results support an expanded regulatory network for early anther development, providing a series of hypotheses for future experimentation.

AMS-dependent and independent regulation of anther transcriptome and comparison with those affected by other Arabidopsis anther genes

BACKGROUND

In flowering plants, the development of male reproductive organs is controlled precisely to achieve successful fertilization and reproduction. Despite the increasing knowledge of genes that contribute to anther development, the regulatory mechanisms controlling this process are still unclear.

RESULTS

In this study, we analyzed the transcriptome profiles of early anthers of sterile mutants aborted microspores (ams) and found that 1,368 genes were differentially expressed in ams compared to wild type anthers, affecting metabolism, transportation, ubiquitination and stress response. Moreover, the lack of significant enrichment of potential AMS binding sites (E-box) in the promoters of differentially expressed genes suggests both direct and indirect regulation for AMS-dependent regulation of anther transcriptome involving other transcription factors. Combining ams transcriptome profiles with those of two other sterile mutants, spl/nzz and ems1/exs, expression of 3,058 genes were altered in at least one mutant. Our investigation of expression patterns of major transcription factor families, such as bHLH, MYB and MADS, suggested that some closely related homologs of known anther developmental genes might also have similar functions. Additionally, comparison of expression levels of genes in different organs suggested that anther-preferential genes could play important roles in anther development.

CONCLUSION

Analysis of ams anther transcriptome and its comparison with those of spl/nzz and ems1/exs anthers uncovered overlapping and distinct sets of regulated genes, including those encoding transcription factors and other proteins. These results support an expanded regulatory network for early anther development, providing a series of hypotheses for future experimentation.

An ungrouped plant kinesin accumulates at the preprophase band in a cell cycle-dependent manner

Past phylogenic studies have identified a plant-specific, ungrouped family of kinesins in which the motor domain does not group to one of the fourteen recognized families. Members of this family contain an N-terminal motor domain, a C-terminal armadillo repeat domain and a conserved destruction box (D-BOX) motif. This domain architecture is unique to plants and to a subset of protists. Further characterization of one representative member from Arabidopsis, Arabidopsis thaliana KINESIN ungrouped clade, gene A (AtKINUa), was completed to ascertain its functional role in plants. Fluorescence confocal microscopy revealed an accumulation of ATKINUA:GFP at the preprophase band (PPB) in a cell cycle-dependent manner in Arabidopsis epidermal cells and tobacco BY-2 cells. Fluorescence accumulation was highest during prophase and decreased after nuclear envelope breakdown. A conserved D-BOX motif was identified through alignment of AtKINU homologous sequences. Mutagenesis work with D-BOX revealed that conserved residues were necessary for the observed degradation pattern of ATKINUA:GFP, as well as the targeted accumulation at the PPB. Overall results suggest that AtKINUa is necessary for normal plant growth and/or development and is likely involved with PPB organization through microtubule association and specific cell cycle regulation. The D-BOX motif may function to bridge microtubule organization with changes that occur during progression through mitosis and may represent a novel regulatory motif in plant microtubule motor proteins.

Sex chromosomes in land plants

Sex chromosomes in land plants can evolve as a consequence of close linkage between the two sex determination genes with complementary dominance required to establish stable dioecious populations, and they are found in at least 48 species across 20 families. The sex chromosomes in hepatics, mosses, and gymnosperms are morphologically heteromorphic. In angiosperms, heteromorphic sex chromosomes are found in at least 19 species from 4 families, while homomorphic sex chromosomes occur in 20 species from 13 families. The prevalence of the XY system found in 44 out of 48 species may reflect the predominance of the evolutionary pathway from gynodioecy towards dioecy. All dioecious species have the potential to evolve sex chromosomes, and reversions back from dioecy to various forms of monoecy, gynodioecy, or androdioecy have also occurred. Such reversals may occur especially during the early stages of sex chromosome evolution before the lethality of the YY (or WW) genotype is established.

Evolutionary trends in the floral transcriptome: insights from one of the basalmost angiosperms, the water lily Nuphar advena (Nymphaeaceae)

Current understanding of floral developmental genetics comes primarily from the core eudicot model Arabidopsis thaliana. Here, we explore the floral transcriptome of the basal angiosperm, Nuphar advena (water lily), for insights into the ancestral developmental program of flowers. We identify several thousand Nuphar genes with significantly upregulated floral expression, including homologs of the well-known ABCE floral regulators, deployed in broadly overlapping transcriptional programs across floral organ categories. Strong similarities in the expression profiles of different organ categories in Nuphar flowers are shared with the magnoliid Persea americana (avocado), in contrast to the largely organ-specific transcriptional cascades evident in Arabidopsis, supporting the inference that this is the ancestral condition in angiosperms. In contrast to most eudicots, floral organs are weakly differentiated in Nuphar and Persea, with staminodial intermediates between stamens and perianth in Nuphar, and between stamens and carpels in Persea. Consequently, the predominantly organ-specific transcriptional programs that characterize Arabidopsis flowers (and perhaps other eudicots) are derived, and correlate with a shift towards morphologically distinct floral organs, including differentiated sepals and petals, and a perianth distinct from stamens and carpels. Our findings suggest that the genetic regulation of more spatially discrete transcriptional programs underlies the evolution of floral morphology.

Comparative transcriptomics among floral organs of the basal eudicot Eschscholzia californica as reference for floral evolutionary developmental studies

BACKGROUND

Molecular genetic studies of floral development have concentrated on several core eudicots and grasses (monocots), which have canalized floral forms. Basal eudicots possess a wider range of floral morphologies than the core eudicots and grasses and can serve as an evolutionary link between core eudicots and monocots, and provide a reference for studies of other basal angiosperms. Recent advances in genomics have enabled researchers to profile gene activities during floral development, primarily in the eudicot Arabidopsis thaliana and the monocots rice and maize. However, our understanding of floral developmental processes among the basal eudicots remains limited.

RESULTS

Using a recently generated expressed sequence tag (EST) set, we have designed an oligonucleotide microarray for the basal eudicot Eschscholzia californica (California poppy). We performed microarray experiments with an interwoven-loop design in order to characterize the E. californica floral transcriptome and to identify differentially expressed genes in flower buds with pre-meiotic and meiotic cells, four floral organs at preanthesis stages (sepals, petals, stamens and carpels), developing fruits, and leaves.

CONCLUSIONS

Our results provide a foundation for comparative gene expression studies between eudicots and basal angiosperms. We identified whorl-specific gene expression patterns in E. californica and examined the floral expression of several gene families. Interestingly, most E. californica homologs of Arabidopsis genes important for flower development, except for genes encoding MADS-box transcription factors, show different expression patterns between the two species. Our comparative transcriptomics study highlights the unique evolutionary position of E. californica compared with basal angiosperms and core eudicots.

Transcriptome analysis of various flower and silique development stages indicates a set of class III peroxidase genes potentially involved in pod shattering in Arabidopsis thaliana

BACKGROUND

Plant class III peroxidases exist as a large multigenic family involved in numerous functions suggesting a functional specialization of each gene. However, few genes have been linked with a specific function. Consequently total peroxidase activity is still used in numerous studies although its relevance is questionable. Transcriptome analysis seems to be a promising tool to overcome the difficulties associated with the study of this family. Nevertheless available microarrays are not completely reliable for this purpose. We therefore used a macroarray dedicated to the 73 class III peroxidase genes of A. thaliana to identify genes potentially involved in flower and fruit development.

RESULTS

The observed increase of total peroxidase activity during development was actually correlated with the induction of only a few class III peroxidase genes which supports the existence of a functional specialization of these proteins. We identified peroxidase genes that are predominantly expressed in one development stage and are probable components of the complex gene networks involved in the reproductive phase. An attempt has been made to gain insight into plausible functions of these genes by collecting and analyzing the expression data of different studies in plants. Peroxidase activity was additionally observed in situ in the silique dehiscence zone known to be involved in pod shattering. Because treatment with a peroxidase inhibitor delayed pod shattering, we subsequently studied mutants of transcription factors (TF) controlling this mechanism. Three peroxidases genes -AtPrx13, AtPrx30 and AtPrx55- were altered by the TFs involved in pod shatter.

CONCLUSIONS

Our data illustrated the problems caused by linking only an increase in total peroxidase activity to any specific development stage or function. The activity or involvement of specific class III peroxidase genes needs to be assessed. Several genes identified in our study had not been linked to any particular development stage or function until now. Notably AtPrx13, which is one of the peroxidase genes not present on commercially available microarrays. A systematic survey of class III peroxidase genes expression is necessary to reveal specific class III peroxidase gene functions and the regulation and evolution of this key multifunctional enzyme family. The approach used in this study highlights key individual genes that merit further investigation.

Transcription profiling of fertilization and early seed development events in a solanaceous species using a 7.7 K cDNA microarray from Solanum chacoense ovules

BACKGROUND

To provide a broad analysis of gene expression changes in developing embryos from a solanaceous species, we produced amplicon-derived microarrays with 7741 ESTs isolated from Solanum chacoense ovules bearing embryos from all developmental stages. Our aims were to: 1) identify genes expressed in a tissue-specific and temporal-specific manner; 2) define clusters of genes showing similar patterns of spatial and temporal expression; and 3) identify stage-specific or transition-specific candidate genes for further functional genomic analyses.

RESULTS

We analyzed gene expression during S. chacoense embryogenesis in a series of experiments with probes derived from ovules isolated before and after fertilization (from 0 to 22 days after pollination), and from leaves, anthers, and styles. From the 6374 unigenes present in our array, 1024 genes were differentially expressed (>or= +/- 2 fold change, p value <or= 0.01) in fertilized ovules compared to unfertilized ovules and only limited expression overlap was observed between these genes and the genes expressed in the other tissues tested, with the vast majority of the fertilization-regulated genes specifically or predominantly expressed in ovules (955 genes). During embryogenesis three major expression profiles corresponding to early, middle and late stages of embryo development were identified. From the early and middle stages, a large number of genes corresponding to cell cycle, DNA processing, signal transduction, and transcriptional regulation were found. Defense and stress response-related genes were found in all stages of embryo development. Protein biosynthesis genes, genes coding for ribosomal proteins and other components of the translation machinery were highly expressed in embryos during the early stage. Genes for protein degradation were overrepresented later in the middle and late stages of embryo development. As expected, storage protein transcripts accumulated predominantly in the late stage of embryo development.

CONCLUSION

Our analysis provides the first study in a solanaceous species of the transcriptional program that takes place during the early phases of plant reproductive development, including all embryogenesis steps during a comprehensive time-course. Our comparative expression profiling strategy between fertilized and unfertilized ovules identified a subset of genes specifically or predominantly expressed in ovules while a closer analysis between each consecutive time point allowed the identification of a subset of stage-specific and transition-specific genes.

Identification of candidate genes for grain number in rice (Oryza sativa L.)

Large number of well-filled grains per panicle is an important yield component trait in rice. A combination of QTL mapping and transcriptome profiling was used to identify candidate genes for grain number. A framework linkage map was constructed using 166 SSR markers evenly distributed over the 12 rice chromosomes. QTL mapping using 3 years phenotyping data on a set of recombinant inbred lines derived from a cross between Pusa 1266 (high grain number) and Pusa Basmati 1 (low grain number) identified one consistent QTL qGN4-1 on the long arm of chromosome 4 with major effect on grain number. This QTL was co-localized with major QTLs for primary and secondary branches per panicle, and number of panicles per plant. The QTL interval was narrowed down to 11.1 cM (0.78 Mbp) by targeted enrichment of the region with six additional markers. Microarray transcriptome profiling revealed eight genes in the qGN4-1 region differentially expressed between the two parents during early panicle development. Synteny of this QTL and potential candidates was examined in wheat, barley, maize, sorghum, and Brachypodium to further validate the association.

A dynamic gene expression atlas covering the entire life cycle of rice

Growth and development of a plant are controlled by programmed expression of suits of genes at the appropriate time, tissue and abundance. Although genomic resources have been developed rapidly in recent years in rice, a model plant for cereal genome research, data of gene expression profiling are still insufficient to relate the developmental processes to transcriptomes, leaving a large gap between the genome sequence and phenotype. In this study, we generated genome-wide expression data by hybridizing 190 Affymetrix GeneChip Rice Genome Arrays with RNA from 39 tissues collected throughout the life cycle of the rice plant from two varieties, Zhenshan 97 and Minghui 63. Analyses of the global transcriptomes revealed many interesting features of dynamic patterns of gene expression across the tissues and stages. In total, 38 793 probe sets were detected as expressed and 69% of the expressed transcripts showed significantly variable expression levels among tissues/organs. We found that similarity of transcriptomes among organs corresponded well to their developmental relatedness. About 5.2% of the expressed transcripts showed tissue-specific expression in one or both varieties and 22.7% of the transcripts exhibited constitutive expression including 19 genes with high and stable expression in all the tissues. This dataset provided a versatile resource for plant genomic research, which can be used for associating the transcriptomes to the developmental processes, understanding the regulatory network of these processes, tracing the expression profile of individual genes and identifying reference genes for quantitative expression analyses.

Genes and functions controlled by floral organ identity genes

Floral organ identity genes specify the identity of floral organs in a manner analogous to the specification of body segments by Hox genes in animals. Different combinations of organ identity genes co-ordinate the expression of genes required for the development of each type of floral organ, from organ initiation until final differentiation. Here, I review what is known about the genes and functions subordinate to the organ identity genes. The sets of target genes change as organ development progresses and ultimately organ identity genes modify the expression of thousands of genes with a multitude of predicted functions, particularly in reproductive organs. However, genes involved in transcriptional control and hormone functions feature prominently among the early and direct targets. Functional analysis showed that control of organ-specific tissues and structures can be delegated to specialised intermediate regulators, but organ identity genes also fine-tune genes with general roles in shoot organ development, consistent with the notion that organ identity genes modify a core leaf-like developmental program. Future challenges include obtaining data with cellular resolution, predictive modelling of the regulatory network, and quantitative analysis of how organ identity genes and their targets control cell behaviour and ultimately organ shape.

Comparison of next generation sequencing technologies for transcriptome characterization

Background

We have developed a simulation approach to help determine the optimal mixture of sequencing methods for most complete and cost effective transcriptome sequencing. We compared simulation results for traditional capillary sequencing with "Next Generation" (NG) ultra high-throughput technologies. The simulation model was parameterized using mappings of 130,000 cDNA sequence reads to the Arabidopsis genome (NCBI Accession SRA008180.19). We also generated 454-GS20 sequences and de novo assemblies for the basal eudicot California poppy (Eschscholzia californica) and the magnoliid avocado (Persea americana) using a variety of methods for cDNA synthesis.

Results

The Arabidopsis reads tagged more than 15,000 genes, including new splice variants and extended UTR regions. Of the total 134,791 reads (13.8 MB), 119,518 (88.7%) mapped exactly to known exons, while 1,117 (0.8%) mapped to introns, 11,524 (8.6%) spanned annotated intron/exon boundaries, and 3,066 (2.3%) extended beyond the end of annotated UTRs. Sequence-based inference of relative gene expression levels correlated significantly with microarray data. As expected, NG sequencing of normalized libraries tagged more genes than non-normalized libraries, although non-normalized libraries yielded more full-length cDNA sequences. The Arabidopsis data were used to simulate additional rounds of NG and traditional EST sequencing, and various combinations of each. Our simulations suggest a combination of FLX and Solexa sequencing for optimal transcriptome coverage at modest cost. We have also developed ESTcalc http://fgp.huck.psu.edu/NG_Sims/ngsim.pl, an online webtool, which allows users to explore the results of this study by specifying individualized costs and sequencing characteristics.

Conclusion

NG sequencing technologies are a highly flexible set of platforms that can be scaled to suit different project goals. In terms of sequence coverage alone, the NG sequencing is a dramatic advance over capillary-based sequencing, but NG sequencing also presents significant challenges in assembly and sequence accuracy due to short read lengths, method-specific sequencing errors, and the absence of physical clones. These problems may be overcome by hybrid sequencing strategies using a mixture of sequencing methodologies, by new assemblers, and by sequencing more deeply. Sequencing and microarray outcomes from multiple experiments suggest that our simulator will be useful for guiding NG transcriptome sequencing projects in a wide range of organisms.

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.

Evolutionary conservation of neighbouring gene pairs in plants

Evolutionary conservation of neighbouring gene pairs has been widely explored in many species, but remains poorly understood in plants. The availability of several plant genome sequences allows for an in-depth investigation of this problem in plants. Here, we analyzed the phylogenetic conservation of physically linked gene pairs in nine plant genomes and compared the conservation in different orientations. We also examined several potential determinants to detect whether they affect the conservation of neighbouring gene pairs. Our results suggested that among the three types of neighbouring gene pairs, closely linked parallel pairs might be the least conserved. Intergenic distance was shown to be a major determinant of linkage conservation, suggesting that the conservation of gene order in plants was determined primarily by chance. The enrichment of housekeeping genes was identified to contribute to the conservation of all three types and the enrichment of genes involved in protein metabolism might contribute to the conservation of parallel pairs. Moreover, a co-expressed signal was detected in conserved divergent pairs, which might be determined by intergenic distance.

Molecular processes underlying the floral transition in the soybean shoot apical meristem

The transition to flowering is characterized by a shift of the shoot apical meristem (SAM) from leaf production to the initiation of a floral meristem. The flowering process is of vital importance for agriculture, but the associated events or regulatory pathways in the SAM are not well understood, especially at a system level. To address this issue, we have used a GeneChip containing 37 744 probe sets to generate a temporal profile of gene expression during the floral initiation process in the SAM of the crop legume, soybean (Glycine max). A total of 331 transcripts displayed significant changes in their expression profiles. The in silico and RT-PCR analysis on differentially regulated transcripts implies the intriguing involvement of sugar, auxin or abscisic acid (ABA) in events prior to the induction of floral homeotic transcripts. The novel involvement of ABA in the floral transition is further implicated by immunoassay, suggesting an increase in ABA levels in the SAM during this developmental transition. Furthermore, in situ localization, together with in silico data demonstrating a marked enhancement of abiotic stress-related transcripts, such as trehalose metabolism genes in SAMs, points to an overlap of abiotic stress and floral signalling pathways.

Analysis of the Arabidopsis floral proteome: detection of over 2 000 proteins and evidence for posttranslational modifications

The proteome of the Arabidopsis flower has not been extensively studied previously. Here, we report a proteomic analysis of the wild type Arabidopsis flower. Using both two-dimensional electrophoresis/mass spectrometry (2-DGE/MS) and multi-dimensional protein identification technology (MudPIT) approaches, we identified 2,446 proteins. Although a single experiment or analysis uncovered only a subset of the proteins we identified, a combination of multiple experiments and analyses facilitated the detection of a greater number of proteins. When proteins are grouped according to RNA expression levels revealed by microarray experiments, we found that proteins encoded by genes with relatively high levels of expression were detected with greater frequencies. On the other hand, at the level of the individual gene/protein, there was not a good correlation between protein spot intensity and microarray values. We also obtained strong evidence for post-translational modification from 2-DGE and MudPIT data. We detected proteins that are annotated to function in protein synthesis, folding, modification, and degradation, as well as the presence of regulatory proteins such as transcription factors and protein kinases. Finally, sequence and evolutionary analysis of genes for active methyl group metabolisms suggests that these genes are highly conserved. Our results allow the formulation of hypotheses regarding post-translational regulation of proteins in the flower, providing new understanding about Arabidopsis flower development and physiology.

Evolutionary history of histone demethylase families: distinct evolutionary patterns suggest functional divergence

Background

Histone methylation can dramatically affect chromatin structure and gene expression and was considered irreversible until recent discoveries of two families of histone demethylases, the KDM1 (previously LSD1) and JmjC domain-containing proteins. These two types of proteins have different functional domains and distinct substrate specificities. Although more and more KDM1 and JmjC proteins have been shown to have histone demethylase activity, our knowledge about their evolution history is limited.

Results

We performed systematic phylogenetic analysis of these histone demethylase families and uncovered different evolutionary patterns. The KDM1 genes have been maintained with a stable low copy number in most organisms except for a few duplication events in flowering plants. In contrast, multiple genes for JmjC proteins with distinct domain architectures were present before the split of major eukaryotic groups, and experienced subsequent birth-and-death evolution. In addition, distinct evolutionary patterns can also be observed between animal and plant histone demethylases in both families. Furthermore, our results showed that some JmjC subfamilies contain only animal genes with specific demethylase activities, but do not have plant members.

Conclusion

Our study improves the understanding about the evolutionary history of KDM1 and JmjC genes and provides valuable insights into their functions. Based on the phylogenetic relationship, we discussed possible histone demethylase activities for several plant JmjC proteins. Finally, we proposed that the observed differences in evolutionary pattern imply functional divergence between animal and plant histone demethylases.

Phylogenetic analysis of the plant-specific zinc finger-homeobox and mini zinc finger gene families

Zinc finger-homeodomain proteins (ZHD) are present in many plants; however, the evolutionary history of the ZHD gene family remains largely unknown. We show here that ZHD genes are plant-specific, nearly all intronless, and related to MINI ZINC FINGER (MIF) genes that possess only the zinc finger. Phylogenetic analyses of ZHD genes from representative land plants suggest that non-seed plant ZHD genes occupy basal positions and angiosperm homologs form seven distinct clades. Several clades contain genes from two or more major angiosperm groups, including eudicots, monocots, magnoliids, and other basal angiosperms, indicating that several duplications occurred before the diversification of flowering plants. In addition, specific lineages have experienced more recent duplications. Unlike the ZHD genes, MIFs are found only from seed plants, possibly derived from ZHDs by loss of the homeodomain before the divergence of seed plants. Moreover, the MIF genes have also undergone relatively recent gene duplications. Finally, genome duplication might have contributed substantially to the expansion of family size in angiosperms and caused a high level of functional redundancy/overlap in these genes.

Functional divergence of the duplicated AtKIN14a and AtKIN14b genes: critical roles in Arabidopsis meiosis and gametophyte development

Gene duplication is important for gene family evolution, allowing for functional divergence and innovation. In flowering plants, duplicated genes are widely observed, and functional redundancy of closely related duplicates has been reported, but few cases of functional divergence of close duplicates have been described. Here, we show that the Arabidopsis AtKIN14a and AtKIN14b genes encoding highly similar kinesins are two of the most closely related Arabidopsis paralogs, which were formed by a duplication event that occurred after the split of Arabidopsis and poplar. In addition, AtKIN14a and AtKIN14b exhibit varying degrees of coding sequence divergence. Further genetic studies of plants carrying atkin14a and/or atkin14b mutations indicate that, although these two genes have similar functions, there is clear evidence for functional divergence. Although both genes are important for male and female meiosis, AtKIN14a plays a more critical role in male meiosis than AtKIN14b. Moreover, either one of these two genes is necessary and sufficient for gametophyte development, indicating that they are redundant for this function. Therefore, AtKIN14a and AtKIN14b together play important roles in controlling plant reproductive development. Our results suggest that the AtKIN14a and AtKIN14b genes have retained similar functions in gametophyte development and female meiosis, but have evolved partially distinct functions in male meiosis, with AtKIN14a playing a more substantive role.

Global analysis of gene expression in flower buds of Ms-cd1 Brassica oleracea conferring male sterility by using an Arabidopsis microarray

The dominant male sterility gene Ms-cd1 is identified in Brassica oleracea. Electron microscopical observations revealed that abortion of pollen development starts after tetrad formation. This important male sterility phenotype is characterized by lack of degradation of the primary pollen mother cell (PMC) wall and delayed degradation of callose surrounding the tetrads and thus arrest of microspore release. Gene expression of the male sterile and fertile buds was analyzed by heterologous hybridization of Brassica oleracea cRNA onto an Arabidopsis whole genome oligonucleotide microarray. A total of 277 suppressed genes including 40 kinase-, 32 cell wall modification and 29 transport related genes were found to be significantly down regulated >3-fold in the male sterile mutant. The vast majority of the differentially expressed transcripts are found to present late pollen stage specific genes. Kinase genes, cell wall modification genes and ion transport genes were greatly over-represented when compared to their percentage of all flower bud expressed genes and represent 36.5% of the genes suppressed by Ms-cd1. Our results also suggest that Ms-cd1 may blocks an anther developmental pathway with a small number of genes suppressed in tapetum cells which prevent the degradation of callose and PMC wall, which further leads to the suppression of a large number of genes involved in signaling pathways, cell wall modification and ion transport in pollen grains.

Global expression profiling applied to the analysis of Arabidopsis stamen development

To obtain detailed information about gene expression during stamen development in Arabidopsis (Arabidopsis thaliana), we compared, by microarray analysis, the gene expression profile of wild-type inflorescences to those of the floral mutants apetala3, sporocyteless/nozzle, and male sterile1 (ms1), in which different aspects of stamen formation are disrupted. These experiments led to the identification of groups of genes with predicted expression at early, intermediate, and late stages of stamen development. Validation experiments using in situ hybridization confirmed the predicted expression patterns. Additional experiments aimed at characterizing gene expression specifically during microspore formation. To this end, we compared the gene expression profiles of wild-type flowers of distinct developmental stages to those of the ms1 mutant. Computational analysis of the datasets derived from this experiment led to the identification of genes that are likely involved in the control of key developmental processes during microsporogenesis. We also identified a large number of genes whose expression is prolonged in ms1 mutant flowers compared to the wild type. This result suggests that MS1, which encodes a putative transcriptional regulator, is involved in the stage-specific repression of these genes. Lastly, we applied reverse genetics to characterize several of the genes identified in the microarray experiments and uncovered novel regulators of microsporogenesis, including the transcription factor MYB99 and a putative phosphatidylinositol 4-kinase.

Differential gene expression in Arabidopsis wild-type and mutant anthers: insights into anther cell differentiation and regulatory networks

In flowering plants, the anther contains highly specialized reproductive and somatic cells that are required for male fertility. Genetic studies have uncovered several genes that are important for anther development. However, little information is available regarding most genes active during anther development, including possible relationships between these genes and genetically defined regulators. In Arabidopsis, two previously isolated male-sterile mutants display dramatically altered anther cell differentiation patterns. The sporocyteless (spl)/nozzle (nzz) mutant is defective in the differentiation of primary sporogenous cells into microsporocytes, and does not properly form the anther wall. The excess microsporocytes1 (ems1)/extrasporogenous cells (exs) mutants produce excess microsporocytes at the expense of the tapetum. To gain additional insights into microsporocyte and tapetum differentiation and to uncover potential genetic interactions, expression profiles were compared between wild-type anthers (stage 4-6) and those of the spl or ems1 mutants. A total of 1954 genes were found to be differentially expressed in the ems1 and/or spl anthers, and these were grouped into 14 co-expression clusters. The presence of genes with known and predicted functions in specific clusters suggests potential functions for other genes in the same cluster. To obtain clues about possible co-regulation within co-expression clusters, we searched for shared cis-regulatory motifs in putative promoter regions. Our analyses were combined with data from previous studies to develop a model of the anther gene regulatory network. This model includes hypotheses that can be tested experimentally to gain further understanding of the mechanisms controlling anther development.

Molecular evolutionary analyses of the Arabidopsis L7 ribosomal protein gene family

Cytoplasmic ribosomal protein (r-protein) genes in Arabidopsis thaliana are encoded by 80 multigene families that contain between two and seven members. Gene family members are typically similar at the protein sequence level, with the most divergent members of any gene family retaining 94% identity, on average. However, three Arabidopsis r-protein families - S15a, L7 and P2 - contain highly divergent family members. Here, we investigated the organization, structure, expression and molecular evolution of the L7 r-protein family. Phylogenetic analyses showed that L7 r-protein gene family members constitute two distinct phylogenetic groups. The first group including RPL7B, RPL7C and RPL7D has homologs in plants, animals and fungi. The second group represented by RPL7A is found in plants but has no orthologs from other fully-sequenced eukaryotic genomes. These two groups may have derived from a duplication event prior to the divergence of animals and plants. All four L7 r-protein genes are expressed and all exhibit a differential expression in inflorescence and flowers. RPL7A and RPL7B are less expressed than the other genes in all tissues analyzed. Molecular characterization of nucleic and protein sequences of L7 r-protein genes and analysis of their codon usage did not indicate any functional divergence. The probable evolution of an extra-ribosomal function of group 2 genes is discussed.

The floral genome: an evolutionary history of gene duplication and shifting patterns of gene expression

Through multifaceted genome-scale research involving phylogenomics, targeted gene surveys, and gene expression analyses in diverse basal lineages of angiosperms, our studies provide insights into the most recent common ancestor of all extant flowering plants. MADS-box gene duplications have played an important role in the origin and diversification of angiosperms. Furthermore, early angiosperms possessed a diverse tool kit of floral genes and exhibited developmental 'flexibility', with broader patterns of expression of key floral organ identity genes than are found in eudicots. In particular, homologs of B-function MADS-box genes are more broadly expressed across the floral meristem in basal lineages. These results prompted formulation of the 'fading borders' model, which states that the gradual transitions in floral organ morphology observed in some basal angiosperms (e.g. Amborella) result from a gradient in the level of expression of floral organ identity genes across the developing floral meristem.

Patterns of gene duplication and functional diversification during the evolution of the AP1/SQUA subfamily of plant MADS-box genes

Members of the AP1/SQUA subfamily of plant MADS-box genes play broad roles in the regulation of reproductive meristems, the specification of sepal and petal identities, and the development of leaves and fruits. It has been shown that AP1/SQUA-like genes are angiosperm-specific, and have experienced several major duplication events. However, the evolutionary history of this subfamily is still uncertain. Here, we report the isolation of 14 new AP1/SQUA-like genes from seven early-diverging eudicots and the identification of 11 previously uncharacterized ESTs and genomic sequences from public databases. Sequence comparisons of these and other published sequences reveal a conserved C-terminal region, the FUL motif, in addition to the known euAP1/paleoAP1 motif, in AP1/SQUA-like proteins. Phylogenetic analyses further suggest that there are three major lineages (euAP1, euFUL, and AGL79) in core eudicots, likely resulting from two close duplication events that predated the divergence of core eudicots. Among the three lineages, euFUL is structurally very similar to FUL-like genes from early-diverging eudicots and basal angiosperms, whereas euAP1 might have originally been generated through a 1-bp deletion in the exon 8 of an ancestral euFUL- or FUL-like gene. Because euFUL- and FUL-like genes usually have broad expression patterns, we speculate that AP1/SQUA-like genes initially had broad functions. Based on these observations, the evolutionary fates of duplicate genes and the contributions of the frameshift mutation and alternative splicing to functional diversity are discussed.

Transcriptional profiling of wheat caryopsis development using cDNA microarrays

The expression of 7,835 genes in developing wheat caryopses was analyzed using cDNA arrays. Using a mixed model analysis of variance (ANOVA) method, 29% (2,237) of the genes on the array were identified to be differentially expressed at the 6 different time-points examined, which covers the developmental stages from coenocytic endosperm to physiological maturity. Comparison of genes differentially expressed between two time-points revealed a dynamic transcript accumulation profile with major re-programming events that occur at 3-7, 7-14 and 21-28 DPA. A k-means clustering algorithm grouped the differentially expressed genes into 10 clusters, revealing co-expression of genes involved in the same pathway such as carbohydrate and protein synthesis or preparation for desiccation. Functional annotation of genes that show peak expression at specific time-points correlated with the developmental events associated with the respective stages. Results provide information on the temporal expression during caryopsis development for a significant number of differentially expressed genes with unknown function.

EST database for early flower development in California poppy (Eschscholzia californica Cham., Papaveraceae) tags over 6,000 genes from a basal eudicot

The Floral Genome Project (FGP) selected California poppy (Eschscholzia californica Cham. ssp. Californica) to help identify new florally-expressed genes related to floral diversity in basal eudicots. A large, non-normalized cDNA library was constructed from premeiotic and meiotic floral buds and sequenced to generate a database of 9,079 high quality Expressed Sequence Tags (ESTs). These sequences clustered into 5,713 unigenes, including 1,414 contigs and 4,299 singletons. Homologs of genes regulating many aspects of flower development were identified, including those for organ identity and development, cell and tissue differentiation, cell cycle control, and secondary metabolism. Over 5% of the transcriptome consisted of homologs to known floral gene families. Most are the first representatives of their respective gene families in basal eudicots and their conservation suggests they are important for floral development and/or function. App. 10% of the transcripts encoded transcription factors and other regulatory genes, including nine genes from the seven major lineages of the important MADS-box family of developmental regulators. Homologs of alkaloid pathway genes were also recovered, providing opportunities to explore adaptive evolution in secondary products. Furthermore, comparison of the poppy ESTs with the Arabidopsis genome provided support for putative Arabidopsis genes that previously lacked annotation. Finally, over 1,800 unique sequences had no observable homology in the public databases. The California poppy EST database and library will help bridge our understanding of flower initiation and development among higher eudicot and monocot model plants and provide new opportunities for comparative analysis of gene families across angiosperm species.

Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis

The basic/helix-loop-helix (bHLH) transcription factors and their homologs form a large family in plant and animal genomes. They are known to play important roles in the specification of tissue types in animals. On the other hand, few plant bHLH proteins have been studied functionally. Recent completion of whole genome sequences of model plants Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) allows genome-wide analysis and comparison of the bHLH family in flowering plants. We have identified 167 bHLH genes in the rice genome, and their phylogenetic analysis indicates that they form well-supported clades, which are defined as subfamilies. In addition, sequence analysis of potential DNA-binding activity, the sequence motifs outside the bHLH domain, and the conservation of intron/exon structural patterns further support the evolutionary relationships among these proteins. The genome distribution of rice bHLH genes strongly supports the hypothesis that genome-wide and tandem duplication contributed to the expansion of the bHLH gene family, consistent with the birth-and-death theory of gene family evolution. Bioinformatics analysis suggests that rice bHLH proteins can potentially participate in a variety of combinatorial interactions, endowing them with the capacity to regulate a multitude of transcriptional programs. In addition, similar expression patterns suggest functional conservation between some rice bHLH genes and their close Arabidopsis homologs.

Alternative oxidases in Arabidopsis: A comparative analysis of differential expression in the gene family provides new insights into function of non-phosphorylating bypasses

The emergence of Arabidopsis as a model plant provides an opportunity to gain insights into the role of the alternative oxidase that cannot be as readily achieved in other plant species. The analysis of extensive mRNA expression data indicates that all five Aox genes (Aox1a, 1b, 1c, 1d and 2) are expressed, but organ and developmental regulation are evident, suggesting regulatory specialisation of Aox gene members. The stress-induced nature of the alternative pathway in a variety of plants is further supported in Arabidopsis as Aox1a and Aox1d are amongst the most stress responsive genes amongst the hundreds of known genes encoding mitochondrial proteins. Analysis of genes co-expressed with Aoxs from studies of responses to various treatments altering mitochondrial functions and/or from plants with altered Aox levels reveals that: (i) this gene set encodes more functions outside the mitochondrion than functions in mitochondria, (ii) several pathways for induction exist and there is a difference in the magnitude of the induction in each pathway, (iii) the magnitude of induction may depend on the endogenous levels of Aox, and (iv) induction of Aox can be oxidative stress-dependent or -independent depending on the gene member and the tissue analysed. An overall role for Aox in re-programming cellular metabolism in response to the ever changing environment encountered by plants is proposed.

Genome-wide analysis of gene expression during early Arabidopsis flower development

Detailed information about stage-specific changes in gene expression is crucial for the understanding of the gene regulatory networks underlying development. Here, we describe the global gene expression dynamics during early flower development, a key process in the life cycle of a plant, during which floral patterning and the specification of floral organs is established. We used a novel floral induction system in Arabidopsis, which allows the isolation of a large number of synchronized floral buds, in conjunction with whole-genome microarray analysis to identify genes with differential expression at distinct stages of flower development. We found that the onset of flower formation is characterized by a massive downregulation of genes in incipient floral primordia, which is followed by a predominance of gene activation during the differentiation of floral organs. Among the genes we identified as differentially expressed in the experiment, we detected a significant enrichment of closely related members of gene families. The expression profiles of these related genes were often highly correlated, indicating similar temporal expression patterns. Moreover, we found that the majority of these genes is specifically up-regulated during certain developmental stages. Because co-expressed members of gene families in Arabidopsis frequently act in a redundant manner, these results suggest a high degree of functional redundancy during early flower development, but also that its extent may vary in a stage-specific manner.