Correction: EAT1 transcription factor, a non-cell-autonomous regulator of pollen production, activates meiotic small RNA biogenesis in rice anther tapetum.

[This corrects the article DOI: 10.1371/journal.pgen.1007238.].

EAT1 transcription factor, a non-cell-autonomous regulator of pollen production, activates meiotic small RNA biogenesis in rice anther tapetum.

The 24-nucleotides (nt) phased secondary small interfering RNA (phasiRNA) is a unique class of plant small RNAs abundantly expressed in monocot anthers at early meiosis. Previously, 44 intergenic regions were identified as the loci for longer precursor RNAs of 24-nt phasiRNAs (24-PHASs) in the rice genome. However, the regulatory mechanism that determines spatiotemporal expression of these RNAs has remained elusive. ETERNAL TAPETUM1 (EAT1) is a basic-helix-loop-helix (bHLH) transcription factor indispensable for induction of programmed cell death (PCD) in postmeiotic anther tapetum, the somatic nursery for pollen production. In this study, EAT1-dependent non-cell-autonomous regulation of male meiosis was evidenced from microscopic observation of the eat1 mutant, in which meiosis with aberrantly decondensed chromosomes was retarded but accomplished somehow, eventually resulting in abortive microspores due to an aberrant tapetal PCD. EAT1 protein accumulated in tapetal-cell nuclei at early meiosis and postmeiotic microspore stages. Meiotic EAT1 promoted transcription of 24-PHAS RNAs at 101 loci, and importantly, also activated DICER-LIKE5 (DCL5, previous DCL3b in rice) mRNA transcription that is required for processing of double-stranded 24-PHASs into 24-nt lengths. From the results of the chromatin-immunoprecipitation and transient expression analyses, another tapetum-expressing bHLH protein, TDR INTERACTING PROTEIN2 (TIP2), was suggested to be involved in meiotic small-RNA biogenesis. The transient assay also demonstrated that UNDEVELOPED TAPETUM1 (UDT1)/bHLH164 is a potential interacting partner of both EAT1 and TIP2 during early meiosis. This study indicates that EAT1 is one of key regulators triggering meiotic phasiRNA biogenesis in anther tapetum, and that other bHLH proteins, TIP2 and UDT1, also play some important roles in this process. Spatiotemporal expression control of these bHLH proteins is a clue to orchestrate precise meiosis progression and subsequent pollen production non-cell-autonomously.

Detection of novel QTLs qDTH4.5 and qDTH6.3, which confer late heading under short-day conditions, by SSR marker-based and QTL-seq analysis.

Heading date is one of the most important traits in rice breeding. It is governed by multiple genes, including known quantitative trait loci (QTLs). In general, almost all japonica cultivars, including Nipponbare, head early under short-day (SD) conditions, but some indica cultivars, including Kasalath, head late. To explain this difference, we identified QTLs controlling heading date under SD conditions. We used NILs, CSSLs, and BILs from a cross between Nipponbare and Kasalath, and evaluated days to heading (DTH) under SD conditions. No NILs or CSSLs showed late heading, but two BILs (BIL-55 and BIL-78) had almost the same DTH as Kasalath. We developed an F2 population from a cross between BIL-55 and Nipponbare and performed QTL analysis using SSR markers. The late-heading phenotype was controlled by two known genes and at least two novel QTLs on chromosomes 4 and 6, named qDTH4.5 and qDTH6.3. These QTLs were confirmed by QTL-seq. The QTLs and polymorphisms detected here will provide useful information for further genetic studies and breeding under SD conditions at lower latitudes.

A BAC physical map of aus rice cultivar 'Kasalath', and the map-based genomic sequence of 'Kasalath' chromosome 1.

Comparative analysis using available genomic resources within closely related species is an effective way to investigate genomic sequence and structural diversity. Rice (Oryza sativa L.) has undergone significant physiological and morphological changes during its domestication and local adaptation. We present a complete bacterial artificial chromosome (BAC) physical map for the aus rice cultivar 'Kasalath', which covers 90% of the sequence of temperate japonica rice cultivar 'Nipponbare'. Examination of physical distances between computational and experimental measurements of 'Kasalath' BAC insert size revealed the presence of more than 500 genomic regions that appear to have significant chromosome structural changes between the two cultivars. In particular, a genomic region on the long arm of 'Kasalath' chromosome 11 carrying a disease-resistance gene cluster was greatly expanded relative to the 'Nipponbare' genome. We also decoded 41.37 Mb of high-quality genomic sequence from 'Kasalath' chromosome 1. Extensive comparisons of chromosome 1 between 'Kasalath' and 'Nipponbare' led to the discovery of 317,843 single-nucleotide polymorphisms (SNPs) and 66,331 insertion/deletion (indel) sites. Nearly two-thirds of the expressed genes on rice chromosome 1 carried natural variations involving SNPs and/or indels that resulted in substitutions, insertions or deletions of amino acids in one cultivar relative to the other. We also observed gain and loss of genes caused by large indels. This study provides an important framework and an invaluable dataset for further understanding of the molecular mechanisms underlying the evolution and functions of the rice genome.

Mapping of seed shattering loci provides insights into origin of weedy rice and rice domestication.

Seed shattering is an important trait that distinguishes crop cultivars from the wild and weedy species. The genetics of seed shattering was investigated in this study to provide insights into rice domestication and the evolution of weedy rice. Quantitative trait locus (QTL) analysis, conducted in 2 recombinant inbred populations involving 2 rice cultivars and a weedy rice accession of the southern United States, revealed 3-5 QTLs that controlled seed shattering with 38-45% of the total phenotypic variation. Two QTLs on chromosomes 4 and 10 were consistent in both populations. Both cultivar and weedy rice contributed alleles for increased seed shattering. Genetic backgrounds affected both QTL number and the magnitude of QTL effects. The major QTL qSH4 and a minor QTL qSH3 were validated in near-isogenic lines, with the former conferring a significantly higher degree of seed shattering than the latter. Although the major QTL qSH4 overlapped with the sh4, the presence of the nonshattering single nucleotide polymorphism allele in the weedy rice accession suggested involvement of a linked locus or an alternative molecular genetic mechanism. Overlapping of several QTLs with those from earlier studies indicated that weedy rice may have been derived from the wild species Oryza rufipogon. Natural hybridization of rice cultivars with the highly variable O. rufipogon present in different geographic regions might be responsible for the evolution of a wide range of phenotypic and genotypic variabilities seen in weedy rice populations worldwide.

Rice Annotation Project Database (RAP-DB): an integrative and interactive database for rice genomics.

The Rice Annotation Project Database (RAP-DB, http://rapdb.dna.affrc.go.jp/) has been providing a comprehensive set of gene annotations for the genome sequence of rice, Oryza sativa (japonica group) cv. Nipponbare. Since the first release in 2005, RAP-DB has been updated several times along with the genome assembly updates. Here, we present our newest RAP-DB based on the latest genome assembly, Os-Nipponbare-Reference-IRGSP-1.0 (IRGSP-1.0), which was released in 2011. We detected 37,869 loci by mapping transcript and protein sequences of 150 monocot species. To provide plant researchers with highly reliable and up to date rice gene annotations, we have been incorporating literature-based manually curated data, and 1,626 loci currently incorporate literature-based annotation data, including commonly used gene names or gene symbols. Transcriptional activities are shown at the nucleotide level by mapping RNA-Seq reads derived from 27 samples. We also mapped the Illumina reads of a Japanese leading japonica cultivar, Koshihikari, and a Chinese indica cultivar, Guangluai-4, to the genome and show alignments together with the single nucleotide polymorphisms (SNPs) and gene functional annotations through a newly developed browser, Short-Read Assembly Browser (S-RAB). We have developed two satellite databases, Plant Gene Family Database (PGFD) and Integrative Database of Cereal Gene Phylogeny (IDCGP), which display gene family and homologous gene relationships among diverse plant species. RAP-DB and the satellite databases offer simple and user-friendly web interfaces, enabling plant and genome researchers to access the data easily and facilitating a broad range of plant research topics.

Distinct evolutionary patterns of Oryza glaberrima deciphered by genome sequencing and comparative analysis.

Here we present the genomic sequence of the African cultivated rice, Oryza glaberrima, and compare these data with the genome sequence of Asian cultivated rice, Oryza sativa. We obtained gene-enriched sequences of O. glaberrima that correspond to about 25% of the gene regions of the O. sativa (japonica) genome by methylation filtration and subtractive hybridization of repetitive sequences. While patterns of amino acid changes did not differ between the two species in terms of the biochemical properties, genes of O. glaberrima generally showed a larger synonymous-nonsynonymous substitution ratio, suggesting that O. glaberrima has undergone a genome-wide relaxation of purifying selection. We further investigated nucleotide substitutions around splice sites and found that eight genes of O. sativa experienced changes at splice sites after the divergence from O. glaberrima. These changes produced novel introns that partially truncated functional domains, suggesting that these newly emerged introns affect gene function. We also identified 2451 simple sequence repeats (SSRs) from the genomes of O. glaberrima and O. sativa. Although tri-nucleotide repeats were most common among the SSRs and were overrepresented in the protein-coding sequences, we found that selection against indels of tri-nucleotide repeats was relatively weak in both African and Asian rice. Our genome-wide sequencing of O. glaberrima and in-depth analyses provide rice researchers not only with useful genomic resources for future breeding but also with new insights into the genomic evolution of the African and Asian rice species.

DDM1 (decrease in DNA methylation) genes in rice (Oryza sativa).

Regulation of cytosine methylation in the plant genome is of pivotal in determining the epigenetic states of chromosome regions. Relative tolerance of plant to deficiency in cytosine methylation provides unparalleled opportunities to study the mechanism for regulation of cytosine methylation. The Decrease in DNA Methylation 1 (DDM1) of Arabidopsis thaliana is one of the best characterized plant epigenetic regulators that are necessary for maintenance of cytosine methylation in genomic DNA. Although cytosine methylation could affect various aspects of plant growth and development including those related to agricultural importance, orthologs of DDM1 in plants other than Arabidopsis has not been studied in detail. In this study, we identified two rice genes with similarity to Arabidopsis DDM1 and designated them OsDDM1a and OsDDM1b. Both of the rice DDM1 homologs are transcribed during development and their amino acid sequences are 93 % identical to each other. Transgenic rice lines expressing the OsDDM1a cDNA in the antisense orientation exhibited genomic DNA hypomethylation. In those lines, repeated sequences were more severely affected than a single copy sequence as is the case in Arabidopsis ddm1 mutants. Transcripts derived from endogenous transposon-related loci were up-regulated in the antisense OsDDM1 lines, opening a possibility to identify and utilize potentially active transposons for rice functional genomics.

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.

The role of casein kinase II in flowering time regulation has diversified during evolution.

Casein kinase II (CK2) is a protein kinase with an evolutionarily conserved function as a circadian clock component in several organisms, including the long-day plant Arabidopsis (Arabidopsis thaliana). The circadian clock component CIRCADIAN CLOCK ASSOCIATED1 (CCA1) is a CK2 target in Arabidopsis, where it influences photoperiodic flowering. In rice (Oryza sativa), a short-day plant, Heading date6 (Hd6) encodes a CK2alpha subunit that delays flowering time under long-day conditions. Here, we demonstrate that control of flowering time in rice by the Hd6 CK2alpha subunit requires a functional Hd1 gene (an Arabidopsis CONSTANS ortholog) and is independent of the circadian clock mechanism. Our findings from overexpressing the dominant-negative CK2 allele in rice support the independence of CK2 function from the circadian clock. This lack of control of the circadian clock by Hd6 CK2alpha might be due to the presence of glutamate in OsLHY (a CCA1 ortholog in rice) instead of the serine at the corresponding CK2 target site in CCA1. However, this glutamate is critical for the control of the OsPRR1 gene (a rice ortholog of the Arabidopsis TOC1/PRR1 gene) by OsLHY for regulation of the circadian clock. We also demonstrated that the other conserved CK2 target sites in OsLHY conferred robust rhythmic expression of OsLHY-LUC under diurnal conditions. These findings imply that the role of CK2 in flowering-time regulation in higher plants has diversified during evolution.

Comparative analysis of complete orthologous centromeres from two subspecies of rice reveals rapid variation of centromere organization and structure.

Centromeres are sites for assembly of the chromosomal structures that mediate faithful segregation at mitosis and meiosis. This function is conserved across species, but the DNA components that are involved in kinetochore formation differ greatly, even between closely related species. To shed light on the nature, evolutionary timing and evolutionary dynamics of rice centromeres, we decoded a 2.25-Mb DNA sequence covering the centromeric region of chromosome 8 of an indica rice variety, 'Kasalath' (Kas-Cen8). Analysis of repetitive sequences in Kas-Cen8 led to the identification of 222 long terminal repeat (LTR)-retrotransposon elements and 584 CentO satellite monomers, which account for 59.2% of the region. A comparison of the Kas-Cen8 sequence with that of japonica rice 'Nipponbare' (Nip-Cen8) revealed that about 66.8% of the Kas-Cen8 sequence was collinear with that of Nip-Cen8. Although the 27 putative genes are conserved between the two subspecies, only 55.4% of the total LTR-retrotransposon elements in 'Kasalath' had orthologs in 'Nipponbare', thus reflecting recent proliferation of a considerable number of LTR-retrotransposons since the divergence of two rice subspecies of indica and japonica within Oryza sativa. Comparative analysis of the subfamilies, time of insertion, and organization patterns of inserted LTR-retrotransposons between the two Cen8 regions revealed variations between 'Kasalath' and 'Nipponbare' in the preferential accumulation of CRR elements, and the expansion of CentO satellite repeats within the core domain of Cen8. Together, the results provide insights into the recent proliferation of LTR-retrotransposons, and the rapid expansion of CentO satellite repeats, underlying the dynamic variation and plasticity of plant centromeres.

The Rice Annotation Project Database (RAP-DB): 2008 update.

The Rice Annotation Project Database (RAP-DB) was created to provide the genome sequence assembly of the International Rice Genome Sequencing Project (IRGSP), manually curated annotation of the sequence, and other genomics information that could be useful for comprehensive understanding of the rice biology. Since the last publication of the RAP-DB, the IRGSP genome has been revised and reassembled. In addition, a large number of rice-expressed sequence tags have been released, and functional genomics resources have been produced worldwide. Thus, we have thoroughly updated our genome annotation by manual curation of all the functional descriptions of rice genes. The latest version of the RAP-DB contains a variety of annotation data as follows: clone positions, structures and functions of 31 439 genes validated by cDNAs, RNA genes detected by massively parallel signature sequencing (MPSS) technology and sequence similarity, flanking sequences of mutant lines, transposable elements, etc. Other annotation data such as Gnomon can be displayed along with those of RAP for comparison. We have also developed a new keyword search system to allow the user to access useful information. The RAP-DB is available at: http://rapdb.dna.affrc.go.jp/ and http://rapdb.lab.nig.ac.jp/.

Evolutionary dynamics and impacts of chromosome regions carrying R-gene clusters in rice.

To elucidate R-gene evolution, we compared the genomic compositions and structures of chromosome regions carrying R-gene clusters among cultivated and wild rice species. Map-based sequencing and gene annotation of orthologous genomic regions (1.2 to 1.9 Mb) close to the terminal end of the long arm of rice chromosome 11 revealed R-gene clusters within six cultivated and ancestral wild rice accessions. NBS-LRR R-genes were much more abundant in Asian cultivated rice (O. sativa L.) than in its ancestors, indicating that homologs of functional genes involved in the same pathway likely increase in number because of tandem duplication of chromosomal segments and were selected during cultivation. Phylogenetic analysis using amino acid sequences indicated that homologs of paired Pikm1-Pikm2 (NBS-LRR) genes conferring rice-blast resistance were likely conserved among all cultivated and wild rice species we examined, and the homolog of Xa3/Xa26 (LRR-RLK) conferring bacterial blight resistance was lacking only in Kasalath.

Chromosome-specific distribution of nucleotide substitutions in telomeric repeats of rice (Oryza sativa L.).

Examination of the genomic sequence of the telomere region makes it possible to understand the evolution of the structure of chromosomal ends. We compared the genomic sequences of 14 chromosomal ends of rice, Oryza sativa, L., on the basis of the variation in TTTAGGG repeats. In the proximal telomere repeats, nucleotide substitution occurred more frequently than in the more distal repeats. The most significant diversity was observed at the 1st, 2nd, or 3rd position of TTTAGGG, suggesting that T has been a target of mutation preferentially. Copies of ATTAGGG, CTTAGGG, GTTAGGG, TTCAGGG, TTGAGGG, or TATAGGG were arrayed in tandem, or the same subtypes were located close to each other. The substituted variants were accumulated in chromosomes 2L, 3L, 7L, and 10S but not in the ends of the other chromosomes. In contrast, deletion variants, almost all of which were TTTAGGG to TTAGGG, were dispersed over approximately 4.9% of the sequenced telomere repeats. In summary, the rice proximal telomeric arrays were composed of blocks of at least 6 types of substituted variants and the canonical sequence in a chromosome-specific manner. These results suggest that the variants might arise from the rapid expansion of a single mutation rather than from the gradual accumulation of random mutations.

Insights into the Musa genome: syntenic relationships to rice and between Musa species.

BACKGROUND: Musa species (Zingiberaceae, Zingiberales) including bananas and plantains are collectively the fourth most important crop in developing countries. Knowledge concerning Musa genome structure and the origin of distinct cultivars has greatly increased over the last few years. Until now, however, no large-scale analyses of Musa genomic sequence have been conducted. This study compares genomic sequence in two Musa species with orthologous regions in the rice genome. RESULTS: We produced 1.4 Mb of Musa sequence from 13 BAC clones, annotated and analyzed them along with 4 previously sequenced BACs. The 443 predicted genes revealed that Zingiberales genes share GC content and distribution characteristics with eudicot and Poaceae genomes. Comparison with rice revealed microsynteny regions that have persisted since the divergence of the Commelinid orders Poales and Zingiberales at least 117 Mya. The previously hypothesized large-scale duplication event in the common ancestor of major cereal lineages within the Poaceae was verified. The divergence time distributions for Musa-Zingiber (Zingiberaceae, Zingiberales) orthologs and paralogs provide strong evidence for a large-scale duplication event in the Musa lineage after its divergence from the Zingiberaceae approximately 61 Mya. Comparisons of genomic regions from M. acuminata and M. balbisiana revealed highly conserved genome structure, and indicated that these genomes diverged circa 4.6 Mya. CONCLUSION: These results point to the utility of comparative analyses between distantly-related monocot species such as rice and Musa for improving our understanding of monocot genome evolution. Sequencing the genome of M. acuminata would provide a strong foundation for comparative genomics in the monocots. In addition a genome sequence would aid genomic and genetic analyses of cultivated Musa polyploid genotypes in research aimed at localizing and cloning genes controlling important agronomic traits for breeding purposes.

The Rice Annotation Project Database (RAP-DB): hub for Oryza sativa ssp. japonica genome information.

With the completion of the rice genome sequencing, a standardized annotation is necessary so that the information from the genome sequence can be fully utilized in understanding the biology of rice and other cereal crops. An annotation jamboree was held in Japan with the aim of annotating and manually curating all the genes in the rice genome. Here we present the Rice Annotation Project Database (RAP-DB), which has been developed to provide access to the annotation data. The RAP-DB has two different types of annotation viewers, BLAST and BLAT search, and other useful features. By connecting the annotations to other rice genomics data, such as full-length cDNAs and Tos17 mutant lines, the RAP-DB serves as a hub for rice genomics. All of the resources can be accessed through http://rapdb.lab.nig.ac.jp/.

Rice Annotation Database (RAD): a contig-oriented database for map-based rice genomics.

A contig-oriented database for annotation of the rice genome has been constructed to facilitate map-based rice genomics. The Rice Annotation Database has the following functional features: (i) extensive effort of manual annotations of P1-derived artificial chromosome/bacterial artificial chromosome clones can be merged at chromosome and contig-level; (ii) concise visualization of the annotation information such as the predicted genes, results of various prediction programs (RiceHMM, Genscan, Genscan+, Fgenesh, GeneMark, etc.), homology to expressed sequence tag, full-length cDNA and protein; (iii) user-friendly clone / gene query system; (iv) download functions for nucleotide, amino acid and coding sequences; (v) analysis of various features of the genome (GC-content, average value, etc.); and (vi) genome-wide homology search (BLAST) of contig- and chromosome-level genome sequence to allow comparative analysis with the genome sequence of other organisms. As of October 2004, the database contains a total of 215 Mb sequence with relevant annotation results including 30 000 manually curated genes. The database can provide the latest information on manual annotation as well as a comprehensive structural analysis of various features of the rice genome. The database can be accessed at http://rad.dna.affrc.go.jp/.