Cell-type-specific transcriptomics uncovers spatial regulatory networks in bioenergy sorghum stems.

Bioenergy sorghum is a low-input, drought-resilient, deep-rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum's 4-5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about cell-type gene expression and regulation in stems was available to enable engineering. To obtain this information, laser capture microdissection was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell-type-specific and cell-preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell-type-specific gene regulatory networks (GRNs) revealed that unique transcription factor families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell-specific transcriptome data was combined with known secondary cell wall (SCW) networks to identify the GRNs that differentially activate SCW formation in vascular sclerenchyma and epidermal cells. The spatial transcriptomic dataset provides a valuable source of information about the function of different sorghum cell types and GRNs that will enable the engineering of bioenergy sorghum stems, and an interactive web application developed during this project will allow easy access and exploration of the data (https://mc-lab.shinyapps.io/lcm-dataset/).

Intracellular pathways for lignin catabolism in white-rot fungi.

Lignin is a biopolymer found in plant cell walls that accounts for 30% of the organic carbon in the biosphere. White-rot fungi (WRF) are considered the most efficient organisms at degrading lignin in nature. While lignin depolymerization by WRF has been extensively studied, the possibility that WRF are able to utilize lignin as a carbon source is still a matter of controversy. Here, we employ 13C-isotope labeling, systems biology approaches, and in vitro enzyme assays to demonstrate that two WRF, Trametes versicolor and Gelatoporia subvermispora, funnel carbon from lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways. These results provide insights into global carbon cycling in soil ecosystems and furthermore establish a foundation for employing WRF in simultaneous lignin depolymerization and bioconversion to bioproducts-a key step toward enabling a sustainable bioeconomy.

Capturing an Early Gene Induction Event during Wood Decay by the Brown Rot Fungus Rhodonia placenta.

Brown rot fungi dominate wood decomposition in coniferous forests, and their carbohydrate-selective mechanisms are of commercial interest. Brown rot was recently described as a two-step, sequential mechanism orchestrated by fungi using differentially expressed genes (DEGs) and consisting of oxidation via reactive oxygen species (ROS) followed by enzymatic saccharification. There have been indications, however, that the initial oxidation step itself might require induction. To capture this early gene regulation event, here, we integrated fine-scale cryosectioning with whole-transcriptome sequencing to dissect gene expression at the single-hyphal-cell scale (tens of micrometers). This improved the spatial resolution 50-fold, relative to previous work, and we were able to capture the activity of the first 100 µm of hyphal front growth by Rhodonia placenta in aspen wood. This early decay period was dominated by delayed gene expression patterns as the fungus ramped up its mechanism. These delayed DEGs included many genes implicated in ROS pathways (lignocellulose oxidation [LOX]) that were previously and incorrectly assumed to be constitutively expressed. These delayed DEGs, which include those with and without predicted functions, also create a focused subset of target genes for functional genomics. However, this delayed pattern was not universal, with a few genes being upregulated immediately at the hyphal front. Most notably, this included a gene commonly implicated in hydroquinone and iron redox cycling: benzoquinone reductase. IMPORTANCE Earth's aboveground terrestrial biomass is primarily wood, and fungi dominate wood decomposition. Here, we studied these fungal pathways in a common "brown rot"-type fungus, Rhodonia placenta, that selectively extracts sugars from carbohydrates embedded within wood lignin. Using a space-for-time design to map fungal gene expression at the extreme hyphal front in wood, we made two discoveries. First, we found that many genes long assumed to be "on" (constitutively expressed) from the very beginning of decay were instead "off" before being upregulated, when mapped (via transcriptome sequencing [RNA-seq]) at a high resolution. Second, we found that the gene encoding benzoquinone reductase was "on" in incipient decay and quickly downregulated, implying a key role in "kick-starting" brown rot.

Colonies of the fungus Aspergillus niger are highly differentiated to adapt to local carbon source variation.

Saprobic fungi, such as Aspergillus niger, grow as colonies consisting of a network of branching and fusing hyphae that are often considered to be relatively uniform entities in which nutrients can freely move through the hyphae. In nature, different parts of a colony are often exposed to different nutrients. We have investigated, using a multi-omics approach, adaptation of A. niger colonies to spatially separated and compositionally different plant biomass substrates. This demonstrated a high level of intra-colony differentiation, which closely matched the locally available substrate. The part of the colony exposed to pectin-rich sugar beet pulp and to xylan-rich wheat bran showed high pectinolytic and high xylanolytic transcript and protein levels respectively. This study therefore exemplifies the high ability of fungal colonies to differentiate and adapt to local conditions, ensuring efficient use of the available nutrients, rather than maintaining a uniform physiology throughout the colony.

Fluctuation localization imaging-based fluorescence in situ hybridization (fliFISH) for accurate detection and counting of RNA copies in single cells.

Quantitative gene expression analysis in intact single cells can be achieved using single molecule-based fluorescence in situ hybridization (smFISH). This approach relies on fluorescence intensity to distinguish between true signals, emitted from an RNA copy hybridized with multiple oligonucleotide probes, and background noise. Thus, the precision in smFISH is often compromised by partial or nonspecific probe binding and tissue autofluorescence, especially when only a small number of probes can be fitted to the target transcript. Here we provide an accurate approach for setting quantitative thresholds between true and false signals, which relies on on-off duty cycles of photoswitchable dyes. This fluctuation localization imaging-based FISH (fliFISH) uses on-time fractions (measured over a series of exposures) collected from transcripts bound to as low as 8 probes, which are distinct from on-time fractions collected from nonspecifically bound probes or autofluorescence. Using multicolor fliFISH, we identified radial gene expression patterns in mouse pancreatic islets for insulin, the transcription factor, NKX2-2 and their ratio (Nkx2-2/Ins2). These radial patterns, showing higher values in ß cells at the islet core and lower values in peripheral cells, were lost in diabetic mouse islets. In summary, fliFISH provides an accurate, quantitative approach for detecting and counting true RNA copies and rejecting false signals by their distinct on-time fractions, laying the foundation for reliable single-cell transcriptomics.

Reference genes for accurate normalization of gene expression in wood-decomposing fungi.

Wood-decomposing fungi efficiently decompose plant lignocellulose, and there is increasing interest in characterizing and perhaps harnessing the fungal gene regulation strategies that enable wood decomposition. Proper interpretation of these fungal mechanisms relies on accurate quantification of gene expression, demanding reliable internal control genes (ICGs) as references. Commonly used ICGs such as actin, however, fluctuate among wood-decomposing fungi under defined conditions. In this study, by mining RNA-seq data in silico and validating ICGs in vitro using qRT-PCR, we targeted more reliable ICGs for studying transcriptional responses in wood-decomposing fungi, particularly responses to changing environments (e.g., carbon sources, decomposition stages) in various culture conditions. Using the model brown rot fungus Postia placenta in a first-pass study, our mining efforts yielded 15 constitutively-expressed genes robust in variable carbon sources (e.g., no carbon, glucose, cellobiose, aspen) and cultivation stages (e.g., 15 h, 72 h) in submerged cultures. Of these, we found 7 genes as most suitable ICGs. Expression stabilities of these newly selected ICGs were better than commonly used ICGs, analyzed by NormFinder algorithm and qRT-PCR. In a second-pass, multi-species study in solid wood, our RNA-seq mining efforts revealed hundreds of highly constitutively expressed genes among four wood-decomposing fungi with varying nutritional modes (brown rot, white rot), including a shared core set of ICGs numbering 11 genes. Together, the newly selected ICGs highlighted here will increase reliability when studying gene regulatory mechanisms of wood-decomposing fungi.

Rosette core fungal resistance in Arabidopsis thaliana.

MAIN CONCLUSION: Unlike rosette leaves, the mature Arabidopsis rosette core can display full resistance to Botrytis cinerea revealing the importance for spatial and developmental aspects of plant fungal resistance. Arabidopsis thaliana is a model host to investigate plant defense against fungi. However, many of the reports investigating Arabidopsis fungal defense against the necrotrophic fungus, Botrytis cinerea, utilize rosette leaves as host tissue. Here we report organ-dependent differences in B. cinerea resistance of Arabidopsis. Although wild-type Arabidopsis rosette leaves mount a jasmonate-dependent defense that slows fungal growth, this defense is incapable of resisting fungal devastation. In contrast, as the fungus spreads through infected leaf petioles towards the plant center, or rosette core, there is a jasmonate- and age-dependent fungal penetration blockage into the rosette core. We report evidence for induced and preformed resistance in the rosette core, as direct rosette core inoculation can also result in resistance, but at a lower penetrance relative to infections that approach the core from infected leaf petioles. The Arabidopsis rosette core displays a distinct transcriptome relative to other plant organs, and BLADE ON PETIOLE (BOP) transcripts are abundant in the rosette core. The BOP genes, with known roles in abscission zone formation, are required for full Arabidopsis rosette core B. cinerea resistance, suggesting a possible role for BOP-dependent modifications that may help to restrict fungal susceptibility of the rosette core. Finally, we demonstrate that cabbage and cauliflower, common Brassicaceae crops, also display leaf susceptibility and rosette core resistance to B. cinerea that can involve leaf abscission. Thus, spatial and developmental aspects of plant host resistance play critical roles in resistance to necrotrophic fungal pathogens and are important to our understanding of plant defense mechanisms.

Wetland Sediments Host Diverse Microbial Taxa Capable of Cycling Alcohols.

Alcohols are commonly derived from the degradation of organic matter and yet are rarely measured in environmental samples. Wetlands in the Prairie Pothole Region (PPR) support extremely high methane emissions and the highest sulfate reduction rates reported to date, likely contributing to a significant proportion of organic matter mineralization in this system. While ethanol and isopropanol concentrations up to 4 to 5 mM in PPR wetland pore fluids have been implicated in sustaining these high rates of microbial activity, the mechanisms that support alcohol cycling in this ecosystem are poorly understood. We leveraged metagenomic and transcriptomic tools to identify genes, pathways, and microorganisms potentially accounting for alcohol cycling in PPR wetlands. Phylogenetic analyses revealed diverse alcohol dehydrogenases and putative substrates. Alcohol dehydrogenase and aldehyde dehydrogenase genes were included in 62 metagenome-assembled genomes (MAGs) affiliated with 16 phyla. The most frequently encoded pathway (in 30 MAGs) potentially accounting for alcohol production was a Pyrococcus furiosus-like fermentation which can involve pyruvate:ferredoxin oxidoreductase (PFOR). Transcripts for 93 of 137 PFOR genes in these MAGs were detected, as well as for 158 of 243 alcohol dehydrogenase genes retrieved from these same MAGs. Mixed acid fermentation and heterofermentative lactate fermentation were also frequently encoded. Finally, we identified 19 novel putative isopropanol dehydrogenases in 15 MAGs affiliated with Proteobacteria, Acidobacteria, Chloroflexi, Planctomycetes, Ignavibacteriae, Thaumarchaeota, and the candidate divisions KSB1 and Rokubacteria We conclude that diverse microorganisms may use uncommon and potentially novel pathways to produce ethanol and isopropanol in PPR wetland sediments.IMPORTANCE Understanding patterns of organic matter degradation in wetlands is essential for identifying the substrates and mechanisms supporting greenhouse gas production and emissions from wetlands, the main natural source of methane in the atmosphere. Alcohols are common fermentation products but are poorly studied as key intermediates in organic matter degradation in wetlands. By investigating genes, pathways, and microorganisms potentially accounting for the high concentrations of ethanol and isopropanol measured in Prairie Pothole wetland sediments, this work advanced our understanding of alcohol fermentations in wetlands linked to extremely high greenhouse gas emissions. Moreover, the novel alcohol dehydrogenases and microbial taxa potentially involved in alcohol metabolism may serve biotechnological efforts in bioengineering commercially valuable alcohol production and in the discovery of novel isopropanol producers or isopropanol fermentation pathways.

Network analysis of transcriptomics expands regulatory landscapes in Synechococcus sp. PCC 7002.

Cyanobacterial regulation of gene expression must contend with a genome organization that lacks apparent functional context, as the majority of cellular processes and metabolic pathways are encoded by genes found at disparate locations across the genome and relatively few transcription factors exist. In this study, global transcript abundance data from the model cyanobacterium Synechococcus sp. PCC 7002 grown under 42 different conditions was analyzed using Context-Likelihood of Relatedness (CLR). The resulting network, organized into 11 modules, provided insight into transcriptional network topology as well as grouping genes by function and linking their response to specific environmental variables. When used in conjunction with genome sequences, the network allowed identification and expansion of novel potential targets of both DNA binding proteins and sRNA regulators. These results offer a new perspective into the multi-level regulation that governs cellular adaptations of the fast-growing physiologically robust cyanobacterium Synechococcus sp. PCC 7002 to changing environmental variables. It also provides a methodological high-throughput approach to studying multi-scale regulatory mechanisms that operate in cyanobacteria. Finally, it provides valuable context for integrating systems-level data to enhance gene grouping based on annotated function, especially in organisms where traditional context analyses cannot be implemented due to lack of operon-based functional organization.

The Toxoplasma gondii cyst wall protein CST1 is critical for cyst wall integrity and promotes bradyzoite persistence.

Toxoplasma gondii infects up to one third of the world's population. A key to the success of T. gondii as a parasite is its ability to persist for the life of its host as bradyzoites within tissue cysts. The glycosylated cyst wall is the key structural feature that facilitates persistence and oral transmission of this parasite. Because most of the antibodies and reagents that recognize the cyst wall recognize carbohydrates, identification of the components of the cyst wall has been technically challenging. We have identified CST1 (TGME49_064660) as a 250 kDa SRS (SAG1 related sequence) domain protein with a large mucin-like domain. CST1 is responsible for the Dolichos biflorus Agglutinin (DBA) lectin binding characteristic of T. gondii cysts. Deletion of CST1 results in reduced cyst number and a fragile brain cyst phenotype characterized by a thinning and disruption of the underlying region of the cyst wall. These defects are reversed by complementation of CST1. Additional complementation experiments demonstrate that the CST1-mucin domain is necessary for the formation of a normal cyst wall structure, the ability of the cyst to resist mechanical stress, and binding of DBA to the cyst wall. RNA-seq transcriptome analysis demonstrated dysregulation of bradyzoite genes within the various cst1 mutants. These results indicate that CST1 functions as a key structural component that confers essential sturdiness to the T. gondii tissue cyst critical for persistence of bradyzoite forms.

Proteome-scale tissue mapping using mass spectrometry based on label-free and multiplexed workflows.

Multiplexed bimolecular profiling of tissue microenvironment, or spatial omics, can provide deep insight into cellular compositions and interactions in both normal and diseased tissues. Proteome-scale tissue mapping, which aims to unbiasedly visualize all the proteins in whole tissue section or region of interest, has attracted significant interest because it holds great potential to directly reveal diagnostic biomarkers and therapeutic targets. While many approaches are available, however, proteome mapping still exhibits significant technical challenges in both protein coverage and analytical throughput. Since many of these existing challenges are associated with mass spectrometry-based protein identification and quantification, we performed a detailed benchmarking study of three protein quantification methods for spatial proteome mapping, including label-free, TMT-MS2, and TMT-MS3. Our study indicates label-free method provided the deepest coverages of ~3500 proteins at a spatial resolution of 50 µm and the largest quantification dynamic range, while TMT-MS2 method holds great benefit in mapping throughput at >125 pixels per day. The evaluation also indicates both label-free and TMT-MS2 provide robust protein quantifications in terms of identifying differentially abundant proteins and spatially co-variable clusters. In the study of pancreatic islet microenvironment, we demonstrated deep proteome mapping not only enables to identify protein markers specific to different cell types, but more importantly, it also reveals unknown or hidden protein patterns by spatial co-expression analysis.

A crucial role for dynamic expression of components encoding the negative arm of the circadian clock.

In the Neurospora circadian system, the White Collar Complex (WCC) drives expression of the principal circadian negative arm component frequency ( frq ). FRQ interacts with FRH (FRQ-interacting helicase) and CK-1 forming a stable complex that represses its own expression by inhibiting WCC. In this study, a genetic screen identified a gene, designated as brd-8 , that encodes a conserved auxiliary subunit of the NuA4 histone acetylation complex. Loss of brd-8 reduces H4 acetylation and RNA polymerase (Pol) II occupancy at frq and other known circadian genes, and leads to a long circadian period, delayed phase, and defective overt circadian output at some temperatures. In addition to strongly associating with the NuA4 histone acetyltransferase complex, BRD-8 is also found complexed with the transcription elongation regulator BYE-1. Expression of brd-8, bye-1, histone hH2Az , and several NuA4 subunits is controlled by the circadian clock, indicating that the molecular clock both regulates the basic chromatin status and is regulated by changes in chromatin. Taken together, our data identify new auxiliary elements of the fungal NuA4 complex having homology to mammalian components, which along with conventional NuA4 subunits, are required for timely and dynamic frq expression and thereby a normal and persistent circadian rhythm.

A crucial role for dynamic expression of components encoding the negative arm of the circadian clock.

In the Neurospora circadian system, the White Collar Complex (WCC) drives expression of the principal circadian negative arm component frequency (frq). FRQ interacts with FRH (FRQ-interacting RNA helicase) and CKI, forming a stable complex that represses its own expression by inhibiting WCC. In this study, a genetic screen identified a gene, designated as brd-8, that encodes a conserved auxiliary subunit of the NuA4 histone acetylation complex. Loss of brd-8 reduces H4 acetylation and RNA polymerase (Pol) II occupancy at frq and other known circadian genes, and leads to a long circadian period, delayed phase, and defective overt circadian output at some temperatures. In addition to strongly associating with the NuA4 histone acetyltransferase complex, BRD-8 is also found complexed with the transcription elongation regulator BYE-1. Expression of brd-8, bye-1, histone h2a.z, and several NuA4 subunits is controlled by the circadian clock, indicating that the molecular clock both regulates the basic chromatin status and is regulated by changes in chromatin. Taken together, our data identify auxiliary elements of the fungal NuA4 complex having homology to mammalian components, which along with conventional NuA4 subunits, are required for timely and dynamic frq expression and thereby a normal and persistent circadian rhythm.

Spatial Proteomics toward Subcellular Resolution by Coupling Deep Ultraviolet Laser Ablation with Nanodroplet Sample Preparation.

Multiplexed molecular profiling of tissue microenvironments, or spatial omics, can provide critical insights into cellular functions and disease pathology. The coupling of laser microdissection with mass spectrometry-based proteomics has enabled deep and unbiased mapping of >1000 proteins. However, the throughput of laser microdissection is often limited due to tedious two-step procedures, sequential laser cutting, and sample collection. The two-step procedure also hinders the further improvement of spatial resolution to <10 µm as needed for subcellular proteomics. Herein, we developed a high-throughput and high-resolution spatial proteomics platform by seamlessly coupling deep ultraviolet (DUV) laser ablation (LA) with nanoPOTS (Nanodroplet Processing in One pot for Trace Samples)-based sample preparation. We demonstrated the DUV-LA system can quickly isolate and collect tissue samples at a throughput of ∼30 spots/min and a spatial resolution down to 2 µm from a 10 µm thick human pancreas tissue section. To improve sample recovery, we developed a proximity aerosol collection approach by placing DMSO droplets close to LA spots. We demonstrated the DUV-LA-nanoPOTS platform can detect an average of 1312, 1533, and 1966 proteins from ablation spots with diameters of 7, 13, and 19 µm, respectively. In a proof-of-concept study, we isolated and profiled two distinct subcellular regions of the pancreas tissue revealed by hematoxylin and eosin (H&E) staining. Quantitative proteomics revealed proteins specifically enriched to subcellular compartments.

Oligosaccharide production and signaling correlate with delayed flowering in an Arabidopsis genotype grown and selected in high [CO2].

Since industrialization began, atmospheric CO2 ([CO2]) has increased from 270 to 415 ppm and is projected to reach 800-1000 ppm this century. Some Arabidopsis thaliana (Arabidopsis) genotypes delayed flowering in elevated [CO2] relative to current [CO2], while others showed no change or accelerations. To predict genotype-specific flowering behaviors, we must understand the mechanisms driving flowering response to rising [CO2]. [CO2] changes alter photosynthesis and carbohydrates in plants. Plants sense carbohydrate levels, and exogenous carbohydrate application influences flowering time and flowering transcript levels. We asked how organismal changes in carbohydrates and transcription correlate with changes in flowering time under elevated [CO2]. We used a genotype (SG) of Arabidopsis that was selected for high fitness at elevated [CO2] (700 ppm). SG delays flowering under elevated [CO2] (700 ppm) relative to current [CO2] (400 ppm). We compared SG to a closely related control genotype (CG) that shows no [CO2]-induced flowering change. We compared metabolomic and transcriptomic profiles in these genotypes at current and elevated [CO2] to assess correlations with flowering in these conditions. While both genotypes altered carbohydrates in response to elevated [CO2], SG had higher levels of sucrose than CG and showed a stronger increase in glucose and fructose in elevated [CO2]. Both genotypes demonstrated transcriptional changes, with CG increasing genes related to fructose 1,6-bisphosphate breakdown, amino acid synthesis, and secondary metabolites; and SG decreasing genes related to starch and sugar metabolism, but increasing genes involved in oligosaccharide production and sugar modifications. Genes associated with flowering regulation within the photoperiod, vernalization, and meristem identity pathways were altered in these genotypes. Elevated [CO2] may alter carbohydrates to influence transcription in both genotypes and delayed flowering in SG. Changes in the oligosaccharide pool may contribute to delayed flowering in SG. This work extends the literature exploring genotypic-specific flowering responses to elevated [CO2].

Diverse electron carriers drive syntrophic interactions in an enriched anaerobic acetate-oxidizing consortium.

In many anoxic environments, syntrophic acetate oxidation (SAO) is a key pathway mediating the conversion of acetate into methane through obligate cross-feeding interactions between SAO bacteria (SAOB) and methanogenic archaea. The SAO pathway is particularly important in engineered environments such as anaerobic digestion (AD) systems operating at thermophilic temperatures and/or with high ammonia. Despite the widespread importance of SAOB to the stability of the AD process, little is known about their in situ physiologies due to typically low biomass yields and resistance to isolation. Here, we performed a long-term (300-day) continuous enrichment of a thermophilic (55 °C) SAO community from a municipal AD system using acetate as the sole carbon source. Over 80% of the enriched bioreactor metagenome belonged to a three-member consortium, including an acetate-oxidizing bacterium affiliated with DTU068 encoding for carbon dioxide, hydrogen, and formate production, along with two methanogenic archaea affiliated with Methanothermobacter_A. Stable isotope probing was coupled with metaproteogenomics to quantify carbon flux into each community member during acetate conversion and inform metabolic reconstruction and genome-scale modeling. This effort revealed that the two Methanothermobacter_A species differed in their preferred electron donors, with one possessing the ability to grow on formate and the other only consuming hydrogen. A thermodynamic analysis suggested that the presence of the formate-consuming methanogen broadened the environmental conditions where ATP production from SAO was favorable. Collectively, these results highlight how flexibility in electron partitioning during SAO likely governs community structure and fitness through thermodynamic-driven mutualism, shedding valuable insights into the metabolic underpinnings of this key functional group within methanogenic ecosystems.

Characterization and improvement of RNA-Seq precision in quantitative transcript expression profiling.

MOTIVATION: Measurement precision determines the power of any analysis to reliably identify significant signals, such as in screens for differential expression, independent of whether the experimental design incorporates replicates or not. With the compilation of large-scale RNA-Seq datasets with technical replicate samples, however, we can now, for the first time, perform a systematic analysis of the precision of expression level estimates from massively parallel sequencing technology. This then allows considerations for its improvement by computational or experimental means. RESULTS: We report on a comprehensive study of target identification and measurement precision, including their dependence on transcript expression levels, read depth and other parameters. In particular, an impressive recall of 84% of the estimated true transcript population could be achieved with 331 million 50 bp reads, with diminishing returns from longer read lengths and even less gains from increased sequencing depths. Most of the measurement power (75%) is spent on only 7% of the known transcriptome, however, making less strongly expressed transcripts harder to measure. Consequently, <30% of all transcripts could be quantified reliably with a relative error<20%. Based on established tools, we then introduce a new approach for mapping and analysing sequencing reads that yields substantially improved performance in gene expression profiling, increasing the number of transcripts that can reliably be quantified to over 40%. Extrapolations to higher sequencing depths highlight the need for efficient complementary steps. In discussion we outline possible experimental and computational strategies for further improvements in quantification precision. CONTACT: rnaseq10@boku.ac.at

Quantitative phosphoproteomics identifies filaggrin and other targets of ionizing radiation in a human skin model.

Our objective here was to perform a quantitative phosphoproteomic study on a reconstituted human skin tissue to identify low- and high-dose ionizing radiation-dependent signalling in a complex three-dimensional setting. Application of an isobaric labelling strategy using sham and three radiation doses (3, 10, 200 cGy) resulted in the identification of 1052 unique phosphopeptides. Statistical analyses identified 176 phosphopeptides showing significant changes in response to radiation and radiation dose. Proteins responsible for maintaining skin structural integrity including keratins and desmosomal proteins (desmoglein, desmoplakin, plakophilin 1, 2 and 3) had altered phosphorylation levels following exposure to both low and high doses of radiation. Altered phosphorylation of multiple sites in profilaggrin linker domains coincided with altered profilaggrin processing suggesting a role for linker phosphorylation in human profilaggrin regulation. These studies demonstrate that the reconstituted human skin system undergoes a coordinated response to both low and high doses of ionizing radiation involving multiple layers of the stratified epithelium that serve to maintain tissue integrity and mitigate effects of radiation exposure.

Synthetic Soil Aggregates: Bioprinted Habitats for High-Throughput Microbial Metaphenomics.

The dynamics of microbial processes are difficult to study in natural soil, owing to the small spatial scales on which microorganisms operate and to the opacity and chemical complexity of the soil habitat. To circumvent these challenges, we have created a 3D-bioprinted habitat that mimics aspects of natural soil aggregates while providing a chemically defined and translucent alternative culturing method for soil microorganisms. Our Synthetic Soil Aggregates (SSAs) retain the porosity, permeability, and patchy resource distribution of natural soil aggregates-parameters that are expected to influence emergent microbial community interactions. We demonstrate the printability and viability of several different microorganisms within SSAs and show how the SSAs can be integrated into a multi-omics workflow for single SSA resolution genomics, metabolomics, proteomics, lipidomics, and biogeochemical assays. We study the impact of the structured habitat on the distribution of a model co-culture microbial community and find that it is significantly different from the spatial organization of the same community in liquid culture, indicating a potential for SSAs to reproduce naturally occurring emergent community phenotypes. The SSAs have the potential as a tool to help researchers quantify microbial scale processes in situ and achieve high-resolution data from the interplay between environmental properties and microbial ecology.

Novel metabolic interactions and environmental conditions mediate the boreal peatmoss-cyanobacteria mutualism.

Interactions between Sphagnum (peat moss) and cyanobacteria play critical roles in terrestrial carbon and nitrogen cycling processes. Knowledge of the metabolites exchanged, the physiological processes involved, and the environmental conditions allowing the formation of symbiosis is important for a better understanding of the mechanisms underlying these interactions. In this study, we used a cross-feeding approach with spatially resolved metabolite profiling and metatranscriptomics to characterize the symbiosis between Sphagnum and Nostoc cyanobacteria. A pH gradient study revealed that the Sphagnum-Nostoc symbiosis was driven by pH, with mutualism occurring only at low pH. Metabolic cross-feeding studies along with spatially resolved matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) identified trehalose as the main carbohydrate source released by Sphagnum, which were depleted by Nostoc along with sulfur-containing choline-O-sulfate, taurine and sulfoacetate. In exchange, Nostoc increased exudation of purines and amino acids. Metatranscriptome analysis indicated that Sphagnum host defense was downregulated when in direct contact with the Nostoc symbiont, but not as a result of chemical contact alone. The observations in this study elucidated environmental, metabolic, and physiological underpinnings of the widespread plant-cyanobacterial symbioses with important implications for predicting carbon and nitrogen cycling in peatland ecosystems as well as the basis of general host-microbe interactions.