Search NASA⌕ Search

SEARCH · Search NASA

Results for “data replication”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

Low-Temperature Geothermal Play Fairway Analysis for the Denver Basin

This project is part of a national initiative to showcase the benefits of incorporating low-temperature geothermal resource assessment into the deployment of geothermal heating, combined heat and power (CHP), and geothermal direct-use (GDU) technologies. The initiative was established to accelerate the country's decarbonization efforts by identifying potential for low-temperature geothermal resource utilization (<150 degrees C, e.g., CHP and GDU) in selected sedimentary basins with numerous population centers. The play fairway analysis (PFA) methodologies in this study were adapted from previous PFA investigations of sedimentary basin geothermal play types (SBGPTs) that evaluated the potential for low-temperature resources (<150 degrees C). Workflows, relevant datasets, a new Python library, and common and composite geological criteria maps are utilized to develop low-temperature geothermal resource favorability maps for the Denver Basin, a sedimentary basin spanning Colorado, Nebraska, and Wyoming. The replication of these methodologies in other SBGPTs can evaluate potential for low-temperature resources. To facilitate future assessment of low-temperature geothermal resources in SBGPTs, this project provides PFA workflows, data, tools, and favorability maps that will ultimately support the utilization of low-temperature geothermal resources in sedimentary basins.

combined heat and power↗

Next-Generation Sequencing Data from a CUT&RUN Study of R. toruloides IFO0880 Cse4 and Orc1 Binding Sites

Rhodotorula toruloides has been increasingly explored as a host for bioproduction of lipids, fatty acid derivatives and terpenoids. Various genetic tools have been developed, but neither a centromere nor an autonomously replicating sequence (ARS), both necessary elements for stable episomal plasmid maintenance, has yet been reported. In this study, cleavage under targets and release using nuclease (CUT&RUN), a method used for genome-wide mapping of DNA–protein interactions, was used to identify R. toruloides IFO0880 genomic regions associated with the centromeric histone H3 protein Cse4, a marker of centromeric DNA. Fifteen putative centromeres ranging from 8 to 19 kb in length were identified and analyzed, and four were tested for, but did not show, ARS activity. These centromeric sequences contained below average GC content, corresponded to transcriptional cold spots, were primarily nonrepetitive and shared some vestigial transposon-related sequences but otherwise did not show significant sequence conservation. Future efforts to identify an ARS in this yeast can utilize these centromeric DNA sequences to improve the stability of episomal plasmids derived from putative ARS elements.

Genome Engineering↗

Tulane virus protease as a structural surrogate for inhibitor screening of human norovirus proteases

Human norovirus (HuNoV) is a significant cause of gastroenteritis worldwide, affecting people of all age groups. There are currently no vaccines or drugs available, leaving susceptible populations vulnerable to severe or protracted illness. A HuNoV cultivation system is pivotal for screening norovirus antivirals. While the human intestinal enteroid cultivation system allows robust replication of multiple HuNoV strains, it presents technical and cost barriers. Tulane virus (TV), a surrogate for HuNoV, replicates well in monkey kidney cell lines and is closely related to norovirus in cellular biology. Here, we determined the structures of TV protease (TV-Pro) alone and in complex with rupintrivir, a picornavirus inhibitor that also inhibits HuNoV proteases (HuNoV-Pro). Our data validate TV as an efficient surrogate system for rapid screening of HuNoV protease inhibitors. The TV protease structure exhibits significant backbone similarity to the GI.1 HuNoV protease in the substrate-binding domain, with the BII-CII loop in an open conformation stabilized by hydrogen bonds as present in the GI.1 protease. Structural differences in the S2 pocket and two amino acid changes in the S4 pocket result in slightly altered P2 and P4 substrate and inhibitor conformations. Despite these differences, we confirm previous findings that the TV protease can cleave the GI.1 and GII HuNoV polyprotein substrates with high and moderate efficiency, respectively. We found that rupintrivir efficiently inhibits TV protease in vitro and inhibits TV replication in cell culture with similar efficacy in combination with P-glycoprotein efflux pump inhibitors. We conclude that TV is a valuable surrogate for HuNoV protease inhibitor screening and outline strategies to improve its compatibility as such.

crystal structures↗

Codon usage revisited: Lack of correlation between codon usage and the number of tRNA genes in enterobacteria

Highlights: • Average number of tRNA genes correlate with codon usage in highly expressed genes. • Changes in the number of tRNA genes do not correlate to codon usage. • Altered tRNA gene copy numbers do not correlate to changes of codon translation. • Our data indicates that tRNA gene copies adapt to codon usage, but not the contrary. It is widely believed that if a high number of genes are found for any tRNA in a rapidly replicating bacteria, then the cytoplasmic levels of that tRNA will be high and an open reading frame containing a higher frequency of the complementary codon will be translated faster. This idea is based on correlations between the number of tRNA genes, tRNA concentration and the frequency of codon usage observed in a limited number of strains as well as from the fact that artificially changing the number of tRNA genes alters translation efficiency and consequently the amount of properly folded protein synthesized. tRNA gene number may greatly vary in a genome due to duplications, deletions and lateral transfer which in turn would alter the levels and functionality of many proteins. Such changes are potentially deleterious for fitness and as a result it is expected that changes in tRNA gene numbers should be accompanied by a modification of the frequency of codon usage. In contrast to this model, when comparing the number of tRNA genes and the frequency of codon usage of several Salmonella enterica and Escherichia coli strains we found that changes in the number of tRNA genes are not correlated to changes in codon usage. Furthermore, these changes are not correlated with a change in the efficiency of codon translation. These results suggest that once a genome gains or loses tRNA genes, it responds by modulating the concentrations of tRNAs rather than modifying its frequency of codon usage.

60 APPLIED LIFE SCIENCES↗

The kinetics of SARS-CoV-2 infection based on a human challenge study

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.

59 BASIC BIOLOGICAL SCIENCES↗

Generating Co-expression Networks for Three Cyanobacteria: Synechococcus sp. PCC 7942, Synechococcus sp. PCC 7002, Synechocystis sp. PCC 6803

Cyanobacteria are photosynthetic organisms capable of high growth rate and represent a promising bioplatform for harnessing the sun’s energy to make biofuel. Additionally, the process of photosynthesis absorbs CO2 from the environment. Understanding the metabolic processes involved in photosynthesis could lead to solutions to the recent rise of CO2 concentration in Earth’s atmosphere and the associated climate change. More research on the transcriptional regulation of these cells is needed to learn how to harness the untapped potential of cyanobacteria for these applications. Transcriptional analysis via RNA-seq provides an understanding of how gene expression changes at the mRNA level under diverse growing conditions. I systematically collected and analyzed RNA-Seq data obtained under a variety of conditions and available on the NCBI database for three cyanobacteria model organisms: Synechococcus elongatus sp. PCC 7942, Synechococcus sp. PCC 7002, and Synechocystis sp. PCC 6803. For each organism, the data was mapped to a reference genome to characterize the RNA expression profile. Samples were checked for quality based on the number of reads and the correlation of the expression profile between labeled replicates. All samples were transformed into transcripts per million reads, followed by a log transformation to account for the wide range of sample sizes. Gene co-expression networks were generated and analyzed for each species using Cytoscape. These networks provide a base level of gene expression for each species. The network topology and high-betweeness nodes of these networks need to be analyzed further to provide insight on potential ways to harness cyanobacteria genetics. Additionally, these datasets can be used together to form a core genome network analysis- one that includes only the genes that are homologous between the three species. This project has prepared the way for a more in-depth study on photosynthetic microbes on a genetic level.

59 BASIC BIOLOGICAL SCIENCES↗

TH-E-17A-07: Improved Cine Four-Dimensional Computed Tomography (4D CT) Acquisition and Processing Method

Purpose: Artifacts arising from the 4D CT acquisition and post-processing methods add systematic uncertainty to the treatment planning process. We propose an alternate cine 4D CT acquisition and post-processing method to consistently reduce artifacts, and explore patient parameters indicative of image quality. Methods: In an IRB-approved protocol, 18 patients with primary thoracic malignancies received a standard cine 4D CT acquisition followed by an oversampling 4D CT that doubled the number of images acquired. A second cohort of 10 patients received the clinical 4D CT plus 3 oversampling scans for intra-fraction reproducibility. The clinical acquisitions were processed by the standard phase sorting method. The oversampling acquisitions were processed using Dijkstras algorithm to optimize an artifact metric over available image data. Image quality was evaluated with a one-way mixed ANOVA model using a correlation-based artifact metric calculated from the final 4D CT image sets. Spearman correlations and a linear mixed model tested the association between breathing parameters, patient characteristics, and image quality. Results: The oversampling 4D CT scans reduced artifact presence significantly by 27% and 28%, for the first cohort and second cohort respectively. From cohort 2, the inter-replicate deviation for the oversampling method was within approximately 13% of the cross scan average at the 0.05 significance level. Artifact presence for both clinical and oversampling methods was significantly correlated with breathing period (ρ=0.407, p-value<0.032 clinical, ρ=0.296, p-value<0.041 oversampling). Artifact presence in the oversampling method was significantly correlated with amount of data acquired, (ρ=-0.335, p-value<0.02) indicating decreased artifact presence with increased breathing cycles per scan location. Conclusion: The 4D CT oversampling acquisition with optimized sorting reduced artifact presence significantly and reproducibly compared to the phase-sorted clinical acquisition.

60 APPLIED LIFE SCIENCES↗

Validation Study of the Spacecraft Background Activation using RHESSI Data

This document serves as a brief writeup of the work done to validate CINDER2008 and the Activation Analysis Tools in the Activation in Accelerator Radiation Environments (AARE) distribution that are used to calculate the activation of spacecraft material while in orbit [1{3]. The validation study involves replicating the decay gamma spectrum as observed onboard the RHESSI spacecraft in the high-purity germanium (HPGe) detectors after crossing the South Atlantic Anomaly (SAA) [4]. The tools involved are a suite of Python scripts, the transmutation code CINDER2008 accompanied with a suite of scripts distributed with AARE,MCNPX, and MCNP6.2 [2,5{7]. Figure 1 shows the gamma ray spectrum this study replicates as the spectrum designated \Difference".

42 ENGINEERING↗

Data-Driven Residential Load Modeling and Validation in GridLAB-D

Accurately characterizing the impacts of high penetrations of distributed energy resources (DER) on the electric distribution system has driven modeling methods from traditional static snap shots, often representing a critical point in time (e.g., summer peak load), to quasi-static time series (QSTS) simulations capturing all the effects of variable DER, associated controls and hence, impacts on the distribution system over a given time period. Unfortunately, the high time resolution DER source and load data required for model inputs is often scarce or non-existent. This paper presents work performed within the GridLAB-D model environment to synthesize, calibrate, and validate 1-second residential load models based on measured transformer loads and physics-based models suitable for QSTS electric distribution system modeling. The modeling and validation approach taken was to create a typical GridLAB-D model home that, when replicated to represent multiple diverse houses on a single transformer, creates a statistically similar load to a measured load for a given weather input. The model homes are constructed to represent the range of actual homes on an instrumented transformer: square footage, thermal integrity, heating and cooling system definition as well as realistic occupancy schedules. House model calibration and validation was performed using the distribution transformer load data and corresponding weather. The modeled loads were found to be similar to the measured loads for four evaluation metrics: 1) daily average energy, 2) daily average and standard deviation of power, 3) power spectral density, and 4) load shape.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Role of a nuclear localization signal on the minor capsid Proteins VP2 and VP3 in BKPyV nuclear entry

BK Polyomavirus (BKPyV) is a ubiquitous nonenveloped human virus that can cause severe disease in immunocompromised populations. After internalization into renal proximal tubule epithelial cells, BKPyV traffics through the ER and enters the cytosol. However, it is unclear how the virus enters the nucleus. In this study, we elucidate a role for the nuclear localization signal located on the minor capsid proteins VP2 and VP3 during infection. Site-directed mutagenesis of a single lysine in the basic region of the C-terminus of the minor capsid proteins abrogated their nuclear localization, and the analogous genomic mutation reduced infectivity. Additionally, through use of the inhibitor ivermectin and knockdown of importin β1, we found that the importin α/β pathway is involved during infection. Overall these data are the first to show the significance of the NLS of the BKPyV minor capsid proteins during infection in a natural host cell. - Highlights: • Polyomaviruses must deliver their genome to the nucleus to replicate. • The minor capsid proteins have a well-conserved nuclear localization signal. • Mutation of this NLS diminishes, but does not completely inhibit, infection.

60 APPLIED LIFE SCIENCES↗

Intratracheal exposure of common marmosets to MERS-CoV Jordan-n3/2012 or MERS-CoV EMC/2012 isolates does not result in lethal disease

Middle East Respiratory Syndrome Coronavirus (MERS-CoV) continues to be a threat to human health in the Middle East. Development of countermeasures is ongoing; however, an animal model that faithfully recapitulates human disease has yet to be defined. A recent study indicated that inoculation of common marmosets resulted in inconsistent lethality. Based on these data we sought to compare two isolates of MERS-CoV. We followed disease progression in common marmosets after intratracheal exposure with: MERS-CoV-EMC/2012, MERS-CoV-Jordan-n3/2012, media, or inactivated virus. Our data suggest that common marmosets developed a mild to moderate non-lethal respiratory disease, which was quantifiable by computed tomography (CT), with limited other clinical signs. Based on CT data, clinical data, and virological data, MERS-CoV inoculation of common marmosets results in mild to moderate clinical signs of disease that are likely due to manipulations of the marmoset rather than as a result of robust viral replication. - Highlights: • Common marmosets infected with MERS-EMC and MERS-JOR did not develop lethal disease. • Infected subjects developed transient signs of clinical disease. • CT indicated few differences between the infected and control groups. • Marmosets do not faithfully replicate human MERS pathogenesis.

60 APPLIED LIFE SCIENCES↗

SARS-CoV-2 Diverges from Other Betacoronaviruses in Only Partially Activating the IRE1α/XBP1 Endoplasmic Reticulum Stress Pathway in Human Lung-Derived Cells

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed over 6 million individuals worldwide and continues to spread in countries where vaccines are not yet widely available or its citizens are hesitant to become vaccinated. Therefore, it is critical to unravel the molecular mechanisms that allow SARS-CoV-2 and other coronaviruses to infect and overtake the host machinery of human cells. Coronavirus replication triggers endoplasmic reticulum (ER) stress and activation of the unfolded protein response (UPR), a key host cell pathway widely believed to be essential for viral replication. We examined the master UPR sensor IRE1α kinase/RNase and its downstream transcription factor effector XBP1s, which is processed through an IRE1α-mediated mRNA splicing event, in human lung-derived cells infected with betacoronaviruses. We found that human respiratory coronavirus OC43 (HCoV-OC43), Middle East respiratory syndrome coronavirus (MERS-CoV), and murine coronavirus (MHV) all induce ER stress and strongly trigger the kinase and RNase activities of IRE1α as well as XBP1 splicing. In contrast, SARS-CoV-2 only partially activates IRE1α through autophosphorylation, but its RNase activity fails to splice XBP1. Moreover, while IRE1α was dispensable for replication in human cells for all coronaviruses tested, it was required for maximal expression of genes associated with several key cellular functions, including the interferon signaling pathway, during SARS-CoV-2 infection. Our data suggest that SARS-CoV-2 actively inhibits the RNase of autophosphorylated IRE1α, perhaps as a strategy to eliminate detection by the host immune system.

60 APPLIED LIFE SCIENCES↗

Correlating transcription and protein expression profiles of immune biomarkers following lipopolysaccharide exposure in lung epithelial cells

Universal and early recognition of pathogens occurs through recognition of evolutionarily conserved pathogen associated molecular patterns (PAMPs) by innate immune receptors and the consequent secretion of cytokines and chemokines. The intrinsic complexity of innate immune signaling and associated signal transduction challenges our ability to obtain physiologically relevant, reproducible and accurate data from experimental systems. One of the reasons for the discrepancy in observed data is the choice of measurement strategy. Immune signaling is regulated by the interplay between pathogen-derived molecules with host cells resulting in cellular expression changes. However, these cellular processes are often studied by the independent assessment of either the transcriptome or the proteome. Correlation between transcription and protein analysis is lacking in a variety of studies. In order to methodically evaluate the correlation between transcription and protein expression profiles associated with innate immune signaling, we measured cytokine and chemokine levels following exposure of human cells to the PAMP lipopolysaccharide (LPS) from the Gram-negative pathogen Pseudomonas aeruginosa . Expression of 84 messenger RNA (mRNA) transcripts and 69 proteins, including 35 overlapping targets, were measured in human lung epithelial cells. We evaluated 50 biological replicates to determine reproducibility of outcomes. Following pairwise normalization, 16 mRNA transcripts and 6 proteins were significantly upregulated following LPS exposure, while only five (CCL2, CSF3, CXCL5, CXCL8/IL8, and IL6) were upregulated in both transcriptomic and proteomic analysis. This lack of correlation between transcription and protein expression data may contribute to the discrepancy in the immune profiles reported in various studies. The use of multiomic assessments to achieve a systems-level understanding of immune signaling processes can result in the identification of host biomarker profiles for a variety of infectious diseases and facilitate countermeasure design and development.

59 BASIC BIOLOGICAL SCIENCES↗

Data from: An updated gene atlas for maize reveals organ-specific and stress-induced genes

Maize(Zea mays L.), a model species for genetic studies, is one of the two most important crop species worldwide.The genome sequence of the reference genotype, B73, representative of the stiff stalk heterotic group was recently updated (AGPv4) using long‐read sequencing and optical mapping technology. To facilitate the use of AGPv4 and to enable functional genomic studies and association of genotype with phenotype, we determined expression abundances for replicated mRNA‐sequencing datasets from 79 tissues and five abiotic/biotic stress treatments revealing 36,207 expressed genes. Characterization of the B73 transcriptome across six organs revealed 4,154 organ‐specific and 7,704 differentially expressed genes following stress treatment. Gene co‐expression network analyses revealed 12 modules associated with distinct biological processes containing 13,590 genes providing a resource for further association of gene function based on co‐expression patterns. Presence/absence variants (PAVs) previously identified using whole genome resequencing data from 61 additional inbred lines were enriched in organ‐specific and stress‐induced differentially expressed genes, suggesting PAVs may function in phenological variation and adaptation to environment. Relative to core genes conserved across the 62 profiled inbreds, PAVs have lower expression abundances which are correlated with their frequency of dispersion across inbreds and on average, have significantly fewer co‐expression network connections suggesting that a subset of PAVs may be on an evolutionary path to pseudogenization. To facilitate use by the community, we developed the Maize Genomics Resource website(maize.plantbiology.msu.edu) for viewing and data‐mining these resources and deployed two new views on the Maize electronic Fluorescent Pictograph Browser (bar.utoronto. ca/efp_maize)

59 BASIC BIOLOGICAL SCIENCES↗

EM structure of a helicase-loader complex depicting a 6:2 binding sub-stoichiometry from Geobacillus kaustophilus HTA426

During DNA replication, bacterial helicase is recruited as a complex in association with loader proteins to unwind the parental duplex. Previous structural studies have reported saturated 6:6 helicase-loader complexes with different conformations. However, structural information on the sub-stoichiometric conformations of these previously-documented helicase-loader complexes remains elusive. Here, with the aid of single particle electron-microscopy (EM) image reconstruction, we present the Geobacillus kaustophilus HTA426 helicase-loader (DnaC-DnaI) complex with a 6:2 binding stoichiometry in the presence of ATPγS. In the 19 Å resolution EM map, the undistorted and unopened helicase ring holds a robust loader density above the C-terminal RecA-like domain. Meanwhile, the path of the central DNA binding channel appears to be obstructed by the reconstructed loader density, implying its potential role as a checkpoint conformation to prevent the loading of immature complex onto DNA. Our data also reveals that the bound nucleotides and the consequently induced conformational changes in the helicase hexamer are essential for active association with loader proteins. These observations provide fundamental insights into the formation of the helicase-loader complex in bacteria that regulates the DNA replication process. - Highlights: • Helicase-loader complex structure with 6:2 sub-stoichiometry is resolved by EM. • Helicase hexamer in 6:2 sub-stoichiometry is constricted and un-opened. • 6:2 binding ratio of helicase-loader complex could act as a DNA loading checkpoint. • Nucleotides stabilize helicase-loader complex at low protein concentrations.

60 APPLIED LIFE SCIENCES↗

Single-molecule epitranscriptomic analysis of full-length HIV-1 RNAs reveals functional roles of site-specific m6As

Abstract Although the significance of chemical modifications on RNA is acknowledged, the evolutionary benefits and specific roles in human immunodeficiency virus (HIV-1) replication remain elusive. Most studies have provided only population-averaged values of modifications for fragmented RNAs at low resolution and have relied on indirect analyses of phenotypic effects by perturbing host effectors. Here we analysed chemical modifications on HIV-1 RNAs at the full-length, single RNA level and nucleotide resolution using direct RNA sequencing methods. Our data reveal an unexpectedly simple HIV-1 modification landscape, highlighting three predominant N 6 -methyladenosine (m 6 A) modifications near the 3′ end. More densely installed in spliced viral messenger RNAs than in genomic RNAs, these m 6 As play a crucial role in maintaining normal levels of HIV-1 RNA splicing and translation. HIV-1 generates diverse RNA subspecies with distinct m 6 A ensembles, and maintaining multiple of these m 6 As on its RNAs provides additional stability and resilience to HIV-1 replication, suggesting an unexplored viral RNA-level evolutionary strategy.

60 APPLIED LIFE SCIENCES↗

Results from a multi-laboratory ocean metaproteomic intercomparison: effects of LC-MS acquisition and data analysis procedures

Metaproteomics is an increasingly popular methodology that provides information regarding the metabolic functions of specific microbial taxa and has potential for contributing to ocean ecology and biogeochemical studies. A blinded multi-laboratory intercomparison was conducted to assess comparability and reproducibility of taxonomic and functional results and their sensitivity to methodological variables. Euphotic zone samples from the Bermuda Atlantic Time-series Study (BATS) in the North Atlantic Ocean collected by in situ pumps and the autonomous underwater vehicle (AUV) Clio were distributed with a paired metagenome, and one-dimensional (1D) liquid chromatographic data-dependent acquisition mass spectrometry analysis was stipulated. Analysis of mass spectra from seven laboratories through a common bioinformatic pipeline identified a shared set of 1056 proteins from 1395 shared peptide constituents. Quantitative analyses showed good reproducibility: pairwise regressions of spectral counts between laboratories yielded R 2 values averaged 0.62±0.11, and a Sørensen similarity analysis of the top 1000 proteins revealed 70 %–80 % similarity between laboratory groups. Taxonomic and functional assignments showed good coherence between technical replicates and different laboratories. A bioinformatic intercomparison study, involving 10 laboratories using eight software packages, successfully identified thousands of peptides within the complex metaproteomic datasets, demonstrating the utility of these software tools for ocean metaproteomic research. Lessons learned and potential improvements in methods were described. Future efforts could examine reproducibility in deeper metaproteomes, examine accuracy in targeted absolute quantitation analyses, and develop standards for data output formats to improve data interoperability. Together, these results demonstrate the reproducibility of metaproteomic analyses and their suitability for microbial oceanography research, including integration into global-scale ocean surveys and ocean biogeochemical models.

59 BASIC BIOLOGICAL SCIENCES↗

Phosphate amendment drives bloom of RNA viruses after soil wet-up

Soil rewetting after a dry period results in a surge of activity and succession in both microbial and DNA virus communities. Less is known about the response of RNA viruses to soil rewetting—while they are highly diverse and widely distributed in soil, they remain understudied. We hypothesized that RNA viruses would show temporal succession following rewetting and that phosphate amendment would influence their trajectory, as viral proliferation may cause phosphorus limitation. Using 39 time-resolved metatranscriptomes and amplicon data, 2190 RNA viral populations were identified across five phyla, with 26 % of these predicted to infect bacteria, and 11 % fungi. Only 1.2 % of viral populations had annotated capsid genes, suggesting most persist via intracellular replication without a free virion phase. Phosphate amendment altered RNA viral community composition within the first week and amended vs. unamended communities remained distinguishable for up to three weeks. While the overall host community remained stable, certain bacterial populations showed reduced abundance in phosphate-amended soils, likely due to increased viral lysis, as RNA bacteriophages proliferated significantly. Notably, 60 % of the viruses with increased abundance under phosphate amendment belonged to basal Lenarviricota clades rather than well-known groups like Leviviricetes. We estimate RNA bacteriophage infections may affect 10 7 –10 9 bacteria per gram of soil, aligning with the total bacterial population (10 7 –10 10 g -1 soil), suggesting that RNA phages significantly influence bacterial communities post-wet-up, with phosphorus availability modulating this effect.

59 BASIC BIOLOGICAL SCIENCES↗