Streptomyces umbrella toxin particles block hyphal growth of competing species

Streptomyces are a genus of ubiquitous soil bacteria from which the majority of clinically utilized antibiotics derive1. The production of these antibacterial molecules reflects the relentless competition Streptomyces engage in with other bacteria, including other Streptomyces species1,2. Here we show that in addition to small-molecule antibiotics, Streptomyces produce and secrete antibacterial protein complexes that feature a large, degenerate repeat-containing polymorphic toxin protein. A cryo-electron microscopy structure of these particles reveals an extended stalk topped by a ringed crown comprising the toxin repeats scaffolding five lectin-tipped spokes, which led us to name them umbrella particles. Streptomyces coelicolor encodes three umbrella particles with distinct toxin and lectin composition. Notably, supernatant containing these toxins specifically and potently inhibits the growth of select Streptomyces species from among a diverse collection of bacteria screened. For one target, Streptomyces griseus, inhibition relies on a single toxin and that intoxication manifests as rapid cessation of vegetative hyphal growth. Our data show that Streptomyces umbrella particles mediate competition among vegetative mycelia of related species, a function distinct from small-molecule antibiotics, which are produced at the onset of reproductive growth and act broadly3,4. Sequence analyses suggest that this role of umbrella particles extends beyond Streptomyces, as we identified umbrella loci in nearly 1,000 species across Actinobacteria.

Phylogenetic distribution and experimental characterization of corrinoid production and dependence in soil bacterial isolates

Soil microbial communities impact carbon sequestration and release, biogeochemical cycling, and agricultural yields. These global effects rely on metabolic interactions that modulate community composition and function. However, the physicochemical and taxonomic complexity of soil and the scarcity of available isolates for phenotypic testing are significant barriers to studying soil microbial interactions. Corrinoids-the vitamin B12 family of cofactors-are critical for microbial metabolism, yet they are synthesized by only a subset of microbiome members. Here, we evaluated corrinoid production and dependence in soil bacteria as a model to investigate the ecological roles of microorganisms involved in metabolic interactions. We isolated and characterized a taxonomically diverse collection of 161 soil bacteria from a single study site. Most corrinoid-dependent bacteria in the collection prefer B12 over other corrinoids, while all tested producers synthesize B12, indicating metabolic compatibility between producers and dependents in the collection. Furthermore, a subset of producers release B12 at levels sufficient to support dependent isolates in laboratory culture at estimated ratios of up to 1000 dependents per producer. Within our isolate collection, we did not find strong phylogenetic patterns in corrinoid production or dependence. Upon investigating trends in the phylogenetic dispersion of corrinoid metabolism categories across sequenced bacteria from various environments, we found that these traits are conserved in 47 out of 85 genera. Together, these phenotypic and genomic results provide evidence for corrinoid-based metabolic interactions among bacteria and provide a framework for the study of nutrient-sharing ecological interactions in microbial communities.

High-throughput genetic engineering of nonmodel and undomesticated bacteria via iterative site-specific genome integration

Efficient genome engineering is critical to understand and use microbial functions. Despite recent development of tools such as CRISPR-Cas gene editing, efficient integration of exogenous DNA with well-characterized functions remains limited to model bacteria. Here, we describe serine recombinase-assisted genome engineering, or SAGE, an easy-to-use, highly efficient, and extensible technology that enables selection marker-free, site-specific genome integration of up to 10 DNA constructs, often with efficiency on par with or superior to replicating plasmids. SAGE uses no replicating plasmids and thus lacks the host range limitations of other genome engineering technologies. We demonstrate the value of SAGE by characterizing genome integration efficiency in five bacteria that span multiple taxonomy groups and biotechnology applications and by identifying more than 95 heterologous promoters in each host with consistent transcription across environmental and genetic contexts. We anticipate that SAGE will rapidly expand the number of industrial and environmental bacteria compatible with high-throughput genetics and synthetic biology.

Photosynthesis in rice is increased by CRISPR/Cas9-mediated transformation of two truncated light-harvesting antenna

Plants compete for light partly by over-producing chlorophyll in leaves. The resulting high light absorption is an effective strategy for out competing neighbors in mixed communities, but it prevents light transmission to lower leaves and limits photosynthesis in dense agricultural canopies. We used a CRISPR/Cas9-mediated approach to engineer rice plants with truncated light-harvesting antenna (TLA) via knockout mutations to individual antenna assembly component genes CpSRP43, CpSRP54a, and its paralog, CpSRP54b. We compared the photosynthetic contributions of these components in rice by studying the growth rates of whole plants, quantum yield of photosynthesis, chlorophyll density and distribution, and phenotypic abnormalities. Additionally, we investigated a Poales-specific duplication of CpSRP54. The Poales are an important family that includes staple crops such as rice, wheat, corn, millet, and sorghum. Mutations in any of these three genes involved in antenna assembly decreased chlorophyll content and light absorption and increased photosynthesis per photon absorbed (quantum yield). These results have significant implications for the improvement of high leaf-area-index crop monocultures.

Plant immunity suppression via PHR1-RALF-FERONIA shapes the root microbiome to alleviate phosphate starvation

The microbiome plays an important role in shaping plant growth and immunity, but few plant genes and pathways impacting plant microbiome composition have been reported. In Arabidopsis thaliana, the phosphate starvation response (PSR) was recently found to modulate the root microbiome upon phosphate (Pi) starvation through the transcriptional regulator PHR1. Here, we report that A. thaliana PHR1 directly binds to the promoters of rapid alkalinization factor (RALF) genes, and activates their expression under phosphate-starvation conditions. RALFs in turn suppress complex formation of pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) receptor through FERONIA, a previously-identified PTI modulator that increases resistance to certain detrimental microorganisms. Suppression of immunity via the PHR1-RALF-FERONIA axis allows colonization by specialized root microbiota that help to alleviate phosphate starvation by upregulating the expression of PSR genes. These findings provide a new paradigm for coordination of host-microbe homeostasis through modulating plant innate immunity after environmental perturbations.

Successional adaptive strategies revealed by correlating arbuscular mycorrhizal fungal abundance with host plant gene expression

The shifts in adaptive strategies revealed by ecological succession and the mechanisms that facilitate these shifts are fundamental to ecology. These adaptive strategies could be particularly important in communities of arbuscular mycorrhizal fungi (AMF) mutualistic with sorghum where strong AMF succession replaces initially ruderal species with competitive ones and where the strongest plant response to drought is to manage these AMF. Although most studies of agriculturally important fungi focus on parasites, the mutualistic symbionts, AMF, constitute a research system of human-associated fungi whose relative simplicity and synchrony are conducive to experimental ecology. First, we hypothesize that, when irrigation is stopped to mimic drought, competitive AMF species should be replaced by AMF species tolerant to drought stress. We then, for the first time, correlate AMF abundance and host plant transcription to test two novel hypotheses about the mechanisms behind the shift from ruderal to competitive AMF. Surprisingly, despite imposing drought stress, we found no stress tolerant AMF, likely due to our agricultural system having been irrigated for nearly six decades. Remarkably, we found strong and differential correlation between the successional shift from ruderal to competitive AMF and sorghum genes whose products (i) produce and release strigolactone signals, (ii) perceive mycorrhizal-lipochitinoligosaccharide (Myc-LCO) signals, (iii) provide plant lipid and sugar to AMF and, (iv) import minerals and water provided by AMF. These novel insights frame new hypotheses about AMF adaptive evolution and suggest a rationale for selecting AMF to reduce inputs and maximize yields in commercial agriculture.

Plant Microbiome-Based Genome-Wide Association Studies

Plants form intimate associations with microorganisms, and these associations are directly impacted by the host genotype. However, identifying specific host genetic pathways that influence these microbial interactions has proved challenging. Genome-wide association-based approaches that use features of microbiome composition as a quantitative trait represent a novel and underutilized strategy to identify such pathways. Several recent studies have demonstrated the potential utility of plant microbiome-based genome-wide association studies (GWAS). In this chapter, we describe the process of implementing GWAS using the plant microbiome as the primary quantitative trait, considering experimental design, sample harvest, and processing, but with an emphasis on data filtering, data normalization, and statistical analyses.

Cell Wall Compositions of Sorghum bicolor Leaves and Roots Remain Relatively Constant Under Drought Conditions

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.

Genome wide association study reveals plant loci controlling heritability of the rhizosphere microbiome

Host genetics has recently been shown to be a driver of plant microbiome composition. However, identifying the underlying genetic loci controlling microbial selection remains challenging. Genome-wide association studies (GWAS) represent a potentially powerful, unbiased method to identify microbes sensitive to the host genotype and to connect them with the genetic loci that influence their colonization. Here, we conducted a population-level microbiome analysis of the rhizospheres of 200 sorghum genotypes. Using 16S rRNA amplicon sequencing, we identify rhizosphere-associated bacteria exhibiting heritable associations with plant genotype, and identify significant overlap between these lineages and heritable taxa recently identified in maize. Furthermore, we demonstrate that GWAS can identify host loci that correlate with the abundance of specific subsets of the rhizosphere microbiome. Finally, we demonstrate that these results can be used to predict rhizosphere microbiome structure for an independent panel of sorghum genotypes based solely on knowledge of host genotypic information.

Keep your friends close: Host compartmentalisation of microbial communities facilitates decoupling from effects of habitat fragmentation

Root-associated fungal communities modify the climatic niches and even the competitive ability of their hosts, yet how the different components of the root microbiome are modified by habitat loss remains a key knowledge gap. Using principles of landscape ecology, we tested how free-living versus host-associated microbes differ in their response to landscape heterogeneity. Further, we explore how compartmentalisation of microbes into specialised root structures filters for key fungal symbionts. Our study demonstrates that free-living fungal community structure correlates with landscape heterogeneity, but that host-associated fungal communities depart from these patterns. Specifically, biotic filtering in roots, especially via compartmentalisation within specialised root structures, decouples the biogeographic patterns of host-associated fungal communities from the soil community. In this way, even as habitat loss and fragmentation threaten fungal diversity in the soils, plant hosts exert biotic controls to ensure associations with critical mutualists, helping to preserve the root mycobiome.

Plant single-cell solutions for energy and the environment

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.

The Plastid-Localized AtFtsHi3 Pseudo-Protease of Arabidopsis thaliana Has an Impact on Plant Growth and Drought Tolerance

While drought severely affects plant growth and crop production, the molecular mechanisms of the drought response of plants remain unclear. In this study, we demonstrated for the first time the effect of the pseudo-protease AtFtsHi3 of Arabidopsis thaliana on overall plant growth and in drought tolerance. An AtFTSHi3 knock-down mutant [ftshi3-1(kd)] displayed a pale-green phenotype with lower photosynthetic efficiency and Darwinian fitness compared to wild type (Wt). An observed delay in seed germination of ftshi3-1(kd) was attributed to overaccumulation of abscisic acid (ABA); ftshi3-1(kd) seedlings showed partial sensitivity to exogenous ABA. Being exposed to similar severity of soil drying, ftshi3-1(kd) was drought-tolerant up to 20 days after the last irrigation, while wild type plants wilted after 12 days. Leaves of ftshi3-1(kd) contained reduced stomata size, density, and a smaller stomatic aperture. During drought stress, ftshi3-1(kd) showed lowered stomatal conductance, increased intrinsic water-use efficiency (WUEi), and slower stress acclimation. Expression levels of ABA-responsive genes were higher in leaves of ftshi3-1(kd) than Wt; DREB1A, but not DREB2A, was significantly upregulated during drought. However, although ftshi3-1(kd) displayed a drought-tolerant phenotype in aboveground tissue, the root-associated bacterial community responded to drought.

Genome-resolved metagenomics reveals role of iron metabolism in drought-induced rhizosphere microbiome dynamics

Recent studies have demonstrated that drought leads to dramatic, highly conserved shifts in the root microbiome. At present, the molecular mechanisms underlying these responses remain largely uncharacterized. Here we employ genome-resolved metagenomics and comparative genomics to demonstrate that carbohydrate and secondary metabolite transport functionalities are overrepresented within drought-enriched taxa. These data also reveal that bacterial iron transport and metabolism functionality is highly correlated with drought enrichment. Using time-series root RNA-Seq data, we demonstrate that iron homeostasis within the root is impacted by drought stress, and that loss of a plant phytosiderophore iron transporter impacts microbial community composition, leading to significant increases in the drought-enriched lineage, Actinobacteria. Finally, we show that exogenous application of iron disrupts the drought-induced enrichment of Actinobacteria, as well as their improvement in host phenotype during drought stress. Collectively, our findings implicate iron metabolism in the root microbiome's response to drought and may inform efforts to improve plant drought tolerance to increase food security.

Holo-omics for deciphering plant-microbiome interactions

Host-microbiome interactions are recognized for their importance to host health. An improved understanding of the molecular underpinnings of host-microbiome relationships will advance our capacity to accurately predict host fitness and manipulate interaction outcomes. Within the plant microbiome research field, unlocking the functional relationships between plants and their microbial partners is the next step to effectively using the microbiome to improve plant fitness. We propose that strategies that pair host and microbial datasets-referred to here as holo-omics-provide a powerful approach for hypothesis development and advancement in this area. We discuss several experimental design considerations and present a case study to highlight the potential for holo-omics to generate a more holistic perspective of molecular networks within the plant microbiome system. In addition, we discuss the biggest challenges for conducting holo-omics studies; specifically, the lack of vetted analytical frameworks, publicly available tools, and required technical expertise to process and integrate heterogeneous data. Finally, we conclude with a perspective on appropriate use-cases for holo-omics studies, the need for downstream validation, and new experimental techniques that hold promise for the plant microbiome research field. We argue that utilizing a holo-omics approach to characterize host-microbiome interactions can provide important opportunities for broadening system-level understandings and significantly inform microbial approaches to improving host health and fitness. Video abstract.

Evaluating domestication and ploidy effects on the assembly of the wheat bacterial microbiome

While numerous studies implicate the microbiome in host fitness, contributions of host evolution to microbial recruitment remain largely uncharacterized. Past work has shown that plant polyploidy and domestication can influence plant biotic and abiotic interactions, yet impacts on broader microbiome assembly are still unknown for many crop species. In this study, we utilized three approaches-two field studies and one greenhouse-based experiment-to determine the degree to which patterns in bacterial community assembly in wheat (Triticum sp.) roots and rhizospheres are attributable to the host factors of ploidy level (2n, 4n, 6n) and domestication status (cultivated vs. wild). Profiling belowground bacterial communities with 16S rRNA gene amplicon sequencing, we analyzed patterns in diversity and composition. From our initial analyses of a subsetted dataset, we observed that host ploidy level was statistically significant in explaining variation in alpha and beta diversity for rhizosphere microbiomes, as well as correlated with distinct phylum-level shifts in composition, in the field. Using a reduced complexity field soil inoculum and controlled greenhouse conditions, we found some evidence suggesting that genomic lineage and ploidy level influence root alpha and beta diversity (p-value<0.05). However, in a follow-up field experiment using an expanded set of Triticum genomes that included both wild and domesticated varieties, we did not find a strong signal for either diploid genome lineages, domestication status, or ploidy level in shaping rhizosphere bacterial communities. Taken together, these results suggest that while host ploidy and domestication may have some minor influence on microbial assembly, these impacts are subtle and difficult to assess in belowground compartments for wheat varieties. By improving our understanding of the degree to which host ploidy and cultivation factors shape the plant microbiome, this research informs perspectives on what key driving forces may underlie microbiome structuring, as well as where future efforts may be best directed towards fortifying plant growth by microbial means. The greatest influence of the host on the wheat microbiome appeared to occur in the rhizosphere compartment, and we suggest that future work focuses on this environment to further characterize how host genomic and phenotypic changes influence plant-microbe communications.

The generalizability of water-deficit on bacterial community composition; Site-specific water-availability predicts the bacterial community associated with coast redwood roots

Experimental drought has been shown to delay the development of the root microbiome and increase the relative abundance of Actinobacteria, however, the generalizability of these findings to natural systems or other diverse plant hosts remains unknown. Bacterial cell wall thickness and growth morphology (e.g., filamentous or unicellular) have been proposed as traits that may mediate bacterial responses to environmental drivers. Leveraging a natural gradient of water-availability across the coast redwood (Sequoia sempervirens) range, we tested three hypotheses: (a) that site-specific water-availability is an important predictor of bacterial community composition for redwood roots and rhizosphere soils; (b) that there is relative enrichment of Actinobacteria and other monoderm bacterial groups within the redwood microbiome in response to drier conditions; and (c) that bacterial growth morphology is an important predictor of bacteria response to water-availability, where filamentous taxa will become more dominant at drier sites compared to unicellular bacteria. We find that both α- and β-diversity of redwood bacterial communities is partially explained by water-availability and that Actinobacterial enrichment is a conserved response of land plants to water-deficit. Further, we highlight how the trend of Actinobacterial enrichment in the redwood system is largely driven by the Actinomycetales. We propose bacterial growth morphology (filamentous vs. unicellular) as an additional mechanism behind the increase in Actinomycetales with increasing aridity. A trait-based approach including cell-wall thickness and growth morphology may explain the distribution of bacterial taxa across environmental gradients and help to predict patterns of bacterial community composition for a wide range of host plants.

Nematode-Encoded RALF Peptide Mimics Facilitate Parasitism of Plants through the FERONIA Receptor Kinase

The molecular mechanism by which plants defend against plant root-knot nematodes (RKNs) is largely unknown. The plant receptor kinase FERONIA and its peptide ligands, rapid alkalinization factors (RALFs), regulate plant immune responses and cell expansion, which are two important factors for successful RKN parasitism. In this study, we found that mutation of FERONIA in Arabidopsis thaliana resulted in plants showing low susceptibility to the RKN Meloidogyne incognita. To identify the underlying mechanisms associated with this phenomenon, we identified 18 novel RALF-likes from multiple species of RKNs and showed that two RALF-likes (i.e., MiRALF1 and MiRALF3) from M. incognita were expressed in the esophageal gland with high expression during the parasitic stages of nematode development. These nematode RALF-likes also possess the typical activities of plant RALFs and can directly bind to the extracellular domain of FERONIA to modulate specific steps of nematode parasitism-related immune responses and cell expansion. Genetically, both MiRALF1/3 and FERONIA are required for RKN parasitism in Arabidopsis and rice. Collectively, our study suggests that nematode-encoded RALFs facilitate parasitism via plant-encoded FERONIA and provides a novel paradigm for studying host-pathogen interactions.

Drought Drives Spatial Variation in the Millet Root Microbiome

Efforts to boost crop yield and meet global food demands while striving to reach sustainability goals are hindered by the increasingly severe impacts of abiotic stress, such as drought. One strategy for alleviating drought stress in crops is to utilize root-associated bacteria, yet knowledge concerning the relationship between plant hosts and their microbiomes during drought remain under-studied. One broad pattern that has recently been reported in a variety of monocot and dicot species from both native and agricultural environments, is the enrichment of Actinobacteria within the drought-stressed root microbiome. In order to better understand the causes of this phenomenon, we performed a series of experiments in millet plants to explore the roles of drought severity, drought localization, and root development in provoking Actinobacteria enrichment within the root endosphere. Through 16S rRNA amplicon-based sequencing, we demonstrate that the degree of drought is correlated with levels of Actinobacterial enrichment in four species of millet. Additionally, we demonstrate that the observed drought-induced enrichment of Actinobacteria occurs along the length of the root, but the response is localized to portions of the root experiencing drought. Finally, we demonstrate that Actinobacteria are depleted in the dead root tissue of Japanese millet, suggesting saprophytic activity is not the main cause of observed shifts in drought-treated root microbiome structure. Collectively, these results help narrow the list of potential causes of drought-induced Actinobacterial enrichment in plant roots by showing that enrichment is dependent upon localized drought responses but not root developmental stage or root death.

A Plant Growth-Promoting Microbial Soil Amendment Dynamically Alters the Strawberry Root Bacterial Microbiome

Despite growing interest in utilizing microbial-based methods for improving crop growth, much work still remains in elucidating how beneficial plant-microbe associations are established, and what role soil amendments play in shaping these interactions. Here, we describe a set of experiments that test the effect of a commercially available soil amendment, VESTA, on the soil and strawberry (Fragaria x ananassa Monterey) root bacterial microbiome. The bacterial communities of the soil, rhizosphere, and root from amendment-treated and untreated fields were profiled at four time points across the strawberry growing season using 16S rRNA gene amplicon sequencing on the Illumina MiSeq platform. In all sample types, bacterial community composition and relative abundance were significantly altered with amendment application. Importantly, time point effects on composition are more pronounced in the root and rhizosphere, suggesting an interaction between plant development and treatment effect. Surprisingly, there was slight overlap between the taxa within the amendment and those enriched in plant and soil following treatment, suggesting that VESTA may act to rewire existing networks of organisms through an, as of yet, uncharacterized mechanism. These findings demonstrate that a commercial microbial soil amendment can impact the bacterial community structure of both roots and the surrounding environment.

OrthoVenn2: a web server for whole-genome comparison and annotation of orthologous clusters across multiple species

OrthoVenn is a powerful web platform for the comparison and analysis of whole-genome orthologous clusters. Here we present an updated version, OrthoVenn2, which provides new features that facilitate the comparative analysis of orthologous clusters among up to 12 species. Additionally, this update offers improvements to data visualization and interpretation, including an occurrence pattern table for interrogating the overlap of each orthologous group for the queried species. Within the occurrence table, the functional annotations and summaries of the disjunctions and intersections of clusters between the chosen species can be displayed through an interactive Venn diagram. To facilitate a broader range of comparisons, a larger number of species, including vertebrates, metazoa, protists, fungi, plants and bacteria, have been added in OrthoVenn2. Finally, a stand-alone version is available to perform large dataset comparisons and to visualize results locally without limitation of species number. In summary, OrthoVenn2 is an efficient and user-friendly web server freely accessible at https://orthovenn2.bioinfotoolkits.net.

Strong succession in arbuscular mycorrhizal fungal communities

The ecology of fungi lags behind that of plants and animals because most fungi are microscopic and hidden in their substrates. Here, we address the basic ecological process of fungal succession in nature using the microscopic, arbuscular mycorrhizal fungi (AMF) that form essential mutualisms with 70-90% of plants. We find a signal for temporal change in AMF community similarity that is 40-fold stronger than seen in the most recent studies, likely due to weekly samplings of roots, rhizosphere and soil throughout the 17 weeks from seedling to fruit maturity and the use of the fungal DNA barcode to recognize species in a simple, agricultural environment. We demonstrate the patterns of nestedness and turnover and the microbial equivalents of the processes of immigration and extinction, that is, appearance and disappearance. We also provide the first evidence that AMF species co-exist rather than simply co-occur by demonstrating negative, density-dependent population growth for multiple species. Our study shows the advantages of using fungi to test basic ecological hypotheses (e.g., nestedness v. turnover, immigration v. extinction, and coexistence theory) over periods as short as one season.

Exploring the Root Microbiome: Extracting Bacterial Community Data from the Soil, Rhizosphere, and Root Endosphere

The intimate interaction between plant host and associated microorganisms is crucial in determining plant fitness, and can foster improved tolerance to abiotic stresses and diseases. As the plant microbiome can be highly complex, low-cost, high-throughput methods such as amplicon-based sequencing of the 16S rRNA gene are often preferred for characterizing its microbial composition and diversity. However, the selection of appropriate methodology when conducting such experiments is critical for reducing biases that can make analysis and comparisons between samples and studies difficult. This protocol describes in detail a standardized methodology for the collection and extraction of DNA from soil, rhizosphere, and root samples. Additionally, we highlight a well-established 16S rRNA amplicon sequencing pipeline that allows for the exploration of the composition of bacterial communities in these samples, and can easily be adapted for other marker genes. This pipeline has been validated for a variety of plant species, including sorghum, maize, wheat, strawberry, and agave, and can help overcome issues associated with the contamination from plant organelles.

Drought and host selection influence bacterial community dynamics in the grass root microbiome

Root endophytes have been shown to have important roles in determining host fitness under periods of drought stress, and yet the effect of drought on the broader root endosphere bacterial community remains largely uncharacterized. In this study, we present phylogenetic profiles of bacterial communities associated with drought-treated root and rhizosphere tissues of 18 species of plants with varying degrees of drought tolerance belonging to the Poaceae family, including important crop plants. Through 16S rRNA gene profiling across two distinct watering regimes and two developmental time points, we demonstrate that there is a strong correlation between host phylogenetic distance and the microbiome dissimilarity within root tissues, and that drought weakens this correlation by inducing conserved shifts in bacterial community composition. We identify a significant enrichment in a wide variety of Actinobacteria during drought within the roots of all hosts, and demonstrate that this enrichment is higher within the root than it is in the surrounding environments. Furthermore, we show that this observed enrichment is the result of an absolute increase in Actinobacterial abundance and that previously hypothesized mechanisms for observed enrichments in Actinobacteria in drought-treated soils are unlikely to fully account for the phenomena observed here within the plant root.

Drought Stress and Root-Associated Bacterial Communities

Root-associated bacterial communities play a vital role in maintaining health of the plant host. These communities exist in complex relationships, where composition and abundance of community members is dependent on a number of factors such as local soil chemistry, plant genotype and phenotype, and perturbations in the surrounding abiotic environment. One common perturbation, drought, has been shown to have drastic effects on bacterial communities, yet little is understood about the underlying causes behind observed shifts in microbial abundance. As drought may affect root bacterial communities both directly by modulating moisture availability, as well as indirectly by altering soil chemistry and plant phenotypes, we provide a synthesis of observed trends in recent studies and discuss possible directions for future research that we hope will provide for more knowledgeable predictions about community responses to future drought events.

PIECE 2.0: an update for the plant gene structure comparison and evolution database

PIECE (Plant Intron Exon Comparison and Evolution) is a web-accessible database that houses intron and exon information of plant genes. PIECE serves as a resource for biologists interested in comparing intron–exon organization and provides valuable insights into the evolution of gene structure in plant genomes. Recently, we updated PIECE to a new version, PIECE 2.0 (http://probes.pw.usda.gov/piece or http://aegilops.wheat.ucdavis.edu/piece). PIECE 2.0 contains annotated genes from 49 sequenced plant species as compared to 25 species in the previous version. In the current version, we also added several new features: (i) a new viewer was developed to show phylogenetic trees displayed along with the structure of individual genes; (ii) genes in the phylogenetic tree can now be also grouped according to KOG (The annotation of Eukaryotic Orthologous Groups) and KO (KEGG Orthology) in addition to Pfam domains; (iii) information on intronless genes are now included in the database; (iv) a statistical summary of global gene structure information for each species and its comparison with other species was added; and (v) an improved GSDraw tool was implemented in the web server to enhance the analysis and display of gene structure. The updated PIECE 2.0 database will be a valuable resource for the plant research community for the study of gene structure and evolution.

GSP: a web-based platform for designing genome-specific primers in polyploids

MOTIVATION

The sequences among subgenomes in a polyploid species have high similarity, making it difficult to design genome-specific primers for sequence analysis.

RESULTS

We present GSP, a web-based platform to design genome-specific primers that distinguish subgenome sequences in a polyploid genome. GSP uses BLAST to extract homeologous sequences of the subgenomes in existing databases, performs a multiple sequence alignment, and design primers based on sequence variants in the alignment. An interactive primers diagram, a sequence alignment viewer and a virtual electrophoresis are displayed as parts of the primer design result. GSP also designs specific primers from multiple sequences uploaded by users.

AVAILABILITY AND IMPLEMENTATION

GSP is a user-friendly and efficient web platform freely accessible at http://probes.pw.usda.gov/GSP Source code and command-line application are available at https://github.com/bioinfogenome/GSP CONTACTS: yong.gu@ars.usda.gov or devin.coleman-derr@ars.usda.gov

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Cacti Microbiome: Interplay between Habitat-Filtering and Host-Specificity

Cactaceae represents one of the most species-rich families of succulent plants native to arid and semi-arid ecosystems, yet the associations Cacti establish with microorganisms and the rules governing microbial community assembly remain poorly understood. We analyzed the composition, diversity, and factors influencing above- and below-ground bacterial, archaeal, and fungal communities associated with two native and sympatric Cacti species: Myrtillocactus geometrizans and Opuntia robusta. Phylogenetic profiling showed that the composition and assembly of microbial communities associated with Cacti were primarily influenced by the plant compartment; plant species, site, and season played only a minor role. Remarkably, bacterial, and archaeal diversity was higher in the phyllosphere than in the rhizosphere of Cacti, while the opposite was true for fungi. Semi-arid soils exhibited the highest levels of microbial diversity whereas the stem endosphere the lowest. Despite their taxonomic distance, M. geometrizans and O. robusta shared most microbial taxa in all analyzed compartments. Influence of the plant host did only play a larger role in the fungal communities of the stem endosphere. These results suggest that fungi establish specific interactions with their host plant inside the stem, whereas microbial communities in the other plant compartments may play similar functional roles in these two species. Biochemical and molecular characterization of seed-borne bacteria of Cacti supports the idea that these microbial symbionts may be vertically inherited and could promote plant growth and drought tolerance for the fitness of the Cacti holobiont. We envision this knowledge will help improve and sustain agriculture in arid and semi-arid regions of the world.

High-resolution phylogenetic microbial community profiling

Over the past decade, high-throughput short-read 16S rRNA gene amplicon sequencing has eclipsed clone-dependent long-read Sanger sequencing for microbial community profiling. The transition to new technologies has provided more quantitative information at the expense of taxonomic resolution with implications for inferring metabolic traits in various ecosystems. We applied single-molecule real-time sequencing for microbial community profiling, generating full-length 16S rRNA gene sequences at high throughput, which we propose to name PhyloTags. We benchmarked and validated this approach using a defined microbial community. When further applied to samples from the water column of meromictic Sakinaw Lake, we show that while community structures at the phylum level are comparable between PhyloTags and Illumina V4 16S rRNA gene sequences (iTags), variance increases with community complexity at greater water depths. PhyloTags moreover allowed less ambiguous classification. Last, a platform-independent comparison of PhyloTags and in silico generated partial 16S rRNA gene sequences demonstrated significant differences in community structure and phylogenetic resolution across multiple taxonomic levels, including a severe underestimation in the abundance of specific microbial genera involved in nitrogen and methane cycling across the Lake's water column. Thus, PhyloTags provide a reliable adjunct or alternative to cost-effective iTags, enabling more accurate phylogenetic resolution of microbial communities and predictions on their metabolic potential.

OrthoVenn: a web server for genome wide comparison and annotation of orthologous clusters across multiple species

Genome wide analysis of orthologous clusters is an important component of comparative genomics studies. Identifying the overlap among orthologous clusters can enable us to elucidate the function and evolution of proteins across multiple species. Here, we report a web platform named OrthoVenn that is useful for genome wide comparisons and visualization of orthologous clusters. OrthoVenn provides coverage of vertebrates, metazoa, protists, fungi, plants and bacteria for the comparison of orthologous clusters and also supports uploading of customized protein sequences from user-defined species. An interactive Venn diagram, summary counts, and functional summaries of the disjunction and intersection of clusters shared between species are displayed as part of the OrthoVenn result. OrthoVenn also includes in-depth views of the clusters using various sequence analysis tools. Furthermore, OrthoVenn identifies orthologous clusters of single copy genes and allows for a customized search of clusters of specific genes through key words or BLAST. OrthoVenn is an efficient and user-friendly web server freely accessible at http://probes.pw.usda.gov/OrthoVenn or http://aegilops.wheat.ucdavis.edu/OrthoVenn.

OrthoVenn: a web server for genome wide comparison and annotation of orthologous clusters across multiple species

Genome wide analysis of orthologous clusters is an important component of comparative genomics studies. Identifying the overlap among orthologous clusters can enable us to elucidate the function and evolution of proteins across multiple species. Here, we report a web platform named OrthoVenn that is useful for genome wide comparisons and visualization of orthologous clusters. OrthoVenn provides coverage of vertebrates, metazoa, protists, fungi, plants and bacteria for the comparison of orthologous clusters and also supports uploading of customized protein sequences from user-defined species. An interactive Venn diagram, summary counts, and functional summaries of the disjunction and intersection of clusters shared between species are displayed as part of the OrthoVenn result. OrthoVenn also includes in-depth views of the clusters using various sequence analysis tools. Furthermore, OrthoVenn identifies orthologous clusters of single copy genes and allows for a customized search of clusters of specific genes through key words or BLAST. OrthoVenn is an efficient and user-friendly web server freely accessible at http://probes.pw.usda.gov/OrthoVenn or http://aegilops.wheat.ucdavis.edu/OrthoVenn.

Building the crops of tomorrow: advantages of symbiont-based approaches to improving abiotic stress tolerance

The exponential growth in world population is feeding a steadily increasing global need for arable farmland, a resource that is already in high demand. This trend has led to increased farming on subprime arid and semi-arid lands, where limited availability of water and a host of environmental stresses often severely reduce crop productivity. The conventional approach to mitigating the abiotic stresses associated with arid climes is to breed for stress-tolerant cultivars, a time and labor intensive venture that often neglects the complex ecological context of the soil environment in which the crop is grown. In recent years, studies have attempted to identify microbial symbionts capable of conferring the same stress-tolerance to their plant hosts, and new developments in genomic technologies have greatly facilitated such research. Here, we highlight many of the advantages of these symbiont-based approaches and argue in favor of the broader recognition of crop species as ecological niches for a diverse community of microorganisms that function in concert with their plant hosts and each other to thrive under fluctuating environmental conditions.

Natural soil microbes alter flowering phenology and the intensity of selection on flowering time in a wild Arabidopsis relative

Plant phenology is known to depend on many different environmental variables, but soil microbial communities have rarely been acknowledged as possible drivers of flowering time. Here, we tested separately the effects of four naturally occurring soil microbiomes and their constituent soil chemistries on flowering phenology and reproductive fitness of Boechera stricta, a wild relative of Arabidopsis. Flowering time was sensitive to both microbes and the abiotic properties of different soils; varying soil microbiota also altered patterns of selection on flowering time. Thus, soil microbes potentially contribute to phenotypic plasticity of flowering time and to differential selection observed between habitats. We also describe a method to dissect the microbiome into single axes of variation that can help identify candidate organisms whose abundance in soil correlates with flowering time. This approach is broadly applicable to search for microbial community members that alter biological characteristics of interest.

A near complete snapshot of the Zea mays seedling transcriptome revealed from ultra-deep sequencing

RNA-sequencing (RNA-seq) enables in-depth exploration of transcriptomes, but typical sequencing depth often limits its comprehensiveness. In this study, we generated nearly 3 billion RNA-Seq reads, totaling 341 Gb of sequence, from a Zea mays seedling sample. At this depth, a near complete snapshot of the transcriptome was observed consisting of over 90% of the annotated transcripts, including lowly expressed transcription factors. A novel hybrid strategy combining de novo and reference-based assemblies yielded a transcriptome consisting of 126,708 transcripts with 88% of expressed known genes assembled to full-length. We improved current annotations by adding 4,842 previously unannotated transcript variants and many new features, including 212 maize transcripts, 201 genes, 10 genes with undocumented potential roles in seedlings as well as maize lineage specific gene fusion events. We demonstrated the power of deep sequencing for large transcriptome studies by generating a high quality transcriptome, which provides a rich resource for the research community.

DNA methylation, H2A.Z, and the regulation of constitutive expression

The most well-studied function of DNA methylation in eukaryotic cells is the transcriptional silencing of genes and transposons. More recent results showed that many eukaryotes methylate the bodies of genes as well and that this methylation correlates with transcriptional activity rather than repression. The purpose of gene body methylation remains mysterious, but is potentially related to the histone variant H2A.Z. Studies in plants and animals have shown that the genome-wide distributions of H2A.Z and DNA methylation are strikingly anticorrelated. Furthermore, we and other investigators have shown that this relationship is likely to be the result of an ancient but unknown mechanism by which DNA methylation prevents the incorporation of H2A.Z. Recently, we discovered strong correlations between the presence of H2A.Z within gene bodies, the degree to which a gene's expression varies across tissue types or environmental conditions, and transcriptional misregulation in an h2a.z mutant. We propose that one basal function of gene body methylation is the establishment of constitutive expression patterns within housekeeping genes by excluding H2A.Z from their bodies.

Deposition of histone variant H2A.Z within gene bodies regulates responsive genes

The regulation of eukaryotic chromatin relies on interactions between many epigenetic factors, including histone modifications, DNA methylation, and the incorporation of histone variants. H2A.Z, one of the most conserved but enigmatic histone variants that is enriched at the transcriptional start sites of genes, has been implicated in a variety of chromosomal processes. Recently, we reported a genome-wide anticorrelation between H2A.Z and DNA methylation, an epigenetic hallmark of heterochromatin that has also been found in the bodies of active genes in plants and animals. Here, we investigate the basis of this anticorrelation using a novel h2a.z loss-of-function line in Arabidopsis thaliana. Through genome-wide bisulfite sequencing, we demonstrate that loss of H2A.Z in Arabidopsis has only a minor effect on the level or profile of DNA methylation in genes, and we propose that the global anticorrelation between DNA methylation and H2A.Z is primarily caused by the exclusion of H2A.Z from methylated DNA. RNA sequencing and genomic mapping of H2A.Z show that H2A.Z enrichment across gene bodies, rather than at the TSS, is correlated with lower transcription levels and higher measures of gene responsiveness. Loss of H2A.Z causes misregulation of many genes that are disproportionately associated with response to environmental and developmental stimuli. We propose that H2A.Z deposition in gene bodies promotes variability in levels and patterns of gene expression, and that a major function of genic DNA methylation is to exclude H2A.Z from constitutively expressed genes.

Eukaryotes package their DNA to fit within the nucleus using well-conserved proteins, called histones, that form the building blocks of nucleosomes, the fundamental units of chromatin. Histone variants are specialized versions of these proteins that change the chromatin landscape by altering the biochemical properties and interacting partners of the nucleosome. H2A.Z, a conserved eukaryotic histone variant, is preferentially enriched at the beginnings of genes, though the significance of this enrichment remains unknown. We and others have shown that H2A.Z is conspicuously absent from methylated DNA across the genome in plants and animals. Typically considered a mark of epigenetic silencing, DNA methylation has more recently been discovered in the bodies of many genes. Here, we present evidence that the genome-wide anticorrelation between DNA methylation and H2A.Z enrichment in Arabidopsis is the result of DNA methylation acting to prevent H2A.Z incorporation. We demonstrate that the presence of H2A.Z within gene bodies is correlated with lower transcription levels and higher variability in expression patterns across tissue types and environmental conditions, and we propose that a major function of gene-body DNA methylation is to exclude H2A.Z from the bodies of highly and constitutively expressed genes.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Structural and transcriptional analysis of plant genes encoding the bifunctional lysine ketoglutarate reductase saccharopine dehydrogenase enzyme

BACKGROUND

Among the dietary essential amino acids, the most severely limiting in the cereals is lysine. Since cereals make up half of the human diet, lysine limitation has quality/nutritional consequences. The breakdown of lysine is controlled mainly by the catabolic bifunctional enzyme lysine ketoglutarate reductase - saccharopine dehydrogenase (LKR/SDH). The LKR/SDH gene has been reported to produce transcripts for the bifunctional enzyme and separate monofunctional transcripts. In addition to lysine metabolism, this gene has been implicated in a number of metabolic and developmental pathways, which along with its production of multiple transcript types and complex exon/intron structure suggest an important node in plant metabolism. Understanding more about the LKR/SDH gene is thus interesting both from applied standpoint and for basic plant metabolism.

RESULTS

The current report describes a wheat genomic fragment containing an LKR/SDH gene and adjacent genes. The wheat LKR/SDH genomic segment was found to originate from the A-genome of wheat, and EST analysis indicates all three LKR/SDH genes in hexaploid wheat are transcriptionally active. A comparison of a set of plant LKR/SDH genes suggests regions of greater sequence conservation likely related to critical enzymatic functions and metabolic controls. Although most plants contain only a single LKR/SDH gene per genome, poplar contains at least two functional bifunctional genes in addition to a monofunctional LKR gene. Analysis of ESTs finds evidence for monofunctional LKR transcripts in switchgrass, and monofunctional SDH transcripts in wheat, Brachypodium, and poplar.

CONCLUSIONS

The analysis of a wheat LKR/SDH gene and comparative structural and functional analyses among available plant genes provides new information on this important gene. Both the structure of the LKR/SDH gene and the immediately adjacent genes show lineage-specific differences between monocots and dicots, and findings suggest variation in activity of LKR/SDH genes among plants. Although most plant genomes seem to contain a single conserved LKR/SDH gene per genome, poplar possesses multiple contiguous genes. A preponderance of SDH transcripts suggests the LKR region may be more rate-limiting. Only switchgrass has EST evidence for LKR monofunctional transcripts. Evidence for monofunctional SDH transcripts shows a novel intron in wheat, Brachypodium, and poplar.

Rapid development of PCR-based genome-specific repetitive DNA junction markers in wheat

In hexaploid wheat (Triticum aestivum L.) (AABBDD, C=17 000 Mb), repeat DNA accounts for approximately 90% of the genome, of which transposable elements (TEs) constitute 60%-80%. Despite the dynamic evolution of TEs, our previous study indicated that the majority of TEs are conserved and collinear between the homologous wheat genomes, based on identical insertion patterns. In this study, we exploited the unique and abundant TE insertion junction regions identified from diploid Aegilops tauschii to develop genome-specific repeat DNA junction markers (RJM) for use in hexaploid wheat. In this study, both BAC end and random shotgun sequences were used to search for RJM. Of the 300 RJM primer pairs tested, 269 (90%) amplified single bands from diploid Ae. tauschii. Of these 269 primer pairs, 260 (97%) amplified hexaploid wheat and 9 (3%) amplified Ae. tauschii only. Among the RJM primers that amplified hexaploid wheat, 88% were successfully assigned to individual chromosomes of the hexaploid D genome. Among the 38 RJM primers mapped on chromosome 6D, 31 (82%) were unambiguously mapped to delineated bins of the chromosome using various wheat deletion lines. Our results suggest that the unique RJM derived from the diploid D genome could facilitate genetic, physical, and radiation mapping of the hexaploid wheat D genome.

Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.

Rapid development of PCR-based genome-specific repetitive DNA junction markers in wheat

In hexaploid wheat (Triticum aestivum L.) (AABBDD, C=17 000 Mb), repeat DNA accounts for approximately 90% of the genome, of which transposable elements (TEs) constitute 60%-80%. Despite the dynamic evolution of TEs, our previous study indicated that the majority of TEs are conserved and collinear between the homologous wheat genomes, based on identical insertion patterns. In this study, we exploited the unique and abundant TE insertion junction regions identified from diploid Aegilops tauschii to develop genome-specific repeat DNA junction markers (RJM) for use in hexaploid wheat. In this study, both BAC end and random shotgun sequences were used to search for RJM. Of the 300 RJM primer pairs tested, 269 (90%) amplified single bands from diploid Ae. tauschii. Of these 269 primer pairs, 260 (97%) amplified hexaploid wheat and 9 (3%) amplified Ae. tauschii only. Among the RJM primers that amplified hexaploid wheat, 88% were successfully assigned to individual chromosomes of the hexaploid D genome. Among the 38 RJM primers mapped on chromosome 6D, 31 (82%) were unambiguously mapped to delineated bins of the chromosome using various wheat deletion lines. Our results suggest that the unique RJM derived from the diploid D genome could facilitate genetic, physical, and radiation mapping of the hexaploid wheat D genome.

Histone H2A.Z and DNA methylation are mutually antagonistic chromatin marks

Eukaryotic chromatin is separated into functional domains differentiated by post-translational histone modifications, histone variants and DNA methylation. Methylation is associated with repression of transcriptional initiation in plants and animals, and is frequently found in transposable elements. Proper methylation patterns are crucial for eukaryotic development, and aberrant methylation-induced silencing of tumour suppressor genes is a common feature of human cancer. In contrast to methylation, the histone variant H2A.Z is preferentially deposited by the Swr1 ATPase complex near 5' ends of genes where it promotes transcriptional competence. How DNA methylation and H2A.Z influence transcription remains largely unknown. Here we show that in the plant Arabidopsis thaliana regions of DNA methylation are quantitatively deficient in H2A.Z. Exclusion of H2A.Z is seen at sites of DNA methylation in the bodies of actively transcribed genes and in methylated transposons. Mutation of the MET1 DNA methyltransferase, which causes both losses and gains of DNA methylation, engenders opposite changes (gains and losses) in H2A.Z deposition, whereas mutation of the PIE1 subunit of the Swr1 complex that deposits H2A.Z leads to genome-wide hypermethylation. Our findings indicate that DNA methylation can influence chromatin structure and effect gene silencing by excluding H2A.Z, and that H2A.Z protects genes from DNA methylation.

The complete chloroplast genome sequence of Brachypodium distachyon: sequence comparison and phylogenetic analysis of eight grass plastomes

BACKGROUND

Wheat, barley, and rye, of tribe Triticeae in the Poaceae, are among the most important crops worldwide but they present many challenges to genomics-aided crop improvement. Brachypodium distachyon, a close relative of those cereals has recently emerged as a model for grass functional genomics. Sequencing of the nuclear and organelle genomes of Brachypodium is one of the first steps towards making this species available as a tool for researchers interested in cereals biology.

FINDINGS

The chloroplast genome of Brachypodium distachyon was sequenced by a combinational approach using BAC end and shotgun sequences derived from a selected BAC containing the entire chloroplast genome. Comparative analysis indicated that the chloroplast genome is conserved in gene number and organization with respect to those of other cereals. However, several Brachypodium genes evolve at a faster rate than those in other grasses. Sequence analysis reveals that rice and wheat have a ~2.1 kb deletion in their plastid genomes and this deletion must have occurred independently in both species.

CONCLUSION

We demonstrate that BAC libraries can be used to sequence plastid, and likely other organellar, genomes. As expected, the Brachypodium chloroplast genome is very similar to those of other sequenced grasses. The phylogenetic analyses and the pattern of insertions and deletions in the chloroplast genome confirmed that Brachypodium is a close relative of the tribe Triticeae. Nevertheless, we show that some large indels can arise multiple times and may confound phylogenetic reconstruction.

The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences

Due in part to its small genome (approximately 350 Mb), Brachypodium distachyon is emerging as a model system for temperate grasses, including important crops like wheat and barley. We present the analysis of 10.9% of the Brachypodium genome based on 64,696 bacterial artificial chromosome (BAC) end sequences (BES). Analysis of repeat DNA content in BES revealed that approximately 11.0% of the genome consists of known repetitive DNA. The vast majority of the Brachypodium repetitive elements are LTR retrotransposons. While Bare-1 retrotransposons are common to wheat and barley, Brachypodium repetitive element sequence-1 (BRES-1), closely related to Bare-1, is also abundant in Brachypodium. Moreover, unique Brachypodium repetitive element sequences identified constitute approximately 7.4% of its genome. Simple sequence repeats from BES were analyzed, and flanking primer sequences for SSR detection potentially useful for genetic mapping are available at http://brachypodium.pw.usda.gov . Sequence analyses of BES indicated that approximately 21.2% of the Brachypodium genome represents coding sequence. Furthermore, Brachypodium BES have more significant matches to ESTs from wheat than rice or maize, although these species have similar sizes of EST collections. A phylogenetic analysis based on 335 sequences shared among seven grass species further revealed a closer relationship between Brachypodium and Triticeae than Brachypodium and rice or maize.

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5' coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the omega-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region.

New insights into Oryza genome evolution: high gene colinearity and differential retrotransposon amplification

An approximately 247-kb genomic region from FF genome of wild rice Oryza brachyantha, possessing the smallest Oryza genome, was compared to the orthologous approximately 450-kb region from AA genome, O. sativa L. ssp. japonica. 37 of 38 genes in the orthologous regions are shared between japonica and O. brachyantha. Analyses of nucleotide substitution in coding regions suggest the two genomes diverged approximately 10 million years ago. Comparisons of transposable elements (TEs) reveal that the density of DNA TEs in O. brachyantha is comparable to O. sativa; however, the density of RNA TEs is dramatically lower. The genomic fraction of RNA TEs in japonica is two times greater than in O. brachyantha. Differences, particularly in RNA TEs, in this region and in BAC end sequences from five wild and two cultivated Oryza species explain major genome size differences between sativa and brachyantha. Gene expression analyses of three ObDREB1 genes in the sequenced region indicate orthologous genes retain similar expression patterns following cold stress. Our results demonstrate that size and number of RNA TEs play a major role in genomic differentiation and evolution in Oryza. Additionally, distantly related O. brachyantha shares colinearity with O. sativa, offering opportunities to use comparative genomics to explore the genetic diversity of wild species to improve cultivated rice.

Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics

Brachypodium is well suited as a model system for temperate grasses because of its compact genome and a range of biological features. In an effort to develop resources for genome research in this emerging model species, we constructed 2 bacterial artificial chromosome (BAC) libraries from an inbred diploid Brachypodium distachyon line, Bd21, using restriction enzymes HindIII and BamHI. A total of 73,728 clones (36,864 per BAC library) were picked and arrayed in 192,384-well plates. The average insert size for the BamHI and HindIII libraries is estimated to be 100 and 105 kb, respectively, and inserts of chloroplast origin account for 4.4% and 2.4%, respectively. The libraries individually represent 9.4- and 9.9-fold haploid genome equivalents with combined 19.3-fold genome coverage, based on a genome size of 355 Mb reported for the diploid Brachypodium, implying a 99.99% probability that any given specific sequence will be present in each library. Hybridization of the libraries with 8 starch biosynthesis genes was used to empirically evaluate this theoretical genome coverage; the frequency at which these genes were present in the library clones gave an estimated coverage of 11.6- and 19.6-fold genome equivalents. To obtain a first view of the sequence composition of the Brachypodium genome, 2185 BAC end sequences (BES) representing 1.3 Mb of random genomic sequence were compared with the NCBI GenBank database and the GIRI repeat database. Using a cutoff expectation value of E<10-10, only 3.3% of the BESs showed similarity to repetitive sequences in the existing database, whereas 40.0% had matches to the sequences in the EST database, suggesting that a considerable portion of the Brachypodium genome is likely transcribed. When the BESs were compared with individual EST databases, more matches hit wheat than maize, although their EST collections are of a similar size, further supporting the close relationship between Brachypodium and the Triticeae. Moreover, 122 BESs have significant matches to wheat ESTs mapped to individual chromosome bin positions. These BACs represent colinear regions containing the mapped wheat ESTs and would be useful in identifying additional markers for specific wheat chromosome regions.