Search

Third-harmonic generation microscopy is a powerful label-free nonlinear imaging technique, providing essential information about structural characteristics of cells and tissues without requiring external labelling agents. In this work, we integrated a recently developed compact adaptive optics module into a third-harmonic generation microscope, to measure and correct for optical aberrations in complex tissues. Taking advantage of the high sensitivity of the third-harmonic generation process to material interfaces and thin membranes, along with the 1,300-nm excitation wavelength used here, our adaptive optical third-harmonic generation microscope enabled high-resolution in vivo imaging within highly scattering biological model systems. Examples include imaging of myelinated axons and vascular structures within the mouse spinal cord and deep cortical layers of the mouse brain, along with imaging of key anatomical features in the roots of the model plant Brachypodium distachyon. In all instances, aberration correction led to significant enhancements in image quality.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Plant growth requires the integration of internal and external cues, perceived and transduced into a developmental programme of cell division, elongation and wall thickening. Mechanical forces contribute to this regulation, and thigmomorphogenesis typically includes reducing stem height, increasing stem diameter, and a canonical transcriptomic response. We present data on a bZIP transcription factor involved in this process in grasses. Brachypodium distachyon SECONDARY WALL INTERACTING bZIP (SWIZ) protein translocated into the nucleus following mechanostimulation. Classical touch-responsive genes were upregulated in B. distachyon roots following touch, including significant induction of the glycoside hydrolase 17 family, which may be unique to grass thigmomorphogenesis. SWIZ protein binding to an E-box variant in exons and introns was associated with immediate activation followed by repression of gene expression. SWIZ overexpression resulted in plants with reduced stem and root elongation. These data further define plant touch-responsive transcriptomics and physiology, offering insights into grass mechanotranduction dynamics.

Perennial grasses are important forage crops and emerging biomass crops and have the potential to be more sustainable grain crops. However, most perennial grass crops are difficult experimental subjects due to their large size, difficult genetics, and/or their recalcitrance to transformation. Thus, a tractable model perennial grass could be used to rapidly make discoveries that can be translated to perennial grass crops. Brachypodium sylvaticum has the potential to serve as such a model because of its small size, rapid generation time, simple genetics, and transformability. Here, we provide a high-quality genome assembly and annotation for B. sylvaticum, an essential resource for a modern model system. In addition, we conducted transcriptomic studies under 4 abiotic stresses (water, heat, salt, and freezing). Our results indicate that crowns are more responsive to freezing than leaves which may help them overwinter. We observed extensive transcriptional responses with varying temporal dynamics to all abiotic stresses, including classic heat-responsive genes. These results can be used to form testable hypotheses about how perennial grasses respond to these stresses. Taken together, these results will allow B. sylvaticum to serve as a truly tractable perennial model system.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

The "genomic shock" hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession "Bhyb26." We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species' similarity in transposable element load may account for the subtlety of the observed genome dominance.

Daylength sensing in many plants is critical for coordinating the timing of flowering with the appropriate season. Temperate climate-adapted grasses such as Brachypodium distachyon flower during the spring when days are becoming longer. The photoreceptor PHYTOCHROME C is essential for long-day (LD) flowering in B. distachyon. PHYC is required for the LD activation of a suite of genes in the photoperiod pathway including PHOTOPERIOD1 (PPD1) that, in turn, result in the activation of FLOWERING LOCUS T (FT1)/FLORIGEN, which causes flowering. Thus, B. distachyon phyC mutants are extremely delayed in flowering. Here we show that PHYC-mediated activation of PPD1 occurs via EARLY FLOWERING 3 (ELF3), a component of the evening complex in the circadian clock. The extreme delay of flowering of the phyC mutant disappears when combined with an elf3 loss-of-function mutation. Moreover, the dampened PPD1 expression in phyC mutant plants is elevated in phyC/elf3 mutant plants consistent with the rapid flowering of the double mutant. We show that loss of PPD1 function also results in reduced FT1 expression and extremely delayed flowering consistent with results from wheat and barley. Additionally, elf3 mutant plants have elevated expression levels of PPD1, and we show that overexpression of ELF3 results in delayed flowering associated with a reduction of PPD1 and FT1 expression, indicating that ELF3 represses PPD1 transcription consistent with previous studies showing that ELF3 binds to the PPD1 promoter. Indeed, PPD1 is the main target of ELF3-mediated flowering as elf3/ppd1 double mutant plants are delayed flowering. Our results indicate that ELF3 operates downstream from PHYC and acts as a repressor of PPD1 in the photoperiod flowering pathway of B. distachyon.

Nitrogen (N) fixation in cereals by root-associated bacteria is a promising solution for reducing use of chemical N fertilizers in agriculture. However, plant and bacterial responses are unpredictable across environments. We hypothesized that cereal responses to N-fixing bacteria are dynamic, depending on N supply and time. To quantify the dynamics, a gnotobiotic, fabricated ecosystem (EcoFAB) was adapted to analyse N mass balance, to image shoot and root growth, and to measure gene expression of Brachypodium distachyon inoculated with the N-fixing bacterium Herbaspirillum seropedicae. Phenotyping throughput of EcoFAB-N was 25-30 plants h-1 with open software and imaging systems. Herbaspirillum seropedicae inoculation of B. distachyon shifted root and shoot growth, nitrate versus ammonium uptake, and gene expression with time; directions and magnitude depended on N availability. Primary roots were longer and root hairs shorter regardless of N, with stronger changes at low N. At higher N, H. seropedicae provided 11% of the total plant N that came from sources other than the seed or the nutrient solution. The time-resolved phenotypic and molecular data point to distinct modes of action: at 5 mM NH4NO3 the benefit appears through N fixation, while at 0.5 mM NH4NO3 the mechanism appears to be plant physiological, with H. seropedicae promoting uptake of N from the root medium.Future work could fine-tune plant and root-associated microorganisms to growth and nutrient dynamics.

Plant response to drought stress involves fungi and bacteria that live on and in plants and in the rhizosphere, yet the stability of these myco- and micro-biomes remains poorly understood. We investigate the resistance and resilience of fungi and bacteria to drought in an agricultural system using both community composition and microbial associations. Here we show that tests of the fundamental hypotheses that fungi, as compared to bacteria, are (i) more resistant to drought stress but (ii) less resilient when rewetting relieves the stress, found robust support at the level of community composition. Results were more complex using all-correlations and co-occurrence networks. In general, drought disrupts microbial networks based on significant positive correlations among bacteria, among fungi, and between bacteria and fungi. Surprisingly, co-occurrence networks among functional guilds of rhizosphere fungi and leaf bacteria were strengthened by drought, and the same was seen for networks involving arbuscular mycorrhizal fungi in the rhizosphere. We also found support for the stress gradient hypothesis because drought increased the relative frequency of positive correlations.

Rapid environmental change can lead to population extinction or evolutionary rescue. The global staple crop sorghum (Sorghum bicolor) has recently been threatened by a global outbreak of an aggressive new biotype of sugarcane aphid (SCA; Melanaphis sacchari). We characterized genomic signatures of adaptation in a Haitian breeding population that had rapidly adapted to SCA infestation, conducting evolutionary population genomics analyses on 296 Haitian lines versus 767 global accessions. Genome scans and geographic analyses suggest that SCA adaptation has been conferred by a globally rare East African allele of RMES1, which spread to breeding programs in Africa, Asia, and the Americas. De novo genome sequencing revealed potential causative variants at RMES1. Markers developed from the RMES1 sweep predicted resistance in eight independent commercial and public breeding programs. These findings demonstrate the value of evolutionary genomics to develop adaptive trait technology and highlight the benefits of global germplasm exchange to facilitate evolutionary rescue.

Domestication and breeding have affected interactions between plants and their microbiomes in ways that are only beginning to be understood but may have important implications for recruitment of rhizosphere microorganisms, particularly under stress conditions. We investigated the responses of a modern maize (Zea mays subsp. mays) cultivar and its wild relative, teosinte (Z. mays subsp. parviglumis), to different phosphate availabilities. We appraised responses of the plant-microbial holobiont to phosphate stresses by profiling root exudate metabolomes, and microbial communities in the root endosphere and rhizosphere. We also performed plate assays to quantify phosphate-solubilizing microorganisms from the rhizosphere. Although root exudate metabolite profiles were distinct between the teosinte and modern maize under high phosphate, both plants shifted exudate compositions in response to phosphate stress toward a common metabolite profile. Root and rhizosphere microbial communities also responded significantly to both plant type and the phosphate availability. A subset of bacterial and fungal taxa were differentially abundant under the different phosphate conditions, with each of the three conditions favoring different taxa. Both teosinte and maize rhizospheres harbored phosphate-solubilizing microorganisms under all growth conditions. These results suggest that the root exudation response to phosphate stress was conserved through the domestication of maize from teosinte, shifting exudation levels of specific metabolites. Although microbial communities also shifted, plate-based assays did not detect selective recruitment of phosphate solubilizers in response to phosphate availability.

BackgroundThe vascular system of plants consists of two main tissue types, xylem and phloem. These tissues are organized into vascular bundles that are arranged into a complex network running through the plant that is essential for the viability of land plants. Despite their obvious importance, the genes involved in the organization of vascular tissues remain poorly understood in grasses.ResultsWe studied in detail the vascular network in stems from the model grass Brachypodium distachyon (Brachypodium) and identified a large set of genes differentially expressed in vascular bundles versus parenchyma tissues. To decipher the underlying molecular mechanisms of vascularization in grasses, we conducted a forward genetic screen for abnormal vasculature. We identified a mutation that severely affected the organization of vascular tissues. This mutant displayed defects in anastomosis of the vascular network and uncommon amphivasal vascular bundles. The causal mutation is a premature stop codon in ERECTA, a LRR receptor-like serine/threonine-protein kinase. Mutations in this gene are pleiotropic indicating that it serves multiple roles during plant development. This mutant also displayed changes in cell wall composition, gene expression and hormone homeostasis.ConclusionIn summary, ERECTA has a pleiotropic role in Brachypodium. We propose a major role of ERECTA in vasculature anastomosis and vascular tissue organization in Brachypodium.

Progress in sequencing, microfluidics, and analysis strategies has revolutionized the granularity at which multicellular organisms can be studied. In particular, single-cell transcriptomics has led to fundamental new insights into animal biology, such as the discovery of new cell types and cell type-specific disease processes. However, the application of single-cell approaches to plants, fungi, algae, or bacteria (environmental organisms) has been far more limited, largely due to the challenges posed by polysaccharide walls surrounding these species' cells. In this perspective, we discuss opportunities afforded by single-cell technologies for energy and environmental science and grand challenges that must be tackled to apply these approaches to plants, fungi and algae. We highlight the need to develop better and more comprehensive single-cell technologies, analysis and visualization tools, and tissue preparation methods. We advocate for the creation of a centralized, open-access database to house plant single-cell data. Finally, we consider how such efforts should balance the need for deep characterization of select model species while still capturing the diversity in the plant kingdom. Investments into the development of methods, their application to relevant species, and the creation of resources to support data dissemination will enable groundbreaking insights to propel energy and environmental science forward.

Arbuscular mycorrhizal fungi (AMF), the mutualistic symbionts with mostcrops, constitute a research system of human-associated fungi whoserelative simplicity and synchrony are conducive to experimental ecology.However, little is known about the shifts in adaptive strategies ofsorghum associated AMFs where strong AMF succession replaces initiallyruderal species with competitive ones and where the strongest plantresponse to drought is to manage these AMF. First, we hypothesize that,when irrigation is stopped to mimic drought, competitive AMF speciesshould be replaced by AMF species tolerant to drought stress. We then,for the first time, correlate AMF abundance and host plant transcriptionto test two novel hypotheses about the mechanisms behind the shift fromruderal to competitive AMF. Surprisingly, despite imposing droughtstress, we found no stress tolerant AMF. Remarkably, we found strong anddifferential correlation between the successional shift from ruderal tocompetitive AMF and sorghum genes whose products (i) produce and releasestrigolactone signals, (ii) perceivemycorrhizal-lipochitinoligosaccharide (Myc-LCO) signals, (iii) provideplant lipid and sugar to AMF and, (iv) import minerals and waterprovided by AMF. These novel insights into host gene expression andsuccession of AMF show adaptive strategies evolved by AMF and theirhosts and provide a rationale for selecting AMF to reduce inputs andmaximize yield in commercial agriculture. Future research opportunitiesinclude testing the specifics and generality of our hypotheses byemploying genetically modified host plants, and exploring additionalgenes underlying the adaptive strategies in natural succession.

Wheat stem rust disease caused by Puccinia graminis f. sp. tritici (Pgt) is a global threat to wheat production. Fast evolving populations of Pgt limit the efficacy of plant genetic resistance and constrain disease management strategies. Understanding molecular mechanisms that lead to rust infection and disease susceptibility could deliver novel strategies to deploy crop resistance through genetic loss of disease susceptibility. We used comparative transcriptome-based and orthology-guided approaches to characterize gene expression changes associated with Pgt infection in susceptible and resistant Triticum aestivum genotypes as well as the non-host Brachypodium distachyon. We targeted our analysis to genes with differential expression in T. aestivum and genes suppressed or not affected in B. distachyon and report several processes potentially linked to susceptibility to Pgt, such as cell death suppression and impairment of photosynthesis. We complemented our approach with a gene co-expression network analysis to identify wheat targets to deliver resistance to Pgt through removal or modification of putative susceptibility genes.

Renewable fuels are needed to replace fossil fuels in the immediate future. Lignocellulosic bioenergy crops provide a renewable alternative that sequesters atmospheric carbon. To prevent displacement of food crops, it would be advantageous to grow biofuel crops on marginal lands. These lands will likely face more frequent and extreme drought conditions than conventional agricultural land, so it is crucial to see how proposed bioenergy crops fare under these conditions and how that may affect lignocellulosic biomass composition and saccharification properties. We found that while drought impacts the plant cell wall of Sorghum bicolor differently according to tissue and timing of drought induction, drought-induced cell wall compositional modifications are relatively minor and produce no negative effect on biomass conversion. This contrasts with the cell wall-related transcriptome, which had a varied range of highly variable genes (HVGs) within four cell wall-related GO categories, depending on the tissues surveyed and time of drought induction. Further, many HVGs had expression changes in which putative impacts were not seen in the physical cell wall or which were in opposition to their putative impacts. Interestingly, most pre-flowering drought-induced cell wall changes occurred in the leaf, with matrix and lignin compositional changes that did not persist after recovery from drought. Most measurable physical post-flowering cell wall changes occurred in the root, affecting mainly polysaccharide composition and cross-linking. This study couples transcriptomics to cell wall chemical analyses of a C4 grass experiencing progressive and differing drought stresses in the field. As such, we can analyze the cell wall-specific response to agriculturally relevant drought stresses on the transcriptomic level and see whether those changes translate to compositional or biomass conversion differences. Our results bolster the conclusion that drought stress does not substantially affect the cell wall composition of specific aerial and subterranean biomass nor impede enzymatic hydrolysis of leaf biomass, a positive result for biorefinery processes. Coupled with previously reported results on the root microbiome and rhizosphere and whole transcriptome analyses of this study, we can formulate and test hypotheses on individual gene candidates' function in mediating drought stress in the grass cell wall, as demonstrated in sorghum.

In recent years, the root microbiome (i.e., microorganisms growing inside, on, or in close proximity to plant roots) has been shown to play an important role in plant health and productivity. Despite its importance, the root microbiome is challenging to study because of its complexity, heterogeneity, and subterranean location. Fortunately, root microbiome research has seen a tremendous influx of novel technologies (e.g., imaging tools, robotics, and molecular analyses), experimental platforms (e.g., micro- and mesocosms), and data integration, modeling, and prediction tools in the past decade that have greatly increased our ability to dissect the complex network of interactions between above- and belowground environmental parameters, plants, bacteria, and fungi that dictate soil and broader ecosystem health. Herein, we discuss methods that are currently used in root microbiome research and that can be expanded to phytobiome research in general ranging from laboratory studies to mesocosm-scale studies and, finally, to field studies; evaluate their relevance to ecosystem studies; and discuss future root microbiome research directions.

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

Perennial grasses will account for approximately 16 billion gallons of renewable fuels by the year 2022, contributing significantly to carbon and nitrogen sequestration. However, perennial grasses productivity can be limited by severe freezing conditions in some geographical areas, although these risks could decrease with the advance of climate warming, the possibility of unpredictable early cold events cannot be discarded. We conducted a study on the model perennial grass Brachypodium sylvaticum to investigate the molecular mechanisms that contribute to cold and freezing adaption. The study was performed on two different B. sylvaticum accessions, Ain1 and Osl1, typical to warm and cold climates, respectively. Both accessions were grown under controlled conditions with subsequent cold acclimation followed by freezing stress. For each treatment a set of morphological parameters, transcription, metabolite, and lipid profiles were measured. State-of-the-art algorithms were employed to analyze cross-component relationships. Phenotypic analysis revealed higher adaption of Osl1 to freezing stress. Our analysis highlighted the differential regulation of the TCA cycle and the GABA shunt between Ain1 and Osl1. Osl1 adapted to freezing stress by repressing the GABA shunt activity, avoiding the detrimental reduction in fatty acid biosynthesis and the concomitant detrimental effects on membrane integrity.

Community assembly of crop-associated fungi is thought to be strongly influenced by deterministic selection exerted by the plant host, rather than stochastic processes. Here we use a simple, sorghum system with abundant sampling to show that stochastic forces (drift or stochastic dispersal) act on fungal community assembly in leaves and roots early in host development and when sorghum is drought stressed, conditions when mycobiomes are small. Unexpectedly, we find no signal for stochasticity when drought stress is relieved, likely due to renewed selection by the host. In our experimental system, the host compartment exerts the strongest effects on mycobiome assembly, followed by the timing of plant development and lastly by plant genotype. Using a dissimilarity-overlap approach, we find a universality in the forces of community assembly of the mycobiomes of the different sorghum compartments and in functional guilds of fungi.

Drought is the most important environmental stress limiting crop yields. The C4 cereal sorghum [Sorghum bicolor (L.) Moench] is a critical food, forage, and emerging bioenergy crop that is notably drought-tolerant. We conducted a large-scale field experiment, imposing preflowering and postflowering drought stress on 2 genotypes of sorghum across a tightly resolved time series, from plant emergence to postanthesis, resulting in a dataset of nearly 400 transcriptomes. We observed a fast and global transcriptomic response in leaf and root tissues with clear temporal patterns, including modulation of well-known drought pathways. We also identified genotypic differences in core photosynthesis and reactive oxygen species scavenging pathways, highlighting possible mechanisms of drought tolerance and of the delayed senescence, characteristic of the stay-green phenotype. Finally, we discovered a large-scale depletion in the expression of genes critical to arbuscular mycorrhizal (AM) symbiosis, with a corresponding drop in AM fungal mass in the plants' roots.

Powdery mildew (Golovinomyces cichoracearum), one of the most prolific obligate biotrophic fungal pathogens worldwide, infects its host by penetrating the plant cell wall without activating the plant's innate immune system. The Arabidopsis mutant powdery mildew resistant 5 (pmr5) carries a mutation in a putative pectin acetyltransferase gene that confers enhanced resistance to powdery mildew. Here, we show that heterologously expressed PMR5 protein transfers acetyl groups from [14 C]-acetyl-CoA to oligogalacturonides. Through site-directed mutagenesis, we show that three amino acids within a highly conserved esterase domain in putative PMR5 orthologs are necessary for PMR5 function. A suppressor screen of mutagenized pmr5 seed selecting for increased powdery mildew susceptibility identified two previously characterized genes affecting the acetylation of plant cell wall polysaccharides, RWA2 and TBR. The rwa2 and tbr mutants also suppress powdery mildew disease resistance in pmr6, a mutant defective in a putative pectate lyase gene. Cell wall analysis of pmr5 and pmr6, and their rwa2 and tbr suppressor mutants, demonstrates minor shifts in cellulose and pectin composition. In direct contrast to their increased powdery mildew resistance, both pmr5 and pmr6 plants are highly susceptibile to multiple strains of the generalist necrotroph Botrytis cinerea, and have decreased camalexin production upon infection with B. cinerea. These results illustrate that cell wall composition is intimately connected to fungal disease resistance and outline a potential route for engineering powdery mildew resistance into susceptible crop species.

BackgroundOur understanding of polyploid genomes is limited by our inability to definitively assign sequences to a specific subgenome without extensive prior knowledge like high resolution genetic maps or genome sequences of diploid progenitors. In theory, existing methods for assigning sequences to individual species from metagenome samples could be used to separate subgenomes in polyploid organisms, however, these methods rely on differences in coarse genome properties like GC content or sequences from related species. Thus, these approaches do not work for subgenomes where gross features are indistinguishable and related genomes are lacking. Here we describe a method that uses rapidly evolving repetitive DNA to circumvent these limitations.ResultsBy using short, repetitive, DNA sequences as species-specific signals we separated closely related genomes from test datasets and subgenomes from two polyploid plants, tobacco and wheat, without any prior knowledge.ConclusionThis approach is ideal for separating the subgenomes of polyploid species with unsequenced or unknown progenitor genomes.

Maize domestication and breeding have resulted in drastic and well documented changes in aboveground traits, but belowground effects on root system functioning and rhizosphere microbial communities remain poorly understood, despite their critical importance for nutrient and water acquisition. We investigated the rhizosphere microbial community composition and structure of ten Zea mays accessions along an evolutionary transect (two teosinte, three inbred maize lines, and five modern maize hybrids) grown in nutrient depleted soil from a low input agricultural system. Microbial community analysis revealed significant differences in community composition between soil compartments (proximal vs. distal rhizosphere) and between plant genetic groups (teosinte, inbred, and modern hybrid). Only a small portion of the microbial community was differentially selected across plant genetic groups: 3.7% of prokaryotic community members and 4.9% of fungal community members were significantly associated with a specific plant genetic group. Indicator species analysis showed the greatest differentiation between modern hybrids and the other two plant genetic groups. Co-occurrence network analysis revealed that microbial co-occurrence patterns of the inbred maize lines' rhizosphere were significantly more similar to those of the teosintes than to the modern hybrids. Our results suggest that advances in hybrid development significantly impacted rhizosphere microbial communities and network assembly.

There is a dynamic reciprocity between plants and their environment: soil physiochemical properties influence plant morphology and metabolism, and root morphology and exudates shape the environment surrounding roots. Here, we investigate the reproducibility of plant trait changes in response to three growth environments. We utilized fabricated ecosystem (EcoFAB) devices to grow the model grass Brachypodium distachyon in three distinct media across four laboratories: phosphate-sufficient and -deficient mineral media allowed assessment of the effects of phosphate starvation, and a complex, sterile soil extract represented a more natural environment with yet uncharacterized effects on plant growth and metabolism. Tissue weight and phosphate content, total root length, and root tissue and exudate metabolic profiles were consistent across laboratories and distinct between experimental treatments. Plants grown in soil extract were morphologically and metabolically distinct, with root hairs four times longer than with other growth conditions. Further, plants depleted half of the metabolites investigated from the soil extract. To interact with their environment, plants not only adapt morphology and release complex metabolite mixtures, but also selectively deplete a range of soil-derived metabolites. The EcoFABs utilized here generated high interlaboratory reproducibility, demonstrating their value in standardized investigations of plant traits.

Predicting phenotypic expression from genomic and environmental information is arguably the greatest challenge in today's biology. Being able to survey genomic content, e.g., as single-nucleotide polymorphism data, within a diverse population and predict the phenotypes of external traits, represents the holy grail across genome-informed disciplines, from personal medicine and nutrition to plant breeding. In the present study, we propose a two-step procedure in bridging the genome to phenome gap where external phenotypes are viewed as emergent properties of internal phenotypes, such as molecular profiles, in interaction with the environment. Using biomass accumulation and shoot-root allometry as external traits in diverse genotypes of the model grass Brachypodium distachyon, we established correlative models between genotypes and metabolite profiles (metabotypes) as internal phenotypes, and between metabotypes and external phenotypes under two contrasting watering regimes. Our results demonstrate the potential for employing metabotypes as an integrator in predicting external phenotypes from genomic information.

The development of model systems requires a detailed assessment of standing genetic variation across natural populations. The Brachypodium species complex has been promoted as a plant model for grass genomics with translation to small grain and biomass crops. To capture the genetic diversity within this species complex, thousands of Brachypodium accessions from around the globe were collected and genotyped by sequencing. Overall, 1897 samples were classified into two diploid or allopolyploid species, and then further grouped into distinct inbred genotypes. A core set of diverse B. distachyon diploid lines was selected for whole genome sequencing and high resolution phenotyping. Genome-wide association studies across simulated seasonal environments was used to identify candidate genes and pathways tied to key life history and agronomic traits under current and future climatic conditions. A total of 8, 22, and 47 QTL were identified for flowering time, early vigor, and energy traits, respectively. The results highlight the genomic structure of the Brachypodium species complex, and the diploid lines provided a resource that allows complex trait dissection within this grass model species.

Amplicon sequencing of 16S, ITS, and 18S regions of microbial genomes is a commonly used first step toward understanding microbial communities of interest for human health, agriculture, and the environment. Correlation network analysis is an emerging tool for investigating the interactions within these microbial communities. However, when data from different habitats (e.g., sampling sites, host genotype, etc.) are combined into one analysis, habitat filtering (co-occurrence of microbes due to habitat sampled rather than biological interactions) can induce apparent correlations, resulting in a network dominated by habitat effects and masking correlations of biological interest. We developed an algorithm to correct for habitat filtering effects in microbial correlation network analysis in order to reveal the true underlying microbial correlations. This algorithm was tested on simulated data that was constructed to exhibit habitat filtering. Our algorithm significantly improved correlation detection accuracy for these data compared to Spearman and Pearson correlations. We then used our algorithm to analyze a two real data sets of 16S variable region amplicon sequences that were expected to exhibit habitat filtering. Our algorithm was found to effectively reduce habitat effects, enabling the construction of consensus correlation networks from data sets combining multiple related sample habitats.

BackgroundSorghum bicolor is the fifth most commonly grown cereal worldwide and is remarkable for its drought and abiotic stress tolerance. For these reasons and the large size of biomass varieties, it has been proposed as a bioenergy crop. However, little is known about the genes underlying sorghum's abiotic stress tolerance and biomass yield.ResultsTo uncover the genetic basis of drought tolerance in sorghum at a genome-wide level, we undertook a high-density phenomics genome wide association study (GWAS) in which 648 diverse sorghum lines were phenotyped at two locations in California once per week by drone over the course of a growing season. Biomass, height, and leaf area were measured by drone for individual field plots, subjected to two drought treatments and a well-watered control. The resulting dataset of ~ 171,000 phenotypic data-points was analyzed along with 183,989 genotype by sequence markers to reveal 213 high-quality, replicated, and conserved GWAS associations.ConclusionsThe genomic intervals defined by the associations include many strong candidate genes, including those encoding heat shock proteins, antifreeze proteins, and other domains recognized as important to plant stress responses. The markers identified by our study can be used for marker assisted selection for drought tolerance and biomass. In addition, our results are a significant step toward identifying specific sorghum genes controlling drought tolerance and biomass yield.

There is a dynamic reciprocity between plants and their environment: On one hand, the physiochemical properties of soil influence plant morphology and metabolism, while on the other, root morphology and exudates shape the environment surrounding roots. Here, we investigate both of these aspects as well as the reproducibility of these responses across laboratories. The model grass Brachypodium distachyon was grown in phosphate-sufficient and phosphate-deficient mineral media, as well as in sterile soil extract, within fabricated ecosystem (EcoFAB) devices across four laboratories. Tissue weight and phosphate content, total root length, root tissue and exudate metabolic profiles were found to be consistent across laboratories and distinct between experimental treatments. Plants grown in soil extract were morphologically and metabolically distinct in all laboratories, with root hairs four times longer compared to other growth conditions. Further, plants depleted half of the investigated metabolites from the soil extract. To interact with their environment, plants not only adapt morphology and release complex metabolite mixtures; they also selectively deplete a range of soil-derived metabolites. The EcoFABs utilized here generated high inter-laboratory reproducibility, demonstrating that their value in standardized investigations of plant traits.

Targeted cellular auxin distribution is required for morphogenesis and adaptive responses of plant organs. In Arabidopsis thaliana (Arabidopsis), this involves the prototypical auxin influx facilitator AUX1 and its LIKE-AUX1 (LAX) homologs, which act partially redundantly in various developmental processes. Interestingly, AUX1 and its homologs are not strictly essential for the Arabidopsis life cycle. Indeed, aux1 lax1 lax2 lax3 quadruple knock-outs are mostly viable and fertile, and strong phenotypes are only observed at low penetrance. Here we investigated the Brachypodium distachyon (Brachypodium) AUX1 homolog BdAUX1 by genetic, cell biological and physiological analyses. We report that BdAUX1 is essential for Brachypodium development. Bdaux1 loss-of-function mutants are dwarfs with aberrant flower development, and consequently infertile. Moreover, they display a counter-intuitive root phenotype. Although Bdaux1 roots are agravitropic as expected, in contrast to Arabidopsis aux1 mutants they are dramatically longer than wild type roots because of exaggerated cell elongation. Interestingly, this correlates with higher free auxin content in Bdaux1 roots. Consistently, their cell wall characteristics and transcriptome signature largely phenocopy other Brachypodium mutants with increased root auxin content. Our results imply fundamentally different wiring of auxin transport in Brachypodium roots and reveal an essential role of BdAUX1 in a broad spectrum of developmental processes, suggesting a central role for AUX1 in pooideae.

Beneficial plant-microbe interactions offer a sustainable biological solution with the potential to boost low-input food and bioenergy production. A better mechanistic understanding of these complex plant-microbe interactions will be crucial to improving plant production as well as performing basic ecological studies investigating plant-soil-microbe interactions. Here, a detailed description for ecosystem fabrication is presented, using widely available 3D printing technologies, to create controlled laboratory habitats (EcoFABs) for mechanistic studies of plant-microbe interactions within specific environmental conditions. Two sizes of EcoFABs are described that are suited for the investigation of microbial interactions with various plant species, including Arabidopsis thaliana, Brachypodium distachyon, and Panicum virgatum. These flow-through devices allow for controlled manipulation and sampling of root microbiomes, root chemistry as well as imaging of root morphology and microbial localization. This protocol includes the details for maintaining sterile conditions inside EcoFABs and mounting independent LED light systems onto EcoFABs. Detailed methods for addition of different forms of media, including soils, sand, and liquid growth media coupled to the characterization of these systems using imaging and metabolomics are described. Together, these systems enable dynamic and detailed investigation of plant and plant-microbial consortia including the manipulation of microbiome composition (including mutants), the monitoring of plant growth, root morphology, exudate composition, and microbial localization under controlled environmental conditions. We anticipate that these detailed protocols will serve as an important starting point for other researchers, ideally helping create standardized experimental systems for investigating plant-microbe interactions.

Transposable element (TE) activity has emerged as a major cause of variation in genome size and structure among species. To what extent TEs contribute to genetic variation and divergence within species, however, is much less clear, mainly because population genomic data have so far only been available for the classical model organisms. In this study, we use the annual Mediterranean grass Brachypodium distachyon to investigate TE dynamics in natural populations. Using whole-genome sequencing data for 53 natural accessions, we identified more than 5,400 TE polymorphisms across the studied genomes. We found, first, that while population bottlenecks and expansions have shaped genetic diversity in B. distachyon, these events did not lead to lineage-specific activations of TE families, as observed in other species. Instead, the same families have been active across the species range and TE activity is homogeneous across populations, indicating the presence of conserved regulatory mechanisms. Second, almost half of the TE insertion polymorphisms are accession-specific, most likely because of recent activity in expanding populations and the action of purifying selection. And finally, although TE insertion polymorphisms are underrepresented in and around genes, more than 1,000 of them occur in genic regions and could thus contribute to functional divergence. Our study shows that while TEs in B. distachyon are "well-behaved" compared with TEs in other species with larger genomes, they are an abundant source of lineage-specific genetic variation and may play an important role in population divergence and adaptation.

While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.

Fusarium crown rot (FCR) of wheat and barley, predominantly caused by the fungal pathogen Fusarium pseudograminearum, is a disease of economic significance. The quantitative nature of FCR resistance within cultivated wheat germplasm has significantly limited breeding efforts to enhanced FCR resistance in wheat. In this study, we characterized the molecular responses of Brachypodium distachyon (Brachypodium hereafter) to F. pseudograminearum infection using RNA-seq to determine whether Brachypodium can be exploited as a model system towards better understanding of F. pseudograminearum-wheat interaction. The transcriptional response to infection in Brachypodium was strikingly similar to that previously reported in wheat, both in shared expression patterns of wheat homologs of Brachypodium genes and functional overlap revealed through comparative gene ontology analysis in both species. Metabolites produced by various biosynthetic pathways induced in both wheat and Brachypodium were quantified, revealing a high degree of overlap between these two species in metabolic response to infection but also showed Brachypodium does not produce certain defence-related metabolites found in wheat. Functional analyses of candidate genes identified in this study will improve our understanding of resistance mechanisms and may lead to the development of new strategies to protect cereal crops from pathogen infection.

In Arabidopsis, development during flowering is coordinated by transport of the hormone auxin mediated by polar-localized PIN-FORMED1 (AtPIN1). However Arabidopsis has lost a PIN clade sister to AtPIN1, Sister-of-PIN1 (SoPIN1), which is conserved in flowering plants. We previously proposed that the AtPIN1 organ initiation and vein patterning functions are split between the SoPIN1 and PIN1 clades in grasses. Here we show that in the grass Brachypodium sopin1 mutants have organ initiation defects similar to Arabidopsis atpin1, while loss of PIN1 function in Brachypodium has little effect on organ initiation but alters stem growth. Heterologous expression of Brachypodium SoPIN1 and PIN1b in Arabidopsis provides further evidence of functional specificity. SoPIN1 but not PIN1b can mediate flower formation in null atpin1 mutants, although both can complement a missense allele. The behavior of SoPIN1 and PIN1b in Arabidopsis illustrates how membrane and tissue-level accumulation, transport activity, and interaction contribute to PIN functional specificity.

Due to a large and growing collection of genomic and experimental resources, Brachypodium distachyon has emerged as a powerful experimental model for the grasses. To add to these resources we sequenced 21 165 T-DNA lines, 15 569 of which were produced in this study. This increased the number of unique insertion sites in the T-DNA collection by 21 078, bringing the overall total to 26 112. Thirty-seven per cent (9754) of these insertion sites are within genes (including untranslated regions and introns) and 28% (7217) are within 500 bp of a gene. Approximately 31% of the genes in the v.2.1 annotation have been tagged in this population. To demonstrate the utility of this collection, we phenotypically characterized six T-DNA lines with insertions in genes previously shown in other systems to be involved in cellulose biosynthesis, hemicellulose biosynthesis, secondary cell wall development, DNA damage repair, wax biosynthesis and chloroplast synthesis. In all cases, the phenotypes observed supported previous studies, demonstrating the utility of this collection for plant functional genomics. The Brachypodium T-DNA collection can be accessed at http://jgi.doe.gov/our-science/science-programs/plant-genomics/brachypodium/brachypodium-t-dna-collection/.

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

The transition to reproductive development is a crucial step in the plant life cycle, and the timing of this transition is an important factor in crop yields. Here, we report new insights into the genetic control of natural variation in flowering time in Brachypodium distachyon, a nondomesticated pooid grass closely related to cereals such as wheat (Triticum spp.) and barley (Hordeum vulgare L.). A recombinant inbred line population derived from a cross between the rapid-flowering accession Bd21 and the delayed-flowering accession Bd1-1 were grown in a variety of environmental conditions to enable exploration of the genetic architecture of flowering time. A genotyping-by-sequencing approach was used to develop SNP markers for genetic map construction, and quantitative trait loci (QTLs) that control differences in flowering time were identified. Many of the flowering-time QTLs are detected across a range of photoperiod and vernalization conditions, suggesting that the genetic control of flowering within this population is robust. The two major QTLs identified in undomesticated B. distachyon colocalize with VERNALIZATION1/PHYTOCHROME C and VERNALIZATION2, loci identified as flowering regulators in the domesticated crops wheat and barley. This suggests that variation in flowering time is controlled in part by a set of genes broadly conserved within pooid grasses.

Due to their economic relevance, the study of plant pathogen interactions is of importance. However, elucidating these interactions and their underlying molecular mechanisms remains challenging since both host and pathogen need to be fully genetically accessible organisms. Here we present milestones in the establishment of a new biotrophic model pathosystem: Ustilago bromivora and Brachypodium sp. We provide a complete toolset, including an annotated fungal genome and methods for genetic manipulation of the fungus and its host plant. This toolset will enable researchers to easily study biotrophic interactions at the molecular level on both the pathogen and the host side. Moreover, our research on the fungal life cycle revealed a mating type bias phenomenon. U. bromivora harbors a haplo-lethal allele that is linked to one mating type region. As a result, the identified mating type bias strongly promotes inbreeding, which we consider to be a potential speciation driver.

Stomata, epidermal valves facilitating plant-atmosphere gas exchange, represent a powerful model for understanding cell fate and pattern in plants. Core basic helix-loop-helix (bHLH) transcription factors regulating stomatal development were identified in Arabidopsis, but this dicot's developmental pattern and stomatal morphology represent only one of many possibilities in nature. Here, using unbiased forward genetic screens, followed by analysis of reporters and engineered mutants, we show that stomatal initiation in the grass Brachypodium distachyon uses orthologs of stomatal regulators known from Arabidopsis but that the function and behavior of individual genes, the relationships among genes, and the regulation of their protein products have diverged. Our results highlight ways in which a kernel of conserved genes may be alternatively wired to produce diversity in patterning and morphology and suggest that the stomatal transcription factor module is a prime target for breeding or genome modification to improve plant productivity.

The genetic transformation of monocot grasses is a resource intensive process, the quality and efficiency of which is dependent in part upon the method of DNA introduction, as well as the ability to effectively separate transformed from wildtype tissue. Agrobacterium-mediated transformation of Brachypodium has relied mainly on Agrobacterium tumefaciens strain AGL1. Currently the antibiotic hygromycin B has been the selective agent of choice for robust identification of transgenic calli in Brachypodium distachyon and Brachypodium sylvaticum but few other chemicals have been shown to work as well for selection of transgenic Brachypodium cells in tissue culture. This study demonstrates that Agrobacterium rhizogenes strain 18r12v and paromomycin selection can be successfully used for the efficient generation of transgenic B. distachyon and B. sylvaticum. Additionally we observed that the transformation rates were similar to or higher than those obtained with A. tumefaciens strain AGL1 and hygromycin selection. The A. rhizogenes strain 18r12v harboring the pARS1 binary vector and paromomycin selection is an effective means of generating transgenic Brachypodium plants. This novel approach will facilitate the transgenic complementation of T-DNA knockout mutants of B. distachyon which were created using hygromycin selection, as well as aid the implementation of more complex genome manipulation strategies which require multiple rounds of transformation.

Utility vectors with promoters that confer desired spatial and temporal expression patterns are useful tools for studying gene and cellular function and for industrial applications. To target the expression of DNA sequences of interest to cells forming plant secondary cell walls, which generate most of the vegetative biomass, upstream regulatory sequences of the Brachypodium distachyon lignin biosynthetic gene BdPMT and the cellulose synthase genes BdCESA7 and BdCESA8 were isolated and cloned into binary vectors designed for Agrobacterium-mediated transformation of monocots. Expression patterns were assessed using the β-glucuronidase gene GUSPlus and X-glucuronide staining. All three promoters showed strong expression levels in stem tissue at the base of internodes where cell wall deposition is most active, in both vascular bundle xylem vessels and tracheids, and in interfascicular tissues, with expression less pronounced in developmentally older tissues. In leaves, BdCESA7 and BdCESA8 promoter-driven expression was strongest in leaf veins, leaf margins, and trichomes; relatively weaker and patchy expression was observed in the epidermis. BdPMT promoter-driven expression was similar to the BdCESA promoters expression patterns, including strong expression in trichomes. The intensity and extent of GUS staining varied considerably between transgenic lines, suggesting that positional effects influenced promoter activity. Introducing the BdPMT and BdCESA8 Open Reading Frames into BdPMT and BdCESA8 utility promoter binary vectors, respectively, and transforming those constructs into Brachypodium pmt and cesa8 loss-of-function mutants resulted in rescue of the corresponding mutant phenotypes. This work therefore validates the functionality of these utility promoter binary vectors for use in Brachypodium and likely other grass species. The identification, in Bdcesa8-1 T-DNA mutant stems, of an 80% reduction in crystalline cellulose levels confirms that the BdCESA8 gene is a secondary-cell-wall-forming cellulose synthase.

Brachypodium hybridum (2n = 30) is a natural allopolyploid with highly divergent sub-genomes derived from two extant diploid species, B. distachyon (2n = 10) and B. stacei (2n = 20) that differ in chromosome evolution and number. We created synthetic B. hybridum allotetraploids by hybridizing various lines of B. distachyon and B. stacei. The initial amphihaploid F1 interspecific hybrids were obtained at low frequencies when B. distachyon was used as the maternal parent (0.15% or 0.245% depending on the line used) and were sterile. No hybrids were obtained from reciprocal crosses or when autotetraploids of the parental species were crossed. Colchicine treatment was used to double the genome of the F1 amphihaploid lines leading to allotetraploids. The genome-doubled F1 plants produced a few S1 (first selfed generation) seeds after self-pollination. S1 plants from one parental combination (Bd3-1×Bsta5) were fertile and gave rise to further generations whereas those of another parental combination (Bd21×ABR114) were sterile, illustrating the importance of the parental lineages crossed. The synthetic allotetraploids were stable and resembled the natural B. hybridum at the phenotypic, cytogenetic and genomic levels. The successful creation of synthetic B. hybridum offers the possibility to study changes in genome structure and regulation at the earliest stages of allopolyploid formation in comparison with the parental species and natural B. hybridum.

Brachypodium hybridum (2n = 30) is a natural allopolyploid with highly divergent sub-genomes derived from two extant diploid species, B. distachyon (2n = 10) and B. stacei (2n = 20) that differ in chromosome evolution and number. We created synthetic B. hybridum allotetraploids by hybridizing various lines of B. distachyon and B. stacei. The initial amphihaploid F1 interspecific hybrids were obtained at low frequencies when B. distachyon was used as the maternal parent (0.15% or 0.245% depending on the line used) and were sterile. No hybrids were obtained from reciprocal crosses or when autotetraploids of the parental species were crossed. Colchicine treatment was used to double the genome of the F1 amphihaploid lines leading to allotetraploids. The genome-doubled F1 plants produced a few S1 (first selfed generation) seeds after self-pollination. S1 plants from one parental combination (Bd3-1×Bsta5) were fertile and gave rise to further generations whereas those of another parental combination (Bd21×ABR114) were sterile, illustrating the importance of the parental lineages crossed. The synthetic allotetraploids were stable and resembled the natural B. hybridum at the phenotypic, cytogenetic and genomic levels. The successful creation of synthetic B. hybridum offers the possibility to study changes in genome structure and regulation at the earliest stages of allopolyploid formation in comparison with the parental species and natural B. hybridum.

We are using the perennial model grass Brachypodium sylvaticum to identify combinations of transgenes that enhance tolerance to multiple, simultaneous abiotic stresses. The most successful transgene combinations will ultimately be used to create improved switchgrass (Panicum virgatum L.) cultivars. To further develop B. sylvaticum as a perennial model grass, and facilitate our planned transcriptional profiling, we are sequencing and annotating the genome. We have generated ~40x genome coverage using PacBio sequencing of the largest possible size selected libraries (18, 22, 25 kb). Our initial assembly using only long-read sequence contained 320 Mb of sequence with an N50 contig length of 315 kb and an N95 contig length of 40 kb. This assembly consists of 2,430 contigs, the largest of which was 1.6 Mb. The estimated genome size based on c-values is 340 Mb indicating that about 20 Mb of presumably repetitive DNA remains yet unassembled. Significantly, this assembly is far superior to an assembly created from paired-end short-read sequence, ~100x genome coverage. The short-read-only assembly contained only 226 Mb of sequence in 19k contigs. To aid the assembly of the scaffolds into chromosome-scale assemblies we produced an F2 mapping population and have genotyped 480 individuals using a genotype by sequence approach. One of the reasons for using B. sylvaticum as a model system is to determine if the transgenes adversely affect perenniality and winter hardiness. Toward this goal, we examined the freezing tolerance of wild type B. sylvaticum lines to determine the optimal conditions for testing the freezing tolerance of the transgenics. A survey of seven accessions noted significant natural variation in freezing tolerance. Seedling or adult Ain-1 plants, the line used for transformation, survived an 8 hour challenge down to -6 oC and 50% survived a challenge down to -9 oC. Thus, we will be able to easily determine if the transgenes compromise freezing tolerance. In the effort to develop biotechnological tools for perennial grass improvement, we have completed the transformation of B. sylvaticum with constructs containing 20 genes shown to be associated with enhanced abiotic stress tolerance in monocots. In addition, we have transformed plants with constructs containing a combination of genes (i.e. SARK::IPT- Ubi::HSR1::Ubi::NHX1) in order to simultaneously overexpress genes associated with drought + heat tolerance + salt tolerance. We generated single copy insert T1 lines for all constructs and the generation and bulking of homozygous T2 lines is well underway. In addition to our B. sylvaticum transgenics, we transformed B. distachyon with many of the same genes. Some of the transgenic B. distachyon plants subjected to a combined stress of both drought and salinity were able to produced higher yields than wild type plants. Our results indicate a great potential for the development of grasses with improved performance and yield in water-limited areas.

Root systems develop different root types that individually sense cues from their local environment and integrate this information with systemic signals. This complex multi-dimensional amalgam of inputs enables continuous adjustment of root growth rates, direction, and metabolic activity that define a dynamic physical network. Current methods for analyzing root biology balance physiological relevance with imaging capability. To bridge this divide, we developed an integrated-imaging system called Growth and Luminescence Observatory for Roots (GLO-Roots) that uses luminescence-based reporters to enable studies of root architecture and gene expression patterns in soil-grown, light-shielded roots. We have developed image analysis algorithms that allow the spatial integration of soil properties, gene expression, and root system architecture traits. We propose GLO-Roots as a system that has great utility in presenting environmental stimuli to roots in ways that evoke natural adaptive responses and in providing tools for studying the multi-dimensional nature of such processes.