Exome-wide variation in a diverse barley panel reveals genetic associations with ten agronomic traits in Eastern landraces.

Barley (Hordeum vulgare ssp. vulgare) is one of the first crops to be domesticated and is adapted to a wide range of environments. Worldwide barley germplasm collections possess valuable allelic variations that could further improve barley productivity. Although barley genomics has offered a global picture of allelic variation among varieties and its association with various agronomic traits, polymorphisms from East Asian varieties remain scarce. In this study, we analyze exome polymorphisms in a panel of 274 barley varieties collected worldwide, including 137 varieties from East Asian countries and Ethiopia. We reveal the underlying population structure and conduct genome-wide association studies for 10 agronomic traits. Moreover, we examin genome-wide associations for traits related to grain size such as awn length and glume length. Our results demonstrate the value of diverse barley germplasm panels containing Eastern varieties, highlighting their distinct genomic signatures relative to Western subpopulations.

Cellular and transcriptomic analyses reveal two-staged chloroplast biogenesis underpinning photosynthesis build-up in the wheat leaf.

BACKGROUND: The developmental gradient in monocot leaves has been exploited to uncover leaf developmental gene expression programs and chloroplast biogenesis processes. However, the relationship between the two is barely understood, which limits the value of transcriptome data to understand the process of chloroplast development. RESULTS: Taking advantage of the developmental gradient in the bread wheat leaf, we provide a simultaneous quantitative analysis for the development of mesophyll cells and of chloroplasts as a cellular compartment. This allows us to generate the first biologically-informed gene expression map of this leaf, with the entire developmental gradient from meristematic to fully differentiated cells captured. We show that the first phase of plastid development begins with organelle proliferation, which extends well beyond cell proliferation, and continues with the establishment and then the build-up of the plastid genetic machinery. The second phase is marked by the development of photosynthetic chloroplasts which occupy the available cellular space. Using a network reconstruction algorithm, we predict that known chloroplast gene expression regulators are differentially involved across those developmental stages. CONCLUSIONS: Our analysis generates both the first wheat leaf transcriptional map and one of the most comprehensive descriptions to date of the developmental history of chloroplasts in higher plants. It reveals functionally distinct plastid and chloroplast development stages, identifies processes occurring in each of them, and highlights our very limited knowledge of the earliest drivers of plastid biogenesis, while providing a basis for their future identification.

Genetic Elucidation for Response of Flowering Time to Ambient Temperatures in Asian Rice Cultivars.

Climate resilience of crops is critical for global food security. Understanding the genetic basis of plant responses to ambient environmental changes is key to developing resilient crops. To detect genetic factors that set flowering time according to seasonal temperature conditions, we evaluated differences of flowering time over years by using chromosome segment substitution lines (CSSLs) derived from japonica rice cultivars "Koshihikari" × "Khao Nam Jen", each with different robustness of flowering time to environmental fluctuations. The difference of flowering times in 9 years' field tests was large in "Khao Nam Jen" (36.7 days) but small in "Koshihikari" (9.9 days). Part of this difference was explained by two QTLs. A CSSL with a "Khao Nam Jen" segment on chromosome 11 showed 28.0 days' difference; this QTL would encode a novel flowering-time gene. Another CSSL with a segment from "Khao Nam Jen" in the region around Hd16 on chromosome 3 showed 23.4 days" difference. A near-isogenic line (NIL) for Hd16 showed 21.6 days' difference, suggesting Hd16 as a candidate for this QTL. RNA-seq analysis showed differential expression of several flowering-time genes between early and late flowering seasons. Low-temperature treatment at panicle initiation stage significantly delayed flowering in the CSSL and NIL compared with "Koshihikari". Our results unravel the molecular control of flowering time under ambient temperature fluctuations.

BdWRKY38 is required for the incompatible interaction of Brachypodium distachyon with the necrotrophic fungus Rhizoctonia solani.

Rhizoctonia solani is a soil-borne necrotrophic fungus that causes sheath blight in grasses. The basal resistance of compatible interactions between R. solani and rice is known to be modulated by some WRKY transcription factors (TFs). However, genes and defense responses involved in incompatible interaction with R. solani remain unexplored, because no such interactions are known in any host plants. Recently, we demonstrated that Bd3-1, an accession of the model grass Brachypodium distachyon, is resistant to R. solani and, upon inoculation with the fungus, undergoes rapid induction of genes responsive to the phytohormone salicylic acid (SA) that encode the WRKY TFs BdWRKY38 and BdWRKY44. Here, we show that endogenous SA and these WRKY TFs positively regulate this accession-specific R. solani resistance. In contrast to a susceptible accession (Bd21), the infection process in the resistant accessions Bd3-1 and Tek-3 was suppressed at early stages before the development of fungal biomass and infection machinery. A comparative transcriptome analysis during pathogen infection revealed that putative WRKY-dependent defense genes were induced faster in the resistant accessions than in Bd21. A gene regulatory network (GRN) analysis based on the transcriptome dataset demonstrated that BdWRKY38 was a GRN hub connected to many target genes specifically in resistant accessions, whereas BdWRKY44 was shared in the GRNs of all three accessions. Moreover, overexpression of BdWRKY38 increased R. solani resistance in Bd21. Our findings demonstrate that these resistant accessions can activate an incompatible host response to R. solani, and BdWRKY38 regulates this response by mediating SA signaling.

Genetic Factors Associated with Heading Responses Revealed by Field Evaluation of 274 Barley Accessions for 20 Seasons.

Heading time is a key trait in cereals affecting the maturation period for optimal grain filling before harvest. Here, we aimed to understand the factors controlling heading time in barley (Hordeum vulgare). We characterized a set of 274 barley accessions collected worldwide by planting them for 20 seasons under different environmental conditions at the same location in Kurashiki, Japan. We examined interactions among accessions, known genetic factors, and an environmental factor to determine the factors controlling heading response. Locally adapted accessions have been selected for genetic factors that stabilize heading responses appropriate for barley cultivation, and these accessions show stable heading responses even under varying environmental conditions. We identified vernalization requirement and PPD-H1 haplotype as major stabilizing mechanisms of the heading response for regional adaptation in Kurashiki.

Parental legacy and regulatory novelty in Brachypodium diurnal transcriptomes accompanying their polyploidy.

Polyploidy is a widespread phenomenon in eukaryotes that can lead to phenotypic novelty and has important implications for evolution and diversification. The modification of phenotypes in polyploids relative to their diploid progenitors may be associated with altered gene expression. However, it is largely unknown how interactions between duplicated genes affect their diurnal expression in allopolyploid species. In this study, we explored parental legacy and hybrid novelty in the transcriptomes of an allopolyploid species and its diploid progenitors. We compared the diurnal transcriptomes of representative Brachypodium cytotypes, including the allotetraploid Brachypodium hybridum and its diploid progenitors Brachypodium distachyon and Brachypodium stacei. We also artificially induced an autotetraploid B. distachyon. We identified patterns of homoeolog expression bias (HEB) across Brachypodium cytotypes and time-dependent gain and loss of HEB in B. hybridum. Furthermore, we established that many genes with diurnal expression experienced HEB, while their expression patterns and peak times were correlated between homoeologs in B. hybridum relative to B. distachyon and B. stacei, suggesting diurnal synchronization of homoeolog expression in B. hybridum. Our findings provide insight into the parental legacy and hybrid novelty associated with polyploidy in Brachypodium, and highlight the evolutionary consequences of diurnal transcriptional regulation that accompanied allopolyploidy.

Genetic Variation for Seed Metabolite Levels in Brachypodium distachyon.

Metabolite composition and concentrations in seed grains are important traits of cereals. To identify the variation in the seed metabolotypes of a model grass, namely Brachypodium distachyon, we applied a widely targeted metabolome analysis to forty inbred lines of B. distachyon and examined the accumulation patterns of 183 compounds in the seeds. By comparing the metabolotypes with the population structure of these lines, we found signature metabolites that represent different accumulation patterns for each of the three B. distachyon subpopulations. Moreover, we found that thirty-seven metabolites exhibited significant differences in their accumulation between the lines Bd21 and Bd3-1. Using a recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21, we identified the quantitative trait loci (QTLs) linked with this variation in the accumulation of thirteen metabolites. Our metabolite QTL analysis illustrated that different genetic factors may presumably regulate the accumulation of 4-pyridoxate and pyridoxamine in vitamin B6 metabolism. Moreover, we found two QTLs on chromosomes 1 and 4 that affect the accumulation of an anthocyanin, chrysanthemin. These QTLs genetically interacted to regulate the accumulation of this compound. This study demonstrates the potential for metabolite QTL mapping in B. distachyon and provides new insights into the genetic dissection of metabolomic traits in temperate grasses.

Transcriptional alterations during proliferation and lignification in Phyllostachys nigra cells.

Highly-lignified culms of bamboo show distinctive anatomical and mechanical properties compared with the culms of other grass species. A cell culture system for Phyllostachys nigra has enabled investigating the alterations in cellular states associated with secondary cell wall formation during its proliferation and lignification in woody bamboos. To reveal transcriptional changes related to lignification in bamboo, we analyzed transcriptome in P. nigra cells treated with the synthetic auxin 2,4-dichlorophenoxyacetic acid (2,4-D) and the synthetic cytokinin benzylaminopurine (BA) by RNA-seq analysis. We found that some genes putatively involved in cell wall biogenesis and cell division were up-regulated in response to the 2,4-D treatment, and the induction of lignification by the BA treatment was correlated with up-regulation of genes involved in the shikimate pathway. We also found that genes encoding MYB transcription factors (TFs) show correlated expression patterns with those encoding cinnamyl alcohol dehydrogenase (CAD), suggesting that MYB TFs presumably regulate secondary cell wall formation in the bamboo cells. These findings suggest that cytokinin signaling may regulate lignification in P. nigra cells through coordinated transcriptional regulation and metabolic alterations. Our results have also produced a useful resource for better understanding of secondary cell wall formation in bamboo plants.

Homoeolog-specific activation of genes for heat acclimation in the allopolyploid grass Brachypodium hybridum.

Background: Allopolyploid plants often show wider environmental tolerances than their ancestors; this is expected to be due to the merger of multiple distinct genomes with a fixed heterozygosity. The complex homoeologous gene expression could have been evolutionarily advantageous for the adaptation of allopolyploid plants. Despite multiple previous studies reporting homoeolog-specific gene expression in allopolyploid species, there are no clear examples of homoeolog-specific function in acclimation to a long-term stress condition. Results: We found that the allopolyploid grass Brachypodium hybridum and its ancestor Brachypodium stacei show long-term heat stress tolerance, unlike its other ancestor, Brachypodium distachyon. To understand the physiological traits of B. hybridum, we compared the transcriptome of the 3 Brachypodium species grown under normal and heat stress conditions. We found that the expression patterns of approximately 26% and approximately 38% of the homoeolog groups in B. hybridum changed toward nonadditive expression and nonancestral expression, respectively, under normal condition. Moreover, we found that B. distachyon showed similar expression patterns between normal and heat stress conditions, whereas B. hybridum and B. stacei significantly altered their transcriptome in response to heat after 3 days of stress exposure, and homoeologs that were inherited from B. stacei may have contributed to the transcriptional stress response to heat in B. hybridum. After 15 days of heat exposure, B. hybridum and B. stacei maintained transcriptional states similar to those under normal conditions. These results suggest that an earlier response to heat that was specific to homoeologs originating from B. stacei contributed to cellular homeostasis under long-term heat stress in B. hybridum. Conclusions: Our results provide insights into different regulatory events of the homoeo-transcriptome that are associated with stress acclimation in allopolyploid plants.

Analysis of single nucleotide polymorphisms based on RNA sequencing data of diverse bio-geographical accessions in barley.

Barley is one of the founder crops of Old world agriculture and has become the fourth most important cereal worldwide. Information on genome-scale DNA polymorphisms allows elucidating the evolutionary history behind domestication, as well as discovering and isolating useful genes for molecular breeding. Deep transcriptome sequencing enables the exploration of sequence variations in transcribed sequences; such analysis is particularly useful for species with large and complex genomes, such as barley. In this study, we performed RNA sequencing of 20 barley accessions, comprising representatives of several biogeographic regions and a wild ancestor. We identified 38,729 to 79,949 SNPs in the 19 domesticated accessions and 55,403 SNPs in the wild barley and revealed their genome-wide distribution using a reference genome. Genome-scale comparisons among accessions showed a clear differentiation between oriental and occidental barley populations. The results based on population structure analyses provide genome-scale properties of sub-populations grouped to oriental, occidental and marginal groups in barley. Our findings suggest that the oriental population of domesticated barley has genomic variations distinct from those in occidental groups, which might have contributed to barley's domestication.

Large-scale collection and analysis of full-length cDNAs from Brachypodium distachyon and integration with Pooideae sequence resources.

A comprehensive collection of full-length cDNAs is essential for correct structural gene annotation and functional analyses of genes. We constructed a mixed full-length cDNA library from 21 different tissues of Brachypodium distachyon Bd21, and obtained 78,163 high quality expressed sequence tags (ESTs) from both ends of ca. 40,000 clones (including 16,079 contigs). We updated gene structure annotations of Brachypodium genes based on full-length cDNA sequences in comparison with the latest publicly available annotations. About 10,000 non-redundant gene models were supported by full-length cDNAs; ca. 6,000 showed some transcription unit modifications. We also found ca. 580 novel gene models, including 362 newly identified in Bd21. Using the updated transcription start sites, we searched a total of 580 plant cis-motifs in the -3 kb promoter regions and determined a genome-wide Brachypodium promoter architecture. Furthermore, we integrated the Brachypodium full-length cDNAs and updated gene structures with available sequence resources in wheat and barley in a web-accessible database, the RIKEN Brachypodium FL cDNA database. The database represents a "one-stop" information resource for all genomic information in the Pooideae, facilitating functional analysis of genes in this model grass plant and seamless knowledge transfer to the Triticeae crops.

Global landscape of a co-expressed gene network in barley and its application to gene discovery in Triticeae crops.

Accumulated transcriptome data can be used to investigate regulatory networks of genes involved in various biological systems. Co-expression analysis data sets generated from comprehensively collected transcriptome data sets now represent efficient resources that are capable of facilitating the discovery of genes with closely correlated expression patterns. In order to construct a co-expression network for barley, we analyzed 45 publicly available experimental series, which are composed of 1,347 sets of GeneChip data for barley. On the basis of a gene-to-gene weighted correlation coefficient, we constructed a global barley co-expression network and classified it into clusters of subnetwork modules. The resulting clusters are candidates for functional regulatory modules in the barley transcriptome. To annotate each of the modules, we performed comparative annotation using genes in Arabidopsis and Brachypodium distachyon. On the basis of a comparative analysis between barley and two model species, we investigated functional properties from the representative distributions of the gene ontology (GO) terms. Modules putatively involved in drought stress response and cellulose biogenesis have been identified. These modules are discussed to demonstrate the effectiveness of the co-expression analysis. Furthermore, we applied the data set of co-expressed genes coupled with comparative analysis in attempts to discover potentially Triticeae-specific network modules. These results demonstrate that analysis of the co-expression network of the barley transcriptome together with comparative analysis should promote the process of gene discovery in barley. Furthermore, the insights obtained should be transferable to investigations of Triticeae plants. The associated data set generated in this analysis is publicly accessible at http://coexpression.psc.riken.jp/barley/.