Dominance in self-compatibility between subgenomes of allopolyploid Arabidopsis kamchatica shown by transgenic restoration of self-incompatibility.

The evolutionary transition to self-compatibility facilitates polyploid speciation. In Arabidopsis relatives, the self-incompatibility system is characterized by epigenetic dominance modifiers, among which small RNAs suppress the expression of a recessive SCR/SP11 haplogroup. Although the contribution of dominance to polyploid self-compatibility is speculated, little functional evidence has been reported. Here we employ transgenic techniques to the allotetraploid plant A. kamchatica. We find that when the dominant SCR-B is repaired by removing a transposable element insertion, self-incompatibility is restored. This suggests that SCR was responsible for the evolution of self-compatibility. By contrast, the reconstruction of recessive SCR-D cannot restore self-incompatibility. These data indicate that the insertion in SCR-B conferred dominant self-compatibility to A. kamchatica. Dominant self-compatibility supports the prediction that dominant mutations increasing selfing rate can pass through Haldane's sieve against recessive mutations. The dominance regulation between subgenomes inherited from progenitors contrasts with previous studies on novel epigenetic mutations at polyploidization termed genome shock.

Two aquaporins, SIP1;1 and PIP1;2, mediate water transport for pollen hydration in the Arabidopsis pistil.

Pollination is the crucial initial step that brings together the male and female gametophytes, and occurs at the surface of the stigmatic papilla cell in Arabidopsis thaliana. After pollen recognition, pollen hydration is initiated as a second critical step to activate desiccated mature pollen grains for germination, and thus water transport from pistil to pollen is essential for this process. In this study, we report a novel aquaporin-mediated water transport process in the papilla cell as a control mechanism for pollen hydration. Coupled with a time-series imaging analysis of pollination and a reverse genetic analysis using T-DNA insertion Arabidopsis mutants, we found that two aquaporins, the ER-bound SIP1;1 and the plasma membrane-bound PIP1;2, are key players in water transport from papilla cell to pollen during pollination. In wild type plant, hydration speed reached its maximal value within 5 min after pollination, remained high until 10-15 min. In contrast, sip1;1 and pip1;2 mutants showed no rapid increase of hydration speed, but instead a moderate increase during ∼25 min after pollination. Pollen of sip1;1 and pip1;2 mutants had normal viability without any functional defects for pollination, indicating that decelerated pollen hydration is due to a functional defect on the female side in sip1;1 and pip1;2 mutants. In addition, sip1;1 pip1;2 double knockout mutant showed a similar impairment of pollen hydration to individual single mutants, suggesting that their coordinated regulation is critical for proper water transport, in terms of speed and amount, in the pistil to accomplish successful pollen hydration.

Spatiogenetic characterization of S receptor kinase (SRK) alleles in naturalized populations of Raphanus sativus L. var. raphanistroides on Yakushima island.

In various coastal areas of Japan, naturalized radish populations are observed. Radish is a cruciferous plant and exhibits self-incompatibility, involving a system controlled by a single locus with multiple S alleles. Although the S allele diversity of radish cultivars and wild radishes has been characterized, the S allele distribution in naturalized populations has not yet been analyzed in relation to the positions of the plants in situ. Here, we show the S allele distribution in naturalized radish populations of Yakushima, a small island in the East China Sea, with positions of the plants. Radish plants were sampled in coastal areas in Yakushima, and their S alleles were detected and characterized. Most of the S alleles had been previously identified in radish cultivars. However, four novel S alleles, which may be unique to Yakushima, were also found. Moreover, seeds in siliques from plants growing in the study areas were sampled, and S allele determination in DNA extracted from these seeds suggested that the plants had exchanged their pollen among their close neighbors. There was also a problem in that the PCR amplification of some SRK alleles was difficult because of their sequence diversity in the naturalized populations, as occurs in cultivars. Our results suggest that the exchange of S alleles between cultivars and naturalized populations occurs and that S alleles in naturalized populations are highly diverse. The methodology established in our study should be applicable to other self-incompatible species to dissect the diversity of S allele distribution in naturalized populations.

Double-Locking Mechanism of Self-Compatibility in Arabidopsis thaliana: The Synergistic Effect of Transcriptional Depression and Disruption of Coding Region in the Male Specificity Gene.

Self-compatibility in Arabidopsis thaliana represents the relatively recent disruption of ancestral obligate cross pollination, recognized as one of the prevalent evolutionary pathways in flowering plants, as noted by Darwin. Our previous study found that inversion of the male specificity gene (SP11/SCR) disrupted self-incompatibility, which was restored by overexpressing the SCR with the reversed inversion. However, SCR in A. thaliana has other mutations aside from the pivotal inversion, in both promoter and coding regions, with probable effects on transcriptional regulation. To examine the functional consequences of these mutations, we conducted reciprocal introductions of native promoters and downstream sequences from orthologous loci of self-compatible A. thaliana and self-incompatible A. halleri. Use of this inter-species pair enabled us to expand the scope of the analysis to transcriptional regulation and deletion in the intron, in addition to inversion in the native genomic background. Initial analysis revealed that A. thaliana has a significantly lower basal expression level of SCR transcripts in the critical reproductive stage compared to that of A. halleri, suggesting that the promoter was attenuated in inducing transcription in A. thaliana. However, in reciprocal transgenic experiments, this A. thaliana promoter was able to restore partial function if coupled with the functional A. halleri coding sequence, despite extensive alterations due to the self-compatible mode of reproduction in A. thaliana. This represents a synergistic effect of the promoter and the inversion resulting in fixation of self-compatibility, primarily enforced by disruption of SCR. Our findings elucidate the functional and evolutionary context of the historical transition in A. thaliana thus contributing to the understanding of the molecular events leading to development of self-compatibility.

Genetic and tissue-specific RNA-sequencing analysis of self-compatible mutant TSC28 in Brassica rapa L. toward identification of a novel self-incompatibility factor.

Self-incompatibility (SI) is a sophisticated system for pollen selectivity to prevent pollination by genetically identical pollen. In Brassica, it is genetically controlled by a single, highly polymorphic S-locus, and the male and female S-determinant factors have been identified as S-locus protein 11 (SP11)/S-locus cysteine-rich protein (SCR) and S-locus receptor kinase (SRK), respectively. However, the overall molecular system and identity of factors in the downstream cascade of the SI reaction remain unclear. Previously, we identified a self-compatible B. rapa mutant line, TSC28, which has a disruption in an unidentified novel factor of the SI signaling cascade. Here, in a genetic analysis of TSC28, using an F2 population from a cross with the reference B. rapa SI line Chiifu-401, the causal gene was mapped to a genetic region of DNA containing markers BrSA64 and ACMP297 in B. rapa chromosome A1. By fine mapping using an F2 population of 1,034 plants, it was narrowed down to a genetic region between DNA markers ACMP297 and BrgMS4028, with physical length approximately 1.01 Mbp. In this genomic region, 113 genes are known to be located and, among these, we identified 55 genes that were expressed in the papilla cells. These are candidates for the gene responsible for the disruption of SI in TSC28. This list of candidate genes will contribute to the discovery of a novel downstream factor in the SP11-SRK signaling cascade in the Brassica SI system.

Abscisic acid-mediated developmental flexibility of stigmatic papillae in response to ambient humidity in Arabidopsis thaliana.

Stigmatic papillae develop at the apex of the gynoecium and play an important role as a site of pollination. The papillae in Brassicaceae are of the dry and unicellular type, and more than 15,000 genes are expressed in the papillae; however, the molecular and physiological mechanisms of their development remain unknown. We found that the papillae in Arabidopsis thaliana change their length in response to altered ambient humidity: papillae of flowers incubated under high humidity elongated more than those under normal humidity conditions. Genetic analysis and transcriptome data suggest that an abscisic acid-mediated abiotic stress response mechanism regulates papilla length. Our data suggest a flexible regulation of papilla elongation at the post-anthesis stage, in response to abiotic stress, as an adaptation to environmental conditions.

Duplicated pollen-pistil recognition loci control intraspecific unilateral incompatibility in Brassica rapa.

In plants, cell-cell recognition is a crucial step in the selection of optimal pairs of gametes to achieve successful propagation of progeny. Flowering plants have evolved various genetic mechanisms, mediated by cell-cell recognition, to enable their pistils to reject self-pollen, thus preventing inbreeding and the consequent reduced fitness of progeny (self-incompatibility, SI), and to reject foreign pollen from other species, thus maintaining species identity (interspecific incompatibility)1. In the genus Brassica, the SI system is regulated by an S-haplotype-specific interaction between a stigma-expressed female receptor (S receptor kinase, SRK) and a tapetum cell-expressed male ligand (S locus protein 11, SP11), encoded by their respective polymorphic genes at the S locus2-6. However, the molecular mechanism for recognition of foreign pollen, leading to reproductive incompatibility, has not yet been identified. Here, we show that recognition between a novel pair of proteins, a pistil receptor SUI1 (STIGMATIC UNILATERAL INCOMPATIBILITY 1) and a pollen ligand PUI1 (POLLEN UNILATERAL INCOMPATIBILITY 1), triggers unilateral reproductive incompatibility between plants of two geographically distant self-incompatible Brassica rapa lines, even though crosses would be predicted to be compatible based on the S haplotypes of pollen and stigma. Interestingly, SUI1 and PUI1 are similar to the SI genes, SRK and SP11, respectively, and are maintained as cryptic incompatibility genes in these two populations. The duplication of the SRK and SP11 followed by reciprocal loss in different populations would provide a molecular mechanism of the emergence of a reproductive barrier in allopatry.

Identification of reference genes for quantitative expression analysis using large-scale RNA-seq data of Arabidopsis thaliana and model crop plants.

In quantitative gene expression analysis, normalization using a reference gene as an internal control is frequently performed for appropriate interpretation of the results. Efforts have been devoted to exploring superior novel reference genes using microarray transcriptomic data and to evaluating commonly used reference genes by targeting analysis. However, because the number of specifically detectable genes is totally dependent on probe design in the microarray analysis, exploration using microarray data may miss some of the best choices for the reference genes. Recently emerging RNA sequencing (RNA-seq) provides an ideal resource for comprehensive exploration of reference genes since this method is capable of detecting all expressed genes, in principle including even unknown genes. We report the results of a comprehensive exploration of reference genes using public RNA-seq data from plants such as Arabidopsis thaliana (Arabidopsis), Glycine max (soybean), Solanum lycopersicum (tomato) and Oryza sativa (rice). To select reference genes suitable for the broadest experimental conditions possible, candidates were surveyed by the following four steps: (1) evaluation of the basal expression level of each gene in each experiment; (2) evaluation of the expression stability of each gene in each experiment; (3) evaluation of the expression stability of each gene across the experiments; and (4) selection of top-ranked genes, after ranking according to the number of experiments in which the gene was expressed stably. Employing this procedure, 13, 10, 12 and 21 top candidates for reference genes were proposed in Arabidopsis, soybean, tomato and rice, respectively. Microarray expression data confirmed that the expression of the proposed reference genes under broad experimental conditions was more stable than that of commonly used reference genes. These novel reference genes will be useful for analyzing gene expression profiles across experiments carried out under various experimental conditions.

Comparative analysis of microRNA profiles of rice anthers between cool-sensitive and cool-tolerant cultivars under cool-temperature stress.

Plants subjected to abiotic stress can regulate gene expression post-transcriptionally by means of small RNAs such as microRNAs. Cool-temperature stress causes abnormal tapetum hypertrophy in rice anthers, leading to pollen sterility. As a first step toward understanding the molecular mechanisms of cool tolerance in developing anthers of rice, we report here a comprehensive comparative analysis of microRNAs between cool-sensitive Sasanishiki and cool-tolerant Hitomebore cultivars. High-throughput Illumina sequencing revealed 241 known and 46 novel microRNAs. Interestingly, 15 of these microRNAs accumulated differentially in the two cultivars at the uninucleate microspore stage under cool conditions. Inverse correlations between expression patterns of microRNAs and their target genes were confirmed by quantitative RT-PCR analysis, and cleavage sites of some of the target genes were determined by 5' RNA ligase-mediated RACE experiments. Thus, our data are useful resources to elucidate microRNA-mediated mechanism(s) of cool tolerance in rice anthers at the booting stage.

Plant Omics Data Center: an integrated web repository for interspecies gene expression networks with NLP-based curation.

Comprehensive integration of large-scale omics resources such as genomes, transcriptomes and metabolomes will provide deeper insights into broader aspects of molecular biology. For better understanding of plant biology, we aim to construct a next-generation sequencing (NGS)-derived gene expression network (GEN) repository for a broad range of plant species. So far we have incorporated information about 745 high-quality mRNA sequencing (mRNA-Seq) samples from eight plant species (Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, Sorghum bicolor, Vitis vinifera, Solanum tuberosum, Medicago truncatula and Glycine max) from the public short read archive, digitally profiled the entire set of gene expression profiles, and drawn GENs by using correspondence analysis (CA) to take advantage of gene expression similarities. In order to understand the evolutionary significance of the GENs from multiple species, they were linked according to the orthology of each node (gene) among species. In addition to other gene expression information, functional annotation of the genes will facilitate biological comprehension. Currently we are improving the given gene annotations with natural language processing (NLP) techniques and manual curation. Here we introduce the current status of our analyses and the web database, PODC (Plant Omics Data Center; http://bioinf.mind.meiji.ac.jp/podc/), now open to the public, providing GENs, functional annotations and additional comprehensive omics resources.

Transcriptional characteristics and differences in Arabidopsis stigmatic papilla cells pre- and post-pollination.

Pollination is an important early step in sexual plant reproduction. In Arabidopsis thaliana, sequential pollination events, from pollen adhesion onto the stigma surface to pollen tube germination and elongation, occur on the stigmatic papilla cells. Following successful completion of these events, the pollen tube penetrates the stigma and finally fertilizes a female gametophyte. The pollination events are thought to be initiated and regulated by interactions between papilla cells and pollen. Here, we report the characterization of gene expression profiles of unpollinated (UP), compatible pollinated (CP) and incompatible pollinated (IP) papilla cells in A. thaliana. Based on cell type-specific transcriptome analysis from a combination of laser microdissection and RNA sequencing, 15,475, 17,360 and 16,918 genes were identified as expressed in UP, CP and IP papilla cells, respectively, and, of these, 14,392 genes were present in all three data sets. Differentially expressed gene (DEG) analyses identified 147 and 71 genes up-regulated in CP and IP papilla cells, respectively, and 115 and 46 genes down-regulated. Gene Ontology and metabolic pathway analyses revealed that papilla cells play an active role as the female reproductive component in pollination, particularly in information exchange, signal transduction, internal physiological changes and external morphological modification. This study provides fundamental information on the molecular mechanisms involved in pollination in papilla cells, furthering our understanding of the reproductive role of papilla cells.

Physiological and genetic analysis of CO2-induced breakdown of self-incompatibility in Brassica rapa.

Self-incompatibility (SI) of the Brassicaceae family can be overcome by CO2 gas treatment. This method has been used for decades as an effective means to obtain a large amount of inbred seeds which can then be used for F1 hybrid seed production; however, the molecular mechanism by which CO2 alters the SI pathway has not been elucidated. In this study, to obtain new insights into the mechanism of CO2-induced SI breakdown, the focus was on two inbred lines of Brassica rapa (syn. campestris) with different CO2 sensitivity. Physiological examination using X-ray microanalysis suggested that SI breakdown in the CO2-sensitive line was accompanied by a significant accumulation of calcium at the pollen-stigma interface. Pre-treatment of pollen or pistil with CO2 gas before pollination showed no effect on the SI reaction, suggesting that some physiological process after pollination is necessary for SI to be overcome. Genetic analyses using F1 progeny of a CO2-sensitive × CO2-insensitive cross suggested that CO2 sensitivity is a semi-dominant trait in these lines. Analysis of F2 progeny suggested that CO2 sensitivity could be a quantitative trait, which is controlled by more than one gene. Quantitative trait locus (QTL) analyses identified two major loci, BrSIO1 and BrSIO2, which work additively in overcoming SI during CO2 treatment. No QTL was detected at the loci previously shown to affect SI stability, suggesting that CO2 sensitivity is determined by novel genes. The QTL data presented here should be useful for determining the responsible genes, and for the marker-assisted selection of desirable parental lines with stable but CO2-sensitive SI in F1 hybrid breeding.

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.

Time-lapse imaging of self- and cross-pollinations in Brassica rapa.

BACKGROUND AND AIMS: Pollination is an important process in the life cycle of plants and is the first step in bringing together the male and female gametophytes for plant reproduction. While pollination has been studied for many years, accurate knowledge of the morphological aspects of this process is still far from complete. This study therefore focuses on a morphological characterization of pollination, using time-series image analysis of self- and cross-pollinations in Brassica rapa. METHODS: Time-lapse imaging of pollen behaviour during self- and cross-pollinations was recorded for 90 min, at 1 min intervals, using a stereoscopic microscope. Using time-series digital images of pollination, characteristic features of pollen behaviours during self- and cross-pollinations were studied. KEY RESULTS: Pollen exhibited various behaviours in both self- and cross-pollinations, and these were classified into six representative patterns: germination, expansion, contraction, sudden contraction, pulsation and no change. It is noteworthy that in 'contraction' pollen grains shrunk within a short period of 30-50 min, and in 'pulsation' repeated expansion and contraction occurred with an interval of 10 min, suggesting that a dehydration system is operating in pollination. All of the six patterns were observed on an individual stigma with both self- and cross-pollinations, and the difference between self- and cross-pollinations was in the ratios of the different behaviours. With regard to water transport to and from pollen grains, this occurred in multiple steps, before, during and after hydration. Thus, pollination is regulated by a combination of multiple components of hydration, rehydration and dehydration systems. CONCLUSIONS: Regulated hydration of pollen is a key process for both pollination and self-incompatibility, and this is achieved by a balanced complex of hydration, dehydration and nutrient supply to pollen grains from stigmatic papilla cells.

Identification and characterization of Crr1a, a gene for resistance to clubroot disease (Plasmodiophora brassicae Woronin) in Brassica rapa L.

Clubroot disease, caused by the obligate biotrophic protist Plasmodiophora brassicae Woronin, is one of the most economically important diseases of Brassica crops in the world. Although many clubroot resistance (CR) loci have been identified through genetic analysis and QTL mapping, the molecular mechanisms of defense responses against P. brassicae remain unknown. Fine mapping of the Crr1 locus, which was originally identified as a single locus, revealed that it comprises two gene loci, Crr1a and Crr1b. Here we report the map-based cloning and characterization of Crr1a, which confers resistance to clubroot in Brassica rapa. Crr1a(G004), cloned from the resistant line G004, encodes a Toll-Interleukin-1 receptor/nucleotide-binding site/leucine-rich repeat (TIR-NB-LRR) protein expressed in the stele and cortex of hypocotyl and roots, where secondary infection of the pathogen occurs, but not in root hairs, where primary infection occurs. Gain-of-function analysis proved that Crr1a(G004) alone conferred resistance to isolate Ano-01 in susceptible Arabidopsis and B. rapa. In comparison, the susceptible allele Crr1a(A9709) encodes a truncated NB-LRR protein, which lacked more than half of the TIR domain on account of the insertion of a solo-long terminal repeat (LTR) in exon 1 and included several substitutions and insertion-deletions in the LRR domain. This study provides a basis for further molecular analysis of defense mechanisms against P. brassicae and will contribute to the breeding of resistant cultivars of Brassica vegetables by marker-assisted selection.Data deposition The sequence reported in this paper has been deposited in the GenBank database (accession no. AB605024).

Demonstration in vivo of the role of Arabidopsis PLIM2 actin-binding proteins during pollination.

In plant reproduction, pollination is the initial key process in bringing together the male and female gametophytes. When a pollen grain lands on the surface of the stigma, information is exchanged between the pollen and stigmatic cell to determine whether the pollen grain will be accepted or rejected. If it is accepted, the stigmatic papilla cell supplies water and other resources to the pollen for germination and pollen tube elongation. Cellular processes involving actin are essential for pollen germination and tube growth, and actin-binding proteins regulate these processes by interacting with actin filaments to assemble cytoskeletal structures and actin networks. LIM proteins, which belong to a subfamily of cysteine-rich proteins, are a family of actin-binding proteins in plants, and are considered to be important for formation of the actin cytoskeleton and maintenance of its dynamics. Although the physiological and biochemical characteristics of LIMs have been elucidated in vitro in a variety of cell types, their exact role in pollen germination and pollen tube growth during pollination remained unclear. In this manuscript, we focus on the pollen-specific LIM proteins, AtPLIM2a and AtPLIM2c, and define their biological function during pollination in Arabidopsis thaliana. The atplim2a/atplim2c double knockdown RNAi plants showed a reduced pollen germination, approximately one-fifth of wild type, and slower pollen tube growth in the pistil, that is 80.4 µm/hr compared to 140.8 µm/hr in wild type. These defects led to an occasional unfertilized ovule at the bottom of the silique in RNAi plants. Our data provide direct evidence of the biological function of LIM proteins during pollination as actin-binding proteins, modulating cytoskeletal structures and actin networks, and their consequent importance in seed production.

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.

Microstructure of a Brassica rapa genome segment homoeologous to the resistance gene cluster on Arabidopsis chromosome 4.

Genome evolution is a continuous process and genomic rearrangement occurs both within and between species. With the sequencing of the Arabidopsis thaliana genome, comparative genetics and genomics offer new insights into plant biology. The genus Brassica offers excellent opportunities with which to compare genomic synteny so as to reveal genome evolution. During a previous genetic analysis of clubroot resistance in Brassica rapa, we identified a genetic region that is highly collinear with Arabidopsis chromosome 4. This region corresponds to a disease resistance gene cluster in the A. thaliana genome. Relying on synteny with Arabidopsis, we fine-mapped the region and found that the location and order of the markers showed good correspondence with those in Arabidopsis. Microsynteny on a physical map indicated an almost parallel correspondence, with a few rearrangements such as inversions and insertions. The results show that this genomic region of Brassica is conserved extensively with that of Arabidopsis and has potential as a disease resistance gene cluster, although the genera diverged 20 million years ago.

Comprehensive network analysis of anther-expressed genes in rice by the combination of 33 laser microdissection and 143 spatiotemporal microarrays.

Co-expression networks systematically constructed from large-scale transcriptome data reflect the interactions and functions of genes with similar expression patterns and are a powerful tool for the comprehensive understanding of biological events and mining of novel genes. In Arabidopsis (a model dicot plant), high-resolution co-expression networks have been constructed from very large microarray datasets and these are publicly available as online information resources. However, the available transcriptome data of rice (a model monocot plant) have been limited so far, making it difficult for rice researchers to achieve reliable co-expression analysis. In this study, we performed co-expression network analysis by using combined 44 K agilent microarray datasets of rice, which consisted of 33 laser microdissection (LM)-microarray datasets of anthers, and 143 spatiotemporal transcriptome datasets deposited in RicexPro. The entire data of the rice co-expression network, which was generated from the 176 microarray datasets by the Pearson correlation coefficient (PCC) method with the mutual rank (MR)-based cut-off, contained 24,258 genes and 60,441 genes pairs. Using these datasets, we constructed high-resolution co-expression subnetworks of two specific biological events in the anther, "meiosis" and "pollen wall synthesis". The meiosis network contained many known or putative meiotic genes, including genes related to meiosis initiation and recombination. In the pollen wall synthesis network, several candidate genes involved in the sporopollenin biosynthesis pathway were efficiently identified. Hence, these two subnetworks are important demonstrations of the efficiency of co-expression network analysis in rice. Our co-expression analysis included the separated transcriptomes of pollen and tapetum cells in the anther, which are able to provide precise information on transcriptional regulation during male gametophyte development in rice. The co-expression network data presented here is a useful resource for rice researchers to elucidate important and complex biological events.

OryzaExpress: an integrated database of gene expression networks and omics annotations in rice.

Similarity of gene expression profiles provides important clues for understanding the biological functions of genes, biological processes and metabolic pathways related to genes. A gene expression network (GEN) is an ideal choice to grasp such expression profile similarities among genes simultaneously. For GEN construction, the Pearson correlation coefficient (PCC) has been widely used as an index to evaluate the similarities of expression profiles for gene pairs. However, calculation of PCCs for all gene pairs requires large amounts of both time and computer resources. Based on correspondence analysis, we developed a new method for GEN construction, which takes minimal time even for large-scale expression data with general computational circumstances. Moreover, our method requires no prior parameters to remove sample redundancies in the data set. Using the new method, we constructed rice GENs from large-scale microarray data stored in a public database. We then collected and integrated various principal rice omics annotations in public and distinct databases. The integrated information contains annotations of genome, transcriptome and metabolic pathways. We thus developed the integrated database OryzaExpress for browsing GENs with an interactive and graphical viewer and principal omics annotations (http://riceball.lab.nig.ac.jp/oryzaexpress/). With integration of Arabidopsis GEN data from ATTED-II, OryzaExpress also allows us to compare GENs between rice and Arabidopsis. Thus, OryzaExpress is a comprehensive rice database that exploits powerful omics approaches from all perspectives in plant science and leads to systems biology.