Evolution and Diversity of the Wild Rice Oryza officinalis Complex, across Continents, Genome Types, and Ploidy Levels.

The Oryza officinalis complex is the largest species group in Oryza, with more than nine species from four continents, and is a tertiary gene pool that can be exploited in breeding programs for the improvement of cultivated rice. Most diploid and tetraploid members of this group have a C genome. Using a new reference C genome for the diploid species O. officinalis, and draft genomes for two other C genome diploid species Oryza eichingeri and Oryza rhizomatis, we examine the influence of transposable elements on genome structure and provide a detailed phylogeny and evolutionary history of the Oryza C genomes. The O. officinalis genome is 1.6 times larger than the A genome of cultivated Oryza sativa, mostly due to proliferation of Gypsy type long-terminal repeat transposable elements, but overall syntenic relationships are maintained with other Oryza genomes (A, B, and F). Draft genome assemblies of the two other C genome diploid species, Oryza eichingeri and Oryza rhizomatis, and short-read resequencing of a series of other C genome species and accessions reveal that after the divergence of the C genome progenitor, there was still a substantial degree of variation within the C genome species through proliferation and loss of both DNA and long-terminal repeat transposable elements. We provide a detailed phylogeny and evolutionary history of the Oryza C genomes and a genomic resource for the exploitation of the Oryza tertiary gene pool.

Functional screening of salt tolerance genes from a halophyte Sporobolus virginicus and transcriptomic and metabolomic analysis of salt tolerant plants expressing glycine-rich RNA-binding protein.

Sporobolus virginicus is a halophytic C4 grass found worldwide, from tropical to warm temperate regions. One Japanese genotype showed a salinity tolerance up to 1.5 M NaCl, a three-fold higher concentration than the salinity of sea water. To identify the key genes involved in the regulation of salt tolerance in S. virginicus, we produced 3500 independent transgenic Arabidopsis lines expressing random cDNA from S. virginicus and screened 10 lines which showed enhanced salt tolerance compared with the wild type in a medium containing 150 mM NaCl. Among the selected lines, two contained cDNA coding glycine-rich RNA-binding proteins (SvGRP1 and SvGRP2). This is the first reports on the function of GRPs from halophytes in salt tolerance though reports have shown GRPs are involved in diverse biological and biochemical processes including salt tolerance in Arabidopsis and some other glycophytes. Transcriptomic analysis and GO enrichment analysis of SvGRP1-expressing Arabidopsis under salt stress revealed upregulation of polyol and downregulation of glucosinolate and indole acetic acid biosynthesis/metabolic pathways. Metabolomic analysis of the SvGRP1-transformant suggested that the increase in 3-aminoppropanoic acid, citramalic acid, and isocitric acid content was associated with enhanced salt tolerance. These findings could provide novel insight into the roles of GRPs in plant salt tolerance.

TodoFirGene: Developing Transcriptome Resources for Genetic Analysis of Abies sachalinensis.

Todo-matsu (Abies sachalinensis) is one of the most important forestry species in Hokkaido, Japan and is distributed from near sea level to the alpine zone. Due to its wide spatial distribution, the species adapts to its environment, displaying phenotypes of ecological relevance. In order to identify candidate genes under natural selection, we collected the transcriptome from the female and male flower, leaf and inner bark. De novo assembly with 34.7 Gb of sequencing reads produced 158,542 transcripts from 69,618 loci, whose estimated coverage reached 95.6% of conserved eukaryotic genes. Homology searches against publicly available databases identified 134,190 (84.6%) transcripts with at least one hit. In total, 28,944 simple sequence repeats (SSRs) and 80,758 single nucleotide variants (SNVs) were detected from 23,570 (14.9%) and 25,366 (16.0%) transcripts, which were valuable for use in genetic analysis of the species. All the annotations were included in a relational database, TodoFirGene, which provides an interface for various queries and homology search, and can be accessed at http://plantomics.mind.meiji.ac.jp/todomatsu/. This database hosts not only the A. sachalinensis transcriptome but also links to the proteomes of 13 other species, allowing a comparative genomic study of plant species.

PlantExpress: A Database Integrating OryzaExpress and ArthaExpress for Single-species and Cross-species Gene Expression Network Analyses with Microarray-Based Transcriptome Data.

Publicly available microarray-based transcriptome data on plants are remarkably valuable in terms of abundance and variation of samples, particularly for Oryza sativa (rice) and Arabidopsis thaliana (Arabidopsis). Here, we introduce the web database PlantExpress (http://plantomics.mind.meiji.ac.jp/PlantExpress/) as a platform for gene expression network (GEN) analysis with the public microarray data of rice and Arabidopsis. PlantExpress has two functional modes. The single-species mode is specialized for GEN analysis within one of the species, while the cross-species mode is optimized for comparative GEN analysis between the species. The single-species mode for rice is the new version of OryzaExpress, which we have maintained since 2006. The single-species mode for Arabidopsis, named ArthaExpress, was newly developed. PlantExpress stores data obtained from three microarrays, the Affymetrix Rice Genome Array, the Agilent Rice Gene Expression 4x44K Microarray, and the Affymetrix Arabidopsis ATH1 Genome Array, with respective totals of 2,678, 1,206, and 10,940 samples. This database employs a 'MyList' function with which users may save lists of arbitrary genes and samples (experimental conditions) to use in analyses. In cross-species mode, the MyList function allows performing comparative GEN analysis between rice and Arabidopsis. In addition, the gene lists saved in MyList can be directly exported to the PODC database, which provides information and a platform for comparative GEN analysis based on RNA-seq data and knowledge-based functional annotation of plant genes. PlantExpress will facilitate understanding the biological functions of plant genes.

TOMATOMICS: A Web Database for Integrated Omics Information in Tomato.

Solanum lycopersicum (tomato) is an important agronomic crop and a major model fruit-producing plant. To facilitate basic and applied research, comprehensive experimental resources and omics information on tomato are available following their development. Mutant lines and cDNA clones from a dwarf cultivar, Micro-Tom, are two of these genetic resources. Large-scale sequencing data for ESTs and full-length cDNAs from Micro-Tom continue to be gathered. In conjunction with information on the reference genome sequence of another cultivar, Heinz 1706, the Micro-Tom experimental resources have facilitated comprehensive functional analyses. To enhance the efficiency of acquiring omics information for tomato biology, we have integrated the information on the Micro-Tom experimental resources and the Heinz 1706 genome sequence. We have also inferred gene structure by comparison of sequences between the genome of Heinz 1706 and the transcriptome, which are comprised of Micro-Tom full-length cDNAs and Heinz 1706 RNA-seq data stored in the KaFTom and Sequence Read Archive databases. In order to provide large-scale omics information with streamlined connectivity we have developed and maintain a web database TOMATOMICS (http://bioinf.mind.meiji.ac.jp/tomatomics/). In TOMATOMICS, access to the information on the cDNA clone resources, full-length mRNA sequences, gene structures, expression profiles and functional annotations of genes is available through search functions and the genome browser, which has an intuitive graphical interface.

Practical Utilization of OryzaExpress and Plant Omics Data Center Databases to Explore Gene Expression Networks in Oryza Sativa and Other Plant Species.

Analysis of a gene expression network (GEN), which is constructed based on similarity of gene expression profiles, is a widely used approach to gain clues for new biological insights. The recent abundant availability of transcriptome data in public databases is enabling GEN analysis under various experimental conditions, and even comparative GEN analysis across species. To provide a platform to gain biological insights from public transcriptome data, valuable databases have been created and maintained. This chapter introduces the web database OryzaExpress, providing omics information on Oryza sativa (rice). The integrated database Plant Omics Data Center, supporting a wide variety of plant species, is also described to compare omics information among multiple plant species.

CATchUP: A Web Database for Spatiotemporally Regulated Genes.

For proper control of biological activity, some key genes are highly expressed in a particular spatiotemporal domain. Mining of such spatiotemporally expressed genes using large-scale gene expression data derived from a broad range of experimental sources facilitates our understanding of genome-scale functional gene networks. However, comprehensive information on spatiotemporally expressed genes is lacking in plants. To collect such information, we devised a new index, Δdmax, which is the maximum difference in relative gene expression levels between sample runs which are neighboring when sorted by the levels. Employing this index, we comprehensively evaluated transcripts using large-scale RNA sequencing (RNA-Seq) data stored in the Sequence Read Archive for eight plant species: Arabidopsis thaliana (Arabidopsis), Solanum lycopersicum (tomato), Solanum tuberosum (potato), Oryza sativa (rice), Sorghum bicolor (sorghum), Vitis vinifera (grape), Medicago truncatula (Medicago), and Glycine max (soybean). Based on the frequency distribution of the Δdmax values, approximately 70,000 transcripts showing 0.3 or larger Δdmax values were extracted for the eight species. Information on these genes including the Δdmax values, functional annotations, conservation among species, and experimental conditions where the genes show high expression levels is provided in a new database, CATchUP (http://plantomics.mind.meiji.ac.jp/CATchUP). The CATchUP database assists in identifying genes specifically expressed under particular conditions with powerful search functions and an intuitive graphical user interface.

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.

Comprehensive analysis of transcriptome response to salinity stress in the halophytic turf grass Sporobolus virginicus.

The turf grass Sporobolus virginicus is halophyte and has high salinity tolerance. To investigate the molecular basis of its remarkable tolerance, we performed Illumina high-throughput RNA sequencing on roots and shoots of a S. virginicus genotype under normal and saline conditions. The 130 million short reads were assembled into 444,242 unigenes. A comparative analysis of the transcriptome with rice and Arabidopsis transcriptome revealed six turf grass-specific unigenes encoding transcription factors. Interestingly, all of them showed root specific expression and five of them encode bZIP type transcription factors. Another remarkable transcriptional feature of S. virginicus was activation of specific pathways under salinity stress. Pathway enrichment analysis suggested transcriptional activation of amino acid, pyruvate, and phospholipid metabolism. Up-regulation of several unigenes, previously shown to respond to salt stress in other halophytes was also observed. Gene Ontology enrichment analysis revealed that unigenes assigned as proteins in response to water stress, such as dehydrin and aquaporin, and transporters such as cation, amino acid, and citrate transporters, and H(+)-ATPase, were up-regulated in both shoots and roots under salinity. A correspondence analysis of the enriched pathways in turf grass cells, but not in rice cells, revealed two groups of unigenes similarly up-regulated in the turf grass in response to salt stress; one of the groups, showing excessive up-regulation under salinity, included unigenes homologos to salinity responsive genes in other halophytes. Thus, the present study identified candidate genes involved in salt tolerance of S. virginicus. This genetic resource should be valuable for understanding the mechanisms underlying high salt tolerance in S. virginicus. This information can also provide insight into salt tolerance in other halophytes.

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.

Immunohistochemical study of endothelial nitric oxide synthase in the trigeminal ganglia of a crotaline snake Trimeresurus flavoviridis.

The immunoreactivity of constitutive endothelial nitric oxide synthase (eNOS) was studied in the trigeminal ganglia (TG) of a crotaline snake, Trimeresurus flavoviridis. eNOS immunoreactivity was found in TG neurons of different sizes. The percentage of eNOS-positive TG neurons was significantly higher in the mandibular division than in the infrared-related divisions, the maxillary division and ophthalmic ganglion (p<0.001). These findings suggest that eNOS in the TG of crotaline snakes is involved in constitutive neurotransmission in the TG, and is minimally involved in processing in the infrared-sensory system.

Snake infrared receptors respond to dimethylsulfoxide in the blood stream.

1. We used extracellular recording of the infrared (IR)-sensitive trigeminal ganglion (TG) neurons (primary neurons) of a crotaline snake, Trimeresurus flavoviridis, which has very sensitive thermoreceptors, to examine changes in the IR response induced by dimethylsulfoxide (DMSO), in vivo. 2. The responses in the TG were recorded after each concentration of DMSO (1, 10, and 25%) was administered in the bloodstream. 3. At a constant temperature, DMSO dose-dependently potentiated the IR-triggered discharges of IR-sensitive TG neurons in this snake. 4. It is suggested that the increased IR response to DMSO is due to its chemical effect, or to an indirect effect via its vasoactive role in the thermoreceptors of IR-sensitive snakes.

Immunohistochemical localization of the delta subspecies of protein kinase C in the trigeminal sensory system of Trimeresurus flavoviridis, an infrared-sensitive snake.

We examined the expression of the protein kinase C (PKC) delta subspecies in the trigeminal sensory system of the infrared-sensitive snake Trimeresurus flavoviridis. In the trigeminal ganglion (TG), diffuse low-intensity PKC delta immunoreactivity was found in TG neurons and fibers, while intense reactions were observed mainly in medium-sized neurons, which include most of the infrared-sensitive neurons. In the brainstem, intense PKC delta immunoreactivity was present in the intermediate layer of the optic tectum of the midbrain and in the nucleus descendens lateralis n. trigemini of the medulla oblongata; these areas are related to the infrared sensory pathway. In the pit organ (the infrared receptor), PKC delta immunoreactivity was present in terminal nerve masses in the pit membrane. These findings suggest that the PKC delta subspecies is involved in the infrared sensory pathway in the trigeminal sensory system of the infrared-sensitive snake.

Degeneration of infrared receptor terminals of snakes caused by capsaicin.

Capsaicin, the main pungent ingredient in hot peppers (genus Capsicum), caused degeneration of the infrared receptor terminals in infrared sensitive snakes, Trimeresurus flavoviridis, when it was applied perineurally to a branch of the trigeminal nerve. The degeneration of the terminals was found 6 h after the application. This finding suggests that capsaicin stimulates this infrared receptor terminal, a kind of warm receptor terminal.

Response of the infrared receptors of a crotaline snake to ethanol.

The pit organs of crotaline snakes can sense infrared (IR). The pit membrane has a finer, flatter, more convoluted vasculature than other sensory organs. Using extracellular recording from IR-sensitive trigeminal ganglion (TG) neurons (primary neurons) and tectal (OT) neurons of the crotaline snake Trimeresurus flavoviridis, we examined the IR response to ethanol (EtOH) in vivo. The response to EtOH was recorded in the TG and OT 20-80 s after 10% EtOH in Ringer's solution (100 microl/ 500 g body weight) was injected via the heart. The responses to EtOH and those to lower or higher temperature stimulation were additive. At a constant temperature (25 degrees C), EtOH significantly potentiated the IR-triggered discharges of IR-sensory pathways in this snake. These results suggest that the IR response to EtOH is due to either its vasodilatory effect on the abundant vasculature of the pit membrane or its chemical effect on temperature-sensitive receptors.

Immunohistochemical analysis of neuronal nitric oxide synthase in the trigeminal ganglia of the crotaline snake Trimeresurus flavoviridis.

The expression of neuronal nitric oxide synthase (nNOS) in the trigeminal ganglia (TG) of the infrared-sensitive crotaline snake Trimeresurus flavoviridis was studied immunohistochemically. The percentage of nNOS-positive (+) neurons in the TG was significantly higher (about 3.5-fold, P<0.001) in the mandibular division than in the infrared-sensory processing area (the maxillary division and ophthalmic ganglion). nNOS was found in varying sizes of TG neurons. However, nNOS (+) neurons were more abundant in small and large neurons than in medium-sized neurons, which include most of the infrared-sensitive neurons of the TG. These findings suggest that nNOS may be involved in normal physiological functions, such as the transmission of tactile, vibrotactile, and nociceptive sensations in the TG, rather than in infrared sensory processing in this species.