Effective use of legacy data in a genome-wide association studies improves the credibility of quantitative trait loci detection in rice.

Genome-wide association studies (GWASs) are used to detect quantitative trait loci (QTL) using genomic and phenotypic data as inputs. While genomic data are obtained with high throughput and low cost, obtaining phenotypic data requires a large amount of effort and time. In past breeding programs, researchers and breeders have conducted a large number of phenotypic surveys and accumulated results as legacy data. In this study, we conducted a GWAS using phenotypic data of temperate japonica rice (Oryza sativa) varieties from a public database. The GWAS using the legacy data detected several known agriculturally important genes, indicating reliability of the legacy data for GWAS. By comparing the GWAS using legacy data (L-GWAS) and a GWAS using phenotypic data that we measured (M-GWAS), we detected reliable QTL for agronomically important traits. These results suggest that an L-GWAS is a strong alternative to replicate tests to confirm the reproducibility of QTL detected by an M-GWAS. In addition, because legacy data have often been accumulated for many traits, it is possible to evaluate the pleiotropic effect of the QTL identified for the specific trait that we focused on with respect to various other traits. This study demonstrates the effectiveness of using legacy data for GWASs and proposes the use of legacy data to accelerate genomic breeding.

Genome-wide association study of leaf photosynthesis using a high-throughput gas exchange system in rice.

Enhancing leaf photosynthetic capacity is essential for improving the yield of rice (Oryza sativa L.). Although the exploitation of natural genetic resources is considered a promising approach to enhance photosynthetic capacity, genomic factors related to the genetic diversity of leaf photosynthetic capacity have yet to be fully elucidated due to the limitation of measurement efficiency. In this study, we aimed to identify novel genomic regions for the net CO2 assimilation rate (A) by combining genome-wide association study (GWAS) and the newly developed rapid closed gas exchange system MIC-100. Using three MIC-100 systems in the field at the vegetative stage, we measured A of 168 temperate japonica rice varieties with six replicates for three years. We found that the modern varieties exhibited higher A than the landraces, while there was no significant relationship between the release year and A among the modern varieties. Our GWAS scan revealed two major peaks located on chromosomes 4 and 8, which were repeatedly detected in the different experiments and in the generalized linear modelling approach. We suggest that high-throughput gas exchange measurements combined with GWAS is a reliable approach for understanding the genetic mechanisms underlying photosynthetic diversities in crop species.

How have breeders adapted rice flowering to the growing region?

Flowering time is a crucial rice trait that influences its adaptation to various environments, cropping schedules, and agronomic characteristics. Rice breeders have exploited spontaneous mutations in heading date genes to regulate the flowering time. In the present study, we investigated how breeders in Fukui Prefecture regulated days to heading while developing promising rice varieties. Genome-wide association studies (GWAS) identified Hd1, Hd16, and Hd18 as the major genes controlling days to heading in the population. However, we suspected that this highly bred population might exhibit genomic stratification, which could lead to spurious or false correlations in the GWAS. Thus, we also conducted correlation and partial correlation analyses, which uncovered another key heading date gene, Hd17, that GWAS failed to detect because of its linkage disequilibrium with the major effect gene Hd16. Examination of haplotype frequencies across different breeding periods revealed that the early-heading Hd16 (Hd16(E)) and late-heading Hd17 (Hd17(L)) were increasingly co-selected in the Hd1 functional population. Varieties carrying this Hd16(E)/Hd17(L) combination exhibited days to heading in the range of 70-80, which corresponds to the peak temperature and sunshine period and is also optimal for grain quality and yield components in the Fukui environment. The present study highlights that it is imperative to remain vigilant for Type I (false positives) and Type II (false negatives) errors when performing GWAS on highly bred populations and to implement appropriate countermeasures by accounting for gene-by-gene interactions established through the breeding process. We also discuss the effectiveness of Hd16(E), which is not used outside Japan for subtle days to heading control but is widely used in Japan at certain latitudes.

GWAS with principal component analysis identifies a gene comprehensively controlling rice architecture.

Elucidation of the genetic control of rice architecture is crucial due to the global demand for high crop yields. Rice architecture is a complex trait affected by plant height, tillering, and panicle morphology. In this study, principal component analysis (PCA) on 8 typical traits related to plant architecture revealed that the first principal component (PC), PC1, provided the most information on traits that determine rice architecture. A genome-wide association study (GWAS) using PC1 as a dependent variable was used to isolate a gene encoding rice, SPINDLY (OsSPY), that activates the gibberellin (GA) signal suppression protein SLR1. The effect of GA signaling on the regulation of rice architecture was confirmed in 9 types of isogenic plant having different levels of GA responsiveness. Further population genetics analysis demonstrated that the functional allele of OsSPY associated with semidwarfism and small panicles was selected in the process of rice breeding. In summary, the use of PCA in GWAS will aid in uncovering genes involved in traits with complex characteristics.

Genome-wide expression quantitative trait locus studies facilitate isolation of causal genes controlling panicle structure.

The morphology of rice (Oryza sativa L.) panicles is an important determinant of grain yield, and elucidation of the genetic control of panicle structure is very important for fulfilling the demand for high yield in breeding programs. In a quantitative trait locus (QTL) study using 82 backcross inbred lines (BILs) derived from Koshihikari and Habataki, 68 QTLs for 25 panicle morphological traits were identified. Gene expression profiling from inflorescence meristems of BILs was obtained. A combination of phenotypic QTL (pQTL) and expression QTL (eQTL) analysis revealed co-localization between pQTLs and eQTLs, consistent with significant correlations between phenotypic traits and gene expression levels. By combining pQTL and eQTL data, two genes were identified as controlling panicle structure: OsMADS18 modulates the average length of the primary rachis and OsFTL1 has pleiotropic effects on the total number of secondary rachides, number of grains per panicle, plant height and the length of flag leaves. Phenotypes were confirmed in RNA interference knocked-down plants and overexpressor lines. The combination of pQTL and eQTL analysis could facilitate identification of genes involved in rice panicle formation.

Making the 'Green Revolution' Truly Green: Improving Crop Nitrogen Use Efficiency.

Traditional breeding for high-yielding crops has mainly relied on the widespread cultivation of gibberellin (GA)-deficient semi-dwarf varieties, as dwarfism increases lodging resistance and allows for high nitrogen use, resulting in high grain yield. Although the adoption of semi-dwarf varieties in rice and wheat breeding brought big success to the 'Green Revolution' in the 20th century, it consequently increased the demand for nitrogen-based fertilizer, which causes severe threat to ecosystems and sustainable agriculture. To make the 'Green Revolution' truly green, it is necessary to develop new varieties with high nitrogen use efficiency (NUE). Under this demand, research on NUE, mainly for rice, has made great strides in the last decade. This mini-review focuses on three aspects of recent epoch-making findings on rice breeding for high NUE. The first one on 'NUE genes related to GA signaling' shows how promising it is to improve NUE in semi-dwarf Green Revolution varieties. The second aspect centers around the nitrate transporter1.1B, NRT1.1B; studies have revealed a nutrient signaling pathway through the discovery of the nitrate-NRT1.1B-SPX4-NLP3 cascade. The last one is based on the recent finding that the teosinte branched1, cycloidea, proliferating cell factor (TCP)-domain protein 19 underlies the genomic basis of geographical adaptation to soil nitrogen; OsTCP19 regulates the expression of a key transacting factor, DLT/SMOS2, which participates in the signaling of four different phytohormones, GA, auxin, brassinosteroid and strigolactone. Collectively, these breakthrough findings represent a significant step toward breeding high-NUE rice in the future.

Evolution and diversification of the plant gibberellin receptor GID1.

The plant gibberellin (GA) receptor GID1 shows sequence similarity to carboxylesterase (CXE). Here, we report the molecular evolution of GID1 from establishment to functionally diverse forms in eudicots. By introducing 18 mutagenized rice GID1s into a rice gid1 null mutant, we identified the amino acids crucial for GID1 activity in planta. We focused on two amino acids facing the C2/C3 positions of ent-gibberellane, not shared by lycophytes and euphyllophytes, and found that adjustment of these residues resulted in increased GID1 affinity toward GA4, new acceptance of GA1 and GA3 carrying C13-OH as bioactive ligands, and elimination of inactive GAs. These residues rendered the GA perception system more sophisticated. We conducted phylogenetic analysis of 169 GID1s from 66 plant species and found that, unlike other taxa, nearly all eudicots contain two types of GID1, named A- and B-type. Certain B-type GID1s showed a unique evolutionary characteristic of significantly higher nonsynonymous-to-synonymous divergence in the region determining GA4 affinity. Furthermore, these B-type GID1s were preferentially expressed in the roots of Arabidopsis, soybean, and lettuce and might be involved in root elongation without shoot elongation for adaptive growth under low-temperature stress. Based on these observations, we discuss the establishment and adaption of GID1s during plant evolution.

Evolution of GA Metabolic Enzymes in Land Plants.

Gibberellins (GAs) play key roles in various developmental processes in land plants. We studied the evolutionary trends of GA metabolic enzymes through a comprehensive homology search and phylogenetic analyses from bryophytes to angiosperms. Our analyses suggest that, in the process of evolution, plants were able to acquire GA metabolic enzymes in a stepwise manner and that the enzymes had rapidly diversified in angiosperms. As a good example of their rapid diversification, we focused on the GA-deactivating enzyme, GA 2-oxidase (GA2ox). Although the establishment of a GA system first occurred in lycophytes, its inactivation system mediated by GA2oxs was established at a much later time: the rise of gymnosperms and the rise of angiosperms through C19-GA2ox and C20-GA2ox development, respectively, as supported by the results of our direct examination of their enzymatic activities in vitro. Based on these comprehensive studies of GA metabolic enzymes, we discuss here that angiosperms rapidly developed a sophisticated system to delicately control the level of active GAs by increasing their copy numbers for their survival under different challenging environments.

Arabidopsis Group IIId ERF proteins positively regulate primary cell wall-type CESA genes.

The cell wall determines morphology and the environmental responses of plant cells. The primary cell wall (PCW) is produced during cell division and expansion, determining the cell shape and volume. After cell expansion, specific types of plant cells produce a lignified wall, known as a secondary cell wall (SCW). We functionally analyzed Group IIId Arabidopsis AP2/EREBP genes, namely ERF34, ERF35, ERF38, and ERF39, which are homologs of a rice ERF gene previously proposed to be related to SCW biosynthesis. Expression analysis revealed that these four genes are expressed in regions related to cell division and/or cell differentiation in seedlings (i.e., shoot apical meristems, the primordia of leaves and lateral roots, trichomes, and central cylinder of primary roots) and flowers (i.e., vascular tissues of floral organs and replums and/or valve margins of pistils). Overexpression of ERF genes significantly upregulated PCW-type, but not SCW-type, CESA genes encoding cellulose synthase catalytic subunits in Arabidopsis seedlings. Transient co-expression reporter analysis indicated that ERF35, ERF38, and ERF39 possess transcriptional activator activity, and that ERF34, ERF35, ERF38, and ERF39 upregulated the promoter activity of CESA1, a PCW-type CESA gene, through the DRECRTCOREAT elements, the core cis-acting elements known to be recognized by AP2/ERF proteins. Together, our findings show that Group IIId ERF genes are positive transcriptional regulators of PCW-type CESA genes in Arabidopsis and are possibly involved in modulating cellulose biosynthesis in response to developmental requirements and environmental stimuli.

Regulation of Strigolactone Biosynthesis by Gibberellin Signaling.

Strigolactones (SLs) are a class of plant hormones that regulate diverse physiological processes, including shoot branching and root development. They also act as rhizosphere signaling molecules to stimulate the germination of root parasitic weeds and the branching of arbuscular mycorrhizal fungi. Although various types of cross talk between SLs and other hormones have been reported in physiological analyses, the cross talk between gibberellin (GA) and SLs is poorly understood. We screened for chemicals that regulate the level of SLs in rice (Oryza sativa) and identified GA as, to our knowledge, a novel SL-regulating molecule. The regulation of SL biosynthesis by GA is dependent on the GA receptor GID1 and F-box protein GID2. GA treatment also reduced the infection of rice plants by the parasitic plant witchers weed (Striga hermonthica). These data not only demonstrate, to our knowledge, the novel plant hormone cross talk between SL and GA, but also suggest that GA can be used to control parasitic weed infections.

OsIDD2, a zinc finger and INDETERMINATE DOMAIN protein, regulates secondary cell wall formation.

Previously, we found 123 transcription factors (TFs) as candidate regulators of secondary cell wall (SCW) formation in rice by using phylogenetic and co-expression network analyses. Among them, we examined in this work the role of OsIDD2, a zinc finger and indeterminate domain (IDD) family TF. Its overexpressors showed dwarfism, fragile leaves, and decreased lignin content, which are typical phenotypes of plants defective in SCW formation, whereas its knockout plants showed slightly increased lignin content. The RNA-seq and quantitative reverse transcription polymerase chain reaction analyses confirmed that some lignin biosynthetic genes were downregulated in the OsIDD2-overexpressing plants, and revealed the same case for other genes involved in cellulose synthesis and sucrose metabolism. The transient expression assay using rice protoplasts revealed that OsIDD2 negatively regulates the transcription of genes involved in lignin biosynthesis, cinnamyl alcohol dehydrogenase 2 and 3 (CAD2 and 3), and sucrose metabolism, sucrose synthase 5 (SUS5), whereas an AlphaScreen assay, which can detect the interaction between TFs and their target DNA sequences, directly confirmed the interaction between OsIDD2 and the target sequences located in the promoter regions of CAD2 and CAD3. Based on these observations, we conclude that OsIDD2 is negatively involved in SCW formation and other biological events by downregulating its target genes.

DELLA protein functions as a transcriptional activator through the DNA binding of the indeterminate domain family proteins.

DELLA protein is a key negative regulator of gibberellin (GA) signaling. Although how DELLA regulates downstream gene expression remains unclear, DELLA has been proposed to function as a transcriptional activator. However, because DELLA lacks a DNA-binding domain, intermediate protein(s) mediating the DELLA/DNA interaction are believed to be necessary for activating DELLA target genes. Here, using yeast hybrid screenings, we identified five members of indeterminate domain (IDD) protein family which bind physically to both DELLA and the promoter sequence of the GA-positive regulator SCARECROW-LIKE 3 (SCL3), which previously was characterized as a DELLA direct target gene. Transient assays using Arabidopsis protoplasts demonstrated that a luciferase reporter controlled by the SCL3 promoter was additively transactivated by REPRESSOR of ga1-3 (RGA) and IDDs. Phenotypic analysis of transgenic plants expressing AtIDD3 (one of the 16 IDDs in the Arabidopsis genome) fused with the plant-specific repression domain (SRDX) supported the possibility that AtIDD3 is positively involved in GA signaling. In addition, we found that SCL3 protein also interacts with IDDs, resulting in the suppression of its target gene expression. In this context, DELLA and SCL3 interact competitively with IDD proteins to regulate downstream gene expression. These results suggest that the coregulators DELLA and SCL3, using IDDs as transcriptional scaffolds for DNA binding, antagonistically regulate the expression of their downstream targets to control the GA signaling pathway.

Engineering the lodging resistance mechanism of post-Green Revolution rice to meet future demands.

Traditional breeding for high-yielding rice has been dependent on the widespread cultivation of gibberellin (GA)-deficient semi-dwarf varieties. Dwarfism lowers the "center of gravity" of the plant body, which increases resistance against lodging and enables plants to support high grain yield. Although this approach was successful in latter half of the 20th century in rice and wheat breeding, this may no longer be enough to sustain rice with even higher yields. This is because relying solely on the semi-dwarf trait is subject to certain limitations, making it necessary to use other important traits to reinforce it. In this review, we present an alternative approach to increase lodging resistance by improving the quality of the culm by identifying genes related to culm quality and introducing these genes into high-yielding rice cultivars through molecular breeding technique.

Bar the windows: an optimized strategy to survive drought and salt adversities.

Hydrogen peroxide (H(2)O(2)) is a central modulator of stomatal closure. It remains unknown, however, how the upstream regulation of H(2)O(2) homeostasis operates. In this issue of Genes & Development, Huang and colleagues (pp. 1805-1817) report that a novel C(2)H(2)-type transcription factor, drought and salt tolerance (DST), mediates H(2)O(2)-induced stomatal closure and abiotic stress tolerance.

Comprehensive gene expression analysis of rice aleurone cells: probing the existence of an alternative gibberellin receptor.

Current gibberellin (GA) research indicates that GA must be perceived in plant nuclei by its cognate receptor, GIBBERELLIN INSENSITIVE DWARF1 (GID1). Recognition of GA by GID1 relieves the repression mediated by the DELLA protein, a model known as the GID1-DELLA GA perception system. There have been reports of potential GA-binding proteins in the plasma membrane that perceive GA and induce α-amylase expression in cereal aleurone cells, which is mechanistically different from the GID1-DELLA system. Therefore, we examined the expression of the rice (Oryza sativa) α-amylase genes in rice mutants impaired in the GA receptor (gid1) and the DELLA repressor (slender rice1; slr1) and confirmed their lack of response to GA in gid1 mutants and constitutive expression in slr1 mutants. We also examined the expression of GA-regulated genes by genome-wide microarray and quantitative reverse transcription-polymerase chain reaction analyses and confirmed that all GA-regulated genes are modulated by the GID1-DELLA system. Furthermore, we studied the regulatory network involved in GA signaling by using a set of mutants defective in genes involved in GA perception and gene expression, namely gid1, slr1, gid2 (a GA-related F-box protein mutant), and gamyb (a GA-related trans-acting factor mutant). Almost all GA up-regulated genes were regulated by the four named GA-signaling components. On the other hand, GA down-regulated genes showed different expression patterns with respect to GID2 and GAMYB (e.g. a considerable number of genes are not controlled by GAMYB or GID2 and GAMYB). Based on these observations, we present a comprehensive discussion of the intricate network of GA-regulated genes in rice aleurone cells.

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.