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.

The ectomycorrhizal fungus Pisolithus microcarpus encodes a microRNA involved in cross-kingdom gene silencing during symbiosis.

Small RNAs (sRNAs) are known to regulate pathogenic plant-microbe interactions. Emerging evidence from the study of these model systems suggests that microRNAs (miRNAs) can be translocated between microbes and plants to facilitate symbiosis. The roles of sRNAs in mutualistic mycorrhizal fungal interactions, however, are largely unknown. In this study, we characterized miRNAs encoded by the ectomycorrhizal fungus Pisolithus microcarpus and investigated their expression during mutualistic interaction with Eucalyptus grandis. Using sRNA sequencing data and in situ miRNA detection, a novel fungal miRNA, Pmic_miR-8, was found to be transported into E. grandis roots after interaction with P. microcarpus Further characterization experiments demonstrate that inhibition of Pmic_miR-8 negatively impacts the maintenance of mycorrhizal roots in E. grandis, while supplementation of Pmic_miR-8 led to deeper integration of the fungus into plant tissues. Target prediction and experimental testing suggest that Pmic_miR-8 may target the host NB-ARC domain containing transcripts, suggesting a potential role for this miRNA in subverting host signaling to stabilize the symbiotic interaction. Altogether, we provide evidence of previously undescribed cross-kingdom sRNA transfer from ectomycorrhizal fungi to plant roots, shedding light onto the involvement of miRNAs during the developmental process of mutualistic symbioses.

The regulatory and transcriptional landscape associated with carbon utilization in a filamentous fungus.

Filamentous fungi, such as Neurospora crassa, are very efficient in deconstructing plant biomass by the secretion of an arsenal of plant cell wall-degrading enzymes, by remodeling metabolism to accommodate production of secreted enzymes, and by enabling transport and intracellular utilization of plant biomass components. Although a number of enzymes and transcriptional regulators involved in plant biomass utilization have been identified, how filamentous fungi sense and integrate nutritional information encoded in the plant cell wall into a regulatory hierarchy for optimal utilization of complex carbon sources is not understood. Here, we performed transcriptional profiling of N. crassa on 40 different carbon sources, including plant biomass, to provide data on how fungi sense simple to complex carbohydrates. From these data, we identified regulatory factors in N. crassa and characterized one (PDR-2) associated with pectin utilization and one with pectin/hemicellulose utilization (ARA-1). Using in vitro DNA affinity purification sequencing (DAP-seq), we identified direct targets of transcription factors involved in regulating genes encoding plant cell wall-degrading enzymes. In particular, our data clarified the role of the transcription factor VIB-1 in the regulation of genes encoding plant cell wall-degrading enzymes and nutrient scavenging and revealed a major role of the carbon catabolite repressor CRE-1 in regulating the expression of major facilitator transporter genes. These data contribute to a more complete understanding of cross talk between transcription factors and their target genes, which are involved in regulating nutrient sensing and plant biomass utilization on a global level.

Transcriptomic analysis of field-droughted sorghum from seedling to maturity reveals biotic and metabolic responses.

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

Serendipita Fungi Modulate the Switchgrass Root Transcriptome to Circumvent Host Defenses and Establish a Symbiotic Relationship.

The fungal family Serendipitaceae encompasses root-associated lineages with endophytic, ericoid, orchid, and ectomycorrhizal lifestyles. Switchgrass is an important bioenergy crop for cellulosic ethanol production owing to high biomass production on marginal soils otherwise unfit for food crop cultivation. The aim of this study was to investigate the host plant responses to Serendipita spp. colonization by characterizing the switchgrass root transcriptome during different stages of symbiosis in vitro. For this, we included a native switchgrass strain, Serendipita bescii, and a related strain, S. vermifera, isolated from Australian orchids. Serendipita colonization progresses from thin hyphae that grow between root cells to, finally, the production of large, bulbous hyphae that fill root cells during the later stages of colonization. We report that switchgrass seems to perceive both fungi prior to physical contact, leading to the activation of chemical and structural defense responses and putative host disease resistance genes. Subsequently, the host defense system appears to be quenched and carbohydrate metabolism adjusted, potentially to accommodate the fungal symbiont. In addition, prior to contact, switchgrass exhibited significant increases in root hair density and root surface area. Furthermore, genes involved in phytohormone metabolism such as gibberellin, jasmonic acid, and salicylic acid were activated during different stages of colonization. Both fungal strains induced plant gene expression in a similar manner, indicating a conserved plant response to members of this fungal order. Understanding plant responsiveness to Serendipita spp. will inform our efforts to integrate them into forages and row crops for optimal plant-microbe functioning, thus facilitating low-input, sustainable agricultural practices.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Intra-species genetic variability drives carbon metabolism and symbiotic host interactions in the ectomycorrhizal fungus Pisolithus microcarpus.

Ectomycorrhizal (ECM) fungi are integral to boreal and temperate forest ecosystem functioning and nutrient cycling. ECM fungi, however, originate from diverse saprotrophic lineages and the impacts of genetic variation across species, and especially within a given ECM species, on function and interactions with the environment is not well understood. Here, we explore the extent of intra-species variation between four isolates of the ECM fungus Pisolithus microcarpus, in terms of gene regulation, carbon metabolism and growth, and interactions with a host, Eucalyptus grandis. We demonstrate that, while a core response to the host is maintained by all of the isolates tested, they have distinct patterns of gene expression and carbon metabolism, resulting in the differential expression of isolate-specific response pathways in the host plant. Together, these results highlight the importance of using a wider range of individuals within a species to understand the broader ecological roles of ECM fungi and their host interactions.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus.

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Inorganic nitrogen availability alters Eucalyptus grandis receptivity to the ectomycorrhizal fungus Pisolithus albus but not symbiotic nitrogen transfer.

Forest trees are able to thrive in nutrient-poor soils in part because they obtain growth-limiting nutrients, especially nitrogen (N), through mutualistic symbiosis with ectomycorrhizal (ECM) fungi. Addition of inorganic N into these soils is known to disrupt this mutualism and reduce the diversity of ECM fungi. Despite its ecological impact, the mechanisms governing the observed effects of elevated inorganic N on mycorrhizal communities remain unknown. We address this by using a compartmentalized in vitro system to independently alter nutrients to each symbiont. Using stable isotopes, we traced the nutrient flux under different nutrient regimes between Eucalyptus grandis and its ectomycorrhizal symbiont, Pisolithus albus. We demonstrate that giving E. grandis independent access to N causes a significant reduction in root colonization by P. albus. Transcriptional analysis suggests that the observed reduction in colonization may be caused, in part, by altered transcription of microbe perception genes and defence genes. We show that delivery of N to host leaves is not increased by host nutrient deficiency but by fungal nutrient availability instead. Overall, this advances our understanding of the effects of N fertilization on ECM fungi and the factors governing nutrient transfer in the E. grandis-P. microcarpus interaction.

A fungal transcription factor essential for starch degradation affects integration of carbon and nitrogen metabolism.

In Neurospora crassa, the transcription factor COL-26 functions as a regulator of glucose signaling and metabolism. Its loss leads to resistance to carbon catabolite repression. Here, we report that COL-26 is necessary for the expression of amylolytic genes in N. crassa and is required for the utilization of maltose and starch. Additionally, the Δcol-26 mutant shows growth defects on preferred carbon sources, such as glucose, an effect that was alleviated if glutamine replaced ammonium as the primary nitrogen source. This rescue did not occur when maltose was used as a sole carbon source. Transcriptome and metabolic analyses of the Δcol-26 mutant relative to its wild type parental strain revealed that amino acid and nitrogen metabolism, the TCA cycle and GABA shunt were adversely affected. Phylogenetic analysis showed a single col-26 homolog in Sordariales, Ophilostomatales, and the Magnaporthales, but an expanded number of col-26 homologs in other filamentous fungal species. Deletion of the closest homolog of col-26 in Trichoderma reesei, bglR, resulted in a mutant with similar preferred carbon source growth deficiency, and which was alleviated if glutamine was the sole nitrogen source, suggesting conservation of COL-26 and BglR function. Our finding provides novel insight into the role of COL-26 for utilization of starch and in integrating carbon and nitrogen metabolism for balanced metabolic activities for optimal carbon and nitrogen distribution.

Dichomitus squalens partially tailors its molecular responses to the composition of solid wood.

White-rot fungi, such as Dichomitus squalens, degrade all wood components and inhabit mixed-wood forests containing both soft- and hardwood species. In this study, we evaluated how D. squalens responded to the compositional differences in softwood [guaiacyl (G) lignin and higher mannan content] and hardwood [syringyl/guaiacyl (S/G) lignin and higher xylan content] using semi-natural solid cultures. Spruce (softwood) and birch (hardwood) sticks were degraded by D. squalens as measured by oxidation of the lignins using 2D-NMR. The fungal response as measured by transcriptomics, proteomics and enzyme activities showed a partial tailoring to wood composition. Mannanolytic transcripts and proteins were more abundant in spruce cultures, while a proportionally higher xylanolytic activity was detected in birch cultures. Both wood types induced manganese peroxidases to a much higher level than laccases, but higher transcript and protein levels of the manganese peroxidases were observed on the G-lignin rich spruce. Overall, the molecular responses demonstrated a stronger adaptation to the spruce rather than birch composition, possibly because D. squalens is mainly found degrading softwoods in nature, which supports the ability of the solid wood cultures to reflect the natural environment.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus.

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody 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.].

Fungal and plant gene expression in the Tulasnella calospora-Serapias vomeracea symbiosis provides clues about nitrogen pathways in orchid mycorrhizas.

Orchids are highly dependent on their mycorrhizal fungal partners for nutrient supply, especially during early developmental stages. In addition to organic carbon, nitrogen (N) is probably a major nutrient transferred to the plant because orchid tissues are highly N-enriched. We know almost nothing about the N form preferentially transferred to the plant or about the key molecular determinants required for N uptake and transfer. We identified, in the genome of the orchid mycorrhizal fungus Tulasnella calospora, two functional ammonium transporters and several amino acid transporters but found no evidence of a nitrate assimilation system, in agreement with the N preference of the free-living mycelium grown on different N sources. Differential expression in symbiosis of a repertoire of fungal and plant genes involved in the transport and metabolism of N compounds suggested that organic N may be the main form transferred to the orchid host and that ammonium is taken up by the intracellular fungus from the apoplatic symbiotic interface. This is the first study addressing the genetic determinants of N uptake and transport in orchid mycorrhizas, and provides a model for nutrient exchanges at the symbiotic interface, which may guide future experiments.

Immunofluorescence and fluorescent-protein tagging show high correlation for protein localization in mammalian cells.

Imaging techniques such as immunofluorescence (IF) and the expression of fluorescent protein (FP) fusions are widely used to investigate the subcellular distribution of proteins. Here we report a systematic analysis of >500 human proteins comparing the localizations obtained in live versus fixed cells using FPs and IF, respectively. We identify systematic discrepancies between IF and FPs as well as between FP tagging at the N and C termini. The analysis shows that for 80% of the proteins, IF and FPs yield the same subcellular distribution, and the locations of 250 previously unlocalized proteins were determined by the overlap between the two methods. Approximately 60% of proteins localize to multiple organelles for both methods, indicating a complex subcellular protein organization. These results show that both IF and FP tagging are reliable techniques and demonstrate the usefulness of an integrative approach for a complete investigation of the subcellular human proteome.

Quantitative image analysis approaches for probing Rab GTPase localization and function in mammalian cells.

Membrane traffic pathways play an essential role in cells, providing a mechanism for organelles of the endomembrane system to communicate and exchange material between each other. A significant number of infections and diseases are associated with trafficking pathways, and as such gaining a greater understanding of their regulation is essential. Fluorescence-based imaging techniques are widely used to probe the trafficking machinery within cells, and many of these methods have the potential to be applied in a quantitative manner. In the present mini-review, we highlight several recent examples of how image intensity, kinetic measurements, co-localization and texture feature analysis have been used to study the function of one key family of membrane traffic regulators, the Rab GTPases. We give specific emphasis to the importance of the quantitative nature of these recent studies and comment on their potential applicability to a high-throughput format.

Dual channel rank-based intensity weighting for quantitative co-localization of microscopy images.

BACKGROUND: Accurate quantitative co-localization is a key parameter in the context of understanding the spatial co-ordination of molecules and therefore their function in cells. Existing co-localization algorithms consider either the presence of co-occurring pixels or correlations of intensity in regions of interest. Depending on the image source, and the algorithm selected, the co-localization coefficients determined can be highly variable, and often inaccurate. Furthermore, this choice of whether co-occurrence or correlation is the best approach for quantifying co-localization remains controversial. RESULTS: We have developed a novel algorithm to quantify co-localization that improves on and addresses the major shortcomings of existing co-localization measures. This algorithm uses a non-parametric ranking of pixel intensities in each channel, and the difference in ranks of co-localizing pixel positions in the two channels is used to weight the coefficient. This weighting is applied to co-occurring pixels thereby efficiently combining both co-occurrence and correlation. Tests with synthetic data sets show that the algorithm is sensitive to both co-occurrence and correlation at varying levels of intensity. Analysis of biological data sets demonstrate that this new algorithm offers high sensitivity, and that it is capable of detecting subtle changes in co-localization, exemplified by studies on a well characterized cargo protein that moves through the secretory pathway of cells. CONCLUSIONS: This algorithm provides a novel way to efficiently combine co-occurrence and correlation components in biological images, thereby generating an accurate measure of co-localization. This approach of rank weighting of intensities also eliminates the need for manual thresholding of the image, which is often a cause of error in co-localization quantification. We envisage that this tool will facilitate the quantitative analysis of a wide range of biological data sets, including high resolution confocal images, live cell time-lapse recordings, and high-throughput screening data sets.

Genomics Characterization of an Engineered Corynebacterium glutamicum in Bioreactor Cultivation Under Ionic Liquid Stress.

Corynebacterium glutamicum is an ideal microbial chassis for production of valuable bioproducts including amino acids and next generation biofuels. Here we resequence engineered isopentenol (IP) producing C. glutamicum BRC-JBEI 1.1.2 strain and assess differential transcriptional profiles using RNA sequencing under industrially relevant conditions including scale transition and compare the presence vs absence of an ionic liquid, cholinium lysinate ([Ch][Lys]). Analysis of the scale transition from shake flask to bioreactor with transcriptomics identified a distinct pattern of metabolic and regulatory responses needed for growth in this industrial format. These differential changes in gene expression corroborate altered accumulation of organic acids and bioproducts, including succinate, acetate, and acetoin that occur when cells are grown in the presence of 50 mM [Ch][Lys] in the stirred-tank reactor. This new genome assembly and differential expression analysis of cells grown in a stirred tank bioreactor clarify the cell response of an C. glutamicum strain engineered to produce IP.

Genomic variation within the maize stiff-stalk heterotic germplasm pool.

The stiff-stalk heterotic group in Maize (Zea mays L.) is an important source of inbreds used in U.S. commercial hybrid production. Founder inbreds B14, B37, B73, and, to a lesser extent, B84, are found in the pedigrees of a majority of commercial seed parent inbred lines. We created high-quality genome assemblies of B84 and four expired Plant Variety Protection (ex-PVP) lines LH145 representing B14, NKH8431 of mixed descent, PHB47 representing B37, and PHJ40, which is a Pioneer Hi-Bred International (PHI) early stiff-stalk type. Sequence was generated using long-read sequencing achieving highly contiguous assemblies of 2.13-2.18 Gbp with N50 scaffold lengths >200 Mbp. Inbred-specific gene annotations were generated using a core five-tissue gene expression atlas, whereas transposable element (TE) annotation was conducted using de novo and homology-directed methodologies. Compared with the reference inbred B73, synteny analyses revealed extensive collinearity across the five stiff-stalk genomes, although unique components of the maize pangenome were detected. Comparison of this set of stiff-stalk inbreds with the original Iowa Stiff Stalk Synthetic breeding population revealed that these inbreds represent only a proportion of variation in the original stiff-stalk pool and there are highly conserved haplotypes in released public and ex-Plant Variety Protection inbreds. Despite the reduction in variation from the original stiff-stalk population, substantial genetic and genomic variation was identified supporting the potential for continued breeding success in this pool. The assemblies described here represent stiff-stalk inbreds that have historical and commercial relevance and provide further insight into the emerging maize pangenome.

Mixtures of aromatic compounds induce ligninolytic gene expression in the wood-rotting fungus Dichomitus squalens.

Heterologous production of fungal ligninolytic cocktails is challenging due to the low yields of catalytically active lignin modifying peroxidases. Production using a natural system, such as a wood-rotting fungus, is a promising alternative if specific or preferential induction of the ligninolytic activities could be achieved. Using transcriptomics, gene expression of the white-rot Dichomitus squalens during growth on mixtures of aromatic compounds, with ring structures representing the two major lignin sub-units, was compared to a wood substrate. Most of the genes encoding lignin modifying enzymes (laccases and peroxidases) categorised as highly or moderately expressed on wood were expressed similarly on aromatic compounds. Higher expression levels of a subset of manganese and versatile peroxidases was observed on di- compared to mono-methoxylated aromatics. The expression of polysaccharide degrading enzymes was lower on aromatic compounds compared to wood, demonstrating that the induction of lignin modifying enzymes became more specific. This study suggests potential for aromatic waste streams, e.g. from lignocellulose pretreatment, to produce a lignin-specific enzyme cocktail from D. squalens or other white-rot fungi.