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.

Time-series transcriptome of Brachypodium distachyon during bacterial flagellin-induced pattern-triggered immunity.

Plants protect themselves from microorganisms by inducing pattern-triggered immunity (PTI) via recognizing microbe-associated molecular patterns (MAMPs), conserved across many microbes. Although the MAMP perception mechanism and initial events during PTI have been well-characterized, knowledge of the transcriptomic changes in plants, especially monocots, is limited during the intermediate and terminal stages of PTI. Here, we report a time-series high-resolution RNA-sequencing (RNA-seq) analysis during PTI in the leaf disks of Brachypodium distachyon. We identified 6,039 differentially expressed genes (DEGs) in leaves sampled at 0, 0.5, 1, 3, 6, and 12 hours after treatment (hat) with the bacterial flagellin peptide flg22. The k-means clustering method classified these DEGs into 10 clusters (6 upregulated and 4 downregulated). Based on the results, we selected 10 PTI marker genes in B. distachyon. Gene ontology (GO) analysis suggested a tradeoff between defense responses and photosynthesis during PTI. The data indicated the recovery of photosynthesis started at least at 12 hat. Over-representation analysis of transcription factor genes and cis-regulatory elements in DEG promoters implied the contribution of 12 WRKY transcription factors in plant defense at the early stage of PTI induction.

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.

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.

Life-Course Monitoring of Endogenous Phytohormone Levels under Field Conditions Reveals Diversity of Physiological States among Barley Accessions.

Agronomically important traits often develop during the later stages of crop growth as consequences of various plant-environment interactions. Therefore, the temporal physiological states that change and accumulate during the crop's life course can significantly affect the eventual phenotypic differences in agronomic traits among crop varieties. Thus, to improve productivity, it is important to elucidate the associations between temporal physiological responses during the growth of different crop varieties and their agronomic traits. However, data representing the dynamics and diversity of physiological states in plants grown under field conditions are sparse. In this study, we quantified the endogenous levels of five phytohormones - auxin, cytokinins (CKs), ABA, jasmonate and salicylic acid - in the leaves of eight diverse barley (Hordeum vulgare) accessions grown under field conditions sampled weekly over their life course to assess the ongoing fluctuations in hormone levels in the different accessions under field growth conditions. Notably, we observed enormous changes over time in the development-related plant hormones, such as auxin and CKs. Using 3' RNA-seq-based transcriptome data from the same samples, we investigated the expression of barley genes orthologous to known hormone-related genes of Arabidopsis throughout the life course. These data illustrated the dynamics and diversity of the physiological states of these field-grown barley accessions. Together, our findings provide new insights into plant-environment interactions, highlighting that there is cultivar diversity in physiological responses during growth under field conditions.

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.

Gene Co-expression Network Analysis Suggests the Existence of Transcriptional Modules Containing a High Proportion of Transcriptionally Differentiated Homoeologs in Hexaploid Wheat.

Genome duplications aid in the formation of novel molecular networks through regulatory differentiation of the duplicated genes and facilitate adaptation to environmental change. Hexaploid wheat, Triticum aestivum, contains three homoeologous chromosome sets, the A-, B-, and D-subgenomes, which evolved through interspecific hybridization and subsequent whole-genome duplication. The divergent expression patterns of the homoeologs in hexaploid wheat suggest that they have undergone transcriptional and/or functional differentiation during wheat evolution. However, the distribution of transcriptionally differentiated homoeologs in gene regulatory networks and their related biological functions in hexaploid wheat are still largely unexplored. Therefore, we retrieved 727 publicly available wheat RNA-sequencing (RNA-seq) datasets from various tissues, developmental stages, and conditions, and identified 10,415 expressed homoeologous triplets. Examining the co-expression modules in the wheat transcriptome, we found that 66% of the expressed homoeologous triplets possess all three homoeologs grouped in the same co-expression modules. Among these, 15 triplets contain co-expressed homoeologs with differential expression levels between homoeoalleles across ≥ 95% of the 727 RNA-seq datasets, suggesting a consistent trend of homoeolog expression bias. In addition, we identified 2,831 differentiated homoeologs that showed gene expression patterns that deviated from those of the other two homoeologs. We found that seven co-expression modules contained a high proportion of such differentiated homoeologs, which accounted for ≥ 20% of the genes in each module. We also found that five of the co-expression modules are abundantly composed of genes involved in biological processes such as chloroplast biogenesis, RNA metabolism, putative defense response, putative posttranscriptional modification, and lipid metabolism, thereby suggesting that, the differentiated homoeologs might highly contribute to these biological functions in the gene network of hexaploid wheat.

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.

Multiplex PCR Targeted Amplicon Sequencing (MTA-Seq): Simple, Flexible, and Versatile SNP Genotyping by Highly Multiplexed PCR Amplicon Sequencing.

Next-generation sequencing (NGS) technologies have enabled genome re-sequencing for exploring genome-wide polymorphisms among individuals, as well as targeted re-sequencing for the rapid and simultaneous detection of polymorphisms in genes associated with various biological functions. Therefore, a simple and robust method for targeted re-sequencing should facilitate genotyping in a wide range of biological fields. In this study, we developed a simple, custom, targeted re-sequencing method, designated "multiplex PCR targeted amplicon sequencing (MTA-seq)," and applied it to the genotyping of the model grass Brachypodium distachyon. To assess the practical usability of MTA-seq, we applied it to the genotyping of genome-wide single-nucleotide polymorphisms (SNPs) identified in natural accessions (Bd1-1, Bd3-1, Bd21-3, Bd30-1, Koz-1, Koz-3, and Koz-4) by comparing the re-sequencing data with that of reference accession Bd21. Examination of SNP-genotyping accuracy in 443 amplicons from eight parental accessions and an F1 progeny derived by crossing of Bd21 and Bd3-1 revealed that ~95% of the SNPs were correctly called. The assessment suggested that the method provided an efficient framework for accurate and robust SNP genotyping. The method described here enables easy design of custom target SNP-marker panels in various organisms, facilitating a wide range of high-throughput genetic applications, such as genetic mapping, population analysis, molecular breeding, and genomic diagnostics.

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.

Draft genome assembly and annotation of Glycyrrhiza uralensis, a medicinal legume.

Chinese liquorice/licorice (Glycyrrhiza uralensis) is a leguminous plant species whose roots and rhizomes have been widely used as a herbal medicine and natural sweetener. Whole-genome sequencing is essential for gene discovery studies and molecular breeding in liquorice. Here, we report a draft assembly of the approximately 379-Mb whole-genome sequence of strain 308-19 of G. uralensis; this assembly contains 34 445 predicted protein-coding genes. Comparative analyses suggested well-conserved genomic components and collinearity of gene loci (synteny) between the genome of liquorice and those of other legumes such as Medicago and chickpea. We observed that three genes involved in isoflavonoid biosynthesis, namely, 2-hydroxyisoflavanone synthase (CYP93C), 2,7,4'-trihydroxyisoflavanone 4'-O-methyltransferase/isoflavone 4'-O-methyltransferase (HI4OMT) and isoflavone-7-O-methyltransferase (7-IOMT) formed a cluster on the scaffold of the liquorice genome and showed conserved microsynteny with Medicago and chickpea. Based on the liquorice genome annotation, we predicted genes in the P450 and UDP-dependent glycosyltransferase (UGT) superfamilies, some of which are involved in triterpenoid saponin biosynthesis, and characterised their gene expression with the reference genome sequence. The genome sequencing and its annotations provide an essential resource for liquorice improvement through molecular breeding and the discovery of useful genes for engineering bioactive components through synthetic biology approaches.

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.