Transcriptome analysis of the Japanese eel (Anguilla japonica) during larval metamorphosis.

BACKGROUND: Anguillid eels spend their larval period as leptocephalus larvae that have a unique and specialized body form with leaf-like and transparent features, and they undergo drastic metamorphosis to juvenile glass eels. Less is known about the transition of leptocephali to the glass eel stage, because it is difficult to catch the metamorphosing larvae in the open ocean. However, recent advances in rearing techniques for the Japanese eel have made it possible to study the larval metamorphosis of anguillid eels. In the present study, we investigated the dynamics of gene expression during the metamorphosis of Japanese eel leptocephali using RNA sequencing. RESULTS: During metamorphosis, Japanese eels were classified into 7 developmental stages according to their morphological characteristics, and RNA sequencing was used to collect gene expression data from each stage. A total of 354.8 million clean reads were generated from the body and 365.5 million from the head, after the processing of raw reads. For filtering of genes that characterize developmental stages, a classification model created by a Random Forest algorithm was built. Using the importance of explanatory variables feature obtained from the created model, we identified 46 genes selected in the body and 169 genes selected in the head that were defined as the "most characteristic genes" during eel metamorphosis. Next, network analysis and subsequently gene clustering were conducted using the most characteristic genes and their correlated genes, and then 6 clusters in the body and 5 clusters in the head were constructed. Then, the characteristics of the clusters were revealed by Gene Ontology (GO) enrichment analysis. The expression patterns and GO terms of each stage were consistent with previous observations and experiments during the larval metamorphosis of the Japanese eel. CONCLUSION: Genome and transcriptome resources have been generated for metamorphosing Japanese eels. Genes that characterized metamorphosis of the Japanese eel were identified through statistical modeling by a Random Forest algorithm. The functions of these genes were consistent with previous observations and experiments during the metamorphosis of anguillid eels.

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.

Phosphate or nitrate imbalance induces stronger molecular responses than combined nutrient deprivation in roots and leaves of chickpea plants.

The negative effects of phosphate (Pi) and/or nitrate (NO3- ) fertilizers on the environment have raised an urgent need to develop crop varieties with higher Pi and/or nitrogen use efficiencies for cultivation in low-fertility soils. Achieving this goal depends upon research that focuses on the identification of genes involved in plant responses to Pi and/or NO3- starvation. Although plant responses to individual deficiency in either Pi (-Pi/+NO3- ) or NO3- (+Pi/-NO3- ) have been separately studied, our understanding of plant responses to combined Pi and NO3- deficiency (-Pi/-NO3- ) is still very limited. Using RNA-sequencing approach, transcriptome changes in the roots and leaves of chickpea cultivated under -Pi/+NO3- , +Pi/-NO3- or -Pi/-NO3- conditions were investigated in a comparative manner. -Pi/-NO3- treatment displayed lesser effect on expression changes of genes related to Pi or NO3- transport, signalling networks, lipid remodelling, nitrogen and Pi scavenging/remobilization/recycling, carbon metabolism and hormone metabolism than -Pi/+NO3- or +Pi/-NO3- treatments. Therefore, the plant response to -Pi/-NO3- is not simply an additive result of plant responses to -Pi/+NO3- and +Pi/-NO3- treatments. Our results indicate that nutrient imbalance is a stronger stimulus for molecular reprogramming than an overall deficiency.

Plastic transport systems of rice for mineral elements in response to diverse soil environmental changes.

Climate change will increase frequency of drought and flooding, which threaten global crop productivity and food security. Rice (Oryza sativa) is unique in that it is able to grow in both flooded and upland conditions, which have large differences in the concentrations and chemical forms of mineral elements available to plants. To comprehensively understand the mechanisms of rice for coping with different water status, we performed ionomics and transcriptomics analysis of the roots, nodes and leaves of rice grown in flooded and upland conditions. Focusing the analysis on genes encoding proteins involved in transport functions for mineral elements, it was found that, although rice plants maintained similar levels of each element in the shoots for optimal growth, different transporters for mineral elements were utilised for nitrogen, iron, copper and zinc to deal with different soil water conditions. For example, under flooded conditions, rice roots take up nitrogen using transporters for both ammonium (OsAMT1/2) and nitrate (OsNPF2.4, OsNRT1.1A and OsNRT2.3), whereas under upland conditions, nitrogen uptake is mediated by different nitrate transporters (OsNRT1.1B and OsNRT1.5A). This study shows that rice possesses plastic transport systems for mineral elements in response to different water conditions (upland and flooding).

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.

Comparative transcriptome analysis of nodules of two Mesorhizobium-chickpea associations with differential symbiotic efficiency under phosphate deficiency.

Phosphate (Pi) deficiency is known to be a major limitation for symbiotic nitrogen fixation (SNF), and hence legume crop productivity globally. However, very little information is available on the adaptive mechanisms, particularly in the important legume crop chickpea (Cicer arietinum L.), which enable nodules to respond to low-Pi availability. Thus, to elucidate these mechanisms in chickpea nodules at molecular level, we used an RNA sequencing approach to investigate transcriptomes of the nodules in Mesorhizobium mediterraneum SWRI9-(MmSWRI9)-chickpea and M. ciceri CP-31-(McCP-31)-chickpea associations under Pi-sufficient and Pi-deficient conditions, of which the McCP-31-chickpea association has a better SNF capacity than the MmSWRI9-chickpea association during Pi starvation. Our investigation revealed that more genes showed altered expression patterns in MmSWRI9-induced nodules than in McCP-31-induced nodules (540 vs. 225) under Pi deficiency, suggesting that the Pi-starvation-more-sensitive MmSWRI9-induced nodules required expression change in a larger number of genes to cope with low-Pi stress than the Pi-starvation-less-sensitive McCP-31-induced nodules. The functional classification of differentially expressed genes (DEGs) was examined to gain an understanding of how chickpea nodules respond to Pi starvation, caused by soil Pi deficiency. As a result, more DEGs involved in nodulation, detoxification, nutrient/ion transport, transcriptional factors, key metabolic pathways, Pi remobilization and signalling were found in Pi-starved MmSWRI9-induced nodules than in Pi-starved McCP-31-induced nodules. Our findings have enabled the identification of molecular processes that play important roles in the acclimation of nodules to Pi deficiency, ultimately leading to the development of Pi-efficient chickpea symbiotic associations suitable for Pi-deficient soils.

Salicylic acid-dependent immunity contributes to resistance against Rhizoctonia solani, a necrotrophic fungal agent of sheath blight, in rice and Brachypodium distachyon.

Rhizoctonia solani is a soil-borne fungus causing sheath blight. In consistent with its necrotrophic life style, no rice cultivars fully resistant to R. solani are known, and agrochemical plant defense activators used for rice blast, which upregulate a phytohormonal salicylic acid (SA)-dependent pathway, are ineffective towards this pathogen. As a result of the unavailability of genetics, the infection process of R. solani remains unclear. We used the model monocotyledonous plants Brachypodium distachyon and rice, and evaluated the effects of phytohormone-induced resistance to R. solani by pharmacological, genetic and microscopic approaches to understand fungal pathogenicity. Pretreatment with SA, but not with plant defense activators used in agriculture, can unexpectedly induce sheath blight resistance in plants. SA treatment inhibits the advancement of R. solani to the point in the infection process in which fungal biomass shows remarkable expansion and specific infection machinery is developed. The involvement of SA in R. solani resistance is demonstrated by SA-deficient NahG transgenic rice and the sheath blight-resistant B. distachyon accessions, Bd3-1 and Gaz-4, which activate SA-dependent signaling on inoculation. Our findings suggest a hemi-biotrophic nature of R. solani, which can be targeted by SA-dependent plant immunity. Furthermore, B. distachyon provides a genetic resource that can confer disease resistance against R. solani to plants.

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.

Transcriptome Analysis of Chloris virgata, Which Shows the Fastest Germination and Growth in the Major Mongolian Grassland Plant.

Plants in Mongolian grasslands are exposed to short, dry summers and long, cold winters. These plants should be prepared for fast germination and growth activity in response to the limited summer rainfall. The wild plant species adapted to the Mongolian grassland environment may allow us to explore useful genes, as a source of unique genetic codes for crop improvement. Here, we identified the Chloris virgata Dornogovi accession as the fastest germinating plant in major Mongolian grassland plants. It germinated just 5 h after treatment for germination initiation and showed rapid growth, especially in its early and young development stages. This indicates its high growth potential compared to grass crops such as rice and wheat. By assessing growth recovery after animal bite treatment (mimicked by cutting the leaves with scissors), we found that C. virgata could rapidly regenerate leaves after being damaged, suggesting high regeneration potential against grazing. To analyze the regulatory mechanism involved in the high growth potential of C. virgata, we performed RNA-seq-based transcriptome analysis and illustrated a comprehensive gene expression map of the species. Through de novo transcriptome assembly with the RNA-seq reads from whole organ samples of C. virgata at the germination stage (2 days after germination, DAG), early young development stage (8 DAG), young development stage (17 DAG), and adult development stage (28 DAG), we identified 21,589 unified transcripts (contigs) and found that 19,346 and 18,156 protein-coding transcripts were homologous to those in rice and Arabidopsis, respectively. The best-aligned sequences were annotated with gene ontology groups. When comparing the transcriptomes across developmental stages, we found an over-representation of genes involved in growth regulation in the early development stage in C. virgata. Plant development is tightly regulated by phytohormones such as brassinosteroids, gibberellic acid, abscisic acid, and strigolactones. Moreover, our transcriptome map demonstrated the expression profiles of orthologs involved in the biosynthesis of these phytohormones and their signaling networks. We discuss the possibility that C. virgata phytohormone signaling and biosynthesis genes regulate early germination and growth advantages. Comprehensive transcriptome information will provide a useful resource for gene discovery and facilitate a deeper understanding of the diversity of the regulatory systems that have evolved in C. virgata while adapting to severe environmental conditions.

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.

Identification of effector candidate genes of Rhizoctonia solani AG-1 IA expressed during infection in Brachypodium distachyon.

Rhizoctonia solani is a necrotrophic phytopathogen belonging to basidiomycetes. It causes rice sheath blight which inflicts serious damage in rice production. The infection strategy of this pathogen remains unclear. We previously demonstrated that salicylic acid-induced immunity could block R. solani AG-1 IA infection in both rice and Brachypodium distachyon. R. solani may undergo biotrophic process using effector proteins to suppress host immunity before necrotrophic stage. To identify pathogen genes expressed at the early infection process, here we developed an inoculation method using B. distachyon which enables to sample an increased amount of semi-synchronous infection hyphae. Sixty-one R. solani secretory effector-like protein genes (RsSEPGs) were identified using in silico approach with the publicly available gene annotation of R. solani AG-1 IA genome and our RNA-sequencing results obtained from hyphae grown on agar medium. Expression of RsSEPGs was analyzed at 6, 10, 16, 24, and 32 h after inoculation by a quantitative reverse transcription-polymerase chain reaction and 52 genes could be detected at least on a single time point tested. Their expressions showed phase-specific patterns which were classified into 6 clusters. The 23 RsSEPGs in the cluster 1-3 and 29 RsSEPGs in the cluster 4-6 are expected to be involved in biotrophic and necrotrophic interactions, respectively.

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.

De novo transcriptome analysis reveals an unperturbed transcriptome under high cadmium conditions in the Cd-hypertolerant fern Athyrium yokoscense.

Athyrium yokoscense shows strong tolerance to cadmium exposure, even at levels that are many times greater than the toxic levels in ordinary plants. To determine the mechanism of Cd tolerance in A. yokoscense, we grew these plants under high Cd conditions and observed the tissue-specific accumulation of Cd and generation of reactive oxygen species, which is one of the major physiological responses to Cd stress. Fuchsin staining indicated the existence of a casparian strip in A. yokoscense roots, which may participate in Cd hypertolerance in A. yokoscense. Moreover, we performed RNA-seq of RNA samples from A. yokoscense plants treated with or without Cd exposure and obtained comprehensive RNA sequences as well as the Cd-responsive expression patterns of individual genes. Through de novo transcriptome assembly and gene expression analysis, we found that A. yokoscense showed normal features with no significant change in the expression levels of any transporter genes, even under high Cd exposure conditions. Our results demonstrate that A. yokoscense has an unusual mechanism that allows the invading Cd to partition into the distal roots, thus avoiding translocation of Cd into the xylem.

Computer vision-based phenotyping for improvement of plant productivity: a machine learning perspective.

Employing computer vision to extract useful information from images and videos is becoming a key technique for identifying phenotypic changes in plants. Here, we review the emerging aspects of computer vision for automated plant phenotyping. Recent advances in image analysis empowered by machine learning-based techniques, including convolutional neural network-based modeling, have expanded their application to assist high-throughput plant phenotyping. Combinatorial use of multiple sensors to acquire various spectra has allowed us to noninvasively obtain a series of datasets, including those related to the development and physiological responses of plants throughout their life. Automated phenotyping platforms accelerate the elucidation of gene functions associated with traits in model plants under controlled conditions. Remote sensing techniques with image collection platforms, such as unmanned vehicles and tractors, are also emerging for large-scale field phenotyping for crop breeding and precision agriculture. Computer vision-based phenotyping will play significant roles in both the nowcasting and forecasting of plant traits through modeling of genotype/phenotype relationships.

Characterization of sulfur-compound metabolism underlying wax-ester fermentation in Euglena gracilis.

Euglena gracilis is a microalga, which has been used as a model organism for decades. Recent technological advances have enabled mass cultivation of this species for industrial applications such as feedstock in nutritional foods and cosmetics. E. gracilis degrades its storage polysaccharide (paramylon) under hypoxic conditions for energy acquisition by an oxygen-independent process and accumulates high amount of wax-ester as a by-product. Using this sequence of reactions referred to as wax-ester fermentation, E. gracilis is studied for its application in biofuel production. Although the wax-ester production pathway is well characterized, little is known regarding the biochemical reactions underlying the main metabolic route, especially, the existence of an unknown sulfur-compound metabolism implied by the nasty odor generation accompanying the wax-ester fermentation. In this study, we show sulfur-metabolomics of E. gracilis in aerobic and hypoxic conditions, to reveal the biochemical reactions that occur during wax-ester synthesis. Our results helped us in identifying hydrogen sulfide (H2S) as the nasty odor-producing component in wax-ester fermentation. In addition, the results indicate that glutathione and protein degrades during hypoxia, whereas cysteine, methionine, and their metabolites increase in the cells. This indicates that this shift of abundance in sulfur compounds is the cause of H2S synthesis.

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.

Statistical and Machine Learning Approaches to Predict Gene Regulatory Networks From Transcriptome Datasets.

Statistical and machine learning (ML)-based methods have recently advanced in construction of gene regulatory network (GRNs) based on high-throughput biological datasets. GRNs underlie almost all cellular phenomena; hence, comprehensive GRN maps are essential tools to elucidate gene function, thereby facilitating the identification and prioritization of candidate genes for functional analysis. High-throughput gene expression datasets have yielded various statistical and ML-based algorithms to infer causal relationship between genes and decipher GRNs. This review summarizes the recent advancements in the computational inference of GRNs, based on large-scale transcriptome sequencing datasets of model plants and crops. We highlight strategies to select contextual genes for GRN inference, and statistical and ML-based methods for inferring GRNs based on transcriptome datasets from plants. Furthermore, we discuss the challenges and opportunities for the elucidation of GRNs based on large-scale datasets obtained from emerging transcriptomic applications, such as from population-scale, single-cell level, and life-course transcriptome analyses.

Benzothiadiazole, a plant defense inducer, negatively regulates sheath blight resistance in Brachypodium distachyon.

Plant defense inducers that mimic functions of the plant immune hormone salicylic acid (SA) often affect plant growth. Although benzothiadiazole (BTH), a synthetic analog of SA, has been widely used to protect crops from diseases by inducing plant defense responses, we recently demonstrated that SA, but not BTH, confers resistance against Rhizoctonia solani, the causal agent of sheath blight disease, in Brachypodium distachyon. Here, we demonstrated that BTH compromised the resistance of Bd3-1 and Gaz4, the two sheath blight-resistant accessions of B. distachyon, which activate SA-dependent signaling following challenge by R. solani. Moreover, upon analyzing our published RNA-seq data from B. distachyon treated with SA or BTH, we found that BTH specifically induces expression of genes related to chloroplast function and jasmonic acid (JA) signaling, suggesting that BTH attenuates R. solani resistance by perturbing growth-defense trade-offs and/or by inducing a JA response that may increase susceptibility to R. solani. Our findings demonstrated that BTH does not work as a simple mimic of SA in B. distachyon, and consequently may presumably cause unfavorable side effects through the transcriptional alteration, particularly with respect to R. solani resistance.