EST sequencing and phylogenetic analysis of the model grass Brachypodium distachyon

Comparative transcriptomic analysis reveals key components controlling spathe color in Anthurium andraeanum (Hort.)

Anthurium andraeanum (Hort.) is an important ornamental in the tropical cut-flower industry. However, there is currently insufficient information to establish a clear connection between the genetic model(s) proposed and the putative genes involved in the differentiation between colors. In this study, 18 cDNA libraries related to the spathe color and developmental stages of A. andraeanum were characterized by transcriptome sequencing (RNA-seq). For the de novo transcriptome, a total of 114,334,082 primary sequence reads were obtained from the Illumina sequencer and were assembled into 151,652 unigenes. Approximately 58,476 transcripts were generated and used for comparative transcriptome analysis between three cultivars that differ in spathe color (‘Sasha’ (white), ‘Honduras’ (red), and ‘Rapido’ (purple)). A large number of differentially expressed genes (8,324), potentially involved in multiple biological and metabolic pathways, were identified, including genes in the flavonoid and anthocyanin biosynthetic pathways. Our results showed that the chalcone isomerase (CHI) gene presented the strongest evidence for an association with differences in color and the highest correlation with other key genes (flavanone 3-hydroxylase (F3H), flavonoid 3’5’ hydroxylase (F3’5’H)/ flavonoid 3’-hydroxylase (F3’H), and leucoanthocyanidin dioxygenase (LDOX)) in the anthocyanin pathway. We also identified a differentially expressed cytochrome P450 gene in the late developmental stage of the purple spathe that appeared to determine the difference between the red- and purple-colored spathes. Furthermore, transcription factors related to putative MYB-domain protein that may control anthocyanin pathway were identified through a weighted gene co-expression network analysis (WGCNA). The results provided basic sequence information for future research on spathe color, which have important implications for this ornamental breeding strategies.

Contrasting effects of Miocene and Anthropocene levels of atmospheric CO2 on silicon accumulation in a model grass

Grasses are hyper-accumulators of silicon (Si), which they acquire from the soil and deposit in tissues to resist environmental stresses. Given the high metabolic costs of herbivore defensive chemicals and structural constituents (e.g. cellulose), grasses may substitute Si for these components when carbon is limited. Indeed, high Si uptake grasses evolved in the Miocene when atmospheric CO₂ concentration was much lower than present levels. It is, however, unknown how pre-industrial CO₂ concentrations affect Si accumulation in grasses. Using Brachypodium distachyon, we hydroponically manipulated Si-supply (0.0, 0.5, 1, 1.5, 2 mM) and grew plants under Miocene (200 ppm) and Anthropocene levels of CO₂ comprising ambient (410 ppm) and elevated (640 ppm) CO₂ concentrations. We showed that regardless of Si treatments, the Miocene CO₂ levels increased foliar Si concentrations by 47% and 56% relative to plants grown under ambient and elevated CO₂, respectively. This is owing to higher accumulation overall, but also the reallocation of Si from the roots into the shoots. Our results suggest that grasses may accumulate high Si concentrations in foliage when carbon is less available (i.e. pre-industrial CO₂ levels) but this is likely to decline under future climate change scenarios, potentially leaving grasses more susceptible to environmental stresses.

Roles of BdUNICULME4 and BdLAXATUM-A in the non-domesticated grass Brachypodium distachyon

In cultivated grasses, tillering, spike architecture and seed shattering represent major agronomical traits. In barley, maize and rice, the NOOT-BOP-COCH-LIKE (NBCL) genes play important roles in development, especially in ligule development, tillering and flower identity. However, compared with dicots, the role of grass NBCL genes is underinvestigated. To better understand the role of grass NBCLs and to overcome any effects of domestication that might conceal their original functions, we studied TILLING nbcl mutants in the non-domesticated grass Brachypodium distachyon. In B. distachyon, the NBCL genes BdUNICULME4 (CUL4) and BdLAXATUM-A (LAXA) are orthologous, respectively, to the barley HvUniculme4 and HvLaxatum-a, to the maize Zmtassels replace upper ears1 and Zmtassels replace upper ears2 and to the rice OsBLADE-ON-PETIOLE1 and OsBLADE-ON-PETIOLE2/3. In B. distachyon, our reverse genetics study shows that CUL4 is not essential for the establishment of the blade-sheath boundary but is necessary for the development of the ligule and auricles. We report that CUL4 also exerts a positive role in tillering and a negative role in spikelet meristem activity. On the other hand, we demonstrate that LAXA plays a negative role in tillering, positively participates in spikelet development and contributes to the control of floral organ number and identity. In this work, we functionally characterized two new NBCL genes in a context of non-domesticated grass and highlighted original roles for grass NBCL genes that are related to important agronomical traits.

Proteogenomic analysis of pitaya reveals cold stress-related molecular signature

Pitayas (Hylocereus spp.) is an attractive, highly nutritious and commercially valuable tropical fruit. However, low-temperature damage limits crop production. Genome of pitaya has not been sequenced yet. In this study, we sequenced the transcriptome of pitaya as the reference and further investigated the proteome under low temperature. By RNAseq technique, approximately 25.3 million reads were obtained, and further trimmed and assembled into 81,252 unigene sequences. The unigenes were searched against UniProt, NR and COGs at NCBI, Pfam, InterPro and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and 57,905 unigenes were retrieved annotations. Among them, 44,337 coding sequences were predicted by Trandecoder (v2.0.1), which served as the reference database for label-free proteomic analysis study of pitaya. Here, we identified 116 Differentially Abundant Proteins (DAPs) associated with the cold stress in pitaya, of which 18 proteins were up-regulated and 98 proteins were down-regulated. KEGG analysis and other results showed that these DAPs mainly related to chloroplasts and mitochondria metabolism. In summary, chloroplasts and mitochondria metabolism-related proteins may play an important role in response to cold stress in pitayas.

Genome-wide identification of the glutathione transferase superfamily in the model organism Brachypodium distachyon

The detoxification of harmful metabolites can determine the effectiveness of plant stress responses. Scavenging some of these toxic stress by-products through the reduced form of glutathione is catalysed by members of the glutathione transferase (GST) enzyme superfamily. The involvement of these enzymes was studied in the model organism Brachypodium distachyon (L.)P.Beauv. Bd21 and in its derivative Bd21-3, a more drought tolerant line. Osmotic stress treatment resulted in a decrease in the water potential of both Brachypodium genotypes, the difference between the control and treated plant's ψw decreased by the last sampling day in Bd21-3, suggesting some degree of adaptation to the applied osmotic stress. Increased GST activity revealed a severe defence reaction against the harmful imbalance of the redox environment. Screening for the gene sequences led to the identification of 91 full-length or partial GST sequences. Although purple false brome has a relatively small genome, the number of identified GST genes was almost as high as the number predicted in wheat. The estimation of GST expression showed stress-induced differences: higher expression levels or the fast induction of BdGSTF8, BdGSTU35 and BdGSTU42 gene products presumably indicate a strong detoxification under osmotic stress.

Isopentenyl Transferase (IPT) Gene Transfer to Perennial Ryegrass Through Sonication-Assisted Agrobacterium-Mediated Transformation (SAAT), Vacuum and Heat Treatment

The successful introduction of isopentenyl transferase (IPT) gene into perennial ryegrass, cultivars Numan and Grassland using Agrobacterium tumefaciens via three explants (callus, seed and meristem tip) under three individual experiment was evaluated. In the first experiment, the calli were inoculated with LBA4404 Agrobacterium strain under vacuum, heat and in combination of both at 42 °C for 5 min followed by vacuum treatment (390 mm Hg pressure) for 15 min. Sonication-assisted Agrobacterium-mediated transformation (SAAT) was applied for seed and meristem tip transformation of perennial ryegrass for the first time. Results showed positive effects of heat treatment on transformation efficiency during Agro-infection in both cultivars. However, heat shock treatment was more effective in 'Grassland' than 'Numan' (14.2% vs 9.2%). In addition, high transformation efficiency of about 46.65% and 29.15% was observed using meristem tip explants of 'Grassland' and 'Numan' based on IPT and RD29A positive PCR results, respectively. Seed transformation efficiency in 'Grassland' and 'Numan' under SAAT method reached to 37.5% and 16.65%, respectively. Results of these experiments revealed that LBA4404 strain was more efficient than GV3101 in transformation of both perennial ryegrass cultivars. The DNA-blot analysis confirmed that a single T-DNA copy of the IPT gene was integrated into the genomic DNA of the positive transgenic T0 plants which obtained from callus and meristem tip explants of 'Grassland' after heat and SAAT treatment, respectively. Because monocots are not the host of Agrobacterium tumefaciens, this novel protocol can be used in further experiments on genetic transformation of perennial ryegrass cultivars.

Application of Brachypodium genotypes to the analysis of type II resistance to Fusarium head blight (FHB)

The resistance to Fusarium head blight (FHB) in wheat is mainly via the restrain of fungal expansion through spike rachis (type II resistance). In order to unravel the resistance mechanisms, Brachypodium distachyon 21 (Bd21), a monocotyledonous model plant, was previously proved to interact with F. graminearum, while the disease development in spike still needs to be explored in detail. Herein, it is found that the fungal spores mainly germinate on pistil of Bd21, then the hyphae rapidly extend to the bottom of floret and enter spike rachis, similar with the infection progress in wheat. However, structural difference of spike rachis was found between Brachypodium and wheat. It was found that the spread of the fungus through the rachis node of inoculated spikelets is an important index for the evaluation of type II FHB resistance in Brachypodium under optimal conditions at 28 °C and 50%-70% humidity. To verify the feasibility of this strategy, the transcription factor TaTGA2 was overexpressed in Bd21, and transgenic plants were found to show improved resistance to F. graminearum in both spikes and detached leaves, which was further supported by the increased disease severity when silencing TaTGA2 in the wheat cultivar "Sumai 3" or in tilling "Kronos" mutants. Except for Bd21, another 49 Brachypodium germplasms were further screened for FHB resistance, and three moderately susceptible germplasms, namely, PI 317418, W6-39284, and PI 254868, feasible for transformation, were determined to be better hosts than Bd21 when evaluating heterologous genes that positively regulate FHB resistance. The present study also observed variations in the levels of FHB resistance between coleoptiles and spikes or transgenic plants and natural germplasms.

Wheat root length and not branching is altered in the presence of neighbours, including blackgrass

The effect of neighbouring plants on crop root system architecture may directly interfere with water and nutrient acquisition, yet this important and interesting aspect of competition remains poorly understood. Here, the effect of the weed blackgrass (Alopecurus myosuroides Huds.) on wheat (Triticum aestivum L.) roots was tested, since a low density of this species (25 plants m-2) can lead to a 10% decrease in wheat yield and herbicide resistance is problematic. We used a simplified growth system based on gelled medium, to grow wheat alongside a neighbour, either another wheat plant, a blackgrass or Brachypodium dystachion individual (a model grass). A detailed analysis of wheat seminal root system architecture showed that the presence of a neighbour principally affected the root length, rather than number or diameter under a high nutrient regime. In particular, the length of first order lateral roots decreased significantly in the presence of blackgrass and Brachypodium. However, this effect was not noted when wheat plants were grown in low nutrient conditions. This suggests that wheat may be less sensitive to the presence of blackgrass when grown in low nutrient conditions. In addition, nutrient availability to the neighbour did not modulate the neighbour effect on wheat root architecture.

Interactive effects of water limitation and elevated temperature on the physiology, development and fitness of diverse accessions of Brachypodium distachyon

An enduring question in plant physiology and evolution is how single genotypes of plants optimize performance in diverse, often highly variable, environments. We grew 35 natural accessions of the grass Brachypodium distachyon in four environments in the glasshouse, contrasting soil water deficit, elevated temperature and their interaction. We modeled treatment, genotype and interactive effects on leaf-level and whole-plant traits, including fecundity. We also assessed the relationship between glasshouse-measured traits and parameters related to climate at the place of origin. We found abundant genetic variation in both constitutive and induced traits related to plant-water relations. Most traits showed strong interaction between temperature and water availability, and we observed genotype-by-environment interaction for several traits. Notably, leaf free proline abundance showed a strong effect of genotype × temperature × water. We found strong associations between phenology, biomass and water use efficiency (WUE) with parameters describing climate of origin. Plants respond to multiple stressors in ways not directly predictable from single stressors, underscoring the complex and trait-specific mechanisms of environmental response. Climate-trait correlations support a role for WUE and phenology in local adaptation to climate in B. distachyon.

Lignification of Sheepgrass Internodes at Different Developmental Stages and Associated Alteration of Cell Wall Saccharification Efficiency

Sheepgrass (Leymus chinensis) is a high-quality cool-season forage crop used as pasture and hay for livestock feeds. The presence of lignin in cell walls, however, impairs forage digestibility of such lignocellulosic feedstock. Here, the structural characterization and cell wall composition of sheepgrass internodes were studied, and a progressive increase in cell wall lignification was observed with internode maturation. Lignin composition analysis further revealed a gradual accumulation of guaiacyl and syringyl lignin units during internode development. Consistently, the transcript abundance of lignin-related genes was upregulated in mature internodes, suggesting their potential roles in lignin biosynthesis. Furthermore, the effects of cell wall composition and lignification extent on biomass saccharification efficiency were examined in sheepgrass. The results showed that lignin content, guaiacyl and syringyl lignin unit levels inversely correlated with cell wall digestibility, indicating that lignin is a crucial obstacle for utilizing sheepgrass feedstock. The baseline information obtained in this work will facilitate establishment, grazing management, harvesting and feedstock utilization of sheepgrass in future.

Characterization of globulin storage proteins of a low prolamin cereal species in relation to celiac disease

Brachypodium distachyon, a small annual grass with seed storage globulins as primary protein reserves was used in our study to analyse the toxic nature of non-prolamin seed storage proteins related to celiac disease. The main storage proteins of B. distachyon are the 7S globulin type proteins and the 11S, 12S seed storage globulins similar to oat and rice. Immunoblot analyses using serum samples from celiac disease patients were carried out followed by the identification of immune-responsive proteins using mass spectrometry. Serum samples from celiac patients on a gluten-free diet, from patients with Crohn’s disease and healthy subjects, were used as controls. The identified proteins with intense serum-IgA reactivity belong to the 7S and 11–12S seed globulin family. Structure prediction and epitope predictions analyses confirmed the presence of celiac disease-related linear B cell epitope homologs and the presence of peptide regions with strong HLA-DQ8 and DQ2 binding capabilities. These results highlight that both MHC-II presentation and B cell response may be developed not only to prolamins but also to seed storage globulins. This is the first study of the non-prolamin type seed storage proteins of Brachypodium from the aspect of the celiac disease.

Genome-Wide Identification, Evolution, and Co-expression Network Analysis of Mitogen-Activated Protein Kinase Kinase Kinases in Brachypodium distachyon

Mitogen-activated protein kinase (MAPK) cascades are the conserved and universal signal transduction modules in all eukaryotes, which play the vital roles in plant growth, development, and in response to multiple stresses. In this study, we used bioinformatics methods to identify 86 MAPKKK protein encoded by 73 MAPKKK genes in Brachypodium. Phylogenetic analysis of MAPKKK family from Arabidopsis, rice, and Brachypodium has classified them into three subfamilies, of which 28 belonged to MEKK, 52 to Raf, and 6 to ZIK subfamily, respectively. Conserved protein motif, exon-intron organization, and splicing intron phase in kinase domains supported the evolutionary relationships inferred from the phylogenetic analysis. And gene duplication analysis suggested the chromosomal segment duplication happened before the divergence of the rice and Brachypodium, while all of three tandem duplicated gene pairs happened after their divergence. We further demonstrated that the MAPKKKs have evolved under strong purifying selection, implying the conservation of them. The splicing transcripts expression analysis showed that the splicesome translating longest protein tended to be adopted. Furthermore, the expression analysis of BdMAPKKKs in different organs and development stages as well as heat, virus and drought stresses revealed that the MAPKKK genes were involved in various signaling pathways. And the circadian analysis suggested there were 41 MAPKKK genes in Brachypodium showing cycled expression in at least one condition, of which seven MAPKKK genes expressed in all conditions and the promoter analysis indicated these genes possessed many cis-acting regulatory elements involved in circadian and light response. Finally, the co-expression network of MAPK, MAPKK, and MAPKKK in Brachypodium was constructed using 144 microarray and RNA-seq datasets, and ten potential MAPK cascades pathway were predicted. To conclude, our study provided the important information for evolutionary and functional characterization of MAPKKK family in Brachypodium, which will facilitate the functional analysis of BdMAPKKK genes, and also will facilitate better understanding the MAPK signal pathway in Brachypodium and beyond.

Mutation in Brachypodium caffeic acid O-methyltransferase 6 alters stem and grain lignins and improves straw saccharification without deteriorating grain quality

Cereal crop by-products are a promising source of renewable raw material for the production of biofuel from lignocellulose. However, their enzymatic conversion to fermentable sugars is detrimentally affected by lignins. Here the characterization of the Brachypodium Bd5139 mutant provided with a single nucleotide mutation in the caffeic acid O-methyltransferase BdCOMT6 gene is reported. This BdCOMT6-deficient mutant displayed a moderately altered lignification in mature stems. The lignin-related BdCOMT6 gene was also found to be expressed in grains, and the alterations of Bd5139 grain lignins were found to mirror nicely those evidenced in stem lignins. The Bd5139 grains displayed similar size and composition to the control. Complementation experiments carried out by introducing the mutated gene into the AtCOMT1-deficient Arabidopsis mutant demonstrated that the mutated BdCOMT6 protein was still functional. Such a moderate down-regulation of lignin-related COMT enzyme reduced the straw recalcitrance to saccharification, without compromising the vegetative or reproductive development of the plant.

Cell Wall Composition and Biomass Recalcitrance Differences Within a Genotypically Diverse Set of Brachypodium distachyon Inbred Lines

Brachypodium distachyon (Brachypodium) has emerged as a useful model system for studying traits unique to graminaceous species including bioenergy crop grasses owing to its amenability to laboratory experimentation and the availability of extensive genetic and germplasm resources. Considerable natural variation has been uncovered for a variety of traits including flowering time, vernalization responsiveness, and above-ground growth characteristics. However, cell wall composition differences remain underexplored. Therefore, we assessed cell wall-related traits relevant to biomass conversion to biofuels in seven Brachypodium inbred lines that were chosen based on their high level of genotypic diversity as well as available genome sequences and recombinant inbred line (RIL) populations. Senesced stems plus leaf sheaths from these lines exhibited significant differences in acetyl bromide soluble lignin (ABSL), cell wall polysaccharide-derived sugars, hydroxycinnamates content, and syringyl:guaiacyl:p-hydroxyphenyl (S:G:H) lignin ratios. Free glucose, sucrose, and starch content also differed significantly in senesced stems, as did the amounts of sugars released from cell wall polysaccharides (digestibility) upon exposure to a panel of thermochemical pretreatments followed by hydrolytic enzymatic digestion. Correlations were identified between inbred line lignin compositions and plant growth characteristics such as biomass accumulation and heading date (HD), and between amounts of cell wall polysaccharides and biomass digestibility. Finally, stem cell wall p-coumarate and ferulate contents and free-sugars content changed significantly with increased duration of vernalization for some inbred lines. Taken together, these results show that Brachypodium displays substantial phenotypic variation with respect to cell wall composition and biomass digestibility, with some compositional differences correlating with growth characteristics. Moreover, besides influencing HD and biomass accumulation, vernalization was found to affect cell wall composition and free sugars accumulation in some Brachypodium inbred lines, suggesting genetic differences in how vernalization affects carbon flux to polysaccharides. The availability of related RIL populations will allow for the genetic and molecular dissection of this natural variation, the knowledge of which may inform ways to genetically improve bioenergy crop grasses.

Molecular characterization and evolutionary origins of farinin genes in Brachypodium distachyon L

Farinins are one of the oldest members of the gluten family in wheat and Aegilops species, and they influence dough properties. Here, we performed the first detailed molecular genetic study on farinin genes in Brachypodium distachyon L., the model species for Triticum aestivum. A total of 51 b-type farinin genes were cloned and characterized, including 27 functional and 24 non-functional pseudogenes from 14 different B. distachyon accessions. All genes were highly similar to those previously reported from wheat and Aegilops species. The identification of deduced amino acid sequences showed that b-type farinins across Triticeae genomes could be classified as b1-, b2-, b3-, and b4-type farinins; however, B. distachyon had only b3- and b4-type farinins. Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) revealed that farinin genes are transcribed into mRNA in B. distachyon at much lower levels than in Triticeae, despite the presence of cis-acting elements in promoter regions. Phylogenetic analysis suggested that Brachypodium farinins may have closer relationships with common wheat and further confirmed four different types of b-type farinins in Triticeae and Brachypodium genomes, corresponding to b1, b2, b3 (group 1), and b4 (group 2). A putative evolutionary origin model of farinin genes in Brachypodium, Triticum, and the related species suggests that all b-type farinins diverged from their common ancestor ~3.2 million years ago (MYA). The b3 and b4 types could be considered older in the farinin family. The results explain the loss of b1- and b2-type farinin alleles in Brachypodium.

Glucanocellulosic ethanol: the undiscovered biofuel potential in energy crops and marine biomass

Converting biomass to biofuels is a key strategy in substituting fossil fuels to mitigate climate change. Conventional strategies to convert lignocellulosic biomass to ethanol address the fermentation of cellulose-derived glucose. Here we used super-resolution fluorescence microscopy to uncover the nanoscale structure of cell walls in the energy crops maize and Miscanthus where the typical polymer cellulose forms an unconventional layered architecture with the atypical (1, 3)-β-glucan polymer callose. This raised the question about an unused potential of (1, 3)-β-glucan in the fermentation of lignocellulosic biomass. Engineering biomass conversion for optimized (1, 3)-β-glucan utilization, we increased the ethanol yield from both energy crops. The generation of transgenic Miscanthus lines with an elevated (1, 3)-β-glucan content further increased ethanol yield providing a new strategy in energy crop breeding. Applying the (1, 3)-β-glucan-optimized conversion method on marine biomass from brown macroalgae with a naturally high (1, 3)-β-glucan content, we not only substantially increased ethanol yield but also demonstrated an effective co-fermentation of plant and marine biomass. This opens new perspectives in combining different kinds of feedstock for sustainable and efficient biofuel production, especially in coastal regions.

Pushing the boundaries of resistance: insights from Brachypodium-rust interactions

The implications of global population growth urge transformation of current food and bioenergy production systems to sustainability. Members of the family Poaceae are of particular importance both in food security and for their applications as biofuel substrates. For centuries, rust fungi have threatened the production of valuable crops such as wheat, barley, oat, and other small grains; similarly, biofuel crops can also be susceptible to these pathogens. Emerging rust pathogenic races with increased virulence and recurrent rust epidemics around the world point out the vulnerability of monocultures. Basic research in plant immunity, especially in model plants, can make contributions to understanding plant resistance mechanisms and improve disease management strategies. The development of the grass Brachypodium distachyon as a genetically tractable model for monocots, especially temperate cereals and grasses, offers the possibility to overcome the experimental challenges presented by the genetic and genomic complexities of economically valuable crop plants. The numerous resources and tools available in Brachypodium have opened new doors to investigate the underlying molecular and genetic bases of plant-microbe interactions in grasses and evidence demonstrating the applicability and advantages of working with B. distachyon is increasing. Importantly, several interactions between B. distachyon and devastating plant pathogens, such rust fungi, have been examined in the context of non-host resistance. Here, we discuss the use of B. distachyon in these various pathosystems. Exploiting B. distachyon to understand the mechanisms underpinning disease resistance to non-adapted rust fungi may provide effective and durable approaches to fend off these pathogens. The close phylogenetic relationship among Brachypodium spp. and grasses with industrial and agronomic value support harnessing this model plant to improve cropping systems and encourage its use in translational research.

Reduced susceptibility to Fusarium head blight in Brachypodium distachyon through priming with the Fusarium mycotoxin deoxynivalenol

The fungal cereal pathogen Fusarium graminearum produces deoxynivalenol (DON) during infection. The mycotoxin DON is associated with Fusarium head blight (FHB), a disease that can cause vast grain losses. Whilst investigating the suitability of Brachypodium distachyon as a model for spreading resistance to F. graminearum, we unexpectedly discovered that DON pretreatment of spikelets could reduce susceptibility to FHB in this model grass. We started to analyse the cell wall changes in spikelets after infection with F. graminearum wild-type and defined mutants: the DON-deficient Δtri5 mutant and the DON-producing lipase disruption mutant Δfgl1, both infecting only directly inoculated florets, and the mitogen-activated protein (MAP) kinase disruption mutant Δgpmk1, with strongly decreased virulence but intact DON production. At 14 days post-inoculation, the glucose amounts in the non-cellulosic cell wall fraction were only increased in spikelets infected with the DON-producing strains wild-type, Δfgl1 and Δgpmk1. Hence, we tested for DON-induced cell wall changes in B. distachyon, which were most prominent at DON concentrations ranging from 1 to 100 ppb. To test the involvement of DON in defence priming, we pretreated spikelets with DON at a concentration of 1 ppm prior to F. graminearum wild-type infection, which significantly reduced FHB disease symptoms. The analysis of cell wall composition and plant defence-related gene expression after DON pretreatment and fungal infection suggested that DON-induced priming of the spikelet tissue contributed to the reduced susceptibility to FHB.

Brachypodium distachyon and Setaria viridis: Model Genetic Systems for the Grasses

The family of grasses encompasses the world's most important food, feed, and bioenergy crops, yet we are only now beginning to develop the genetic resources to explore the diversity of form and function that underlies economically important traits. Two emerging model systems, Brachypodium distachyon and Setaria viridis, promise to greatly accelerate the process of gene discovery in the grasses and to serve as bridges in the exploration of panicoid and pooid grasses, arguably two of the most important clades of plants from a food security perspective. We provide both a historical view of the development of plant model systems and highlight several recent reports that are providing these developing communities with the tools for gene discovery and pathway engineering.

Comparative analysis of antiviral responses in Brachypodium distachyon and Setaria viridis reveals conserved and unique outcomes among C3 and C4 plant defenses

Viral diseases cause significant losses in global agricultural production, yet little is known about grass antiviral defense mechanisms. We previously reported on host immune responses triggered by Panicum mosaic virus (PMV) and its satellite virus (SPMV) in the model C3 grass Brachypodium distachyon. To aid comparative analyses of C3 and C4 grass antiviral defenses, here, we establish B. distachyon and Setaria viridis (a C4 grass) as compatible hosts for seven grass-infecting viruses, including PMV and SPMV, Brome mosaic virus, Barley stripe mosaic virus, Maize mild mottle virus, Sorghum yellow banding virus, Wheat streak mosaic virus (WSMV), and Foxtail mosaic virus (FoMV). Etiological and molecular characterization of the fourteen grass-virus pathosystems showed evidence for conserved crosstalk among salicylic acid (SA), jasmonic acid, and ethylene pathways in B. distachyon and S. viridis. Strikingly, expression of PHYTOALEXIN DEFICIENT4, an upstream modulator of SA signaling, was consistently suppressed during most virus infections in B. distachyon and S. viridis. Hierarchical clustering analyses further identified unique antiviral responses triggered by two morphologically similar viruses, FoMV and WSMV, and uncovered other host-dependent effects. Together, the results of this study establish B. distachyon and S. viridis as models for the analysis of plant-virus interactions and provide the first framework for conserved and unique features of C3 and C4 grass antiviral defenses.

Phylogeny in defining model plants for lignocellulosic ethanol production: a comparative study of Brachypodium distachyon, wheat, maize, and Miscanthus x giganteus leaf and stem biomass

The production of ethanol from pretreated plant biomass during fermentation is a strategy to mitigate climate change by substituting fossil fuels. However, biomass conversion is mainly limited by the recalcitrant nature of the plant cell wall. To overcome recalcitrance, the optimization of the plant cell wall for subsequent processing is a promising approach. Based on their phylogenetic proximity to existing and emerging energy crops, model plants have been proposed to study bioenergy-related cell wall biochemistry. One example is Brachypodium distachyon, which has been considered as a general model plant for cell wall analysis in grasses. To test whether relative phylogenetic proximity would be sufficient to qualify as a model plant not only for cell wall composition but also for the complete process leading to bioethanol production, we compared the processing of leaf and stem biomass from the C3 grasses B. distachyon and Triticum aestivum (wheat) with the C4 grasses Zea mays (maize) and Miscanthus x giganteus, a perennial energy crop. Lambda scanning with a confocal laser-scanning microscope allowed a rapid qualitative analysis of biomass saccharification. A maximum of 108-117 mg ethanol·g(-1) dry biomass was yielded from thermo-chemically and enzymatically pretreated stem biomass of the tested plant species. Principal component analysis revealed that a relatively strong correlation between similarities in lignocellulosic ethanol production and phylogenetic relation was only given for stem and leaf biomass of the two tested C4 grasses. Our results suggest that suitability of B. distachyon as a model plant for biomass conversion of energy crops has to be specifically tested based on applied processing parameters and biomass tissue type.

Digital imaging approaches for phenotyping whole plant nitrogen and phosphorus response in Brachypodium distachyon

This work evaluates the phenotypic response of the model grass (Brachypodium distachyon (L.) P. Beauv.) to nitrogen and phosphorus nutrition using a combination of imaging techniques and destructive harvest of shoots and roots. Reference line Bd21-3 was grown in pots using 11 phosphorus and 11 nitrogen concentrations to establish a dose-response curve. Shoot biovolume and biomass, root length and biomass, and tissue phosphorus and nitrogen concentrations increased with nutrient concentration. Shoot biovolume, estimated by imaging, was highly correlated with dry weight (R(2) > 0.92) and both biovolume and growth rate responded strongly to nutrient availability. Higher nutrient supply increased nodal root length more than other root types. Photochemical efficiency was strongly reduced by low phosphorus concentrations as early as 1 week after germination, suggesting that this measurement may be suitable for high throughput screening of phosphorus response. In contrast, nitrogen concentration had little effect on photochemical efficiency. Changes in biovolume over time were used to compare growth rates of four accessions in response to nitrogen and phosphorus supply. We demonstrate that a time series image-based approach coupled with mathematical modeling provides higher resolution of genotypic response to nutrient supply than traditional destructive techniques and shows promise for high throughput screening and determination of genomic regions associated with superior nutrient use efficiency.

De novo assembly, gene annotation, and simple sequence repeat marker development using Illumina paired-end transcriptome sequences in the pearl oyster Pinctada maxima

We analyzed the mantle transcriptome of pearl oyster Pinctada maxima and developed EST-SSR markers using Illumina HiSeq 2000 paired-end sequencing technology. A total of 49,500,748 raw reads were generated. De novo assembly generated 108,704 unigenes with an average length of 407 bp. Sequence similarity search with known proteins or nucleotides revealed that 30,200 (27.78%) and 25,824 (23.76%) consensus sequences were homologous with the sequences in the non-redundant protein and Swiss-Prot databases, respectively, and that 19,701 (18.12%) of these unigenes were possibly involved in approximately 234 known signaling pathways in the Kyoto Encyclopedia of Genes and Genomes database. Ninety one biomineralization-related unigenes were detected. In a cultured stock, 1764 simple sequence repeats were identified and 56 primer pairs were randomly selected and tested. The rate of successful amplification was 68.3%. The developed molecular markers are helpful for further studies on genetic linkage analysis, gene localization, and quantitative trait loci mapping.

Brachypodium distachyon exhibits compatible interactions with Oculimacula spp. and Ramularia collo-cygni, providing the first pathosystem model to study eyespot and ramularia leaf spot diseases

Brachypodium distachyon (Bd) has established itself as an essential tool for comparative genomic studies in cereals and increasing attention is being paid to its potential as a model pathosystem. Eyespot and ramularia leaf spot (RLS) are important diseases of wheat, barley and other small-grain cereals for which very little is known about the mechanisms of host resistance despite urgent requirements for plant breeders to develop resistant varieties. This work aimed to test the compatibility of interaction of two Bd accessions with the cereal pathogens Oculimacula spp. and Ramularia collo-cygni, the causal agents of eyespot and RLS diseases, respectively. Results showed that both Bd accessions developed symptoms similar to those on the natural host for all pathogen species tested. Microscopy images demonstrated that R. collo-cygni produced secondary conidia and both Oculimacula spp. formed characteristic infection structures on successive tissue layers. Visual disease assessment revealed that quantitative differences in disease severity exist between the two Bd accessions. The results presented here provide the first evidence that Bd is compatible with the main causal agents of eyespot and RLS diseases, and suggest that future functional genetic studies can be undertaken to investigate the mechanisms of eyespot and RLS disease resistance using Bd.

Transcriptome sequencing and gene expression analysis of Trichoderma brevicompactum under different culture conditions

Background

Trichoderma brevicompactum is the Trichoderma species producing simple trichothecenes-trichodermin, a potential antifungal antibiotic and a protein synthesis inhibitor. However, the biosynthetic pathway of trichodermin in Trichoderma is not completely clarified. Therefore, transcriptome and gene expression profiling data for this species are needed as an important resource to better understand the mechanism of the trichodermin biosynthesis and provide a blueprint for further study of T. brevicompactum.

Results

In this study, de novo assembly of the T. brevicompactum transcriptome using the short-read sequencing technology (Illumina) was performed. In addition, two digital gene expression (DGE) libraries of T. brevicompactum under the trichodermin-producing and trichodermin-nonproducing culture conditions, respectively, were constructed to identify the differences in gene expression. A total of 23,351 unique transcripts with a mean length of 856 bp were obtained by a new Trinity de novo assembler. The variations of the gene expression under different culture conditions were also identified. The expression profiling data revealed that 3,282 unique transcripts had a significantly differential expression under the trichodermin-producing condition, as compared to the trichodermin-nonproducing condition. This study provides a large amount of transcript sequence data that will contribute to the study of the trichodermin biosynthesis in T. brevicompactum. Furthermore, quantitative real-time PCR (qRT-PCR) was found to be useful to confirm the differential expression of the unique transcripts.

Conclusion

Our study provides considerable gene expression information of T. brevicompactum at the transcriptional level,which will help accelerate the research on the trichodermin biosynthesis. Additionally, we have demonstrated the feasibility of using the Illumina sequencing based DGE system for gene expression profiling, and have shed new light on functional studies of the genes involved in T. brevicompactum biosynthesis.

p-Coumaroyl-CoA:monolignol transferase (PMT) acts specifically in the lignin biosynthetic pathway in Brachypodium distachyon

Grass lignins contain substantial amounts of p-coumarate (pCA) that acylate the side-chains of the phenylpropanoid polymer backbone. An acyltransferase, named p-coumaroyl-CoA:monolignol transferase (OsPMT), that could acylate monolignols with pCA in vitro was recently identified from rice. In planta, such monolignol-pCA conjugates become incorporated into lignin via oxidative radical coupling, thereby generating the observed pCA appendages; however p-coumarates also acylate arabinoxylans in grasses. To test the authenticity of PMT as a lignin biosynthetic pathway enzyme, we examined Brachypodium distachyon plants with altered BdPMT gene function. Using newly developed cell wall analytical methods, we determined that the transferase was involved specifically in monolignol acylation. A sodium azide-generated Bdpmt-1 missense mutant had no (<0.5%) residual pCA on lignin, and BdPMT RNAi plants had levels as low as 10% of wild-type, whereas the amounts of pCA acylating arabinosyl units on arabinoxylans in these PMT mutant plants remained unchanged. pCA acylation of lignin from BdPMT-overexpressing plants was found to be more than three-fold higher than that of wild-type, but again the level on arabinosyl units remained unchanged. Taken together, these data are consistent with a defined role for grass PMT genes in encoding BAHD (BEAT, AHCT, HCBT, and DAT) acyltransferases that specifically acylate monolignols with pCA and produce monolignol p-coumarate conjugates that are used for lignification in planta.

Molecular characterisation and evolution of HMW glutenin subunit genes in Brachypodium distachyon L

Brachypodium distachyon, a small wild grass within the Pooideae family, is a new model organism for exploring the functional genomics of cereal crops. It was shown to have close relationships to wheat, barley and rice. Here, we describe the molecular characterisation and evolutionary relationships of high molecular weight glutenin subunits (HMW-GS) genes from B. distachyon. Sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE), high performance capillary electrophoresis (HPCE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses demonstrated that there was no HMW-GS expression in the Brachypodium grains due to the silencing of their encoding genes. Through allele-specific polymerase chain reaction (AS-PCR) amplification and cloning, a total of 13 HMW-GS encoding genes from diploid, tetraploid and hexaploid Brachypodium species were obtained, and all of them had typical structural features of y-type HMW-GS genes from common wheat and related species, particularly more similar to the 1Dy12 gene. However, the presence of an in-frame premature stop codon (TAG) at position 1521 in the coding region resulted in the conversion of all the genes to pseudogenes. Further, quantitative real-time PCR (qRT-PCR) analysis revealed that HMW-GS genes in B. distachyon displayed a similar trend, but with a low transcriptional expression profile during grain development due to the occurrence of the stop codon. Phylogenetic analysis showed that the highly conserved Glu-1-2 loci were presented in B. distachyon, which displayed close phylogenetic evolutionary relationships with Triticum and related species.

Brachypodium sylvaticum, a model for perennial grasses: transformation and inbred line development

Perennial species offer significant advantages as crops including reduced soil erosion, lower energy inputs after the first year, deeper root systems that access more soil moisture, and decreased fertilizer inputs due to the remobilization of nutrients at the end of the growing season. These advantages are particularly relevant for emerging biomass crops and it is projected that perennial grasses will be among the most important dedicated biomass crops. The advantages offered by perennial crops could also prove favorable for incorporation into annual grain crops like wheat, rice, sorghum and barley, especially under the dryer and more variable climate conditions projected for many grain-producing regions. Thus, it would be useful to have a perennial model system to test biotechnological approaches to crop improvement and for fundamental research. The perennial grass Brachypodiumsylvaticum is a candidate for such a model because it is diploid, has a small genome, is self-fertile, has a modest stature, and short generation time. Its close relationship to the annual model Brachypodiumdistachyon will facilitate comparative studies and allow researchers to leverage the resources developed for B. distachyon. Here we report on the development of two keystone resources that are essential for a model plant: high-efficiency transformation and inbred lines. Using Agrobacterium tumefaciens-mediated transformation we achieved an average transformation efficiency of 67%. We also surveyed the genetic diversity of 19 accessions from the National Plant Germplasm System using SSR markers and created 15 inbred lines.

Transcriptome analysis in sheepgrass (Leymus chinensis): a dominant perennial grass of the Eurasian Steppe

BACKGROUND

Sheepgrass [Leymus chinensis (Trin.) Tzvel.] is an important perennial forage grass across the Eurasian Steppe and is known for its adaptability to various environmental conditions. However, insufficient data resources in public databases for sheepgrass limited our understanding of the mechanism of environmental adaptations, gene discovery and molecular marker development.

RESULTS

The transcriptome of sheepgrass was sequenced using Roche 454 pyrosequencing technology. We assembled 952,328 high-quality reads into 87,214 unigenes, including 32,416 contigs and 54,798 singletons. There were 15,450 contigs over 500 bp in length. BLAST searches of our database against Swiss-Prot and NCBI non-redundant protein sequences (nr) databases resulted in the annotation of 54,584 (62.6%) of the unigenes. Gene Ontology (GO) analysis assigned 89,129 GO term annotations for 17,463 unigenes. We identified 11,675 core Poaceae-specific and 12,811 putative sheepgrass-specific unigenes by BLAST searches against all plant genome and transcriptome databases. A total of 2,979 specific freezing-responsive unigenes were found from this RNAseq dataset. We identified 3,818 EST-SSRs in 3,597 unigenes, and some SSRs contained unigenes that were also candidates for freezing-response genes. Characterizations of nucleotide repeats and dominant motifs of SSRs in sheepgrass were also performed. Similarity and phylogenetic analysis indicated that sheepgrass is closely related to barley and wheat.

CONCLUSIONS

This research has greatly enriched sheepgrass transcriptome resources. The identified stress-related genes will help us to decipher the genetic basis of the environmental and ecological adaptations of this species and will be used to improve wheat and barley crops through hybridization or genetic transformation. The EST-SSRs reported here will be a valuable resource for future gene-phenotype studies and for the molecular breeding of sheepgrass and other Poaceae species.

Genome-wide landscape of alternative splicing events in Brachypodium distachyon

Recently, Brachypodium distachyon has emerged as a model plant for studying monocot grasses and cereal crops. Using assembled expressed transcript sequences and subsequent mapping to the corresponding genome, we identified 1219 alternative splicing (AS) events spanning across 2021 putatively assembled transcripts generated from 941 genes. Approximately, 6.3% of expressed genes are alternatively spliced in B. distachyon. We observed that a majority of the identified AS events were related to retained introns (55.5%), followed by alternative acceptor sites (16.7%). We also observed a low percentage of exon skipping (5.0%) and alternative donor site events (8.8%). The 'complex event' that consists of a combination of two or more basic splicing events accounted for ∼14.0%. Comparative AS transcript analysis revealed 163 and 39 homologous pairs between B. distachyon and Oryza sativa and between B. distachyon and Arabidopsis thaliana, respectively. In all, we found 16 AS transcripts to be conserved in all 3 species. AS events and related putative assembled transcripts annotation can be systematically browsed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/plant/).

MORE SPIKELETS1 is required for spikelet fate in the inflorescence of Brachypodium

Grasses produce florets on a structure called a spikelet, and variation in the number and arrangement of both branches and spikelets contributes to the great diversity of grass inflorescence architecture. In Brachypodium (Brachypodium distachyon), the inflorescence is an unbranched spike with a terminal spikelet and a limited number of lateral spikelets. Spikelets are indeterminate and give rise to a variable number of florets. Here, we provide a detailed description of the stages of inflorescence development in Brachypodium. To gain insight into the genetic regulation of Brachypodium inflorescence development, we generated fast neutron mutant populations and screened for phenotypic mutants. Among the mutants identified, the more spikelets1 (mos1) mutant had an increased number of axillary meristems produced from inflorescence meristem compared with the wild type. These axillary meristems developed as branches with production of higher order spikelets. Using a candidate gene approach, mos1 was found to have a genomic rearrangement disrupting the expression of an ethylene response factor class of APETALA2 transcription factor related to the spikelet meristem identity genes branched silkless1 (bd1) in maize (Zea mays) and FRIZZY PANICLE (FZP) in rice (Oryza sativa). We propose MOS1 likely corresponds to the Brachypodium bd1 and FZP ortholog and that the function of this gene in determining spikelet meristem fate is conserved with distantly related grass species. However, MOS1 also appears to be involved in the timing of initiation of the terminal spikelet. As such, MOS1 may regulate the transition to terminal spikelet development in other closely related and agriculturally important species, particularly wheat (Triticum aestivum).

Disrupting the cinnamyl alcohol dehydrogenase 1 gene (BdCAD1) leads to altered lignification and improved saccharification in Brachypodium distachyon

Brachypodium distachyon (Brachypodium) has been proposed as a model for grasses, but there is limited knowledge regarding its lignins and no data on lignin-related mutants. The cinnamyl alcohol dehydrogenase (CAD) genes involved in lignification are promising targets to improve the cellulose-to-ethanol conversion process. Down-regulation of CAD often induces a reddish coloration of lignified tissues. Based on this observation, we screened a chemically induced population of Brachypodium mutants (Bd21-3 background) for red culm coloration. We identified two mutants (Bd4179 and Bd7591), with mutations in the BdCAD1 gene. The mature stems of these mutants displayed reduced CAD activity and lower lignin content. Their lignins were enriched in 8-O-4- and 4-O-5-coupled sinapaldehyde units, as well as resistant inter-unit bonds and free phenolic groups. By contrast, there was no increase in coniferaldehyde end groups. Moreover, the amount of sinapic acid ester-linked to cell walls was measured for the first time in a lignin-related CAD grass mutant. Functional complementation of the Bd4179 mutant with the wild-type BdCAD1 allele restored the wild-type phenotype and lignification. Saccharification assays revealed that Bd4179 and Bd7591 lines were more susceptible to enzymatic hydrolysis than wild-type plants. Here, we have demonstrated that BdCAD1 is involved in lignification of Brachypodium. We have shown that a single nucleotide change in BdCAD1 reduces the lignin level and increases the degree of branching of lignins through incorporation of sinapaldehyde. These changes make saccharification of cells walls pre-treated with alkaline easier without compromising plant growth.

Tapping the promise of genomics in species with complex, nonmodel genomes

Genomics is enabling a renaissance in all disciplines of plant biology. However, many plant genomes are complex and remain recalcitrant to current genomic technologies. The complexities of these nonmodel plant genomes are attributable to gene and genome duplication, heterozygosity, ploidy, and/or repetitive sequences. Methods are available to simplify the genome and reduce these barriers, including inbreeding and genome reduction, making these species amenable to current sequencing and assembly methods. Some, but not all, of the complexities in nonmodel genomes can be bypassed by sequencing the transcriptome rather than the genome. Additionally, comparative genomics approaches, which leverage phylogenetic relatedness, can aid in the interpretation of complex genomes. Although there are limitations in accessing complex nonmodel plant genomes using current sequencing technologies, genome manipulation and resourceful analyses can allow access to even the most recalcitrant plant genomes.

Molecular characterization of LMW-GS genes in Brachypodium distachyon L. reveals highly conserved Glu-3 loci in Triticum and related species

BACKGROUND

Brachypodium distachyon L. is a newly emerging model plant system for temperate cereal crop species. However, its grain protein compositions are still not clear. In the current study, we carried out a detailed proteomics and molecular genetics study on grain glutenin proteins in B. distachyon.

RESULTS

SDS-PAGE and RP-HPLC analysis of grain proteins showed that Brachypodium has few gliadins and high molecular weight glutenin subunits. In contrast the electrophoretic patterns for the albumin, globulin and low molecular weight glutenin subunit (LMW-GS) fractions of the grain protein were similar to those in wheat. In particular, the LMW-C type subunits in Brachypodium were more abundant than the equivalent proteins in common wheat. Southern blotting analysis confirmed that Brachypodium has 4-5 copies of LMW-GS genes. A total of 18 LMW-GS genes were cloned from Brachypodium by allele specific PCR. LMW-GS and 4 deduced amino acid sequences were further confirmed by using Western-blotting and MALDI-TOF-MS. Phylogenetic analysis indicated that Brachypodium was closer to Ae. markgrafii and Ae. umbellulata than to T. aestivum.

CONCLUSIONS

Brachypodium possessed a highly conserved Glu-3 locus that is closely related to Triticum and related species. The presence of LMW-GS in B. distachyon grains indicates that B. distachyon may be used as a model system for studying wheat quality attributes.

The family of DOF transcription factors in Brachypodium distachyon: phylogenetic comparison with rice and barley DOFs and expression profiling

BACKGROUND

Transcription factors (TFs) are proteins that have played a central role both in evolution and in domestication, and are major regulators of development in living organisms. Plant genome sequences reveal that approximately 7% of all genes encode putative TFs. The DOF (DNA binding with One Finger) TF family has been associated with vital processes exclusive to higher plants and to their close ancestors (algae, mosses and ferns). These are seed maturation and germination, light-mediated regulation, phytohormone and plant responses to biotic and abiotic stresses, etc. In Hordeum vulgare and Oryza sativa, 26 and 30 different Dof genes, respectively, have been annotated. Brachypodium distachyon has been the first Pooideae grass to be sequenced and, due to its genomic, morphological and physiological characteristics, has emerged as the model system for temperate cereals, such as wheat and barley.

RESULTS

Through searches in the B. distachyon genome, 27 Dof genes have been identified and a phylogenetic comparison with the Oryza sativa and the Hordeum vulgare DOFs has been performed. To explore the evolutionary relationship among these DOF proteins, a combined phylogenetic tree has been constructed with the Brachypodium DOFs and those from rice and barley. This phylogenetic analysis has classified the DOF proteins into four Major Cluster of Orthologous Groups (MCOGs). Using RT-qPCR analysis the expression profiles of the annotated BdDof genes across four organs (leaves, roots, spikes and seeds) has been investigated. These results have led to a classification of the BdDof genes into two groups, according to their expression levels. The genes highly or preferentially expressed in seeds have been subjected to a more detailed expression analysis (maturation, dry stage and germination).

CONCLUSIONS

Comparison of the expression profiles of the Brachypodium Dof genes with the published functions of closely related DOF sequences from the cereal species considered here, deduced from the phylogenetic analysis, indicates that although the expression profile has been conserved in many of the putative orthologs, in some cases duplication followed by subsequent divergence may have occurred (neo-functionalization).

Brachypodium distachyon as a model for defining the allergen potential of non-prolamin proteins

Epitope databases and the protein sequences of published plant genomes are suitable to identify some of the proteins causing food allergies and sensitivities. Brachypodium distachyon, a diploid wild grass with a sequenced genome and low prolamin content, is the closest relative of the allergen cereals, such as wheat or barley. Using the Brachypodium genome sequence, a workflow has been developed to identify potentially harmful proteins which may cause either celiac disease or wheat allergy-related symptoms. Seed tissue-specific expression of the potential allergens has been determined, and intact epitopes following an in silico digestion with several endopeptidases have been identified. Molecular function of allergen proteins has been evaluated using Gene Ontology terms. Biologically overrepresented proteins and potentially allergen protein families have been identified.

Natural genetic variation of root system architecture from Arabidopsis to Brachypodium: towards adaptive value

Root system architecture is a trait that displays considerable plasticity because of its sensitivity to environmental stimuli. Nevertheless, to a significant degree it is genetically constrained as suggested by surveys of its natural genetic variation. A few regulators of root system architecture have been isolated as quantitative trait loci through the natural variation approach in the dicotyledon model, Arabidopsis. This provides proof of principle that allelic variation for root system architecture traits exists, is genetically tractable, and might be exploited for crop breeding. Beyond Arabidopsis, Brachypodium could serve as both a credible and experimentally accessible model for root system architecture variation in monocotyledons, as suggested by first glimpses of the different root morphologies of Brachypodium accessions. Whether a direct knowledge transfer gained from molecular model system studies will work in practice remains unclear however, because of a lack of comprehensive understanding of root system physiology in the native context. For instance, apart from a few notable exceptions, the adaptive value of genetic variation in root system modulators is unknown. Future studies should thus aim at comprehensive characterization of the role of genetic players in root system architecture variation by taking into account the native environmental conditions, in particular soil characteristics.

Fine mapping of the Bsr1 barley stripe mosaic virus resistance gene in the model grass Brachypodium distachyon

The ND18 strain of Barley stripe mosaic virus (BSMV) infects several lines of Brachypodium distachyon, a recently developed model system for genomics research in cereals. Among the inbred lines tested, Bd3-1 is highly resistant at 20 to 25°C, whereas Bd21 is susceptible and infection results in an intense mosaic phenotype accompanied by high levels of replicating virus. We generated an F_6∶7 recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21 and used the RILs, and an F₂ population of a second Bd21 × Bd3-1 cross to evaluate the inheritance of resistance. The results indicate that resistance segregates as expected for a single dominant gene, which we have designated Barley stripe mosaic virus resistance 1 (Bsr1). We constructed a genetic linkage map of the RIL population using SNP markers to map this gene to within 705 Kb of the distal end of the top of chromosome 3. Additional CAPS and Indel markers were used to fine map Bsr1 to a 23 Kb interval containing five putative genes. Our study demonstrates the power of using RILs to rapidly map the genetic determinants of BSMV resistance in Brachypodium. Moreover, the RILs and their associated genetic map, when combined with the complete genomic sequence of Brachypodium, provide new resources for genetic analyses of many other traits.

Identification and comparative mapping of a powdery mildew resistance gene derived from wild emmer (Triticum turgidum var. dicoccoides) on chromosome 2BS

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is an important foliar disease of wheat worldwide. Wild emmer (Triticum turgidum var. dicoccoides) is a valuable genetic resource for improving disease resistance in common wheat. A powdery mildew resistance gene conferring resistance to B. graminis f. sp. tritici isolate E09 at the seedling and adult stages was identified in wild emmer accession IW170 introduced from Israel. An incomplete dominant gene, temporarily designated MlIW170, was responsible for the resistance. Through molecular marker and bulked segregant analyses of an F(2) population and F(3) families derived from a cross between susceptible durum wheat line 81086A and IW170, MlIW170 was located in the distal chromosome bin 2BS3-0.84-1.00 and flanked by SSR markers Xcfd238 and Xwmc243. MlIW170 co-segregated with Xcau516, an STS marker developed from RFLP marker Xwg516 that co-segregated with powdery mildew resistance gene Pm26 on 2BS. Four EST-STS markers, BE498358, BF201235, BQ160080, and BF146221, were integrated into the genetic linkage map of MlIW170. Three AFLP markers, XPaacMcac, XPagcMcta, XPaacMcag, and seven AFLP-derived SCAR markers, XcauG2, XcauG3, XcauG6, XcauG8, XcauG10, XcauG20, and XcauG25, were linked to MlIW170. XcauG3, a resistance gene analog (RGA)-like sequence, co-segregated with MlIW170. The non-glaucousness locus Iw1 was 18.77 cM distal to MlIW170. By comparative genomics of wheat-Brachypodium-rice genomic co-linearity, four EST-STS markers, CJ658408, CJ945509, BQ169830, CJ945085, and one STS marker XP2430, were developed and MlIW170 was mapped in an 2.69 cM interval that is co-linear with a 131 kb genomic region in Brachypodium and a 105 kb genomic region in rice. Four RGA-like sequences annotated in the orthologous Brachypodium genomic region could serve as chromosome landing target regions for map-based cloning of MlIW170.

QTLs for resistance to the false brome rust Puccinia brachypodii in the model grass Brachypodium distachyon L

The potential of the model grass Brachypodium distachyon L. (Brachypodium) for studying grass-pathogen interactions is still underexploited. We aimed to identify genomic regions in Brachypodium associated with quantitative resistance to the false brome rust fungus Puccinia brachypodii . The inbred lines Bd3-1 and Bd1-1, differing in their level of resistance to P. brachypodii, were crossed to develop an F(2) population. This was evaluated for reaction to a virulent isolate of P. brachypodii at both the seedling and advanced growth stages. To validate the results obtained on the F(2), resistance was quantified in F(2)-derived F(3) families in two experiments. Disease evaluations showed quantitative and transgressive segregation for resistance. A new AFLP-based Brachypodium linkage map consisting of 203 loci and spanning 812 cM was developed and anchored to the genome sequence with SSR and SNP markers. Three false brome rust resistance QTLs were identified on chromosomes 2, 3, and 4, and they were detected across experiments. This study is the first quantitative trait analysis in Brachypodium. Resistance to P. brachypodii was governed by a few QTLs: two acting at the seedling stage and one acting at both seedling and advanced growth stages. The results obtained offer perspectives to elucidate the molecular basis of quantitative resistance to rust fungi.

Host Status of False Brome Grass to the Leaf Rust Fungus Puccinia brachypodii and the Stripe Rust Fungus P. striiformis

Purple false brome grass (Brachypodium distachyon) has recently emerged as a model system for temperate grasses and is also a potential model plant to investigate plant interactions with economically important pathogens such as rust fungi. We determined the host status of five Brachypodium species to three isolates of Puccinia brachypodii, the prevalent rust species on Brachypodium sylvaticum in nature, and to one isolate each of three formae speciales of the stripe rust fungus P. striiformis. Two P. striiformis isolates produced sporulating lesions, both in only one of the tested interactions, suggesting a marginal host status of B. distachyon. P. brachypodii formed sporulating uredinia on the five Brachypodium species tested, and a range of reactions was observed. Surprisingly, the B. sylvaticum-derived rust isolates were more frequently pathogenic to B. distachyon than to their original host species. The B. distachyon diploid inbred lines, developed and distributed as reference material to the Brachypodium research community, include susceptible and resistant genotypes to at least three of the four P. brachypodii isolates tested. This creates the opportunity to use B. distachyon/P. brachypodii as a model pathosystem. In one B. distachyon accession, heavy infection by the loose smut fungus Ustilago bromivora occurred. That pathogen could also serve as a model pathogen of Brachypodium.

Targeted mining of drought stress-responsive genes from EST resources in Cleistogenes songorica

Cleistogenes songorica is an important perennial grass found in the pastoral steppe of Inner Mongolia. C. songorica flourishes in drought prone environments, and therefore provides an ideal candidate plant system for the identification of drought-tolerance conferring genes. We constructed cDNA libraries from leaves and roots of drought-stressed C. songorica seedlings. Expressed sequence tag (EST) sequencing of 5664 random cDNA clones produced 3579 high quality, trimmed sequences. The average read length of trimmed ESTs was 613bp. Clustering and assembly identified a non-redundant set of 1499 contigs, including 805 singleton unigenes and 694 multi-member unigenes. The resulting unigenes were functionally categorized according to the Gene Ontology (GO) hierarchy using the in house Bioinformatic Advanced Scientific Computing (BASC) annotation pipeline. Among the total 2.2Mbp of EST sequence data, 161 putative SSRs were found, a frequency similar to that previously observed in oat and Arabidopsis ESTs. Sixty-three unigenes were functionally annotated as being stress responsive, of which 22 were similar to genes implicated in drought stress response. Using quantitative real time RT-PCR, transcripts of 13 of these 22 genes were shown to be at least three fold more, or less abundant in drought-stressed leaves or roots, with 8 increased and 5 decreased in relative transcript abundance. The C. songorica EST and cDNA collections generated in this study are a valuable resource for microarray-based expression profiling, and functional genomics in order to elucidate their role, and to understand the underlying mechanisms of drought-tolerance in C. songorica.

Comparison of a high-density genetic linkage map to genome features in the model grass Brachypodium distachyon

The small annual grass Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to study questions unique to the grasses. Many Brachypodium resources have been developed including a whole genome sequence, highly efficient transformation and a large germplasm collection. We developed a genetic linkage map of Brachypodium using single nucleotide polymorphism (SNP) markers and an F(2) mapping population of 476 individuals. SNPs were identified by targeted resequencing of single copy genomic sequences. Using the Illumina GoldenGate Genotyping platform we placed 558 markers into five linkage groups corresponding to the five chromosomes of Brachypodium. The unusually long total genetic map length, 1,598 centiMorgans (cM), indicates that the Brachypodium mapping population has a high recombination rate. By comparing the genetic map to genome features we found that the recombination rate was positively correlated with gene density and negatively correlated with repetitive regions and sites of ancestral chromosome fusions that retained centromeric repeat sequences. A comparison of adjacent genome regions with high versus low recombination rates revealed a positive correlation between interspecific synteny and recombination rate.

Exploiting the Brachypodium Tool Box in cereal and grass research

It is now a decade since Brachypodium distachyon (Brachypodium) was suggested as a model species for temperate grasses and cereals. Since then transformation protocols, large expressed sequence tag (EST) databases, tools for forward and reverse genetic screens, highly refined cytogenetic probes, germplasm collections and, recently, a complete genome sequence have been generated. In this review, we will describe the current status of the Brachypodium Tool Box and how it is beginning to be applied to study a range of biological traits. Further, as genomic analysis of larger cereals and forage grasses genomes are becoming easier, we will re-evaluate Brachypodium as a model species. We suggest that there remains an urgent need to employ reverse genetic and functional genomic approaches to identify the functionality of key genetic elements, which could be employed subsequently in plant breeding programmes; and a requirement for a Pooideae reference genome to aid assembling large pooid genomes. Brachypodium is an ideal system for functional genomic studies, because of its easy growth requirements, small physical stature, and rapid life cycle, coupled with the resources offered by the Brachypodium Tool Box.

Brachypodium distachyon: a new pathosystem to study Fusarium head blight and other Fusarium diseases of wheat

BACKGROUND

Fusarium species cause Fusarium head blight (FHB) and other important diseases of cereals. The causal agents produce trichothecene mycotoxins such as deoxynivalenol (DON). The dicotyledonous model species Arabidopsis thaliana has been used to study Fusarium-host interactions but it is not ideal for model-to-crop translation. Brachypodium distachyon (Bd) has been proposed as a new monocotyledonous model species for functional genomic studies in grass species. This study aims to assess the interaction between the most prevalent FHB-causing Fusarium species and Bd in order to develop and exploit Bd as a genetic model for FHB and other Fusarium diseases of wheat.

RESULTS

The ability of Fusarium graminearum and Fusarium culmorum to infect a range of Bd tissues was examined in various bioassays which showed that both species can infect all Bd tissues examined, including intact foliar tissues. DON accumulated in infected spike tissues at levels similar to those of infected wheat spikes. Histological studies revealed details of infection, colonisation and host response and indicate that hair cells are important sites of infection. Susceptibility to Fusarium and DON was assessed in two Bd ecotypes and revealed variation in resistance between ecotypes.

CONCLUSIONS

Bd exhibits characteristics of susceptibility highly similar to those of wheat, including susceptibility to spread of disease in the spikelets. Bd is the first reported plant species to allow successful infection on intact foliar tissues by FHB-causing Fusarium species. DON appears to function as a virulence factor in Bd as it does in wheat. Bd is proposed as a valuable model for undertaking studies of Fusarium head blight and other Fusarium diseases of wheat.

Construction of high quality Gateway™ entry libraries and their application to yeast two-hybrid for the monocot model plant Brachypodium distachyon

Background

Monocots, especially the temperate grasses, represent some of the most agriculturally important crops for both current food needs and future biofuel development. Because most of the agriculturally important grass species are difficult to study (e.g., they often have large, repetitive genomes and can be difficult to grow in laboratory settings), developing genetically tractable model systems is essential. Brachypodium distachyon (hereafter Brachypodium) is an emerging model system for the temperate grasses. To fully realize the potential of this model system, publicly accessible discovery tools are essential. High quality cDNA libraries that can be readily adapted for multiple downstream purposes are a needed resource. Additionally, yeast two-hybrid (Y2H) libraries are an important discovery tool for protein-protein interactions and are not currently available for Brachypodium.

Results

We describe the creation of two high quality, publicly available Gateway™ cDNA entry libraries and their derived Y2H libraries for Brachypodium. The first entry library represents cloned cDNA populations from both short day (SD, 8/16-h light/dark) and long day (LD, 20/4-h light/dark) grown plants, while the second library was generated from hormone treated tissues. Both libraries have extensive genome coverage (~5 × 10⁷ primary clones each) and average clone lengths of ~1.5 Kb. These entry libraries were then used to create two recombination-derived Y2H libraries. Initial proof-of-concept screens demonstrated that a protein with known interaction partners could readily re-isolate those partners, as well as novel interactors.

Conclusions

Accessible community resources are a hallmark of successful biological model systems. Brachypodium has the potential to be a broadly useful model system for the grasses, but still requires many of these resources. The Gateway™ compatible entry libraries created here will facilitate studies for multiple user-defined purposes and the derived Y2H libraries can be immediately applied to large scale screening and discovery of novel protein-protein interactions. All libraries are freely available for distribution to the research community.

Discovery of barley miRNAs through deep sequencing of short reads

Background

MicroRNAs are important components of the regulatory network of biological systems and thousands have been discovered in both animals and plants. Systematic investigations performed in species with sequenced genomes such as Arabidopsis, rice, poplar and Brachypodium have provided insights into the evolutionary relationships of this class of small RNAs among plants. However, miRNAs from barley, one of the most important cereal crops, remain unknown.

Results

We performed a large scale study of barley miRNAs through deep sequencing of small RNAs extracted from leaves of two barley cultivars. By using the presence of miRNA precursor sequences in related genomes as one of a number of supporting criteria, we identified up to 100 miRNAs in barley. Of these only 56 have orthologs in wheat, rice or Brachypodium that are known to be expressed, while up to 44 appear to be specifically expressed in barley.

Conclusions

Our study, the first large scale investigation of small RNAs in barley, has identified up to 100 miRNAs. We demonstrate that reliable identification of miRNAs via deep sequencing in a species whose genome has not been sequenced requires a more careful analysis of sequencing errors than is commonly performed. We devised a read filtering procedure for dealing with errors. In addition, we found that the use of a large dataset of almost 35 million reads permits the use of read abundance distributions along putative precursor sequences as a practical tool for isolating miRNAs in a large background of reads originating from other non-coding and coding RNAs. This study therefore provides a generic approach for discovering novel miRNAs where no genome sequence is available.