Comprehensive evolutionary analysis and nomenclature of plant G3BPs.

Stress induces extensive reprogramming of mRNA metabolism, which includes the transcription and translation of stress-related genes and the formation of stress granules. RasGAP SH3 domain-binding proteins (G3BPs, also called Rasputins) form a highly conserved family of proteins found throughout eukaryotic evolution, which coordinate signal transduction and posttranscriptional gene regulation and play a key role in the formation of stress granules. G3BPs play a role in osmotic, oxidative, and biotic stress in mammals, and recent results revealed that they play similar functions in higher plants. Although simple eukaryotes such as yeast have only one G3BP gene, higher plants show a massive expansion of their G3BP genes into distinct subfamilies. However, because this family of genes has not been well-characterized in plants, functions that have evolved during this expansion remain unidentified. Therefore, we carried out a phylogenetic analysis of G3BPs in different eukaryotes, particularly focusing on the green lineage. On the basis of this evolutionary analysis of G3BPs in eukaryotes, we propose a uniform nomenclature for plant G3BPs that should help predict the evolutionary and functional diversification in this family.

A cautionary signal from the Red Sea on the impact of increased dust activity on marine microbiota.

BACKGROUND: Global climate change together with growing desertification is leading to increased dust emissions to the atmosphere, drawing attention to possible impacts on marine ecosystems receiving dust deposition. Since microorganisms play important roles in maintaining marine homeostasis through nutrient cycling and carbon flow, detrimental changes in the composition of marine microbiota in response to increased dust input could negatively impact marine health, particularly so in seas located within the Global Dust Belt. Due to its strategic location between two deserts and unique characteristics, the Red Sea provides an attractive semi-enclosed "megacosm" to examine the impacts of large dust deposition on the vastly diverse microbiota in its exceptionally warm oligotrophic waters. RESULTS: We used culture-independent metagenomic approaches to assess temporal changes in the Red Sea microbiota in response to two severe sandstorms, one originated in the Nubian Desert in the summer 2016 and a second one originated in the Libyan Desert in the spring 2017. Despite differences in sandstorm origin and meteorological conditions, both sandstorms shifted bacterial and Archaeal groups in a similar mode. In particular, the relative abundance of autotrophic bacteria declined while those of heterotrophic bacteria, particularly Bacteroidetes, and Archaea increased. The changes peaked within six days from the start of sandstorms, and the community recovered the original assemblage within one month. CONCLUSION: Our results suggest that increased dust emission with expanding desertification could lead to undesirable impacts in ocean function, enhancing heterotrophic processes while reducing autotrophic ones, thereby affecting the marine food web in seas receiving dust deposition.

An alternative, zeaxanthin epoxidase-independent abscisic acid biosynthetic pathway in plants.

Abscisic acid (ABA) is an important carotenoid-derived phytohormone that plays essential roles in plant response to biotic and abiotic stresses as well as in various physiological and developmental processes. In Arabidopsis, ABA biosynthesis starts with the epoxidation of zeaxanthin by the ABA DEFICIENT 1 (ABA1) enzyme, leading to epoxycarotenoids; e.g., violaxanthin. The oxidative cleavage of 9-cis-epoxycarotenoids, a key regulatory step catalyzed by 9-CIS-EPOXYCAROTENOID DIOXYGENASE, forms xanthoxin, which is converted in further reactions mediated by ABA DEFICIENT 2 (ABA2), ABA DEFICIENT 3 (ABA3), and ABSCISIC ALDEHYDE OXIDASE 3 (AAO3) into ABA. By combining genetic and biochemical approaches, we unravel here an ABA1-independent ABA biosynthetic pathway starting upstream of zeaxanthin. We identified the carotenoid cleavage products (i.e., apocarotenoids, ß-apo-11-carotenal, 9-cis-ß-apo-11-carotenal, 3-OH-ß-apo-11-carotenal, and 9-cis-3-OH-ß-apo-11-carotenal) as intermediates of this ABA1-independent ABA biosynthetic pathway. Using labeled compounds, we showed that ß-apo-11-carotenal, 9-cis-ß-apo-11-carotenal, and 3-OH-ß-apo-11-carotenal are successively converted into 9-cis-3-OH-ß-apo-11-carotenal, xanthoxin, and finally into ABA in both Arabidopsis and rice. When applied to Arabidopsis, these ß-apo-11-carotenoids exert ABA biological functions, such as maintaining seed dormancy and inducing the expression of ABA-responsive genes. Moreover, the transcriptomic analysis revealed a high overlap of differentially expressed genes regulated by ß-apo-11-carotenoids and ABA, suggesting that ß-apo-11-carotenoids exert ABA-independent regulatory activities. Taken together, our study identifies a biological function for the common plant metabolites, ß-apo-11-carotenoids, extends our knowledge about ABA biosynthesis, and provides new insights into plant apocarotenoid metabolic networks.

OryzaGenome2.1: Database of Diverse Genotypes in Wild Oryza Species.

BACKGROUND: OryzaGenome ( http://viewer.shigen.info/oryzagenome21detail/index.xhtml ), a feature within Oryzabase ( https://shigen.nig.ac.jp/rice/oryzabase/ ), is a genomic database for wild Oryza species that provides comparative and evolutionary genomics approaches for the rice research community. RESULTS: Here we release OryzaGenome2.1, the first major update of OryzaGenome. The main feature in this version is the inclusion of newly sequenced genotypes and their meta-information, giving a total of 217 accessions of 19 wild Oryza species (O. rufipogon, O. barthii, O. longistaminata, O. meridionalis, O. glumaepatula, O. punctata, O. minuta, O. officinalis, O. rhizomatis, O. eichingeri, O. latifolia, O. alta, O. grandiglumis, O. australiensis, O. brachyantha, O. granulata, O. meyeriana, O. ridleyi, and O. longiglumis). These 19 wild species belong to 9 genome types (AA, BB, CC, BBCC, CCDD, EE, FF, GG, and HHJJ), representing wide genomic diversity in the genus. Using the genotype information, we analyzed the genome diversity of Oryza species. Other features of OryzaGenome facilitate the use of information on single nucleotide polymorphisms (SNPs) between O. sativa and its wild progenitor O. rufipogon in rice research, including breeding as well as basic science. For example, we provide Variant Call Format (VCF) files for genome-wide SNPs of 33 O. rufipogon accessions against the O. sativa reference genome, IRGSP1.0. In addition, we provide a new SNP Effect Table function, allowing users to identify SNPs or small insertion/deletion polymorphisms in the 33 O. rufipogon accessions and to search for the effect of these polymorphisms on protein function if they reside in the coding region (e.g., are missense or nonsense mutations). Furthermore, the SNP Viewer for 446 O. rufipogon accessions was updated by implementing new tracks for possible selective sweep regions and highly mutated regions that were potentially exposed to selective pressures during the process of domestication. CONCLUSION: OryzaGenome2.1 focuses on comparative genomic analysis of diverse wild Oryza accessions collected around the world and on the development of resources to speed up the identification of critical trait-related genes, especially from O. rufipogon. It aims to promote the use of genotype information from wild accessions in rice breeding and potential future crop improvements. Diverse genotypes will be a key resource for evolutionary studies in Oryza, including polyploid biology.

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.

TRANSNAP: a web database providing comprehensive information on Japanese pear transcriptome.

Japanese pear (Pyrus pyrifolia) is a major fruit tree in the family Rosaceae and is bred for fruit production. To promote the development of breeding strategies and molecular research for Japanese pear, we sequenced the transcripts of Japanese pear variety 'Hosui'. To exhaustively collect information of total gene expression, RNA samples from various organs and stages of Japanese pear were sequenced by three technologies, single-molecule real-time (SMRT) sequencing, 454 pyrosequencing, and Sanger sequencing. Using all those reads, we determined comprehensive reference sequences of Japanese pear. Then, their protein sequences were predicted, and biological functional annotations were assigned. Finally, we developed a web database, TRANSNAP (http://plantomics.mind.meiji.ac.jp/nashi), which is the first web resource of Japanese pear omics information. This database provides highly reliable information via a user-friendly web interface: the reference sequences, gene functional annotations, and gene expression profiles from microarray experiments. In addition, based on sequence comparisons among Japanese, Chinese and European pears, similar protein sequences among the pears and species-specific proteins in Japanese pear can be quickly and efficiently identified. TRANSNAP will aid molecular research and breeding in Japanese pear, and its information is available for comparative analysis among other pear species and families.

Complete chloroplast genome sequence and phylogenetic analysis of wasabi (Eutrema japonicum) and its relatives.

In Japan, two Eutrema species, wasabi (Eutrema japonicum, the important traditional Japanese condiment) and yuriwasabi (E. tenue), have been recognized as endemic species. We sequenced complete chloroplast (cp) genomes of seven wasabi and yuriwasabi accessions from Japan to study their phylogeny and evolution, using molecular dating of species divergence. Phylogenetic analyses of the complete cp DNA of these two Japanese species and five other Eurasian Eutrema species revealed that wasabi and yuriwasabi did not form a monophyletic group. One yuriwasabi accession (Gifu) formed a clade with E. yunnanense from China, indicating that this accession should be considered as a different species from the other yuriwasabi accessions. We reveal that Japanese Eutrema species diverged from the 'E. yunnanense-yuriwasabi (Gifu)' clade approximately 1.3 million years ago (Mya), suggesting that the connection between Japan and the Eurasian continent has existed more recently than the Quaternary period. The abundance of cp sequence data in this study also allowed the detection of genetic differentiation among wasabi cultivars. The two polymorphic sites detected between 'Fujidaruma' and 'Shimane No.3' were used to develop genotyping markers. The cp genome information provided here will thus inform the evolutionary histories of Japanese Eutrema species and help in genotyping wasabi cultivars.

Finding Nemo's Genes: A chromosome-scale reference assembly of the genome of the orange clownfish Amphiprion percula.

The iconic orange clownfish, Amphiprion percula, is a model organism for studying the ecology and evolution of reef fishes, including patterns of population connectivity, sex change, social organization, habitat selection and adaptation to climate change. Notably, the orange clownfish is the only reef fish for which a complete larval dispersal kernel has been established and was the first fish species for which it was demonstrated that antipredator responses of reef fishes could be impaired by ocean acidification. Despite its importance, molecular resources for this species remain scarce and until now it lacked a reference genome assembly. Here, we present a de novo chromosome-scale assembly of the genome of the orange clownfish Amphiprion percula. We utilized single-molecule real-time sequencing technology from Pacific Biosciences to produce an initial polished assembly comprised of 1,414 contigs, with a contig N50 length of 1.86 Mb. Using Hi-C-based chromatin contact maps, 98% of the genome assembly were placed into 24 chromosomes, resulting in a final assembly of 908.8 Mb in length with contig and scaffold N50s of 3.12 and 38.4 Mb, respectively. This makes it one of the most contiguous and complete fish genome assemblies currently available. The genome was annotated with 26,597 protein-coding genes and contains 96% of the core set of conserved actinopterygian orthologs. The availability of this reference genome assembly as a community resource will further strengthen the role of the orange clownfish as a model species for research on the ecology and evolution of reef fishes.

Potential of Oryza officinalis to augment the cold tolerance genetic mechanisms of Oryza sativa by network complementation.

Oryza officinalis is an accessible alien donor for genetic improvement of rice. Comparison across a representative panel of Oryza species showed that the wild O. officinalis and cultivated O. sativa ssp. japonica have similar cold tolerance potentials. The possibility that either distinct or similar genetic mechanisms are involved in the low temperature responses of each species was addressed by comparing their transcriptional networks. General similarities were supported by shared transcriptomic signatures indicative of equivalent metabolic, hormonal, and defense status. However, O. officinalis has maintained an elaborate cold-responsive brassinosteroid-regulated BES1-network that appeared to have been fragmented in O. sativa. BES1-network is potentially important for integrating growth-related responses with physiological adjustments and defenses through the protection of photosynthetic machinery and maintenance of stomatal aperture, oxidative defenses, and osmotic adjustment. Equivalent physiological processes are functional in O. sativa but their genetic mechanisms are under the direct control of ABA-dependent, DREB-dependent and/or oxidative-mediated networks uncoupled to BES1. While O. officinalis and O. sativa represent long periods of speciation and domestication, their comparable cold tolerance potentials involve equivalent physiological processes but distinct genetic networks. BES1-network represents a novel attribute of O. officinalis with potential applications in diversifying or complementing other mechanisms in the cultivated germplasm.

Integration of gel-based and gel-free proteomic data for functional analysis of proteins through Soybean Proteome Database.

The Soybean Proteome Database (SPD) stores data on soybean proteins obtained with gel-based and gel-free proteomic techniques. The database was constructed to provide information on proteins for functional analyses. The majority of the data is focused on soybean (Glycine max 'Enrei'). The growth and yield of soybean are strongly affected by environmental stresses such as flooding. The database was originally constructed using data on soybean proteins separated by two-dimensional polyacrylamide gel electrophoresis, which is a gel-based proteomic technique. Since 2015, the database has been expanded to incorporate data obtained by label-free mass spectrometry-based quantitative proteomics, which is a gel-free proteomic technique. Here, the portions of the database consisting of gel-free proteomic data are described. The gel-free proteomic database contains 39,212 proteins identified in 63 sample sets, such as temporal and organ-specific samples of soybean plants grown under flooding stress or non-stressed conditions. In addition, data on organellar proteins identified in mitochondria, nuclei, and endoplasmic reticulum are stored. Furthermore, the database integrates multiple omics data such as genomics, transcriptomics, metabolomics, and proteomics. The SPD database is accessible at http://proteome.dc.affrc.go.jp/Soybean/. BIOLOGICAL SIGNIFICANCE: The Soybean Proteome Database stores data obtained from both gel-based and gel-free proteomic techniques. The gel-free proteomic database comprises 39,212 proteins identified in 63 sample sets, such as different organs of soybean plants grown under flooding stress or non-stressed conditions in a time-dependent manner. In addition, organellar proteins identified in mitochondria, nuclei, and endoplasmic reticulum are stored in the gel-free proteomics database. A total of 44,704 proteins, including 5490 proteins identified using a gel-based proteomic technique, are stored in the SPD. It accounts for approximately 80% of all predicted proteins from genome sequences, though there are over lapped proteins. Based on the demonstrated application of data stored in the database for functional analyses, it is suggested that these data will be useful for analyses of biological mechanisms in soybean. Furthermore, coupled with recent advances in information and communication technology, the usefulness of this database would increase in the analyses of biological mechanisms.

Corrigendum: The genome of Chenopodium quinoa.

This corrects the article DOI: 10.1038/nature21370.

Heap: a highly sensitive and accurate SNP detection tool for low-coverage high-throughput sequencing data.

Recent availability of large-scale genomic resources enables us to conduct so called genome-wide association studies (GWAS) and genomic prediction (GP) studies, particularly with next-generation sequencing (NGS) data. The effectiveness of GWAS and GP depends on not only their mathematical models, but the quality and quantity of variants employed in the analysis. In NGS single nucleotide polymorphism (SNP) calling, conventional tools ideally require more reads for higher SNP sensitivity and accuracy. In this study, we aimed to develop a tool, Heap, that enables robustly sensitive and accurate calling of SNPs, particularly with a low coverage NGS data, which must be aligned to the reference genome sequences in advance. To reduce false positive SNPs, Heap determines genotypes and calls SNPs at each site except for sites at the both ends of reads or containing a minor allele supported by only one read. Performance comparison with existing tools showed that Heap achieved the highest F-scores with low coverage (7X) restriction-site associated DNA sequencing reads of sorghum and rice individuals. This will facilitate cost-effective GWAS and GP studies in this NGS era. Code and documentation of Heap are freely available from https://github.com/meiji-bioinf/heap (29 March 2017, date last accessed) and our web site (http://bioinf.mind.meiji.ac.jp/lab/en/tools.html (29 March 2017, date last accessed)).

The genome of Chenopodium quinoa.

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

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.

Loss of variation of state detected in soybean metabolic and human myelomonocytic leukaemia cell transcriptional networks under external stimuli.

Soybean (Glycine max) is sensitive to flooding stress, and flood damage at the seedling stage is a barrier to growth. We constructed two mathematical models of the soybean metabolic network, a control model and a flooded model, from metabolic profiles in soybean plants. We simulated the metabolic profiles with perturbations before and after the flooding stimulus using the two models. We measured the variation of state that the system could maintain from a state-space description of the simulated profiles. The results showed a loss of variation of state during the flooding response in the soybean plants. Loss of variation of state was also observed in a human myelomonocytic leukaemia cell transcriptional network in response to a phorbol-ester stimulus. Thus, we detected a loss of variation of state under external stimuli in two biological systems, regardless of the regulation and stimulus types. Our results suggest that a loss of robustness may occur concurrently with the loss of variation of state in biological systems. We describe the possible applications of the quantity of variation of state in plant genetic engineering and cell biology. Finally, we present a hypothetical "external stimulus-induced information loss" model of biological systems.

OryzaGenome: Genome Diversity Database of Wild Oryza Species.

The species in the genus Oryza, encompassing nine genome types and 23 species, are a rich genetic resource and may have applications in deeper genomic analyses aiming to understand the evolution of plant genomes. With the advancement of next-generation sequencing (NGS) technology, a flood of Oryza species reference genomes and genomic variation information has become available in recent years. This genomic information, combined with the comprehensive phenotypic information that we are accumulating in our Oryzabase, can serve as an excellent genotype-phenotype association resource for analyzing rice functional and structural evolution, and the associated diversity of the Oryza genus. Here we integrate our previous and future phenotypic/habitat information and newly determined genotype information into a united repository, named OryzaGenome, providing the variant information with hyperlinks to Oryzabase. The current version of OryzaGenome includes genotype information of 446 O. rufipogon accessions derived by imputation and of 17 accessions derived by imputation-free deep sequencing. Two variant viewers are implemented: SNP Viewer as a conventional genome browser interface and Variant Table as a text-based browser for precise inspection of each variant one by one. Portable VCF (variant call format) file or tab-delimited file download is also available. Following these SNP (single nucleotide polymorphism) data, reference pseudomolecules/scaffolds/contigs and genome-wide variation information for almost all of the closely and distantly related wild Oryza species from the NIG Wild Rice Collection will be available in future releases. All of the resources can be accessed through http://viewer.shigen.info/oryzagenome/.

Global expression differences and tissue specific expression differences in rice evolution result in two contrasting types of differentially expressed genes.

BACKGROUND: Since the development of transcriptome analysis systems, many expression evolution studies characterized evolutionary forces acting on gene expression, without explicit discrimination between global expression differences and tissue specific expression differences. However, different types of gene expression alteration should have different effects on an organism, the evolutionary forces that act on them might be different, and different types of genes might show different types of differential expression between species. To confirm this, we studied differentially expressed (DE) genes among closely related groups that have extensive gene expression atlases, and clarified characteristics of different types of DE genes including the identification of regulating loci for differential expression using expression quantitative loci (eQTL) analysis data. RESULTS: We detected differentially expressed (DE) genes between rice subspecies in five homologous tissues that were verified using japonica and indica transcriptome atlases in public databases. Using the transcriptome atlases, we classified DE genes into two types, global DE genes and changed-tissues DE genes. Global type DE genes were not expressed in any tissues in the atlas of one subspecies, however changed-tissues type DE genes were expressed in both subspecies with different tissue specificity. For the five tissues in the two japonica-indica combinations, 4.6 ± 0.8 and 5.9 ± 1.5 % of highly expressed genes were global and changed-tissues DE genes, respectively. Changed-tissues DE genes varied in number between tissues, increasing linearly with the abundance of tissue specifically expressed genes in the tissue. Molecular evolution of global DE genes was rapid, unlike that of changed-tissues DE genes. Based on gene ontology, global and changed-tissues DE genes were different, having no common GO terms. Expression differences of most global DE genes were regulated by cis-eQTLs. Expression evolution of changed-tissues DE genes was rapid in tissue specifically expressed genes and those rapidly evolved changed-tissues DE genes were regulated not by cis-eQTLs, but by complicated trans-eQTLs. CONCLUSIONS: Global DE genes and changed-tissues DE genes had contrasting characteristics. The two contrasting types of DE genes provide possible explanations for the previous controversial conclusions about the relationships between molecular evolution and expression evolution of genes in different species, and the relationship between expression breadth and expression conservation in evolution.