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.

Corrigendum: The genome of Chenopodium quinoa.

This corrects the article DOI: 10.1038/nature21370.

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.

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.

Genome-wide study of KNOX regulatory network reveals brassinosteroid catabolic genes important for shoot meristem function in rice.

In flowering plants, knotted1-like homeobox (KNOX) transcription factors play crucial roles in establishment and maintenance of the shoot apical meristem (SAM), from which aerial organs such as leaves, stems, and flowers initiate. We report that a rice (Oryza sativa) KNOX gene Oryza sativa homeobox1 (OSH1) represses the brassinosteroid (BR) phytohormone pathway through activation of BR catabolism genes. Inducible overexpression of OSH1 caused BR insensitivity, whereas loss of function showed a BR-overproduction phenotype. Genome-wide identification of loci bound and regulated by OSH1 revealed hormonal and transcriptional regulation as the major function of OSH1. Among these targets, BR catabolism genes CYP734A2, CYP734A4, and CYP734A6 were rapidly upregulated by OSH1 induction. Furthermore, RNA interference knockdown plants of CYP734A genes arrested growth of the SAM and mimicked some osh1 phenotypes. Thus, we suggest that local control of BR levels by KNOX genes is a key regulatory step in SAM function.

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/.

Rice germline-specific Argonaute MEL1 protein binds to phasiRNAs generated from more than 700 lincRNAs.

Small RNAs that interact with Argonaute (AGO) proteins play central roles in RNA-mediated silencing. MEIOSIS ARRESTED AT LEPTOTENE1 (MEL1), a rice AGO, has specific functions in the development of pre-meiotic germ cells and the progression of meiosis. Here, we show that MEL1, which is located mostly in the cytoplasm of germ cells, associates preferentially with 21-nucleotide phased small interfering RNAs (phasiRNAs) that bear a 5'-terminal cytosine. Most phasiRNAs are derived from 1171 intergenic clusters distributed on all rice chromosomes. From these clusters, over 700 large intergenic, non-coding RNAs (lincRNAs) that contain the consensus sequence complementary to miR2118 are transcribed specifically in inflorescences, and cleaved within the miR2118 site. Cleaved lincRNAs are processed via DICER-LIKE4 (DCL4) protein, resulting in production of phasiRNAs. This study provides the evidence that the miR2118-dependent and the DCL4-dependent pathways are both required for biogenesis of 21-nt phasiRNAs associated with germline-specific MEL1 AGO in rice, and over 700 lincRNAs are key factors for induction of this biogenesis during reproductive-specific stages.

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.

RiTE database: a resource database for genus-wide rice genomics and evolutionary biology.

BACKGROUND: Comparative evolutionary analysis of whole genomes requires not only accurate annotation of gene space, but also proper annotation of the repetitive fraction which is often the largest component of most if not all genomes larger than 50 kb in size. RESULTS: Here we present the Rice TE database (RiTE-db)--a genus-wide collection of transposable elements and repeated sequences across 11 diploid species of the genus Oryza and the closely-related out-group Leersia perrieri. The database consists of more than 170,000 entries divided into three main types: (i) a classified and curated set of publicly-available repeated sequences, (ii) a set of consensus assemblies of highly-repetitive sequences obtained from genome sequencing surveys of 12 species; and (iii) a set of full-length TEs, identified and extracted from 12 whole genome assemblies. CONCLUSIONS: This is the first report of a repeat dataset that spans the majority of repeat variability within an entire genus, and one that includes complete elements as well as unassembled repeats. The database allows sequence browsing, downloading, and similarity searches. Because of the strategy adopted, the RiTE-db opens a new path to unprecedented direct comparative studies that span the entire nuclear repeat content of 15 million years of Oryza diversity.

De novo transcriptome assembly of a fern, Lygodium japonicum, and a web resource database, Ljtrans DB.

During plant evolution, ferns originally evolved as a major vascular plant with a distinctive life cycle in which the haploid and diploid generations are completely separated. However, the low level of genetic resources has limited studies of their physiological events, as well as hindering research on the evolutionary history of land plants. In this study, to identify a comprehensive catalog of transcripts and characterize their expression traits in the fern Lygodium japonicum, nine different RNA samples isolated from prothalli, trophophylls, rhizomes and sporophylls were sequenced using Roche 454 GS-FLX and Illumina HiSeq sequencers. The hybrid assembly of the high-quality 454 GS-FLX and Illumina HiSeq reads generated a set of 37,830 isoforms with an average length of 1,444 bp. Using four open reading frame (ORF) predictors, 38,142 representative ORFs were identified from a total of 37,830 transcript isoforms and 95 contigs, which were annotated by searching against several public databases. Furthermore, an orthoMCL analysis using the protein sequences of L. japonicum and five model plants revealed various sets of lineage-specific genes, including those detected among land plant lineages and those detected in only L. japonicum. We have also examined the expression patterns of all contigs/isoforms, along with the life cycle of L. japonicum, and identified the tissue-specific transcripts using statistical expression analyses. Finally, we developed a public web resource, the L. japonicum transcriptome database at http://bioinf.mind.meiji.ac.jp/kanikusa/, which provides important opportunities to accelerate molecular research in ferns.

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.

Rapid turnover of antimicrobial-type cysteine-rich protein genes in closely related Oryza genomes.

Defensive and reproductive protein genes undergo rapid evolution. Small, cysteine-rich secreted peptides (CRPs) act as antimicrobial agents and function in plant intercellular signaling and are over-represented among reproductively expressed proteins. Because of their roles in defense, reproduction and development and their presence in multigene families, CRP variation can have major consequences for plant phenotypic and functional diversification. We surveyed the CRP genes of six closely related Oryza genomes comprising Oryza sativa ssp. japonica and ssp. indica, Oryza glaberrima and three accessions of Oryza rufipogon to observe patterns of evolution in these gene families and the effects of variation on their gene expression. These Oryza genomes, like other plant genomes, have accumulated large reservoirs of CRP sequences, comprising 26 groups totaling between 676 and 843 genes, in contrast to antimicrobial CRPs in animal genomes. Despite the close evolutionary relationships between the genomes, we observed rapid changes in number and structure among CRP gene families. Many CRP sequences are in gene clusters generated by local duplications, have undergone rapid turnover and are more likely to be silent or specifically expressed. By contrast, conserved CRP genes are more likely to be highly and broadly expressed. Variable CRP genes created by repeated duplication, gene modification and inactivation can gain new functions and expression patterns in newly evolved gene copies. For the CRP proteins, the process of gain/loss by deletion or duplication at gene clusters seems to be an important mechanism in evolution of the gene families, which also contributes to their expression evolution.

RiceXPro version 3.0: expanding the informatics resource for rice transcriptome.

A wide range of resources on gene expression profiling enhance various strategies in plant molecular biology particularly in characterization of gene function. We have updated our gene expression profile database, RiceXPro (http://ricexpro.dna.affrc.go.jp/), to provide more comprehensive information on the transcriptome of rice encompassing the entire growth cycle and various experimental conditions. The gene expression profiles are currently grouped into three categories, namely, 'field/development' with 572 data corresponding to 12 data sets, 'plant hormone' with 143 data corresponding to 13 data sets and 'cell- and tissue-type' comprising of 38 microarray data. In addition to the interface for retrieving expression information of a gene/genes in each data set, we have incorporated an interface for a global approach in searching an overall view of the gene expression profiles from multiple data sets within each category. Furthermore, we have also added a BLAST search function that enables users to explore expression profile of a gene/genes with similarity to a query sequence. Therefore, the updated version of RiceXPro can be used more efficiently to survey the gene expression signature of rice in sufficient depth and may also provide clues on gene function of other cereal crops.

RiceFREND: a platform for retrieving coexpressed gene networks in rice.

Similarity of gene expression across a wide range of biological conditions can be efficiently used in characterization of gene function. We have constructed a rice gene coexpression database, RiceFREND (http://ricefrend.dna.affrc.go.jp/), to identify gene modules with similar expression profiles and provide a platform for more accurate prediction of gene functions. Coexpression analysis of 27 201 genes was performed against 815 microarray data derived from expression profiling of various organs and tissues at different developmental stages, mature organs throughout the growth from transplanting until harvesting in the field and plant hormone treatment conditions, using a single microarray platform. The database is provided with two search options, namely, 'single guide gene search' and 'multiple guide gene search' to efficiently retrieve information on coexpressed genes. A user-friendly web interface facilitates visualization and interpretation of gene coexpression networks in HyperTree, Cytoscape Web and Graphviz formats. In addition, analysis tools for identification of enriched Gene Ontology terms and cis-elements provide clue for better prediction of biological functions associated with the coexpressed genes. These features allow users to clarify gene functions and gene regulatory networks that could lead to a more thorough understanding of many complex agronomic traits.

Development of INDEL markers to discriminate all genome types rapidly in the genus Oryza.

The wild Oryza species are rich in genetic diversity and are good resources for modern breeding of rice varieties. The reliable ex situ conservation of various genetic resources supports both basic and applied rice research. For this purpose, we developed PCR-based and co-dominant insertion/deletion (INDEL) markers which enable the discrimination of the genome types or species in the genus Oryza. First, 12,107 INDEL candidate sequences were found in the BAC end sequences for 12 Oryza species available in public databases. Next, we designed PCR primers for INDEL-flanking sequences to match the characteristics of each INDEL, based on an assessment of their likelihood to give rise to a single or few PCR products in all 102 wild accessions, covering most Oryza genome types. Then, we selected 22 INDEL markers to discriminate all genome types in the genus Oryza. A phylogenetic tree of 102 wild accessions and two cultivars according to amplicon polymorphisms for the 22 INDEL markers corresponded well to those in previous studies, indicating that the INDEL markers developed in this study were a useful tool to improve the reliability of identification of wild Oryza species in the germplasm stocks.

DaizuBase, an integrated soybean genome database including BAC-based physical maps.

Soybean [Glycine max (L) Merrill] is one of the most important leguminous crops and ranks fourth after to rice, wheat and maize in terms of world crop production. Soybean contains abundant protein and oil, which makes it a major source of nutritious food, livestock feed and industrial products. In Japan, soybean is also an important source of traditional staples such as tofu, natto, miso and soy sauce. The soybean genome was determined in 2010. With its enormous size, physical mapping and genome sequencing are the most effective approaches towards understanding the structure and function of the soybean genome. We constructed bacterial artificial chromosome (BAC) libraries from the Japanese soybean cultivar, Enrei. The end-sequences of approximately 100,000 BAC clones were analyzed and used for construction of a BAC-based physical map of the genome. BLAST analysis between Enrei BAC-end sequences and the Williams82 genome was carried out to increase the saturation of the map. This physical map will be used to characterize the genome structure of Japanese soybean cultivars, to develop methods for the isolation of agronomically important genes and to facilitate comparative soybean genome research. The current status of physical mapping of the soybean genome and construction of database are presented.

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.

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.