Engineering Papers⌕ Search
Home/Archive/DOE OSTI

DOE OSTI · 3388428

wavess 1.2: presenting an HLA-aware within-host virus sequence simulation framework

Lapp, Zena Malakie [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000346742176)·Leitner, Thomas Kenneth [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000181602588)
View original source ↗Read source document ↗Export metadata ↓

Abstract

Motivation Understanding how virus sequences are shaped by selection can inform vaccine design and transmission inference. Modeling within-host evolution to interrogate these questions requires a detailed mechanistic framework that accurately captures sequence diversification. The CD8 + cytotoxic T-lymphocyte (CTL) response plays an important role in immune-mediated selection and can leave strong signatures in virus sequences; however, existing sequence-based within-host virus modeling frameworks do not explicitly include a human leukocyte antigen (HLA)-aware CTL response. Results We extended our previously published within-host sequence evolution simulator, wavess, to include an explicit CTL response, and share a method for identifying HLA-specific CTL epitopes given a founder virus sequence. We also updated the model to permit a variable recombination rate, which allows for modeling non-adjacent genes, segmented genomes, and recombination hotspots. These extensions to wavess allow for more accurate simulation of viruses and virus genes, particularly in regions of the genome where the immune response is dominated by CTLs (rather than antibodies). It also provides the foundation for investigations of how these newly-added biological mechanisms influence within-host evolution. Availability and implementation The core of wavess is written in Python 3, with helper functions written in R. It is available at https://github.com/MolEvolEpid/wavess.

Explore related subjects

60 APPLIED LIFE SCIENCES

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Lapp, Zena Malakie [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000346742176), Leitner, Thomas Kenneth [Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)] (ORCID:0000000181602588). 2026-06-18. wavess 1.2: presenting an HLA-aware within-host virus sequence simulation framework. https://doi.org/10.1093/bioadv%2Fvbag174

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

DOE OSTI2026

Metabolic flux and resource balance in the oleaginous yeast Rhodotorula toruloides

The yeast Rhodotorula toruloides is a promising bioproduction organism due to its high lipid yields and ability to grow on cheap and abundant substrates. Quantitative, systems-level assessment of its metabolic activity is accordingly merited. Resource-balance analysis (RBA) models capture not only reaction stoichiometry but also enzyme requirements for catalysis, providing valuable tools for understanding metabolic trade-offs and optimizing metabolic engineering strategies. Here, in this work, we present systems-level measurements of R. toruloides metabolic flux based on isotope tracing and metabolic flux analysis. In combination with new proteomic measurements, these flux data are used to parameterize a genome-scale resource balance model rtRBA. We find that S. cerevisiae and R. toruloides grow at nearly indistinguishable rates using similar biosynthetic but dramatically different central metabolic programs. R. toruloides consumes one-fifth as much glucose, which it metabolizes primarily via the pentose phosphate pathway and TCA cycle unlike primarily glycolysis in S. cerevisiae . Overall, across these two divergent yeasts, protein abundances aligned more closely than metabolic flux. Resource balance modeling of these metabolic programs predicts superior theoretical yields but lower productivities in R. toruloides than S. cerevisiae for industrial chemicals, highlighting the value of rapid glucose uptake for productivity but respiratory metabolism for yields.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Opportunities for bioenergy crops to support transitions from irrigated agriculture and conserve the U.S. High Plains Aquifer

This study investigates potential economic and groundwater driven transitions from irrigated maize production—the dominant irrigated cropping system in the High Plains Aquifer (HPA) region—to alternative crops such as sorghum and switchgrass, two common bioenergy feedstocks. Unsustainable groundwater extraction in the U.S. High Plains presents critical challenges including reduced irrigation capacities, diminished crop yields, lower land values, and escalating energy costs. Using a spatially explicit optimization framework combined with comprehensive economic and hydrogeological data, we evaluate optimal land-use strategies and irrigation system investments over a 30-year planning period. Results indicate significant regional variability in the future economic viability of irrigated agriculture, driven by differences in aquifer recharge rates, groundwater availability, and market conditions. Nebraska and parts of northern Texas can sustain irrigated maize profitability due to relatively favorable groundwater conditions and lower land rents, respectively. In contrast, many portions of Kansas and southern Texas are more likely to transition to dryland agriculture within two decades. Colorado and New Mexico show potential for significant adoption of switchgrass production as an alternative biomass-based energy crop. Overall, over the 30-year horizon, our model implies that approximately 23% of currently irrigated maize area across the HPA region may transition to dryland farming under business-as-usual conditions. This transition is complemented by a threefold increase in non-irrigated sorghum production, from 0.20 to 0.77 Mt yr−1, indicating that groundwater-driven shifts in agricultural production may support a larger regional base for bioenergy feedstocks. The study also reveals opportunities for producers to optimize economic returns and biomass production potential.

60 APPLIED LIFE SCIENCES↗