Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

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 contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible.

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Phyllosticta paracitricarpa is synonymous with the EU quarantine fungus P. citricarpa based on phylogenomic analyses.

Phyllosticta citricarpa is an important citrus-pathogen and a quarantine organism in the European Union. Its recently described relative, P. paracitricarpa, is very closely related and not listed as a quarantine organism. P. paracitricarpa is very difficult to distinguish from P. citricarpa, since its morphological features overlap and the barcoding gene sequences that were originally used to delimit them as distinct species have a low number of species-specific polymorphisms that have subsequently been shown to overlap between the two clades. Therefore, we performed extensive genomic analyses to determine whether the genetic variation between P. citricarpa and P. paracitricarpa strains should be considered to represent infraspecific variation within P. citricarpa, or whether it is indicative of distinct species. Using a phylogenomic analysis with 3,000 single copy ortholog genes and whole-genome comparisons, we determined that the variation between P. citricarpa and P. paracitricarpa can be considered as infraspecies variation within P. citricarpa. We also determined the level of variation in mitochondrial assemblies of several Phyllosticta species and concluded there are only minimal differences between the assemblies of P. citricarpa and P. paracitricarpa. Thus, using several orthogonal approaches, we here demonstrate that variation within the nuclear and mitochondrial genomes of other Phyllosticta species is larger than variation between genomes obtained from P. citricarpa and P. paracitricarpa strains. Thus, P. citricarpa and P. paracitricarpa should be considered as conspecific.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus.

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales.

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens.

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Construction and comparison of three reference-quality genome assemblies for soybean.

We report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

Transgenic expression of the dicotyledonous pattern recognition receptor EFR in rice leads to ligand-dependent activation of defense responses.

Plant plasma membrane localized pattern recognition receptors (PRRs) detect extracellular pathogen-associated molecules. PRRs such as Arabidopsis EFR and rice XA21 are taxonomically restricted and are absent from most plant genomes. Here we show that rice plants expressing EFR or the chimeric receptor EFR::XA21, containing the EFR ectodomain and the XA21 intracellular domain, sense both Escherichia coli- and Xanthomonas oryzae pv. oryzae (Xoo)-derived elf18 peptides at sub-nanomolar concentrations. Treatment of EFR and EFR::XA21 rice leaf tissue with elf18 leads to MAP kinase activation, reactive oxygen production and defense gene expression. Although expression of EFR does not lead to robust enhanced resistance to fully virulent Xoo isolates, it does lead to quantitatively enhanced resistance to weakly virulent Xoo isolates. EFR interacts with OsSERK2 and the XA21 binding protein 24 (XB24), two key components of the rice XA21-mediated immune response. Rice-EFR plants silenced for OsSERK2, or overexpressing rice XB24 are compromised in elf18-induced reactive oxygen production and defense gene expression indicating that these proteins are also important for EFR-mediated signaling in transgenic rice. Taken together, our results demonstrate the potential feasibility of enhancing disease resistance in rice and possibly other monocotyledonous crop species by expression of dicotyledonous PRRs. Our results also suggest that Arabidopsis EFR utilizes at least a subset of the known endogenous rice XA21 signaling components.

Correction: Transgenic expression of the dicotyledonous pattern recognition receptor EFR in rice leads to ligand-dependent activation of defense responses.

[This corrects the article DOI: 10.1371/journal.ppat.1004809.].

Analysis of clock-regulated genes in Neurospora reveals widespread posttranscriptional control of metabolic potential.

Neurospora crassa has been for decades a principal model for filamentous fungal genetics and physiology as well as for understanding the mechanism of circadian clocks. Eukaryotic fungal and animal clocks comprise transcription-translation-based feedback loops that control rhythmic transcription of a substantial fraction of these transcriptomes, yielding the changes in protein abundance that mediate circadian regulation of physiology and metabolism: Understanding circadian control of gene expression is key to understanding eukaryotic, including fungal, physiology. Indeed, the isolation of clock-controlled genes (ccgs) was pioneered in Neurospora where circadian output begins with binding of the core circadian transcription factor WCC to a subset of ccg promoters, including those of many transcription factors. High temporal resolution (2-h) sampling over 48 h using RNA sequencing (RNA-Seq) identified circadianly expressed genes in Neurospora, revealing that from ∼10% to as much 40% of the transcriptome can be expressed under circadian control. Functional classifications of these genes revealed strong enrichment in pathways involving metabolism, protein synthesis, and stress responses; in broad terms, daytime metabolic potential favors catabolism, energy production, and precursor assembly, whereas night activities favor biosynthesis of cellular components and growth. Discriminative regular expression motif elicitation (DREME) identified key promoter motifs highly correlated with the temporal regulation of ccgs. Correlations between ccg abundance from RNA-Seq, the degree of ccg-promoter activation as reported by ccg-promoter-luciferase fusions, and binding of WCC as measured by ChIP-Seq, are not strong. Therefore, although circadian activation is critical to ccg rhythmicity, posttranscriptional regulation plays a major role in determining rhythmicity at the mRNA level.

Neurospora intermedia from a traditional fermented food enables waste-to-food conversion.

Fungal fermentation of food and agricultural by-products holds promise for improving food sustainability and security. However, the molecular basis of fungal waste-to-food upcycling remains poorly understood. Here we use a multi-omics approach to characterize oncom, a fermented food traditionally produced from soymilk by-products in Java, Indonesia. Metagenomic sequencing of samples from small-scale producers in Western Java indicated that the fungus Neurospora intermedia dominates oncom. Further transcriptomic, metabolomic and phylogenomic analysis revealed that oncom-derived N. intermedia utilizes pectin and cellulose degradation during fermentation and belongs to a genetically distinct subpopulation associated with human-generated by-products. Finally, we found that N. intermedia grew on diverse by-products such as fruit and vegetable pomace and plant-based milk waste, did not encode mycotoxins, and could create foods that were positively perceived by consumers outside Indonesia. These results showcase the traditional significance and future potential of fungal fermentation for creating delicious and nutritious foods from readily available by-products.

High-quality RNA extraction and the regulation of genes encoding cellulosomes are correlated with growth stage in anaerobic fungi.

Anaerobic fungi produce biomass-degrading enzymes and natural products that are important to harness for several biotechnology applications. Although progress has been made in the development of methods for extracting nucleic acids for genomic and transcriptomic sequencing of these fungi, most studies are limited in that they do not sample multiple fungal growth phases in batch culture. In this study, we establish a method to harvest RNA from fungal monocultures and fungal-methanogen co-cultures, and also determine an optimal time frame for high-quality RNA extraction from anaerobic fungi. Based on RNA quality and quantity targets, the optimal time frame in which to harvest anaerobic fungal monocultures and fungal-methanogen co-cultures for RNA extraction was 2-5 days of growth post-inoculation. When grown on cellulose, the fungal strain Anaeromyces robustus cocultivated with the methanogen Methanobacterium bryantii upregulated genes encoding fungal carbohydrate-active enzymes and other cellulosome components relative to fungal monocultures during this time frame, but expression patterns changed at 24-hour intervals throughout the fungal growth phase. These results demonstrate the importance of establishing methods to extract high-quality RNA from anaerobic fungi at multiple time points during batch cultivation.

Genomic prediction of switchgrass winter survivorship across diverse lowland populations.

In the North-Central United States, lowland ecotype switchgrass can increase yield by up to 50% compared with locally adapted but early flowering cultivars. However, lowland ecotypes are not winter tolerant. The mechanism for winter damage is unknown but previously has been associated with late flowering time. This study investigated heading date (measured for two years) and winter survivorship (measured for three years) in a multi-generation population generated from two winter-hardy lowland individuals and diverse southern lowland populations. Sequencing data (311,776 markers) from 1,306 individuals were used to evaluate genome-wide trait prediction through cross-validation and progeny prediction (n = 52). Genetic variance for heading date and winter survivorship was additive with high narrow-sense heritability (0.64 and 0.71, respectively) and reliability (0.68 and 0.76, respectively). The initial negative correlation between winter survivorship and heading date degraded across generations (F1r = -0.43, pseudo-F2r = -0.28, pseudo-F2 progeny r = -0.15). Within-family predictive ability was moderately high for heading date and winter survivorship (0.53 and 0.52, respectively). A multi-trait model did not improve predictive ability for either trait. Progeny predictive ability was 0.71 for winter survivorship and 0.53 for heading date. These results suggest that lowland ecotype populations can obtain sufficient survival rates in the northern United States with two or three cycles of effective selection. Despite accurate genomic prediction, naturally occurring winter mortality successfully isolated winter tolerant genotypes and appears to be an efficient method to develop high-yielding, cold-tolerant switchgrass cultivars.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa.

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi.

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Comparative genomics analysis of drought response between obligate CAM and C3 photosynthesis plants.

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C3 and C4 plants through an inverted pattern of day/night stomatal closure and opening for CO2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C3 species (Arabidopsis thaliana, Populus deltoides 'WV94') share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

Genome of Paspalum vaginatum and the role of trehalose mediated autophagy in increasing maize biomass.

A number of crop wild relatives can tolerate extreme stress to a degree outside the range observed in their domesticated relatives. However, it is unclear whether or how the molecular mechanisms employed by these species can be translated to domesticated crops. Paspalum (Paspalum vaginatum) is a self-incompatible and multiply stress-tolerant wild relative of maize and sorghum. Here, we describe the sequencing and pseudomolecule level assembly of a vegetatively propagated accession of P. vaginatum. Phylogenetic analysis based on 6,151 single-copy syntenic orthologues conserved in 6 related grass species places paspalum as an outgroup of the maize-sorghum clade. In parallel metabolic experiments, paspalum, but neither maize nor sorghum, exhibits a significant increase in trehalose when grown under nutrient-deficit conditions. Inducing trehalose accumulation in maize, imitating the metabolic phenotype of paspalum, results in autophagy dependent increases in biomass accumulation.

A multi-omic characterization of temperature stress in a halotolerant Scenedesmus strain for algal biotechnology.

Microalgae efficiently convert sunlight into lipids and carbohydrates, offering bio-based alternatives for energy and chemical production. Improving algal productivity and robustness against abiotic stress requires a systems level characterization enabled by functional genomics. Here, we characterize a halotolerant microalga Scenedesmus sp. NREL 46B-D3 demonstrating peak growth near 25 °C that reaches 30 g/m2/day and the highest biomass accumulation capacity post cell division reported to date for a halotolerant strain. Functional genomics analysis revealed that genes involved in lipid production, ion channels and antiporters are expanded and expressed. Exposure to temperature stress shifts fatty acid metabolism and increases amino acids synthesis. Co-expression analysis shows that many fatty acid biosynthesis genes are overexpressed with specific transcription factors under cold stress. These and other genes involved in the metabolic and regulatory response to temperature stress can be further explored for strain improvement.