Sean P. Gordon

Genome diversity in Brachypodium distachyon: deep sequencing of highly diverse inbred lines

Brachypodium distachyon is small annual grass that has been adopted as a model for the grasses. Its small genome, high-quality reference genome, large germplasm collection, and selfing nature make it an excellent subject for studies of natural variation. We sequenced six divergent lines to identify a comprehensive set of polymorphisms and analyze their distribution and concordance with gene expression. Multiple methods and controls were utilized to identify polymorphisms and validate their quality. mRNA-Seq experiments under control and simulated drought-stress conditions, identified 300 genes with a genotype-dependent treatment response. We showed that large-scale sequence variants had extremely high concordance with altered expression of hundreds of genes, including many with genotype-dependent treatment responses. We generated a deep mRNA-Seq dataset for the most divergent line and created a de novo transcriptome assembly. This led to the discovery of >2400 previously unannotated transcripts and hundreds of genes not present in the reference genome. We built a public database for visualization and investigation of sequence variants among these widely used inbred lines.

Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing.

Genes in prokaryotic genomes are often arranged into clusters and co-transcribed into polycistronic RNAs. Isolated examples of polycistronic RNAs were also reported in some higher eukaryotes but their presence was generally considered rare. Here we developed a long-read sequencing strategy to identify polycistronic transcripts in several mushroom forming fungal species including Plicaturopsis crispa, Phanerochaete chrysosporium, Trametes versicolor, and Gloeophyllum trabeum. We found genome-wide prevalence of polycistronic transcription in these Agaricomycetes, involving up to 8% of the transcribed genes. Unlike polycistronic mRNAs in prokaryotes, these co-transcribed genes are also independently transcribed. We show that polycistronic transcription may interfere with expression of the downstream tandem gene. Further comparative genomic analysis indicates that polycistronic transcription is conserved among a wide range of mushroom forming fungi. In summary, our study revealed, for the first time, the genome prevalence of polycistronic transcription in a phylogenetic range of higher fungi. Furthermore, we systematically show that our long-read sequencing approach and combined bioinformatics pipeline is a generic powerful tool for precise characterization of complex transcriptomes that enables identification of mRNA isoforms not recovered via short-read assembly.

Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure

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.The role of differential gene content in the evolution and function of eukaryotic genomes remains poorly explored. Here the authors assemble and annotate the Brachypodium distachyon pan-genome consisting of 54 diverse lines and reveal the differential present genes as a major driver of phenotypic variation.

A complete toolset for the study of Ustilago bromivora and Brachypodium sp. as a fungal-temperate grass pathosystem

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. DOI: http://dx.doi.org/10.7554/eLife.20522.001

Analysis of plant pan-genomes and transcriptomes with GET_HOMOLOGUES-EST, a clustering solution for sequences of the same species

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.

Recreating stable Brachypodium hybridum allotetraploids by uniting the divergent genomes of B. distachyon and B. stacei

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.

Recent Activity in Expanding Populations and Purifying Selection Have Shaped Transposable Element Landscapes across Natural Accessions of the Mediterranean Grass Brachypodium distachyon

Abstract 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.

Widespread polycistronic transcripts in mushroom-forming fungi revealed by single-molecule long-read mRNA sequencing

Genes in prokaryotic genomes are often arranged into clusters and co-transcribed into polycistronic RNAs. Isolated examples of polycistronic RNAs were also reported in some eukaryotes but their presence was generally considered rare. Here we developed a long-read sequencing strategy to identify polycistronic transcripts in several mushroom forming fungal species including Plicaturopsis crispa, Phanerochaete chrysosporium, Trametes versicolor and Gloeophyllum trabeum1. We found genome-wide prevalence of polycistronic transcription in these Agaricomycetes, and it involves up to 8% of the transcribed genes. Unlike polycistronic mRNAs in prokaryotes, these co-transcribed genes are also independently transcribed, and upstream transcription may interfere downstream transcription. Further comparative genomic analysis indicates that polycistronic transcription is likely a feature unique to these fungi. In addition, we also systematically demonstrated that short-read assembly is insufficient for mRNA isoform discovery, especially for isoform-rich loci. In summary, our study revealed, for the first time, the genome prevalence of polycistronic transcription in a subset of fungi. Futhermore, our long-read sequencing approach combined with bioinformatics pipeline is a generic powerful tool for precise characterization of complex transcriptomes.

The Genus Brachypodium as a Model for Perenniality and Polyploidy

The genus Brachypodium contains annual and perennial species with both diploid and polyploid genomes. Like the annual species B. distachyon, some of the perennial and polyploid species have traits compatible with use as a model system (e.g. small genomes, rapid generation time, self-fertile and easy to grow). Thus, there is an opportunity to leverage the resources and knowledge developed for B. distachyon to use other Brachypodium species as models for perenniality and the regulation and evolution of polyploid genomes. There are two factors driving an increased interest in perenniality. First, several perennial grasses are being developed as biomass crops for the sustainable production of biofuel and it would be useful to have a perennial model system to rapidly test biotechnological crop improvement strategies for undesirable impacts on perenniality and winter hardiness. In addition, a deeper understanding of the molecular mechanisms underlying perenniality could be used to design strategies for improving energy crops, for example, by changing resource allocation during growth or by altering the onset of dormancy. The second factor driving increased interest in perenniality is the potential environmental benefits of developing perennial grain crops. B. sylvaticum is a perennial with attributes suitable for use as a perennial model system. A high efficiency transformation system has been developed and a genome sequencing project is underway. Since many important crops, including emerging biomass crops, are polyploid, there is a pressing need to understand the rules governing the evolution and regulation of polyploid genomes. Unfortunately, it is difficult to study polyploid crop genomes because of their size and the difficulty of manipulating those plants in the laboratory. By contrast, B. hybridum has a small polyploid genome and is easy to work with in the laboratory. In addition, analysis of the B. hybridum genome, will be greatly aided by the genome sequences of the two extant diploid species (B. distachyon and B. stacei) that apparently gave rise to B. hybridum. Availability of high quality reference genomes for these three species will be a powerful resource for the study of polyploidy.