Genomic basis of the distinct biosynthesis of ß-glucogallin, a biochemical marker for hydrolyzable tannin production, in three oak species.

Hydrolyzable tannins (HTs), predominant polyphenols in oaks, are widely used in grape wine aging, feed additives, and human healthcare. However, the limited availability of a high-quality reference genome of oaks greatly hampered the recognition of the mechanism of HT biosynthesis. Here, high-quality reference genomes of three Asian oak species (Quercus variabilis, Quercus aliena, and Quercus dentata) that have different HT contents were generated. Multi-omics studies were carried out to identify key genes regulating HT biosynthesis. In vitro enzyme activity assay was also conducted. Dual-luciferase and yeast one-hybrid assays were used to reveal the transcriptional regulation. Our results revealed that ß-glucogallin was a biochemical marker for HT production in the cupules of the three Asian oaks. UGT84A13 was confirmed as the key enzyme for ß-glucogallin biosynthesis. The differential expression of UGT84A13, rather than enzyme activity, was the main reason for different ß-glucogallin and HT accumulation. Notably, sequence variations in UGT84A13 promoters led to different trans-activating activities of WRKY32/59, explaining the different expression patterns of UGT84A13 among the three species. Our findings provide three high-quality new reference genomes for oak trees and give new insights into different transcriptional regulation for understanding ß-glucogallin and HT biosynthesis in closely related oak species.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees.

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar.

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Genomic basis of seed colour in quinoa inferred from variant patterns using extreme gradient boosting.

Quinoa is an agriculturally important crop species originally domesticated in the Andes of central South America. One of its most important phenotypic traits is seed colour. Seed colour variation is determined by contrasting abundance of betalains, a class of strong antioxidant and free radicals scavenging colour pigments only found in plants of the order Caryophyllales. However, the genetic basis for these pigments in seeds remains to be identified. Here we demonstrate the application of machine learning (extreme gradient boosting) to identify genetic variants predictive of seed colour. We show that extreme gradient boosting outperforms the classical genome-wide association approach. We provide re-sequencing and phenotypic data for 156 South American quinoa accessions and identify candidate genes potentially controlling betalain content in quinoa seeds. Genes identified include novel cytochrome P450 genes and known members of the betalain synthesis pathway, as well as genes annotated as being involved in seed development. Our work showcases the power of modern machine learning methods to extract biologically meaningful information from large sequencing data sets.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

Flexure wood formation via growth reprogramming in hybrid aspen involves jasmonates and polyamines and transcriptional changes resembling tension wood development.

Stem bending in trees induces flexure wood but its properties and development are poorly understood. Here, we investigated the effects of low-intensity multidirectional stem flexing on growth and wood properties of hybrid aspen, and on its transcriptomic and hormonal responses. Glasshouse-grown trees were either kept stationary or subjected to several daily shakes for 5 wk, after which the transcriptomes and hormones were analyzed in the cambial region and developing wood tissues, and the wood properties were analyzed by physical, chemical and microscopy techniques. Shaking increased primary and secondary growth and altered wood differentiation by stimulating gelatinous-fiber formation, reducing secondary wall thickness, changing matrix polysaccharides and increasing cellulose, G- and H-lignin contents, cell wall porosity and saccharification yields. Wood-forming tissues exhibited elevated jasmonate, polyamine, ethylene and brassinosteroids and reduced abscisic acid and gibberellin signaling. Transcriptional responses resembled those during tension wood formation but not opposite wood formation and revealed several thigmomorphogenesis-related genes as well as novel gene networks including FLA and XTH genes encoding plasma membrane-bound proteins. Low-intensity stem flexing stimulates growth and induces wood having improved biorefinery properties through molecular and hormonal pathways similar to thigmomorphogenesis in herbaceous plants and largely overlapping with the tension wood program of hardwoods.

Seiðr: Efficient calculation of robust ensemble gene networks.

Gene regulatory and gene co-expression networks are powerful research tools for identifying biological signal within high-dimensional gene expression data. In recent years, research has focused on addressing shortcomings of these techniques with regard to the low signal-to-noise ratio, non-linear interactions and dataset dependent biases of published methods. Furthermore, it has been shown that aggregating networks from multiple methods provides improved results. Despite this, few useable and scalable software tools have been implemented to perform such best-practice analyses. Here, we present Seidr (stylized Seiðr), a software toolkit designed to assist scientists in gene regulatory and gene co-expression network inference. Seidr creates community networks to reduce algorithmic bias and utilizes noise corrected network backboning to prune noisy edges in the networks. Using benchmarks in real-world conditions across three eukaryotic model organisms, Saccharomyces cerevisiae, Drosophila melanogaster, and Arabidopsis thaliana, we show that individual algorithms are biased toward functional evidence for certain gene-gene interactions. We further demonstrate that the community network is less biased, providing robust performance across different standards and comparisons for the model organisms. Finally, we apply Seidr to a network of drought stress in Norway spruce (Picea abies (L.) H. Krast) as an example application in a non-model species. We demonstrate the use of a network inferred using Seidr for identifying key components, communities and suggesting gene function for non-annotated genes.

Genetic markers and tree properties predicting wood biorefining potential in aspen (Populus tremula) bioenergy feedstock.

BACKGROUND: Wood represents the majority of the biomass on land and constitutes a renewable source of biofuels and other bioproducts. However, wood is recalcitrant to bioconversion, raising a need for feedstock improvement in production of, for instance, biofuels. We investigated the properties of wood that affect bioconversion, as well as the underlying genetics, to help identify superior tree feedstocks for biorefining. RESULTS: We recorded 65 wood-related and growth traits in a population of 113 natural aspen genotypes from Sweden ( https://doi.org/10.5061/dryad.gtht76hrd ). These traits included three growth and field performance traits, 20 traits for wood chemical composition, 17 traits for wood anatomy and structure, and 25 wood saccharification traits as indicators of bioconversion potential. Glucose release after saccharification with acidic pretreatment correlated positively with tree stem height and diameter and the carbohydrate content of the wood, and negatively with the content of lignin and the hemicellulose sugar units. Most of these traits displayed extensive natural variation within the aspen population and high broad-sense heritability, supporting their potential in genetic improvement of feedstocks towards improved bioconversion. Finally, a genome-wide association study (GWAS) revealed 13 genetic loci for saccharification yield (on a whole-tree-biomass basis), with six of them intersecting with associations for either height or stem diameter of the trees. CONCLUSIONS: The simple growth traits of stem height and diameter were identified as good predictors of wood saccharification yield in aspen trees. GWAS elucidated the underlying genetics, revealing putative genetic markers for bioconversion of bioenergy tree feedstocks.

Cone-setting in spruce is regulated by conserved elements of the age-dependent flowering pathway.

Reproductive phase change is well characterized in angiosperm model species, but less studied in gymnosperms. We utilize the early cone-setting acrocona mutant to study reproductive phase change in the conifer Picea abies (Norway spruce), a gymnosperm. The acrocona mutant frequently initiates cone-like structures, called transition shoots, in positions where wild-type P. abies always produces vegetative shoots. We collect acrocona and wild-type samples, and RNA-sequence their messenger RNA (mRNA) and microRNA (miRNA) fractions. We establish gene expression patterns and then use allele-specific transcript assembly to identify mutations in acrocona. We genotype a segregating population of inbred acrocona trees. A member of the SQUAMOSA BINDING PROTEIN-LIKE (SPL) gene family, PaSPL1, is active in reproductive meristems, whereas two putative negative regulators of PaSPL1, miRNA156 and the conifer specific miRNA529, are upregulated in vegetative and transition shoot meristems. We identify a mutation in a putative miRNA156/529 binding site of the acrocona PaSPL1 allele and show that the mutation renders the acrocona allele tolerant to these miRNAs. We show co-segregation between the early cone-setting phenotype and trees homozygous for the acrocona mutation. In conclusion, we demonstrate evolutionary conservation of the age-dependent flowering pathway and involvement of this pathway in regulating reproductive phase change in the conifer P. abies.

Identification of growth regulators using cross-species network analysis in plants.

With the need to increase plant productivity, one of the challenges plant scientists are facing is to identify genes that play a role in beneficial plant traits. Moreover, even when such genes are found, it is generally not trivial to transfer this knowledge about gene function across species to identify functional orthologs. Here, we focused on the leaf to study plant growth. First, we built leaf growth transcriptional networks in Arabidopsis (Arabidopsis thaliana), maize (Zea mays), and aspen (Populus tremula). Next, known growth regulators, here defined as genes that when mutated or ectopically expressed alter plant growth, together with cross-species conserved networks, were used as guides to predict novel Arabidopsis growth regulators. Using an in-depth literature screening, 34 out of 100 top predicted growth regulators were confirmed to affect leaf phenotype when mutated or overexpressed and thus represent novel potential growth regulators. Globally, these growth regulators were involved in cell cycle, plant defense responses, gibberellin, auxin, and brassinosteroid signaling. Phenotypic characterization of loss-of-function lines confirmed two predicted growth regulators to be involved in leaf growth (NPF6.4 and LATE MERISTEM IDENTITY2). In conclusion, the presented network approach offers an integrative cross-species strategy to identify genes involved in plant growth and development.

Metatranscriptomics captures dynamic shifts in mycorrhizal coordination in boreal forests.

Carbon storage and cycling in boreal forests-the largest terrestrial carbon store-is moderated by complex interactions between trees and soil microorganisms. However, existing methods limit our ability to predict how changes in environmental conditions will alter these associations and the essential ecosystem services they provide. To address this, we developed a metatranscriptomic approach to analyze the impact of nutrient enrichment on Norway spruce fine roots and the community structure, function, and tree-microbe coordination of over 350 root-associated fungal species. In response to altered nutrient status, host trees redefined their relationship with the fungal community by reducing sugar efflux carriers and enhancing defense processes. This resulted in a profound restructuring of the fungal community and a collapse in functional coordination between the tree and the dominant Basidiomycete species, and an increase in functional coordination with versatile Ascomycete species. As such, there was a functional shift in community dominance from Basidiomycetes species, with important roles in enzymatically cycling recalcitrant carbon, to Ascomycete species that have melanized cell walls that are highly resistant to degradation. These changes were accompanied by prominent shifts in transcriptional coordination between over 60 predicted fungal effectors, with more than 5,000 Norway spruce transcripts, providing mechanistic insight into the complex molecular dialogue coordinating host trees and their fungal partners. The host-microbe dynamics captured by this study functionally inform how these complex and sensitive biological relationships may mediate the carbon storage potential of boreal soils under changing nutrient conditions.

Molecular basis of differential adventitious rooting competence in poplar genotypes.

Recalcitrant adventitious root (AR) development is a major hurdle in propagating commercially important woody plants. Although significant progress has been made to identify genes involved in subsequent steps of AR development, the molecular basis of differences in apparent recalcitrance to form AR between easy-to-root and difficult-to-root genotypes remains unknown. To address this, we generated cambium tissue-specific transcriptomic data from stem cuttings of hybrid aspen, T89 (difficult-to-root) and hybrid poplar OP42 (easy-to-root), and used transgenic approaches to verify the role of several transcription factors in the control of adventitious rooting. Increased peroxidase activity was positively correlated with better rooting. We found differentially expressed genes encoding reactive oxygen species scavenging proteins to be enriched in OP42 compared with T89. A greater number of differentially expressed transcription factors in cambium cells of OP42 compared with T89 was revealed by a more intense transcriptional reprograming in the former. PtMYC2, a potential negative regulator, was less expressed in OP42 compared with T89. Using transgenic approaches, we demonstrated that PttARF17.1 and PttMYC2.1 negatively regulate adventitious rooting. Our results provide insights into the molecular basis of genotypic differences in AR and implicate differential expression of the master regulator MYC2 as a critical player in this process.

Demographic History and Natural Selection Shape Patterns of Deleterious Mutation Load and Barriers to Introgression across Populus Genome.

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F1 hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

Norway spruce deploys tissue-specific responses during acclimation to cold.

Climate change in the conifer-dominated boreal forest is expected to lead to warmer but more dynamic winter air temperatures, reducing the depth and duration of snow cover and lowering winter soil temperatures. To gain insight into the mechanisms that have enabled conifers to dominate extreme cold environments, we performed genome-wide RNA-Seq analysis from needles and roots of non-dormant two-year Norway spruce (Picea abies (L.) H. Karst), and contrasted these response to herbaceous model Arabidopsis We show that the main transcriptional response of Norway spruce needles exposed to cold was delayed relative to Arabidopsis, and this delay was associated with slower development of freezing tolerance. Despite this difference in timing, Norway spruce principally utilizes early response transcription factors (TFs) belonging to the same gene families as Arabidopsis, indicating broad evolutionary conservation of cold response networks. In keeping with their different metabolic and developmental states, needles and root of Norway spruce showed contrasting results. Regulatory network analysis identified both conserved TFs with known roles in cold acclimation (e.g. homologs of ICE1, AKS3, and of the NAC and AP2/ERF superfamilies), but also a root-specific bHLH101 homolog, providing functional insights into cold stress response strategies in Norway spruce.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE).

BACKGROUND: Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking. RESULTS: We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases. CONCLUSIONS: qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance.

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

Adaptive Introgression Facilitates Adaptation to High Latitudes in European Aspen (Populus tremula L.).

Understanding local adaptation has become a key research area given the ongoing climate challenge and the concomitant requirement to conserve genetic resources. Perennial plants, such as forest trees, are good models to study local adaptation given their wide geographic distribution, largely outcrossing mating systems, and demographic histories. We evaluated signatures of local adaptation in European aspen (Populus tremula) across Europe by means of whole-genome resequencing of a collection of 411 individual trees. We dissected admixture patterns between aspen lineages and observed a strong genomic mosaicism in Scandinavian trees, evidencing different colonization trajectories into the peninsula from Russia, Central and Western Europe. As a consequence of the secondary contacts between populations after the last glacial maximum, we detected an adaptive introgression event in a genome region of ∼500 kb in chromosome 10, harboring a large-effect locus that has previously been shown to contribute to adaptation to the short growing seasons characteristic of Northern Scandinavia. Demographic simulations and ancestry inference suggest an Eastern origin-probably Russian-of the adaptive Nordic allele which nowadays is present in a homozygous state at the north of Scandinavia. The strength of introgression and positive selection signatures in this region is a unique feature in the genome. Furthermore, we detected signals of balancing selection, shared across regional populations, that highlight the importance of standing variation as a primary source of alleles that facilitate local adaptation. Our results, therefore, emphasize the importance of migration-selection balance underlying the genetic architecture of key adaptive quantitative traits.

The BOP-type co-transcriptional regulator NODULE ROOT1 promotes stem secondary growth of the tropical Cannabaceae tree Parasponia andersonii.

Tree stems undergo a massive secondary growth in which secondary xylem and phloem tissues arise from the vascular cambium. Vascular cambium activity is driven by endogenous developmental signalling cues and environmental stimuli. Current knowledge regarding the genetic regulation of cambium activity and secondary growth is still far from complete. The tropical Cannabaceae tree Parasponia andersonii is a non-legume research model of nitrogen-fixing root nodulation. Parasponia andersonii can be transformed efficiently, making it amenable for CRISPR-Cas9-mediated reverse genetics. We considered whether P. andersonii also could be used as a complementary research system to investigate tree-related traits, including secondary growth. We established a developmental map of stem secondary growth in P. andersonii plantlets. Subsequently, we showed that the expression of the co-transcriptional regulator PanNODULE ROOT1 (PanNOOT1) is essential for controlling this process. PanNOOT1 is orthologous to Arabidopsis thaliana BLADE-ON-PETIOLE1 (AtBOP1) and AtBOP2, which are involved in the meristem-to-organ-boundary maintenance. Moreover, in species forming nitrogen-fixing root nodules, NOOT1 is known to function as a key nodule identity gene. Parasponia andersonii CRISPR-Cas9 loss-of-function Pannoot1 mutants are altered in the development of the xylem and phloem tissues without apparent disturbance of the cambium organization and size. Transcriptomic analysis showed that the expression of key secondary growth-related genes is significantly down-regulated in Pannoot1 mutants. This allows us to conclude that PanNOOT1 positively contributes to the regulation of stem secondary growth. Our work also demonstrates that P. andersonii can serve as a tree research system.

Comparative Fungal Community Analyses Using Metatranscriptomics and Internal Transcribed Spacer Amplicon Sequencing from Norway Spruce.

The health, growth, and fitness of boreal forest trees are impacted and improved by their associated microbiomes. Microbial gene expression and functional activity can be assayed with RNA sequencing (RNA-Seq) data from host samples. In contrast, phylogenetic marker gene amplicon sequencing data are used to assess taxonomic composition and community structure of the microbiome. Few studies have considered how much of this structural and taxonomic information is included in transcriptomic data from matched samples. Here, we described fungal communities using both host-derived RNA-Seq and fungal ITS1 DNA amplicon sequencing to compare the outcomes between the methods. We used a panel of root and needle samples from the coniferous tree species Picea abies (Norway spruce) growing in untreated (nutrient-deficient) and nutrient-enriched plots at the Flakaliden forest research site in boreal northern Sweden. We show that the relationship between samples and alpha and beta diversity indicated by the fungal transcriptome is in agreement with that generated by the ITS data, while also identifying a lack of taxonomic overlap due to limitations imposed by current database coverage. Furthermore, we demonstrate how metatranscriptomics data additionally provide biologically informative functional insights. At the community level, there were changes in starch and sucrose metabolism, biosynthesis of amino acids, and pentose and glucuronate interconversions, while processing of organic macromolecules, including aromatic and heterocyclic compounds, was enriched in transcripts assigned to the genus Cortinarius IMPORTANCE A deeper understanding of microbial communities associated with plants is revealing their importance for plant health and productivity. RNA extracted from plant field samples represents the host and other organisms present. Typically, gene expression studies focus on the plant component or, in a limited number of studies, expression in one or more associated organisms. However, metatranscriptomic data are rarely used for taxonomic profiling, which is currently performed using amplicon approaches. We created an assembly-based, reproducible, and hardware-agnostic workflow to taxonomically and functionally annotate fungal RNA-Seq data obtained from Norway spruce roots, which we compared to matching ITS amplicon sequencing data. While we identified some limitations and caveats, we show that functional, taxonomic, and compositional insights can all be obtained from RNA-Seq data. These findings highlight the potential of metatranscriptomics to advance our understanding of interaction, response, and effect between host plants and their associated microbial communities.