An integrated metagenomic, metabolomic and transcriptomic survey of Populus across genotypes and environments.

Bridging molecular information to ecosystem-level processes would provide the capacity to understand system vulnerability and, potentially, a means for assessing ecosystem health. Here, we present an integrated dataset containing environmental and metagenomic information from plant-associated microbial communities, plant transcriptomics, plant and soil metabolomics, and soil chemistry and activity characterization measurements derived from the model tree species Populus trichocarpa. Soil, rhizosphere, root endosphere, and leaf samples were collected from 27 different P. trichocarpa genotypes grown in two different environments leading to an integrated dataset of 318 metagenomes, 98 plant transcriptomes, and 314 metabolomic profiles that are supported by diverse soil measurements. This expansive dataset will provide insights into causal linkages that relate genomic features and molecular level events to system-level properties and their environmental influences.

Predicted structural proteome of Sphagnum divinum and proteome-scale annotation.

MOTIVATION: Sphagnum-dominated peatlands store a substantial amount of terrestrial carbon. The genus is undersampled and under-studied. No experimental crystal structure from any Sphagnum species exists in the Protein Data Bank and fewer than 200 Sphagnum-related genes have structural models available in the AlphaFold Protein Structure Database. Tools and resources are needed to help bridge these gaps, and to enable the analysis of other structural proteomes now made possible by accurate structure prediction. RESULTS: We present the predicted structural proteome (25 134 primary transcripts) of Sphagnum divinum computed using AlphaFold, structural alignment results of all high-confidence models against an annotated nonredundant crystallographic database of over 90,000 structures, a structure-based classification of putative Enzyme Commission (EC) numbers across this proteome, and the computational method to perform this proteome-scale structure-based annotation. AVAILABILITY AND IMPLEMENTATION: All data and code are available in public repositories, detailed at https://github.com/BSDExabio/SAFA. The structural models of the S. divinum proteome have been deposited in the ModelArchive repository at https://modelarchive.org/doi/10.5452/ma-ornl-sphdiv.

Clonality, local population structure and gametophyte sex ratios in cryptic species of the Sphagnum magellanicum complex.

BACKGROUND AND AIMS: Sphagnum (peatmoss) comprises a moss (Bryophyta) clade with ~300-500 species. The genus has unparalleled ecological importance because Sphagnum-dominated peatlands store almost a third of the terrestrial carbon pool and peatmosses engineer the formation and microtopography of peatlands. Genomic resources for Sphagnum are being actively expanded, but many aspects of their biology are still poorly known. Among these are the degree to which Sphagnum species reproduce asexually, and the relative frequencies of male and female gametophytes in these haploid-dominant plants. We assess clonality and gametophyte sex ratios and test hypotheses about the local-scale distribution of clones and sexes in four North American species of the S. magellanicum complex. These four species are difficult to distinguish morphologically and are very closely related. We also assess microbial communities associated with Sphagnum host plant clones and sexes at two sites. METHODS: Four hundred and five samples of the four species, representing 57 populations, were subjected to restriction site-associated DNA sequencing (RADseq). Analyses of population structure and clonality based on the molecular data utilized both phylogenetic and phenetic approaches. Multi-locus genotypes (genets) were identified using the RADseq data. Sexes of sampled ramets were determined using a molecular approach that utilized coverage of loci on the sex chromosomes after the method was validated using a sample of plants that expressed sex phenotypically. Sex ratios were estimated for each species, and populations within species. Difference in fitness between genets was estimated as the numbers of ramets each genet comprised. Degrees of clonality [numbers of genets/numbers of ramets (samples)] within species, among sites, and between gametophyte sexes were estimated. Sex ratios were estimated for each species, and populations within species. Sphagnum-associated microbial communities were assessed at two sites in relation to Sphagnum clonality and sex. KEY RESULTS: All four species appear to engage in a mixture of sexual and asexual (clonal) reproduction. A single ramet represents most genets but two to eight ramets were dsumbers ansd text etected for some genets. Only one genet is represented by ramets in multiple populations; all other genets are restricted to a single population. Within populations ramets of individual genets are spatially clustered, suggesting limited dispersal even within peatlands. Sex ratios are male-biased in S. diabolicum but female-biased in the other three species, although significantly so only in S. divinum. Neither species nor males/females differ in levels of clonal propagation. At St Regis Lake (NY) and Franklin Bog (VT), microbial community composition is strongly differentiated between the sites, but differences between species, genets and sexes were not detected. Within S. divinum, however, female gametophytes harboured two to three times the number of microbial taxa as males. CONCLUSIONS: These four Sphagnum species all exhibit similar reproductive patterns that result from a mixture of sexual and asexual reproduction. The spatial patterns of clonally replicated ramets of genets suggest that these species fall between the so-called phalanx patterns, where genets abut one another but do not extensively mix because of limited ramet fragmentation, and the guerrilla patterns, where extensive genet fragmentation and dispersal result in greater mixing of different genets. Although sex ratios in bryophytes are most often female-biased, both male and female biases occur in this complex of closely related species. The association of far greater microbial diversity for female gametophytes in S. divinum, which has a female-biased sex ratio, suggests additional research to determine if levels of microbial diversity are consistently correlated with differing patterns of sex ratio biases.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss).

BACKGROUND AND AIMS: New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. METHODS: We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. KEY RESULTS: Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. CONCLUSIONS: We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

A rapid assay for assessing bacterial effects on Arabidopsis thermotolerance.

BACKGROUND: The role of beneficial microbes in mitigating plant abiotic stress has received considerable attention. However, the lack of a reproducible and relatively high-throughput screen for microbial contributions to plant thermotolerance has greatly limited progress in this area, this slows the discovery of novel beneficial isolates and the processes by which they operate. RESULTS: We designed a rapid phenotyping method to assess the effects of bacteria on plant host thermotolerance. After testing multiple growth conditions, a hydroponic system was selected and used to optimize an Arabidopsis heat shock regime and phenotypic evaluation. Arabidopsis seedlings germinated on a PTFE mesh disc were floated onto a 6-well plate containing liquid MS media, then subjected to heat shock at 45 °C for various duration. To characterize phenotype, plants were harvested after four days of recovery to measure chlorophyll content. The method was extended to include bacterial isolates and to quantify bacterial contributions to host plant thermotolerance. As an exemplar, the method was used to screen 25 strains of the plant growth promoting Variovorax spp. for enhanced plant thermotolerance. A follow-up study demonstrated the reproducibility of this assay and led to the discovery of a novel beneficial interaction. CONCLUSIONS: This method enables rapid screening of individual bacterial strains for beneficial effects on host plant thermotolerance. The throughput and reproducibility of the system is ideal for testing many genetic variants of Arabidopsis and bacterial strains.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics.

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria.

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Draft Metagenome Sequences of the Sphagnum (Peat Moss) Microbiome from Ambient and Warmed Environments across Europe.

We present 49 metagenome assemblies of the microbiome associated with Sphagnum (peat moss) collected from ambient, artificially warmed, and geothermally warmed conditions across Europe. These data will enable further research regarding the impact of climate change on plant-microbe symbiosis, ecology, and ecosystem functioning of northern peatland ecosystems.

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta).

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. We report herein that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum. The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f4 ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.

Natural selection on a carbon cycling trait drives ecosystem engineering by Sphagnum (peat moss).

Sphagnum peat mosses have an extraordinary impact on the global carbon cycle as they control long-term carbon sequestration in boreal peatland ecosystems. Sphagnum species engineer peatlands, which harbour roughly a quarter of all terrestrial carbon, through peat accumulation by constructing their own niche that allows them to outcompete other plants. Interspecific variation in peat production, largely resulting from differences in tissue decomposability, is hypothesized to drive niche differentiation along microhabitat gradients thereby alleviating competitive pressure. However, little empirical evidence exists for the role of selection in the creation and maintenance of such gradients. In order to document how niche construction and differentiation evolved in Sphagnum, we quantified decomposability for 54 species under natural conditions and used phylogenetic comparative methods to model the evolution of this carbon cycling trait. We show that decomposability tracks the phylogenetic diversification of peat mosses, that natural selection favours different levels of decomposability corresponding to optimum niche and that divergence in this trait occurred early in the evolution of the genus prior to the divergence of most extant species. Our results demonstrate the evolution of ecosystem engineering via natural selection on an extended phenotype, of a fundamental ecosystem process, and one of the Earth's largest soil carbon pools.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus.

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Phylogenomics reveals convergent evolution of red-violet coloration in land plants and the origins of the anthocyanin biosynthetic pathway.

The flavonoids, one of the largest classes of plant secondary metabolites, are found in lineages that span the land plant phylogeny and play important roles in stress responses and as pigments. Perhaps the most well-studied flavonoids are the anthocyanins that have human health benefits and help plants attract pollinators, regulate hormone production, and confer resistance to abiotic and biotic stresses. The canonical biochemical pathway responsible for the production of these pigments is well-characterized for flowering plants yet its conservation across deep divergences in land plants remains debated and poorly understood. Many early land plants such as mosses, liverworts, and ferns produce flavonoid pigments, but their biosynthetic origins and homologies to the anthocyanin pathway remain uncertain. We conducted phylogenetic analyses using full genome sequences representing nearly all major green plant lineages to reconstruct the evolutionary history of the anthocyanin biosynthetic pathway then test the hypothesis that genes in this pathway are present in early land plants. We found that the entire pathway was not intact until the most recent common ancestor of seed plants and that orthologs of many downstream enzymes are absent from seedless plants including mosses, liverworts, and ferns. Our results also highlight the utility of phylogenetic inference, as compared to pairwise sequence similarity, in orthology assessment within large gene families that have complex duplication-loss histories. We suggest that the production of red-violet flavonoid pigments widespread in seedless plants, including the 3-deoxyanthocyanins, requires the activity of novel, as-yet discovered enzymes, and represents convergent evolution of red-violet coloration across land plants.

Functional trait evolution in Sphagnum peat mosses and its relationship to niche construction.

Species in the genus Sphagnum create, maintain, and dominate boreal peatlands through 'extended phenotypes' that allow these organisms to engineer peatland ecosystems and thereby impact global biogeochemical cycles. One such phenotype is the production of peat, or incompletely decomposed biomass, that accumulates when rates of growth exceed decomposition. Interspecific variation in peat production is thought to be responsible for the establishment and maintenance of ecological gradients such as the microtopographic hummock-hollow gradient, along which sympatric species sort within communities. This study investigated the mode and tempo of functional trait evolution across 15 species of Sphagnum using data from the most extensive studies of Sphagnum functional traits to date and phylogenetic comparative methods. We found evidence for phylogenetic conservatism of the niche descriptor height-above-water-table and of traits related to growth, decay and litter quality. However, we failed to detect the influence of phylogeny on interspecific variation in other traits such as shoot density and suggest that environmental context can obscure phylogenetic signal. Trait correlations indicate possible adaptive syndromes that may relate to niche and its construction. This study is the first to formally test the extent to which functional trait variation among Sphagnum species is a result of shared evolutionary history.

The Sphagnome Project: enabling ecological and evolutionary insights through a genus-level sequencing project.

Considerable progress has been made in ecological and evolutionary genetics with studies demonstrating how genes underlying plant and microbial traits can influence adaptation and even 'extend' to influence community structure and ecosystem level processes. Progress in this area is limited to model systems with deep genetic and genomic resources that often have negligible ecological impact or interest. Thus, important linkages between genetic adaptations and their consequences at organismal and ecological scales are often lacking. Here we introduce the Sphagnome Project, which incorporates genomics into a long-running history of Sphagnum research that has documented unparalleled contributions to peatland ecology, carbon sequestration, biogeochemistry, microbiome research, niche construction, and ecosystem engineering. The Sphagnome Project encompasses a genus-level sequencing effort that represents a new type of model system driven not only by genetic tractability, but by ecologically relevant questions and hypotheses.

Hornwort Stomata: Architecture and Fate Shared with 400-Million-Year-Old Fossil Plants without Leaves.

As one of the earliest plant groups to evolve stomata, hornworts are key to understanding the origin and function of stomata. Hornwort stomata are large and scattered on sporangia that grow from their bases and release spores at their tips. We present data from development and immunocytochemistry that identify a role for hornwort stomata that is correlated with sporangial and spore maturation. We measured guard cells across the genera with stomata to assess developmental changes in size and to analyze any correlation with genome size. Stomata form at the base of the sporophyte in the green region, where they develop differential wall thickenings, form a pore, and die. Guard cells collapse inwardly, increase in surface area, and remain perched over a substomatal cavity and network of intercellular spaces that is initially fluid filled. Following pore formation, the sporophyte dries from the outside inwardly and continues to do so after guard cells die and collapse. Spore tetrads develop in spore mother cell walls within a mucilaginous matrix, both of which progressively dry before sporophyte dehiscence. A lack of correlation between guard cell size and DNA content, lack of arabinans in cell walls, and perpetually open pores are consistent with the inactivity of hornwort stomata. Stomata are expendable in hornworts, as they have been lost twice in derived taxa. Guard cells and epidermal cells of hornworts show striking similarities with the earliest plant fossils. Our findings identify an architecture and fate of stomata in hornworts that is ancient and common to plants without sporophytic leaves.

Infection with Salmonella enterica Serovar Typhimurium Leads to Increased Proportions of F4/80+ Red Pulp Macrophages and Decreased Proportions of B and T Lymphocytes in the Spleen.

Infection of mice with Salmonella enterica serovar Typhimurium (Salmonella) causes systemic inflammatory disease and enlargement of the spleen (splenomegaly). Splenomegaly has been attributed to a general increase in the numbers of phagocytes, lymphocytes, as well as to the expansion of immature CD71+Ter119+ reticulocytes. The spleen is important for recycling senescent red blood cells (RBCs) and for the capture and eradication of blood-borne pathogens. Conservation of splenic tissue architecture, comprised of the white pulp (WP), marginal zone (MZ), and red pulp (RP) is essential for initiation of adaptive immune responses to captured pathogens. Using flow cytometry and four color immunofluorescence microscopy (IFM), we show that Salmonella-induced splenomegaly is characterized by drastic alterations of the splenic tissue architecture and cell population proportions, as well as in situ cell distributions. A major cause of splenomegaly appears to be the significant increase in immature RBC precursors and F4/80+ macrophages that are important for recycling of heme-associated iron. In contrast, the proportions of B220+, CD4+ and CD8+ lymphocytes, as well as MZ MOMA+ macrophages decrease significantly as infection progresses. Spleen tissue sections show visible tears and significantly altered tissue architecture with F4/80+ macrophages and RBCs expanding beyond the RP and taking over most of the spleen tissue. Additionally, F4/80+ macrophages actively phagocytose not only RBCs, but also lymphocytes, indicating that they may contribute to declining lymphocyte proportions during Salmonella infection. Understanding how these alterations of spleen microarchitecture impact the generation of adaptive immune responses to Salmonella has implications for understanding Salmonella pathogenesis and for the design of more effective Salmonella-based vaccines.