Akifumi Shimizu

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.

Genetic analysis of root elongation induced by phosphorus deficiency in rice (Oryza sativa L.): fine QTL mapping and multivariate analysis of related traits

Root elongation induced by phosphorus deficiency has been reported as one of the adaptive mechanisms in plants. Genetic differences were found in rice for the root elongation under phosphorus deficiency (REP), for which a distinct quantitative trait locus (QTL) was detected on the long arm of chromosome 6. Subsequently, the effect and position of the QTL, designated as qREP-6, were confirmed using chromosome segment substitution lines (CSSLs), in which the background of a japonica cultivar, ‘Nipponbare’ with non-REP, was partially substituted by chromosomal segments from an indica cultivar, ‘Kasalath’ with remarkable REP. Out of 54 CSSLs, two lines (CSSL28 and CSSL29) that retain a common ‘Kasalath’-derived segment on the long arm of chromosome 6 showed a significantly high REP. The high REP lines also showed high adaptabilities such as enhanced tillering ability and shoot phosphorus content. Accordingly, conditional dependencies between the related traits were assessed using a graphical Gaussian model (GGM). Direct interactions between REP and root length, and between root length and tiller number were detected under P deficiency in CSSLs. Furthermore, qREP-6 for REP and qTNP-6 for tiller number under P deficiency were fine-mapped with an F2 population of a cross between Nipponbare and CSSL29. A region containing qREP-6 accounted for more than half of the phenotypic variance, the most plausible interval of which contained 37 candidate genes. The result provides a foundation for cloning of the qREP-6 gene which will be applicable to study P deficiency-dependent response and to improve rice’s adaptability to P deficiency stress.

Genome-wide marker development for the wheat D genome based on single nucleotide polymorphisms identified from transcripts in the wild wheat progenitor Aegilops tauschii

Key message13,347 high-confidence SNPs were discovered through transcriptome sequencing ofAegilops tauschii, which are useful for genomic analysis and molecular breeding of hexaploid wheat.AbstractIn organisms with large and complex genomes, such as wheat, RNA-seq analysis is cost-effective for discovery of genome-wide single nucleotide polymorphisms (SNPs). In this study, deep sequencing of the spike transcriptome from two Aegilops tauschii accessions representing two major lineages led to the discovery of 13,347 high-confidence (HC) SNPs in 4,872 contigs. After removing redundant SNPs detected in the leaf transcriptome from the same accessions in an earlier study, 10,589 new SNPs were discovered. In total, 5,642 out of 5,808 contigs with HC SNPs were assigned to the Ae. tauschii draft genome sequence. On average, 732 HC polymorphic contigs were mapped in silico to each Ae. tauschii chromosome. Based on the polymorphic data, we developed markers to target the short arm of chromosome 2D and validated the polymorphisms using 20 Ae. tauschii accessions. Of the 29 polymorphic markers, 28 were successfully mapped to 2DS in the diploid F2 population of Ae. tauschii. Among ten hexaploid wheat lines, which included wheat synthetics and common wheat cultivars, 25 of the 43 markers were polymorphic. In the hexaploid F2 population between a common wheat cultivar and a synthetic wheat line, 23 of the 25 polymorphic markers between the parents were available for genotyping of the F2 plants and 22 markers mapped to chromosome 2DS. These results indicate that molecular markers that developed from polymorphisms between two distinct lineages of Ae. tauschii might be useful for analysis not only of the diploid, but also of the hexaploid wheat genome.

Discovery of High-Confidence Single Nucleotide Polymorphisms from Large-Scale De Novo Analysis of Leaf Transcripts of Aegilops tauschii, A Wild Wheat Progenitor

Construction of high-resolution genetic maps is important for genetic and genomic research, as well as for molecular breeding. Single nucleotide polymorphisms (SNPs) are the predominant class of genetic variation and can be used as molecular markers. Aegilops tauschii, the D-genome donor of common wheat, is considered a valuable genetic resource for wheat improvement. Our previous study implied that Ae. tauschii accessions can be genealogically divided into two major lineages. In this study, the transcriptome of two Ae. tauschii accessions from each lineage, lineage 1 (L1) and 2 (L2), was sequenced, yielding 9435 SNPs and 739 insertion/deletion polymorphisms (indels) after de novo assembly of the reads. Based on 36 contig sequences, 31 SNPs and six indels were validated on 20 diverse Ae. tauschii accessions. Because almost all of the SNP markers were polymorphic between L1 and L2, and the D-genome donor of common wheat is presumed to belong to L2, these markers are available for D-genome typing in crosses between common wheat varieties and L1-derived synthetic wheat. Due to the conserved synteny between wheat and barley chromosomes, the high-density expressed sequence tag barley map and the hypothetical gene order in barley can be applied to develop markers on target chromosomal regions in wheat.

Boron toxicity in rice (Oryza sativa L.). I. Quantitative trait locus (QTL) analysis of tolerance to boron toxicity

Boron toxicity tolerance of rice plants was studied. Modern japonica subspecies such as Koshihikari, Nipponbare, and Sasanishiki were tolerant, whereas indica subspecies such as Kasalath and IR36 were intolerant to excessive application of boron (B), even though their shoot B contents under B toxicity were not significantly different. Recombinant inbred lines (RILs) of japonica Nekken-1 and indica IR36 were used for quantitative trait locus (QTL) analysis to identify the gene responsible for B toxicity tolerance. A major QTL that could explain 45% of the phenotypic variation was detected in chromosome 4. The QTL was confirmed using a population derived from a recombinant inbred line which is heterogenic at the QTL region. The QTL was also confirmed in other chromosome segment substitution lines (CSSLs).

QTL Analysis of Genetic Tolerance to Iron Toxicity in Rice (Oryza Sativa L.) by Quantification of Bronzing Score

A QTL analysis for iron-toxicity tolerance was conducted in rice. Bronzing scores, which are regarded as an index of sensitivity to iron toxicity, were converted to digital data via a scanner and quantified to parametric data using chromaticity coordinates. On the basis of quantified score, QTL analysis for bronzing tolerance was conducted using F3 lines from a cross between a japonica cultivar ‘Gimbozu’ (tolerant) and an indica cultivar ‘Kasalath’ (susceptible). A single QTL near RM221 marker on chromosome 2 was detected by composite interval mapping and additional five QTL were detected by multiple interval mapping. These QTL explained 99% of the total phenotypic variance and they will be used to develop tolerant varieties.

Changes in the Gene Expression of the White Rot Fungus Phanerochaete chrysosporium Due to the Addition of Atropine

We constructed a LongSAGE (Long Serial Analysis of Gene Expression) library from a 3-d culture of Phanerochaete chrysosporium supplemented with atropine, which inhibits the production of lignin-degrading enzymes. The library (the atropine library) contains 13,108 LongSAGE tags and 6,783 unique tags. The gene expression profile represented by the tags was compared with those of two previously constructed libraries, one of which was constructed using 2-d cultures in which the fungus had not yet produced ligninolytic enzymes (the 2-d library) and the other was constructed using 3-d cultures in which the fungus had just started to produce the enzymes (the 3-d library). We found a total of 595 genes that were at least twice more highly or at least twice less highly expressed in the 3-d library than in the 2-d library or the atropine library, and the fluctuations were statistically significant. The relationships among these 595 genes were considered using cluster analysis. Of the 595 genes, 164 showed expression patterns similar to those of four ligninolytic enzyme genes, which were more expressed on day 3 than under any other conditions. Many of these 164 genes comprised genes possibly involved in lignin degradation, lipid metabolism, xenobiotic degradation, stress response, or signal transduction pathways.

Metal accumulation by Arabidopsis halleri subsp. gemmifera at a limestone mining site

Arabidopsis halleri subsp. gemmifera, a close wild relative of A. thaliana in eastern Asia, is an important model species of Brassicaceae used to study cadmium (Cd) and zinc (Zn) tolerance and hyperaccumulation in plants. To investigate the effects of soil factors on metal accumulation in this plant, we collected plants and rhizosphere soil samples from a limestone mining site on Mount Ibuki, and compared them with those collected from non-calcareous soil in Japan. Irrespective of the sampling site, all the plants efficiently accumulated Cd in shoot tissues. The plants growing on non-calcareous soil also accumulated Zn in shoot tissues, but shoot Zn concentration in plants growing on calcareous soils was below the level required for hyperaccumulators. The pH of calcareous soil samples was between 7.68 and 8.21. Total Zn contents were similar in calcareous and non-calcareous soils, but the amounts of Zn extractable by 0.1 M HCl were lower in calcareous than in non-calcareous soil. These results indicate that the properties of calcareous soil affect metal accumulation in plants.

Rapid construction of a high-density rice linkage map by high efficiency genome scanning (HEGS) system.

High-density linkage maps are essential tools for genome analysis of various biological traits. Our developed compact multi-gel system, HEGS (high efficiency genome scanning) is a high-throughput and high-cost-performance electrophoresis apparatus. Using this system, a high-density (average interval 2.3 cM) map with 1 065 AFLP and 63 SSR markers was constructed from recombinant inbred lines of a japonica and indica hybrid in just two months of electrophoreses by a single person. More than 50% of the mapped AFLP markers were commonly polymorphic for several combinations between japonica and indica rice and 15% were applicable for genetically closer crosses between upland and lowland types of japonica rice. This system can be used for rapid analyses of all kinds of markers.