Separated transcriptomes of male gametophyte and tapetum in rice: validity of a laser microdissection (LM) microarray

Meta-Analysis of Microarray Data and Their Utility in Dissecting the Mapped QTLs for Heat Acclimation in Rice

In the current global warming scenario, it is imperative to develop crops with improved heat tolerance or acclimation, for which knowledge of major heat stress-tolerant genes or genomic regions is a prerequisite. Though several quantitative trait loci (QTLs) for heat tolerance have been mapped in rice, candidate genes from these QTLs have not been reported yet. The meta-analysis of microarray datasets for heat stress in rice can give us a better genomic resource for the dissection of QTLs and the identification of major candidate genes for heat stress tolerance. In the present study, a database, RiceMetaSys-H, comprising 4227 heat stress-responsive genes (HRGs), was created using seven publicly available microarray datasets. This included in-house-generated microarray datasets of Nagina 22 (N22) and IR64 subjected to 8 days of heat stress. The database has provisions for searching the HRGs through genotypes, growth stages, tissues, and physical intervals in the genome, as well as Locus IDs, which provide complete information on the HRGs with their annotations and fold changes, along with the experimental material used for the analysis. The up-regulation of genes involved in hormone biosynthesis and signalling, sugar metabolism, carbon fixation, and the ROS pathway were found to be the key mechanisms of enhanced heat tolerance. Integrating variant and expression analysis, the database was used for the dissection of the major effect of QTLs on chromosomes 4, 5, and 9 from the IR64/N22 mapping population. Out of the 18, 54, and 62 genes in these three QTLs, 5, 15, and 12 genes harboured non-synonymous substitutions. Fifty-seven interacting genes of the selected QTLs were identified by a network analysis of the HRGs in the QTL regions. Variant analysis revealed that the proportion of unique amino acid substitutions (between N22/IR64) in the QTL-specific genes was much higher than the common substitutions, i.e., 2.58:0.88 (2.93-fold), compared to the network genes at a 0.88:0.67 (1.313-fold) ratio. An expression analysis of these 89 genes showed 43 DEGs between IR64/N22. By integrating the expression profiles, allelic variations, and the database, four robust candidates (LOC_Os05g43870, LOC_Os09g27830, LOC_Os09g27650, andLOC_Os09g28000) for enhanced heat stress tolerance were identified. The database thus developed in rice can be used in breeding to combat high-temperature stress.

Comparative transcriptomic analysis of germinating rice seedlings to individual and combined anaerobic and cold stress

BACKGROUND

Rice is one of the most important cereals consumed worldwide. Two major abiotic factors affecting rice plants in different growth stages are flooding stress and cold stress. These abiotic stresses can take place independently or simultaneously and significantly affect rice plants during germination and seedling growth. Fortunately, a wide array of phenotypic responses conferring flooding stress and chilling stress tolerance exist within the rice germplasm, indicating the presence of different molecular mechanisms underlying tolerance to these stresses. Understanding these differences may assist in developing improved rice cultivars having higher tolerance to both stresses. In this study, we conducted a comparative global gene expression analysis of two rice genotypes with contrasting phenotypes under cold stress, anaerobic stress, and combined cold and anaerobic stress during germination.

RESULTS

The differential gene expression analysis revealed that 5571 differentially expressed genes (DEGs), 7206 DEGs, and 13279 DEGs were identified under anaerobic stress, cold stress, and combined stress, respectively. Genes involved in the carbohydrate metabolic process, glucosyltransferase activity, regulation of nitrogen compound metabolic process, protein metabolic process, lipid metabolic process, cellular nitrogen compound biosynthetic process, lipid biosynthetic process, and a microtubule-based process were enriched across all stresses. Notably, the common Gene Ontology (GO) analysis identified three hub genes, namely Os08g0176800 (similar to mRNA-associated protein mrnp 41), Os11g0454200 (dehydrin), and OS10g0505900 (expressed protein).

CONCLUSION

A large number of differentially expressed genes were identified under anaerobic, cold conditions during germination and the combination of the two stress conditions in rice. These results will assist in the identification of promising candidate genes for possible manipulation toward rice crops that are more tolerant under flooding and cold during germination, both independently and concurrently.

RNA-seq analysis of synchronized developing pollen isolated from a single anther

Pollen development, from unicellular microspores to anthesis, is a complex process involving the coordinated specification, differentiation and functions of different cell types. Key to understanding this development is identifying the genes expressed at precise stages of development. However, transcriptomic studies on pollen prior to anthesis are complicated by the inaccessible nature of pollen developing in the anther and the resistant pollen wall. To assist with understanding gene expression during pollen development we have developed a protocol to perform RNA-Seq on pollen isolated from a single anther (SA RNA-Seq). The protocol involves removing pollen from a single anther for analysis and viewing the remaining pollen to determine the developmental stage. The isolated pollen is chemically lysed and mRNA isolated from the lysate using an oligo-dT column before library preparation. Here, we report on the development and testing of our method and the generation of a transcriptome for three stages of pollen development from Arabidopsis (Arabidopsis thaliana) and two stages from male kiwifruit (Actinidia chinensis). This protocol enables the transcriptome of precise developmental stages of pollen to be analyzed, and uses a small number of plants, potentially facilitating studies that require a range of treatments or the analysis of the first generation of transgenic plants.

Use of RNAi With OsMYB76R as a Reporter for Candidate Genes Can Efficiently Create and Verify Gametophytic Male Sterility in Rice

Gametophytic male sterility (GMS) plays an important role in the study of pollen development and seed propagation of recessive nuclear male sterile lines insensitive to the environmental conditions in hybrid rice breeding. Since the inherent phenotypic and genetic characteristics of GMS, it is very difficult to find and identify the GMS mutants. However, due to the abundance of gene transcription data, a large number of pollen-specific genes have been found, and most of them may be associated with GMS. To promote the study of these genes in pollen development and heterosis utilization, in this study, an easy and efficient method of creating and identifying GMS was established using RNAi and OsMYB76R as a reporter. First, the OsC1/OsMYB76 gene involved in anthocyanin synthesis was modified, and we have validated that the modified OsMYB76R is workable as the same as the pre-modified OsMYB76 gene. Then, the ascorbic acid oxidase gene OsPTD1 was downregulated using RNAi, driven by its own promoter that resulted in abnormal pollen tube growth. Finally, the RNAi elements were linked with OsMYB76R and transformed into an osmyb76 mutant, and the distortion of purple color segregation was found in T₁ and F₁ generations. This indicates that the OsPTD1 GMS was prepared successfully. Compared to current methods, there are several advantages to this method. First, time is saved in material preparation, as one generation less needs to be compared than in the conventional method, and mutation screening can be avoided. In addition, for identification, the cost is lower; PCR, electrophoresis, and other processes are not needed; and no expensive chemicals or instruments are required. Finally, the results are more accurate, with much lower background effects, and no damage to the plant. The result is an easy, efficient, low-cost, and accurate method of preparing and identifying GMS genes.

Morpho-Physiological and Transcriptome Changes in Tomato Anthers of Different Developmental Stages under Drought Stress

Drought limits the growth and productivity of plants. Reproductive development is sensitive to drought but the underlying physiological and molecular mechanisms remain unclear in tomatoes. Here, we investigated the effect of drought on tomato floral development using morpho-physiological and transcriptome analyses. Drought-induced male sterility through abnormal anther development includes pollen abortion, inadequate pollen starch accumulation and anther indehiscence which caused floral bud and opened flower abortions and reduced fruit set/yield. Under drought stress (DS), pollen mother cell to meiotic (PMC-MEI) anthers survived whereas tetrad to vacuolated uninucleate microspore (TED-VUM) anthers aborted. PMC-MEI anthers had lower ABA increase, reduced IAA and elevated sugar contents under DS relative to well-watered tomato plants. However, TED-VUM anthers had higher ABA increase and IAA levels, and lower accumulation of soluble sugars, indicating abnormal carbohydrate and hormone metabolisms when exposed to drought-stress conditions. Moreover, RNA-Seq analysis identified altogether >15,000 differentially expressed genes that were assigned to multiple pathways, suggesting that tomato anthers utilize complicated mechanisms to cope with drought. In particular, we found that tapetum development and ABA homeostasis genes were drought-induced while sugar utilization and IAA metabolic genes were drought-repressed in PMC-MEI anthers. Our results suggest an important role of phytohormones metabolisms in anther development under DS and provide novel insight into the molecular mechanism underlying drought resistance in tomatoes.

A Rapid Pipeline for Pollen- and Anther-Specific Gene Discovery Based on Transcriptome Profiling Analysis of Maize Tissues

Recently, crop breeders have widely adopted a new biotechnology-based process, termed Seed Production Technology (SPT), to produce hybrid varieties. The SPT does not produce nuclear male-sterile lines, and instead utilizes transgenic SPT maintainer lines to pollinate male-sterile plants for propagation of nuclear-recessive male-sterile lines. A late-stage pollen-specific promoter is an essential component of the pollen-inactivating cassette used by the SPT maintainers. While a number of plant pollen-specific promoters have been reported so far, their usefulness in SPT has remained limited. To increase the repertoire of pollen-specific promoters for the maize community, we conducted a comprehensive comparative analysis of transcriptome profiles of mature pollen and mature anthers against other tissue types. We found that maize pollen has much less expressed genes (>1 FPKM) than other tissue types, but the pollen grain has a large set of distinct genes, called pollen-specific genes, which are exclusively or much higher (100 folds) expressed in pollen than other tissue types. Utilizing transcript abundance and correlation coefficient analysis, 1215 mature pollen-specific (MPS) genes and 1009 mature anther-specific (MAS) genes were identified in B73 transcriptome. These two gene sets had similar GO term and KEGG pathway enrichment patterns, indicating that their members share similar functions in the maize reproductive process. Of the genes, 623 were shared between the two sets, called mature anther- and pollen-specific (MAPS) genes, which represent the late-stage pollen-specific genes of the maize genome. Functional annotation analysis of MAPS showed that 447 MAPS genes (71.7% of MAPS) belonged to genes encoding pollen allergen protein. Their 2-kb promoters were analyzed for cis-element enrichment and six well-known pollen-specific cis-elements (AGAAA, TCCACCA, TGTGGTT, [TA]AAAG, AAATGA, and TTTCT) were found highly enriched in the promoters of MAPS. Interestingly, JA-responsive cis-element GCC box (GCCGCC) and ABA-responsive cis-element-coupling element1 (ABRE-CE1, CCACC) were also found enriched in the MAPS promoters, indicating that JA and ABA signaling likely regulate pollen-specific MAPS expression. This study describes a robust and straightforward pipeline to discover pollen-specific promotes from publicly available data while providing maize breeders and the maize industry a number of late-stage (mature) pollen-specific promoters for use in SPT for hybrid breeding and seed production.

Deep learning-based six-type classifier for lung cancer and mimics from histopathological whole slide images: a retrospective study

BACKGROUND

Targeted therapy and immunotherapy put forward higher demands for accurate lung cancer classification, as well as benign versus malignant disease discrimination. Digital whole slide images (WSIs) witnessed the transition from traditional histopathology to computational approaches, arousing a hype of deep learning methods for histopathological analysis. We aimed at exploring the potential of deep learning models in the identification of lung cancer subtypes and cancer mimics from WSIs.

METHODS

We initially obtained 741 WSIs from the First Affiliated Hospital of Sun Yat-sen University (SYSUFH) for the deep learning model development, optimization, and verification. Additional 318 WSIs from SYSUFH, 212 from Shenzhen People's Hospital, and 422 from The Cancer Genome Atlas were further collected for multi-centre verification. EfficientNet-B5- and ResNet-50-based deep learning methods were developed and compared using the metrics of recall, precision, F1-score, and areas under the curve (AUCs). A threshold-based tumour-first aggregation approach was proposed and implemented for the label inferencing of WSIs with complex tissue components. Four pathologists of different levels from SYSUFH reviewed all the testing slides blindly, and the diagnosing results were used for quantitative comparisons with the best performing deep learning model.

RESULTS

We developed the first deep learning-based six-type classifier for histopathological WSI classification of lung adenocarcinoma, lung squamous cell carcinoma, small cell lung carcinoma, pulmonary tuberculosis, organizing pneumonia, and normal lung. The EfficientNet-B5-based model outperformed ResNet-50 and was selected as the backbone in the classifier. Tested on 1067 slides from four cohorts of different medical centres, AUCs of 0.970, 0.918, 0.963, and 0.978 were achieved, respectively. The classifier achieved high consistence to the ground truth and attending pathologists with high intraclass correlation coefficients over 0.873.

CONCLUSIONS

Multi-cohort testing demonstrated our six-type classifier achieved consistent and comparable performance to experienced pathologists and gained advantages over other existing computational methods. The visualization of prediction heatmap improved the model interpretability intuitively. The classifier with the threshold-based tumour-first label inferencing method exhibited excellent accuracy and feasibility in classifying lung cancers and confused nonneoplastic tissues, indicating that deep learning can resolve complex multi-class tissue classification that conforms to real-world histopathological scenarios.

Dynamics of cell wall structure and related genomic resources for drought tolerance in rice

Cell wall plasticity plays a very crucial role in vegetative and reproductive development of rice under drought and is a highly potential trait for improving rice yield under drought. Drought is a major constraint in rice (Oryza sativa L.) cultivation severely affecting all developmental stages, with the reproductive stage being the most sensitive. Rice plants employ multiple strategies to cope with drought, in which modification in cell wall dynamics plays a crucial role. Over the years, significant progress has been made in discovering the cell wall-specific genomic resources related to drought tolerance at vegetative and reproductive stages of rice. However, questions remain about how the drought-induced changes in cell wall made by these genomic resources potentially influence the vegetative and reproductive development of rice. The possibly major candidate genes underlying the function of quantitative trait loci directly or indirectly associated with the cell wall plasticization-mediated drought tolerance of rice might have a huge promise in dissecting the putative genomic regions associated with cell wall plasticity under drought. Furthermore, engineering the drought tolerance of rice using cell wall-related genes from resurrection plants may have huge prospects for rice yield improvement. Here, we review the comprehensive multidisciplinary analyses to unravel different components and mechanisms involved in drought-induced cell wall plasticity at vegetative and reproductive stages that could be targeted for improving rice yield under drought.

Deciphering the dynamic gene expression patterns of pollen abortion in a male sterile line of Avena sativa through transcriptome analysis at different developmental stages

BACKGROUND

Male sterility (MS) has important applications in hybrid seed production, and the abortion of anthers has been observed in many plant species. While most studies have focused on the genetic factors affecting male sterility, the dynamic gene expression patterns of pollen abortion in male sterile lines have not been fully elucidated. In addition, there is still no hybrid oat that is commercially planted due to the lack of a suitable system of male sterility for hybrid breeding.

RESULTS

In this study, we cultivated a male sterile oat line and a near-isogenic line by crossbreeding to elucidate the expression patterns of genes that may be involved in sterility. The first reported CA male sterile (CAMS) oat line was used for cross-testing and hybridization experiments and was confirmed to exhibit a type of nuclear sterility controlled by recessive genes. Oat stamens of two lines were sampled at four different developmental stages separately. Paired-end RNA sequencing was performed for each sample and generated 252.84 Gb sequences. There were 295,462 unigenes annotated in public databases in all samples, and we compared the histological characteristics and transcriptomes of oat stamens from the two oat lines at different developmental stages. Our results demonstrate that the sterility of the male sterile oat line occurs in the early stage of stamen development and is primarily attributable to abnormal meiosis and the excessive accumulation of superoxide.

CONCLUSIONS

To the best of our knowledge, this study is the first to decipher the dynamic expression profiles of pollen abortion CAMS and CA male fertile (CAMF) oat lines, which may represent a valuable resource for further studies attempting to understand pollen abortion and anther development in oats.

Cytological and transcriptome analyses reveal OsPUB73 defect affects the gene expression associated with tapetum or pollen exine abnormality in rice

BACKGROUND

As one of the main crops in the world, sterility of rice (Oryza sativa L.) significantly affects the production and leads to yield decrease. Our previous research showed that OsPUB73, which encodes U-box domain-containing protein 73, may be associated with male sterility. However, little information is available on this gene that is required for anther development. In the present study, we knocked out OsPUB73 by using the CRISPR/Cas9 system and studied the cytological and transcriptome of the gene-defect associated with pollen development and sterility in the rice variety (Taichung 65).

RESULTS

The sequence analysis indicated that OsPUB73 was comprised of 3 exons and 2 introns, of which CDS encoded 586 amino acids including a U-box domain. The expression pattern of OsPUB73 showed that it was highly expressed in the anther during meiosis stage. The ospub73 displayed low pollen fertility (19.45%), which was significantly lower than wild type (WT) (85.37%). Cytological observation showed tapetum vacuolated at the meiosis stage and pollen exine was abnormal at the bi-cellular pollen stage of ospub73. RNA-seq analysis detected 2240 down and 571 up-regulated genes in anther of ospub73 compared with WT during meiosis stage. Among of 2240 down-regulated genes, seven known genes were associated with tapetal cell death or pollen exine development, including CYP703A3 (Cytochrome P450 Hydroxylase703A3), CYP704B2 (Cytochrome P450 Hydroxylase704B2), DPW (Defective Pollen Wall), PTC1 (Persistant Tapetal Cell1), UDT1 (Undeveloped Tapetum1), OsAP37 (Aspartic protease37) and OsABCG15 (ATP binding cassette G15), which were validated by quantitative real-time polymerase chain reaction (qRT-PCR). These results suggested OsPUB73 may play an important role in tapetal or pollen exine development and resulted in pollen partial sterility.

CONCLUSION

Our results revealed that OsPUB73 plays an important role in rice male reproductive development, which provides valuable information about the molecular mechanisms of the U-box in rice male reproductive development.

MeioCapture: an efficient method for staging and isolation of meiocytes in the prophase I sub-stages of meiosis in wheat

BACKGROUND

Molecular analysis of meiosis has been hindered by difficulties in isolating high purity subpopulations of sporogenous cells representing the succeeding stages of meiosis. Isolation of purified male meiocytes from defined meiotic stages is crucial in discovering meiosis specific genes and associated regulatory networks.

RESULTS

We describe an optimized method termed MeioCapture for simultaneous isolation of uncontaminated male meiocytes from wheat (Triticum spp.), specifically from the pre-meiotic G2 and the five sub-stages of meiotic prophase I. The MeioCapture protocol builds on the traditional anther squash technique and the capillary collection method, and involves extrusion of intact sporogenous archesporial columns (SACs) containing meiocytes. This improved method exploits the natural meiotic synchrony between anthers of the same floret, the correlation between the length of anthers and meiotic stage, and the occurrence of meiocytes in intact SACs largely free of somatic cells. The main advantage of MeioCapture, compared to previous methods, is that it allows simultaneous collection of meiocytes from different sub-stages of prophase I at a very high level of purity, through correlation of stages with anther sizes. A detailed description is provided for all steps, including the collection of tissue, isolation and size sorting of anthers, extrusion of intact SACs, and staging of meiocytes. Precautions for individual steps throughout the procedure are also provided to facilitate efficient isolation of pure meiocytes. The proof-of-concept was successfully established in wheat, and a light microscopic atlas of meiosis, encompassing all stages from pre-meiosis to telophase II, was developed.

CONCLUSION

The MeioCapture method provides an essential technique to study the molecular basis of chromosome pairing and exchange of genetic information in wheat, leading to strategies for manipulating meiotic recombination frequencies. The method also provides a foundation for similar studies in other crop species.

The coupling of transcriptome and proteome adaptation during development and heat stress response of tomato pollen

BACKGROUND

Pollen development is central for plant reproduction and is assisted by changes of the transcriptome and proteome. At the same time, pollen development and viability is largely sensitive to stress, particularly to elevated temperatures. The transcriptomic and proteomic changes during pollen development and of different stages in response to elevated temperature was targeted to define the underlying molecular principles.

RESULTS

The analysis of the transcriptome and proteome of Solanum lycopersicum pollen at tetrad, post-meiotic and mature stage before and after heat stress yielded a decline of the transcriptome but an increase of the proteome size throughout pollen development. Comparison of the transcriptome and proteome led to the discovery of two modes defined as direct and delayed translation. Here, genes of distinct functional processes are under the control of direct and delayed translation. The response of pollen to elevated temperature occurs rather at proteome, but not as drastic at the transcriptome level. Heat shock proteins, proteasome subunits, ribosomal proteins and eukaryotic initiation factors are most affected. On the example of heat shock proteins we demonstrate a decoupling of transcript and protein levels as well as a distinct regulation between the developmental stages.

CONCLUSIONS

The transcriptome and proteome of developing pollen undergo drastic changes in composition and quantity. Changes at the proteome level are a result of two modes assigned as direct and delayed translation. The response of pollen to elevated temperature is mainly regulated at the proteome level, whereby proteins related to synthesis and degradation of proteins are most responsive and might play a central role in the heat stress response of pollen.

Genome-wide analyses of late pollen-preferred genes conserved in various rice cultivars and functional identification of a gene involved in the key processes of late pollen development

BACKGROUND

Understanding late pollen development, including the maturation and pollination process, is a key component in maintaining crop yields. Transcriptome data obtained through microarray or RNA-seq technologies can provide useful insight into those developmental processes. Six series of microarray data from a public transcriptome database, the Gene Expression Omnibus of the National Center for Biotechnology Information, are related to anther and pollen development.

RESULTS

We performed a systematic and functional study across the rice genome of genes that are preferentially expressed in the late stages of pollen development, including maturation and germination. By comparing the transcriptomes of sporophytes and male gametes over time, we identified 627 late pollen-preferred genes that are conserved among japonica and indica rice cultivars. Functional classification analysis with a MapMan tool kit revealed a significant association between cell wall organization/metabolism and mature pollen grains. Comparative analysis of rice and Arabidopsis demonstrated that genes involved in cell wall modifications and the metabolism of major carbohydrates are unique to rice. We used the GUS reporter system to monitor the expression of eight of those genes. In addition, we evaluated the significance of our candidate genes, using T-DNA insertional mutant population and the CRISPR/Cas9 system. Mutants from T-DNA insertion and CRISPR/Cas9 systems of a rice gene encoding glycerophosphoryl diester phosphodiesterase are defective in their male gamete transfer.

CONCLUSION

Through the global analyses of the late pollen-preferred genes from rice, we found several biological features of these genes. First, biological process related to cell wall organization and modification is over-represented in these genes to support rapid tube growth. Second, comparative analysis of late pollen preferred genes between rice and Arabidopsis provide a significant insight on the evolutional disparateness in cell wall biogenesis and storage reserves of pollen. In addition, these candidates might be useful targets for future examinations of late pollen development, and will be a valuable resource for accelerating the understanding of molecular mechanisms for pollen maturation and germination processes in rice.

Transcriptional profiling of cork oak phellogenic cells isolated by laser microdissection

The phenylpropanoid pathway impacts the cork quality development. In cork of bad quality, the flavonoid route is favored, whereas in good quality, cork lignin and suberin production prevails. Cork oaks develop a thick cork tissue as a protective shield that results of the continuous activity of a secondary meristem, the cork cambium, or phellogen. Most studies applied to developmental processes do not consider the cell types from which the samples were extracted. Here, laser microdissection (LM) coupled with transcript profiling using RNA sequencing (454 pyrosequencing) was applied to phellogen cells of trees producing low- and good quality cork. Functional annotation and functional enrichment analyses showed that stress-related genes are enriched in samples extracted from trees producing good quality cork (GQC). This process is under tight transcriptional (transcription factors, kinases) regulation and also hormonal control involving ABA, ethylene, and auxins. The phellogen cells collected from trees producing bad quality cork (BQC) show a consistent up-regulation of genes belonging to the flavonoid pathway as a response to stress. They also display a different modulation of cell wall genes resulting into a thinner cork layer, i.e., less meristematic activity. Based on the analysis of the phenylpropanoid pathway regulating genes, in GQC, the synthesis of lignin and suberin is promoted, whereas in BQC, the same pathway favors the biosynthesis of free phenolic compounds. This study provided new insights of how cell-specific gene expression can determine tissue and organ morphology and physiology and identified robust candidate genes that can be used in breeding programs aiming at improving cork quality.

Application of rice microspore-preferred promoters to manipulate early pollen development in Arabidopsis: a heterologous system

Rice microspore-promoters. Based on microarray data analyzed for developing anthers and pollen grains, we identified nine rice microspore-preferred (RMP) genes, designated RMP1 through RMP9. To extend their biotechnological applicability, we then investigated the activity of RMP promoters originating from monocotyledonous rice in a heterologous system of dicotyledonous Arabidopsis. Expression of GUS was significantly induced in transgenic plants from the microspore to the mature pollen stages and was driven by the RMP1, RMP3, RMP4, RMP5, and RMP9 promoters. We found it interesting that, whereas RMP2 and RMP6 directed GUS expression in microspore at the early unicellular and bicellular stages, RMP7 and RMP8 seemed to be expressed at the late tricellular and mature pollen stages. Moreover, GUS was expressed in seven promoters, RMP3 through RMP9, during the seedling stage, in immature leaves, cotyledons, and roots. To confirm microspore-specific expression, we used complementation analysis with an Arabidopsis male-specific gametophytic mutant, sidecar pollen-2 (scp-2), to verify the activity of three promoters. That mutant shows defects in microspore development prior to pollen mitosis I. These results provide strong evidence that the SIDECAR POLLEN gene, driven by RMP promoters, successfully complements the scp-2 mutation, and they strongly suggest that these promoters can potentially be applied for manipulating the expression of target genes at the microspore stage in various species.

Genome-wide identification and analysis of rice genes preferentially expressed in pollen at an early developmental stage

Microspore production using endogenous developmental programs has not been well studied. The main limitation is the difficulty in identifying genes preferentially expressed in pollen grains at early stages. To overcome this limitation, we collected transcriptome data from anthers and microspore/pollen and performed meta-expression analysis. Subsequently, we identified 410 genes showing preferential expression patterns in early developing pollen samples of both japonica and indica cultivars. The expression patterns of these genes are distinguishable from genes showing pollen mother cell or tapetum-preferred expression patterns. Gene Ontology enrichment and MapMan analyses indicated that microspores in rice are closely linked with protein degradation, nucleotide metabolism, and DNA biosynthesis and regulation, while the pollen mother cell or tapetum are strongly associated with cell wall metabolism, lipid metabolism, secondary metabolism, and RNA biosynthesis and regulation. We also generated transgenic lines under the control of the promoters of eight microspore-preferred genes and confirmed the preferred expression patterns in plants using the GUS reporting system. Furthermore, cis-regulatory element analysis revealed that pollen specific elements such as POLLEN1LELAT52, and 5659BOXLELAT5659 were commonly identified in the promoter regions of eight rice genes with more frequency than estimation. Our study will provide new sights on early pollen development in rice, a model crop plant.

Characterization of two PEBP genes, SrFT and SrMFT, in thermogenic skunk cabbage (Symplocarpus renifolius)

Floral thermogenesis has been found in dozens of primitive seed plants and the reproductive organs in these plants produce heat during anthesis. Thus, characterization of the molecular mechanisms underlying flowering is required to fully understand the role of thermogenesis, but this aspect of thermogenic plant development is largely unknown. In this study, extensive database searches and cloning experiments suggest that thermogenic skunk cabbage (Symplocarpus renifolius), which is a member of the family Araceae, possesses two genes encoding phosphatidyl ethanolamine-binding proteins (PEBP), FLOWERING LOCUS T (SrFT) and MOTHER OF FT AND TFL1 (SrMFT). Functional analyses of SrFT and SrMFT in Arabidopsis indicate that SrFT promotes flowering, whereas SrMFT does not. In S. renifolius, the stage- and tissue-specific expression of SrFT was more evident than that of SrMFT. SrFT was highly expressed in flowers and leaves and was mainly localized in fibrovascular tissues. In addition, microarray analysis revealed that, within floral tissues, SrFT was co-regulated with the genes associated with cellular respiration and mitochondrial function, including ALTERNATIVE OXIDASE gene proposed to play a major role in floral thermogenesis. Taken together, these data suggest that, among the PEBP genes, SrFT plays a role in flowering and floral development in the thermogenic skunk cabbage.

Bioinformatics resources for pollen

Bioinformatics for Pollen. Pollen plays a key role in crop production, and its development is the most delicate phase in reproduction. Different metabolic pathways are involved in pollen development, and changes in the level of some metabolites, as well as responses to stress, are correlated with the reduction in pollen viability, leading consequently to a decrease in the fruit production. However, studies on pollen may be hard because gamete development and fertilization are complex processes that occur during a short window of time. The rise of the so-called -omics sciences provided key strategies to promote molecular research in pollen tissues, starting from model organisms and moving to increasing number of species. An integrated multi-level approach based on investigations from genomics, transcriptomics, proteomics and metabolomics appears now feasible to clarify key molecular processes in pollen development and viability. To this aim, bioinformatics has a fundamental role for data production and analysis, contributing varied and ad hoc methodologies, endowed with different sensitivity and specificity, necessary for extracting added-value information from the large amount of molecular data achievable. Bioinformatics is also essential for data management, organization, distribution and integration in suitable resources. This is necessary to catch the biological features of the pollen tissues and to design effective approaches to identifying structural or functional properties, enabling the modeling of the major involved processes in normal or in stress conditions. In this review, we provide an overview of the available bioinformatics resources for pollen, ranging from raw data collections to complete databases or platforms, when available, which include data and/or results from -omics efforts on the male gametophyte. Perspectives in the fields will also be described.

Global Dynamic Transcriptome Programming of Rapeseed (Brassica napus L.) Anther at Different Development Stages

Rapeseed (Brassica napus L.) is an important oil crop worldwide and exhibits significant heterosis. Effective pollination control systems, which are closely linked to anther development, are a prerequisite for utilizing heterosis. The anther, which is the male organ in flowering plants, undergoes many metabolic processes during development. Although the gene expression patterns underlying pollen development are well studied in model plant Arabidopsis, the regulatory networks of genome-wide gene expression during rapeseed anther development is poorly understood, especially regarding metabolic regulations. In this study, we systematically analyzed metabolic processes occurring during anther development in rapeseed using ultrastructural observation and global transcriptome analysis. Anther ultrastructure exhibited that numerous cellular organelles abundant with metabolic materials, such as elaioplast, tapetosomes, plastids (containing starch deposits) etc. appeared, accompanied with anther structural alterations during anther development, suggesting many metabolic processes occurring. Global transcriptome analysis revealed dynamic changes in gene expression during anther development that corresponded to dynamic functional alterations between early and late anther developmental stages. The early stage anthers preferentially expressed genes involved in lipid metabolism that are related to pollen extine formation as well as elaioplast and tapetosome biosynthesis, whereas the late stage anthers expressed genes associated with carbohydrate metabolism to form pollen intine and to accumulate starch in mature pollen grains. Finally, a predictive gene regulatory module responsible for early pollen extine formation was generated. Taken together, this analysis provides a comprehensive understanding of dynamic gene expression programming of metabolic processes in the rapeseed anther, especially with respect to lipid and carbohydrate metabolism during pollen development.

Cytoplasmic male sterility (CMS) in hybrid breeding in field crops

A comprehensive understanding of CMS/Rf system enabled by modern omics tools and technologies considerably improves our ability to harness hybrid technology for enhancing the productivity of field crops. Harnessing hybrid vigor or heterosis is a promising approach to tackle the current challenge of sustaining enhanced yield gains of field crops. In the context, cytoplasmic male sterility (CMS) owing to its heritable nature to manifest non-functional male gametophyte remains a cost-effective system to promote efficient hybrid seed production. The phenomenon of CMS stems from a complex interplay between maternally-inherited (mitochondrion) and bi-parental (nucleus) genomic elements. In recent years, attempts aimed to comprehend the sterility-inducing factors (orfs) and corresponding fertility determinants (Rf) in plants have greatly increased our access to candidate genomic segments and the cloned genes. To this end, novel insights obtained by applying state-of-the-art omics platforms have substantially enriched our understanding of cytoplasmic-nuclear communication. Concomitantly, molecular tools including DNA markers have been implicated in crop hybrid breeding in order to greatly expedite the progress. Here, we review the status of diverse sterility-inducing cytoplasms and associated Rf factors reported across different field crops along with exploring opportunities for integrating modern omics tools with CMS-based hybrid breeding.

Pre-Meiotic Anther Development: Cell Fate Specification and Differentiation

Research into anther ontogeny has been an active and developing field, transitioning from a strictly lineage-based view of cellular differentiation events to a more complex understanding of cell fate specification. Here we describe the modern interpretation of pre-meiotic anther development, from the earliest cell specifications within the anther lobes through SPL/NZZ-, MSP1-, and MEL1-dependent pathways as well as the initial setup of the abaxial and adaxial axes and outgrowth of the anther lobes. We then continue with a look at the known information regarding further differentiation of the somatic layers of the anther (the epidermis, endothecium, middle layer, and tapetum), with an emphasis on male-sterile mutants identified as defective in somatic cell specification. We also describe the differences in developmental stages among species and use this information to discuss molecular studies that have analyzed transcriptome, proteome, and small-RNA information in the anther.

Identification of cis-elements and evaluation of upstream regulatory region of a rice anther-specific gene, OSIPP3, conferring pollen-specific expression in Oryza sativa (L.) ssp. indica

Pollen-specific expression. Promoters comprise of various cis-regulatory elements which control development and physiology of plants by regulating gene expression. To understand the promoter specificity and also identification of functional cis-acting elements, progressive 5' deletion analysis of the promoter fragments is widely used. We have evaluated the activity of regulatory elements of 5' promoter deletion sequences of anther-specific gene OSIPP3, viz. OSIPP3-∆1 (1504 bp), OSIPP3-∆2 (968 bp), OSIPP3-∆3 (388 bp) and OSIPP3-∆4 (286 bp) through the expression of transgene GUS in rice. In silico analysis of 1504-bp sequence harboring different copy number of cis-acting regulatory elements such as POLLENLELAT52, GTGANTG10, enhancer element of LAT52 and LAT56 indicated that they were essential for high level of expression in pollen. Histochemical GUS analysis of the transgenic plants revealed that 1504- and 968-bp fragments directed GUS expression in roots and anthers, while the 388- and 286-bp fragments restricted the GUS expression to only pollen, of which 388 bp conferred strong GUS expression. Further, GUS staining analysis of different panicle development stages (P1-P6) confirmed that the GUS gene was preferentially expressed only at P6 stage (late pollen stage). The qRT-PCR analysis of GUS transcript revealed 23-fold higher expression of GUS transcript in OSIPP3-Δ1 followed by OSIPP3-Δ2 (eightfold) and OSIPP3-Δ3 (threefold) when compared to OSIPP3-Δ4. Based on our results, we proposed that among the two smaller fragments, the 388-bp upstream regulatory region could be considered as a promising candidate for pollen-specific expression of agronomically important transgenes in rice.

"Omics" of maize stress response for sustainable food production: opportunities and challenges

Maize originated in the highlands of Mexico approximately 8700 years ago and is one of the most commonly grown cereal crops worldwide, followed by wheat and rice. Abiotic stresses (primarily drought, salinity, and high and low temperatures), together with biotic stresses (primarily fungi, viruses, and pests), negatively affect maize growth, development, and eventually production. To understand the response of maize to abiotic and biotic stresses and its mechanism of stress tolerance, high-throughput omics approaches have been used in maize stress studies. Integrated omics approaches are crucial for dissecting the temporal and spatial system-level changes that occur in maize under various stresses. In this comprehensive analysis, we review the primary types of stresses that threaten sustainable maize production; underscore the recent advances in maize stress omics, especially proteomics; and discuss the opportunities, challenges, and future directions of maize stress omics, with a view to sustainable food production. The knowledge gained from studying maize stress omics is instrumental for improving maize to cope with various stresses and to meet the food demands of the exponentially growing global population. Omics systems science offers actionable potential solutions for sustainable food production, and we present maize as a notable case study.

Flowers and mycorrhizal roots--closer than we think?

Roots and flowers are formed at the extreme ends of plants and they differ in almost every aspect of their development and function; even so, they exhibit surprising molecular commonalities. For example, the calcium and calmodulin-dependent protein kinase (CCaMK) plays a central role in root symbioses with fungi and bacteria, but is also highly expressed in developing anthers. Moreover, independent evidence from transcriptomics, phylogenomics, and genetics reveals common developmental elements in root symbioses and reproductive development. We discuss the significance of these overlaps, and we argue that an integrated comparative view of the two phenomena will stimulate research and provide new insight, not only into shared components, but also into the specific aspects of anther development and root symbioses.

A decade of pollen transcriptomics

KEY MESSAGE

Overview of pollen transcriptome studies. Pollen development is driven by gene expression, and knowledge of the molecular events underlying this process has undergone a quantum leap in the last decade through studies of the transcriptome. Here, we outline historical evidence for male haploid gene expression and review the wealth of pollen transcriptome data now available. Knowledge of the transcriptional capacity of pollen has progressed from genetic studies to the direct analysis of RNA and from gene-by-gene studies to analyses on a genomic scale. Microarray and/or RNA-seq data can now be accessed for all phases and cell types of developing pollen encompassing 10 different angiosperms. These growing resources have accelerated research and will undoubtedly inspire new directions and the application of system-based research into the mechanisms that govern the development, function and evolution of angiosperm pollen.

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.

The mature anther-preferentially expressed genes are associated with pollen fertility, pollen germination and anther dehiscence in rice

Background

The anthers and pollen grains are critical for male fertility and hybrid rice breeding. The development of rice mature anther and pollen consists of multiple continuous stages. However, molecular mechanisms regulating mature anther development were poorly understood.

Results

In this study, we have identified 291 mature anther-preferentially expressed genes (OsSTA) in rice based on Affymetrix microarray data. Gene Ontology (GO) analysis indicated that OsSTA genes mainly participated in metabolic and cellular processes that are likely important for rice anther and pollen development. The expression patterns of OsSTA genes were validated using real-time PCR and mRNA in situ hybridizations. Cis-element identification showed that most of the OsSTA genes had the cis-elements responsive to phytohormone regulation. Co-expression analysis of OsSTA genes showed that genes annotated with pectinesterase and calcium ion binding activities were rich in the network, suggesting that OsSTA genes could be involved in pollen germination and anther dehiscence. Furthermore, OsSTA RNAi transgenic lines showed male-sterility and pollen germination defects.

Conclusions

The results suggested that OsSTA genes function in rice male fertility, pollen germination and anther dehiscence and established molecular regulating networks that lay the foundation for further functional studies.

Electronic supplementary material

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

The meiotic transcriptome architecture of plants

Although a number of genes that play key roles during the meiotic process have been characterized in great detail, the whole process of meiosis is still not completely unraveled. To gain insight into the bigger picture, large-scale approaches like RNA-seq and microarray can help to elucidate the transcriptome landscape during plant meiosis, discover co-regulated genes, enriched processes, and highly expressed known and unknown genes which might be important for meiosis. These high-throughput studies are gaining more and more popularity, but their beginnings in plant systems reach back as far as the 1960's. Frequently, whole anthers or post-meiotic pollen were investigated, while less data is available on isolated cells during meiosis, and only few studies addressed the transcriptome of female meiosis. For this review, we compiled meiotic transcriptome studies covering different plant species, and summarized and compared their key findings. Besides pointing to consistent as well as unique discoveries, we finally draw conclusions what can be learned from these studies so far and what should be addressed next.

Polyamine oxidase 7 is a terminal catabolism-type enzyme in Oryza sativa and is specifically expressed in anthers

Polyamine oxidase (PAO), which requires FAD as a cofactor, functions in polyamine catabolism. Plant PAOs are classified into two groups based on their reaction modes. The terminal catabolism (TC) reaction always produces 1,3-diaminopropane (DAP), H2O2, and the respective aldehydes, while the back-conversion (BC) reaction produces spermidine (Spd) from tetraamines, spermine (Spm) and thermospermine (T-Spm) and/or putrescine from Spd, along with 3-aminopropanal and H2O2. The Oryza sativa genome contains seven PAO-encoded genes termed OsPAO1-OsPAO7. To date, we have characterized four OsPAO genes. The products of these genes, i.e. OsPAO1, OsPAO3, OsPAO4 and OsPAO5, catalyze BC-type reactions. Whereas OsPAO1 remains in the cytoplasm, the other three PAOs localize to peroxisomes. Here, we examined OsPAO7 and its gene product. OsPAO7 shows high identity to maize ZmPAO1, the best characterized plant PAO having TC-type activity. OsPAO7 seems to remain in a peripheral layer of the plant cell with the aid of its predicted signal peptide and transmembrane domain. Recombinant OsPAO7 prefers Spm and Spd as substrates, and it produces DAP from both substrates in a time-dependent manner, indicating that OsPAO7 is the first TC-type enzyme identified in O. sativa. The results clearly show that two types of PAOs co-exist in O. sativa. Furthermore, OsPAO7 is specifically expressed in anthers, with an expressional peak at the bicellular pollen stage. The physiological function of OsPAO7 in anthers is discussed.

Laser microdissection of plant cells

Different plant cell types express unique transcriptomes, proteomes, and metabolomes. Therefore, the isolation of specific cell types prior to molecular analyses is important to understand the specification, differentiation, and function of these cells. Isolation of specific plant cell types from composite organs can be achieved by laser microdissection (LMD). A wide variety of methods to fix and embed tissues prior to LMD and downstream molecular analyses have been developed for different plant species and tissues. The present review summarizes and highlights the most recently applied LMD approaches in plant science.

Comparative transcriptomics of early meiosis in Arabidopsis and maize

Though sexually reproductive plants share the same principle and most processes in meiosis, there are distinct features detectable. To address the similarities and differences of early meiosis transcriptomes from the dicot model system Arabidopsis and monocot model system maize, we performed comparative analyses of RNA-seq data of isolated meiocytes, anthers and seedlings from both species separately and via orthologous genes. Overall gene expression showed similarities, such as an increased number of reads mapping to unannotated features, and differences, such as the amount of differentially expressed genes. We detected major similarities and differences in functional annotations of genes up-regulated in meiocytes, which point to conserved features as well as unique features. Transcriptional regulation seems to be quite similar in Arabidopsis and maize, and we could reveal known and novel transcription factors and cis-regulatory elements acting in early meiosis. Taken together, meiosis between Arabidopsis and maize is conserved in many ways, but displays key distinctions that lie in the patterns of gene expression.

Roles of pollen-specific boron efflux transporter, OsBOR4, in the rice fertilization process

Arabidopsis thaliana BOR1 was the first boron (B) transporter identified in living systems. There are four AtBOR1-like genes, OsBOR1, 2, 3 and 4, present in the rice genome. We characterized the activity, expression and physiological function of OsBOR4. OsBOR4 is an active efflux transporter of B. Quantitative PCR analysis and OsBOR4 promoter-green fluorescent protein (GFP) fusion revealed that OsBOR4 was both highly and specifically expressed in pollen. We obtained five Tos17 insertion mutants of osbor4. The pollen grains were viable and development of floral organs was normal in the homozygous osbor4 mutants. We observed that in all Tos17 insertion lines tested, the frequency of osbor4 homozygous plants was lower than expected in the progeny of self-fertilized heterozygous plants. These results establish that OsBOR4 is essential for normal reproductive processes. Pollen from osbor4 homozygous plants elongated fewer tubes on wild-type stigmas, and tube elongation of mutant pollen was less efficient compared with the wild-type pollen, suggesting reduced competence of osbor4 mutant pollen. The reduced competence of mutant pollen was further supported by the crosses of independent Tos17-inserted alleles of OsBOR4. Our results suggest that OsBOR4, a boron efflux transporter, is required for normal pollen germination and/or tube elongation.

Cell type-specific transcriptome of Brassicaceae stigmatic papilla cells from a combination of laser microdissection and RNA sequencing

Pollination is an early and critical step in plant reproduction, leading to successful fertilization. It consists of many sequential processes, including adhesion of pollen grains onto the surface of stigmatic papilla cells, foot formation to strengthen pollen-stigma interaction, pollen hydration and germination, and pollen tube elongation and penetration. We have focused on an examination of the expressed genes in papilla cells, to increase understanding of the molecular systems of pollination. From three representative species of Brassicaceae (Arabidopsis thaliana, A. halleri and Brassica rapa), stigmatic papilla cells were isolated precisely by laser microdissection, and cell type-specific gene expression in papilla cells was determined by RNA sequencing. As a result, 17,240, 19,260 and 21,026 unigenes were defined in papilla cells of A. thaliana, A. halleri and B. rapa, respectively, and, among these, 12,311 genes were common to all three species. Among the17,240 genes predicted in A. thaliana, one-third were papilla specific while approximately half of the genes were detected in all tissues examined. Bioinformatics analysis revealed that genes related to a wide range of reproduction and development functions are expressed in papilla cells, particularly metabolism, transcription and membrane-mediated information exchange. These results reflect the conserved features of general cellular function and also the specific reproductive role of papilla cells, highlighting a complex cellular system regulated by a diverse range of molecules in these cells. This study provides fundamental biological knowledge to dissect the molecular mechanisms of pollination in papilla cells and will shed light on our understanding of plant reproduction mechanisms.

Rice male development under drought stress: phenotypic changes and stage-dependent transcriptomic reprogramming

Drought affects rice reproduction and results in severe yield loss. The developmental defects and changes of gene regulation network in reproductive tissues under drought stress are largely unknown. In this study, rice plants subjected to reproductive stage drought stress were examined for floral development and transcriptomic changes. The results showed that male fertility was dramatically affected, with differing pollen viability in flowers of the same panicle due to aberrant anther development under water stress. Examination of local starch distribution revealed that starch accumulated abnormally in terms of position and abundance in anthers of water-stressed plants. Microarray analysis using florets of different sizes identified >1000 drought-responsive genes, most of which were specifically regulated in only one or two particular sizes of florets, suggesting developmental stage-dependent responses to drought. Genes known to be involved in tapetum and/or microspore development, cell wall formation or expansion, and starch synthesis were found more frequently among the genes affected by drought than genome average, while meiosis and MADS-box genes were less frequently affected. In addition, pathways related to gibberellin acid signaling and abscisic acid catabolism were reprogrammed by drought. Our results strongly suggest interactions between reproductive development, phytohormone signaling, and carbohydrate metabolism in water-stressed plants.

Collapsed abnormal pollen1 gene encoding the Arabinokinase-like protein is involved in pollen development in rice

We isolated a pollen-defective mutant, collapsed abnormal pollen1 (cap1), from Tos17 insertional mutant lines of rice (Oryza sativa). The cap1 heterozygous plant produced equal numbers of normal and collapsed abnormal grains. The abnormal pollen grains lacked almost all cytoplasmic materials, nuclei, and intine cell walls and did not germinate. Genetic analysis of crosses revealed that the cap1 mutation did not affect female reproduction or vegetative growth. CAP1 encodes a protein consisting of 996 amino acids that showed high similarity to Arabidopsis (Arabidopsis thaliana) l-arabinokinase, which catalyzes the conversion of l-arabinose to l-arabinose 1-phosphate. A wild-type genomic DNA segment containing CAP1 restored mutants to normal pollen grains. During rice pollen development, CAP1 was preferentially expressed in anthers at the bicellular pollen stage, and the effects of the cap1 mutation were mainly detected at this stage. Based on the metabolic pathway of l-arabinose, cap1 pollen phenotype may have been caused by toxic accumulation of l-arabinose or by inhibition of cell wall metabolism due to the lack of UDP-l-arabinose derived from l-arabinose 1-phosphate. The expression pattern of CAP1 was very similar to that of another Arabidopsis homolog that showed 71% amino acid identity with CAP1. Our results suggested that CAP1 and related genes are critical for pollen development in both monocotyledonous and dicotyledonous plants.

Transcriptome analysis of developing ovules in rice isolated by laser microdissection

Comprehensive genome-wide gene expression profiles during plant male gametogenesis have been thoroughly analyzed over the last decade. In contrast, gene expression profiles during female gametogenesis have been studied relatively little, and our knowledge concerning plant female gametogenesis is limited. We determined the genome-wide gene expression profiles of developing ovules containing female gametophytes from the megaspore mother cell at the pre-meiotic stage to the mature embryo sac in rice (Oryza sativa) using microarrays. In order to separate ovules from scutellum, we used a laser microdissection (LM) technique. Dynamic gene expression was revealed in developing ovules, and a major transition of the transcriptome was observed between middle and late meiotic stages, where many genes were down-regulated >10-fold. Many potential players in female gametogenesis, that showed dynamic or enriched expression, were highlighted. We identified the temporal and dramatic up-regulation of a subset of transposable elements during female meiotic stages that were not observed in males. Transcription factor genes enriched in developing ovules were also uncovered, which may play crucial roles during female gametogenesis. This is the first report of comprehensive genome-wide gene expression profiles during female gametogenesis useful for plant reproductive studies. Combined with additional experiments, our data may provide important clues to understand female gametogenesis in plants.

Molecular control of male reproductive development and pollen fertility in rice

Anther development and male fertility are essential biological processes for flowering plants and are important for crop seed production. Genetic manipulation of male fertility/sterility is critical for crop hybrid breeding. Rice (Oryza sativa L.) male sterility phenotypes, including genic male sterility, hybrid male sterility, and cytoplasmic male sterility, are generally caused by mutations of fertility-related genes, by incompatible interactions between divergent allelic or non-allelic genes, or by genetic incompatibilities between cytoplasmic and nuclear genomes. Here, we review the recent advances in the molecular basis of anther development and male fertility-sterility conversion in specific genetic backgrounds, and the interactions with certain environmental factors. The highlighted findings in this review have significant implications in both basic studies and rice genetic improvement.

Transcriptome profiling analysis reveals that flavonoid and ascorbate-glutathione cycle are important during anther development in Upland cotton

BACKGROUND

Previous transcriptome profiling studies have investigated the molecular mechanisms of pollen and anther development, and identified many genes involved in these processes. However, only 51 anther ESTs of Upland cotton (Gossypium hirsutum) were found in NCBI and there have been no reports of transcriptome profiling analyzing anther development in Upland cotton, a major fiber crop in the word.

METHODOLOGY/PRINCIPAL FINDING

Ninety-eight hundred and ninety-six high quality ESTs were sequenced from their 3'-ends and assembled into 6,643 unigenes from a normalized, full-length anther cDNA library of Upland cotton. Combined with previous sequenced anther-related ESTs, 12,244 unigenes were generated as the reference genes for digital gene expression (DGE) analysis. The DGE was conducted on anthers that were isolated at tetrad pollen (TTP), uninucleate pollen (UNP), binucleate pollen (BNP) and mature pollen (MTP) periods along with four other tissues, i.e., roots (RO), stems (ST), leaves (LV) and embryos (EB). Through transcriptome profiling analysis, we identified 1,165 genes that were enriched at certain anther development periods, and many of them were involved in starch and sucrose metabolism, pentose and glucuronate interconversion, flavonoid biosynthesis, and ascorbate and aldarate metabolism.

CONCLUSIONS/SIGNIFICANCE

We first generated a normalized, full-length cDNA library from anthers and performed transcriptome profiling analysis of anther development in Upland cotton. From these results, 10,178 anther expressed genes were identified, among which 1,165 genes were stage-enriched in anthers. And many of these stage-enriched genes were involved in some important processes regulating anther development.

Transcriptome profiling analysis using rice oligonucleotide microarrays

Transcriptome analysis using oligonucleotide microarrays is a powerful tool for detecting changes in genome-wide transcripts under a given biological condition. Although the rice genome sequence is available, the number of functionally characterized genes in rice is still very limited. Genome-wide transcriptome analysis is a useful tool for elucidating the functions of rice genes that have not yet been determined. Currently, more than 3,000 arrays are publicly available. Here, we introduce methods for genome-wide transcriptome analysis in rice.

A novel approach to dissect the abscission process in Arabidopsis

Abscission is the consequence of a specialized layer of cells undergoing a complex series of molecular and biochemical events. Analysis of the specific molecular changes associated with abscission is hampered by contamination from neighboring nonseparating tissues. Moreover, studies of abscission frequently involve the examination of events that take place in isolated segments of tissue exposed to nonphysiological concentrations of ethylene or indole-3-acetic acid for protracted periods (more than 24 h) of time. To resolve these problems, we have adopted the use of a transgenic line of Arabidopsis (Arabidopsis thaliana) where the promoter of an abscission-specific polygalacturonase gene (At2g41850/ARABIDOPSIS DEHISCENCE ZONE POLYGALACTURONASE2) has been fused to a green fluorescent protein reporter. RNA was extracted from green fluorescent protein-tagged cells, released from abscising floral organs, and used to generate a complementary DNA library. This library was used to probe a microarray, and a population of abscission-related transcripts was studied in detail. Seven novel abscission-related genes were identified, four of which encode proteins of unknown function. Reverse transcription-polymerase chain reaction analyses and promoter fusions to the β-glucuronidase reporter gene confirmed the expression of these genes in the abscission zone and revealed other places of expression during seedling development. Three of these genes were studied further by crossing reporter lines to the abscission mutants inflorescence deficient in abscission (ida) and blade-on-petiole1 (bop1)/bop2 and an IDA-overexpressing line. Phenotypic analysis of an At3g14380 transfer DNA insertion line indicates that this gene plays a functional role in floral organ shedding. This strategy has enabled us to uncover new genes involved in abscission, and their possible contribution to the process is discussed.

Isolation and gene expression analysis of a papain-type cysteine protease in thermogenic skunk cabbage (Symplocarpus renifolius)

Skunk cabbage (Symplocarpus renifolius) spadices contain abundant transcripts for cysteine protease (CP). From thermogenic spadices, we isolated SrCPA, a highly expressed CP gene that encoded a papain-type CP. SrCPA is structurally similar to other plant CPs, including the senescence-associated CPs found in aroids. The expression of SrCPA increased during floral development, and was observed in all floral tissues except for the stamens.

Spatial and temporal activity of upstream regulatory regions of rice anther-specific genes in transgenic rice and Arabidopsis

Upstream regulatory regions (URRs) of rice anther-specific genes, namely OSbHLH (coding for basic helix-loop-helix-containing protein) and OSFbox (F-box protein encoding gene), selected from the microarray data have been cloned to control expression of GUS and GFP reporter genes in stably transformed rice. Quantitative real time PCR analysis shows maximum transcript accumulation of these two genes in the meiotic anthers. Analysis of OSbHLH and OSFbox URRs by PLACE database reveal the presence of known pollen-specific cis elements. The URRs of both OSbHLH and OSFbox genes have maximum activity in the meiotic anther stage in rice, but confer constitutive expression in the heterologous dicot system, Arabidopsis, indicative of monocot specificity. Another rice gene (OSIPK; with homology to genes encoding calcium-dependent protein kinases) URR already reported to have anther-specific activity in Arabidopsis and tobacco also confers anther-specific expression in rice and is active in the pollen tubes, suggesting it belongs to the category of late expressed genes. The spatial activity of three URRs has also been analysed by histochemical evaluation of GUS activity in different anther cells/tissues. The activity of OSIPK URR in rice is strongest among the three URRs.

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.

Analysis of plant germline development by high-throughput RNA profiling: technical advances and new insights

Reproduction is a crucial step in the life cycle of plants. The male and female germline lineages develop in the reproductive organs of the flower, which in higher plants are the anthers and ovules, respectively. Development of the germline lineage initiates from a dedicated sporophytic cell that undergoes meiosis to form spores that subsequently give rise to the gametophytes through mitotic cell divisions. The mature male and female gametophytes harbour the male (sperm cells) and female gametes (egg and central cell), respectively. Those unite during double fertilization to initiate embryo and endosperm development in sexually reproducing higher plants. While cytological changes involved in development of the germline lineages have been well characterized in a number of species, investigation of the transcriptional basis underlying their development and the specification of the gametes proved challenging. This is largely due to the inaccessibility of the cells constituting the germline lineages, which are enclosed by sporophytic tissues. Only recently, these technical limitations could be overcome by combining new methods to isolate the relevant cells with powerful transcriptional profiling methods, such as microarrays or high-throughput sequencing of RNA. This review focuses on these technical advances and the new insights gained from them concerning the transcriptional basis and molecular mechanisms underlying germline 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.

The gene expression landscape of thermogenic skunk cabbage suggests critical roles for mitochondrial and vacuolar metabolic pathways in the regulation of thermogenesis

Floral thermogenesis has been described in several plant species. Because of the lack of comprehensive gene expression profiles in thermogenic plants, the molecular mechanisms by which floral thermogenesis is regulated remain to be established. We examined the gene expression landscape of skunk cabbage (Symplocarpus renifolius) during thermogenic and post-thermogenic stages and identified expressed sequence tags from different developmental stages of the inflorescences using super serial analysis of gene expression (SuperSAGE). In-depth analysis suggested that cellular respiration and mitochondrial functions are significantly enhanced during the thermogenic stage. In contrast, genes involved in stress responses and protein degradation were significantly up-regulated during post-thermogenic stages. Quantitative comparisons indicated that the expression levels of genes involved in cellular respiration were higher in thermogenic spadices than in Arabidopsis inflorescences. Thermogenesis-associated genes seemed to be expressed abundantly in the peripheral tissues of the spadix. Our results suggest that cellular respiration and mitochondrial metabolism play key roles in heat production during floral thermogenesis. On the other hand, vacuolar cysteine protease and other degradative enzymes seem to accelerate senescence and terminate thermogenesis in the post-thermogenic stage.

Comparative transcriptome analysis of stylar canal cells identifies novel candidate genes implicated in the self-incompatibility response of Citrus clementina

BACKGROUND

Reproductive biology in citrus is still poorly understood. Although in recent years several efforts have been made to study pollen-pistil interaction and self-incompatibility, little information is available about the molecular mechanisms regulating these processes. Here we report the identification of candidate genes involved in pollen-pistil interaction and self-incompatibility in clementine (Citrus clementina Hort. ex Tan.). These genes have been identified comparing the transcriptomes of laser-microdissected stylar canal cells (SCC) isolated from two genotypes differing for self-incompatibility response ('Comune', a self-incompatible cultivar and 'Monreal', a self- compatible mutation of 'Comune').

RESULTS

The transcriptome profiling of SCC indicated that the differential regulation of few specific, mostly uncharacterized transcripts is associated with the breakdown of self-incompatibility in 'Monreal'. Among them, a novel F-box gene showed a drastic up-regulation both in laser microdissected stylar canal cells and in self-pollinated whole styles with stigmas of 'Comune' in concomitance with the arrest of pollen tube growth. Moreover, we identify a non-characterized gene family as closely associated to the self-incompatibility genetic program activated in 'Comune'. Three different aspartic-acid rich (Asp-rich) protein genes, located in tandem in the clementine genome, were over-represented in the transcriptome of 'Comune'. These genes are tightly linked to a DELLA gene, previously found to be up-regulated in the self-incompatible genotype during pollen-pistil interaction.

CONCLUSION

The highly specific transcriptome survey of the stylar canal cells identified novel genes which have not been previously associated with self-pollen rejection in citrus and in other plant species. Bioinformatic and transcriptional analyses suggested that the mutation leading to self-compatibility in 'Monreal' affected the expression of non-homologous genes located in a restricted genome region. Also, we hypothesize that the Asp-rich protein genes may act as Ca2+ "entrapping" proteins, potentially regulating Ca2+ homeostasis during self-pollen recognition.

Molecular genetics, physiology and biology of self-incompatibility in Brassicaceae

Self-incompatibility (SI) is defined as the inability to produce zygotes after self-pollination in a fertile hermaphrodite plant, which has stamens and pistils in the same flower. This structural organization of the hermaphrodite flower increases the risk of self-pollination, leading to low genetic diversity. To avoid this problem plants have established several pollination systems, among which the most elegant system is surely SI. The SI trait can be observed in Brassica crops, including cabbage, broccoli, turnip and radish. To produce hybrid seed of these crops efficiently, the SI trait has been employed in an agricultural context. From another point of view, the recognition reaction of SI during pollen-stigma interaction is an excellent model system for cell-cell communication and signal transduction in higher plants. In this review, we describe the molecular mechanisms of SI in Brassicaceae, which have been dissected by genetic, physiological, and biological approaches, and we discuss the future prospects in relation to associated scientific fields and new technologies.