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.

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

•

Heading of 274 barley worldwide accessions were evaluated for 20 seasons

•

Locally adapted accessions show stable heading responses

•

Vernalization requirement and PPD-H1 haplotype stabilize the heading response

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.

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.

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.

Diurnal Transcriptome and Gene Network Represented through Sparse Modeling in Brachypodium distachyon

We report the comprehensive identification of periodic genes and their network inference, based on a gene co-expression analysis and an Auto-Regressive eXogenous (ARX) model with a group smoothly clipped absolute deviation (SCAD) method using a time-series transcriptome dataset in a model grass, Brachypodium distachyon. To reveal the diurnal changes in the transcriptome in B. distachyon, we performed RNA-seq analysis of its leaves sampled through a diurnal cycle of over 48 h at 4 h intervals using three biological replications, and identified 3,621 periodic genes through our wavelet analysis. The expression data are feasible to infer network sparsity based on ARX models. We found that genes involved in biological processes such as transcriptional regulation, protein degradation, and post-transcriptional modification and photosynthesis are significantly enriched in the periodic genes, suggesting that these processes might be regulated by circadian rhythm in B. distachyon. On the basis of the time-series expression patterns of the periodic genes, we constructed a chronological gene co-expression network and identified putative transcription factors encoding genes that might be involved in the time-specific regulatory transcriptional network. Moreover, we inferred a transcriptional network composed of the periodic genes in B. distachyon, aiming to identify genes associated with other genes through variable selection by grouping time points for each gene. Based on the ARX model with the group SCAD regularization using our time-series expression datasets of the periodic genes, we constructed gene networks and found that the networks represent typical scale-free structure. Our findings demonstrate that the diurnal changes in the transcriptome in B. distachyon leaves have a sparse network structure, demonstrating the spatiotemporal gene regulatory network over the cyclic phase transitions in B. distachyon diurnal growth.

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