Search NASA⌕ Search

SEARCH · Search NASA

Results for “Functional Traits”

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 451 records · Page 25

Cyber–physical vulnerability and resiliency analysis for DER integration: A review, challenges and research needs

High penetration of renewable and sustainable Distributed Energy Resources (DER) into the traditional distribution system requires a well-coordinated control strategy for the improvement of system-wide reliability and resiliency. Implementation of such a holistic control architecture requires a flexible, near real-time, and bi-directional communication framework for facilitating the participation of various agents in a multi-vendor heterogeneous smart grid. While the sustainability of energy generation is ensured, this exposes the smart grid to extrinsic cyber threats, and appropriate defense mechanism(s) must be deployed to guarantee continued reliability and resiliency of the power grid. Further, the comprehensive literature review presented in this paper discusses the latest trends in the DER control schemes with fast communication requirements and their accompanying cyber–physical vulnerabilities. These control schemes are compared and contrasted for various traits. A three-level DER system architecture has been depicted, facilitating the deployment of these control schemes. The current developments of standard communication protocols, key security mechanisms, and best practices along major standards and guidelines are explored. The impacts of different attack types with miscellaneous DER functions based on various control schemes and associated mitigation solutions are also provided. Finally, challenges and future research directions for limiting cyber-power susceptibility to enhance resiliency are summarized. The work presented here will help us enabling a cyber-resilient and sustainable smart electric grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

State of Wildfires

Fire is an essential component of ecosystems and acts as key driver of biogeochemical cycling with impacts on vegetation structure and composition, soil conditions, and climate feedbacks. Fire behavior and effects vary based on the types of fuel, fuel dryness, and frequency of ignition. Due to changing climate and land use patterns, fire danger is increasing in many regions globally, and fires are having increasingly devastating impacts on human health, infrastructure, and ecosystem services. Recurrent fires help to determine the distribution of trees and grasses, and overall fuel load, which inform the behavior of future fires. Process-based models can be used to capture multi-scale impacts at the forest stand-level up to the landscape level, and across minutes to centuries, but must capture variation in fire behavior and intensity as a function of the fuel and climate, from high intensity forest crown fires within boreal regions to low intensity rapid grass fires of the tropics. Fire model development demonstrates our ability to capture large scale fire influenced biogeography and vegetation distribution at the earth system scale through vegetation traits and fire feedbacks. At the individual stand scale the importance of interactions and feedbacks between above and below ground process is essential to capturing the vegetation dynamics across boreal forest systems. Improving the mechanics of including plant physiology and specifically live fuel moisture content is the next step to advancing the capability of process based models to inform fire research. Advances in the testing and creation of a mechanistic live fuel moisture model demonstrate the foundation for future live fuel dynamics research that can be informed by field and remote sensing information. Improved remote sensing, technological and modeling capabilities support a more comprehensive and cohesive fire response that will be better equipped to overcome current barriers and anticipate the new reality of fires in a warming world.

wildfires↗

The maximum carboxylation rate of Rubisco affects CO2 refixation in temperate broadleaved forest trees

Mesophyll resistance to CO 2 diffusion (rm) and the maximum carboxylation rate of Rubisco (Vcmax) affect photosynthetic rates, and can potentially also influence the percentage of respiratory and photorespiratory CO 2 being refixated (Pr) by mesophyll cells. Here we investigated how various leaf anatomical traits (e.g. leaf mass per area [LMA] and leaf dry matter content [LDMC]) influenced rm in leaves of mature forest trees. We further explored how rm and Vcmax in turn affected Pr, and if these traits varied among species and leaves along a light gradient. Photosynthetic CO 2 response of leaves grown in high-, medium-, and low-light environments was measured, from Pinus sylvestris [Scots pine], Picea abies [Norway spruce], Quercus robur [English oak], and Betula pendula [Silver birch] in southern Sweden. A modified version of the Farquhar-von Caemmerer-Berry model was fitted to the leaf gas exchange data to estimate Vcmax, rm and Pr. We found that of all leaf traits measured, only LMA for Q. robur was significantly higher in leaves from high-light environments. When comparing species, both rm and LMA were significantly higher in the conifers, and rm had a negative correlation with Vcmax. We found that Pr was similar between different species and functional groups, with an average of 73.2% (and SD of ±10.4) across all species. There was a strong, positive correlation between Pr and Vcmax in broadleaves, and we hypothesise that this effect might derive from a higher CO 2 drawdown near Rubisco in leaves with high Vcmax.

59 BASIC BIOLOGICAL SCIENCES↗

A Field-Deployable Magnetic Resonance Imaging Rhizotron for Modeling and Enhancing Root Growth and Biogeochemical Function

A collaborative team from Texas A&M AgriLife Research, ABQMR Inc., the Soil Health Institute, the Athinoula A. Martinos Center for Biomedical Imaging, and NIST developed low-field magnetic resonance imaging (LF-MRI) instrumentation capable of imaging intact soil-root systems. The system measured root biomass, architecture, 3D mass distribution, and growth rates, providing a non-destructive means to evaluate ideal plant characteristics based on root metrics. It also successfully generated three-dimensional images of soil water content, a key property influencing root growth and exploration. Operating much like an MRI used in a medical setting, the system functioned in field conditions without damaging plants, overcoming the limitations of traditional methods such as trenching, soil coring, and root excavation. Over the course of the project, the team designed and built three functional prototype systems. These prototypes provided new insights into root–water–soil interactions that drive processes such as nutrient uptake, water use, and carbon management. This information contributed to efforts to optimize plants for carbon sequestration without sacrificing economic yield. The project also supported the identification of desirable traits for energy sorghum, including high root growth rates, more vertical root angles, and enhanced drought resilience under water-limiting conditions.

09 BIOMASS FUELS↗

Small and mighty: adaptation of superphylum Patescibacteria to groundwater environment drives their genome simplicity

The newly defined superphylum Patescibacteria such as Parcubacteria (OD1) and Microgenomates (OP11) has been found to be prevalent in groundwater, sediment, lake, and other aquifer environments. Recently increasing attention has been paid to this diverse superphylum including > 20 candidate phyla (a large part of the candidate phylum radiation, CPR) because it refreshed our view of the tree of life. However, adaptive traits contributing to its prevalence are still not well known. Results Here, we investigated the genomic features and metabolic pathways of Patescibacteria in groundwater through genome-resolved metagenomics analysis of >600 Gbp sequence data. We observed that, while the members of Patescibacteria have reduced genomes (~ 1 Mbp) exclusively, functions essential to growth and reproduction such as genetic information processing were retained. Surprisingly, they have sharply reduced redundant and nonessential functions, including specific metabolic activities and stress response systems. The Patescibacteria have ultra-small cells and simplified membrane structures, including flagellar assembly, transporters, and two-component systems. Despite the lack of CRISPR viral defense, the bacteria may evade predation through deletion of common membrane phage receptors and other alternative strategies, which may explain the low representation of prophage proteins in their genomes and lack of CRISPR. By establishing the linkages between bacterial features and the groundwater environmental conditions, our results provide important insights into the functions and evolution of this CPR group.

59 BASIC BIOLOGICAL SCIENCES↗

Machine Learning for the Validation of Expert-Elicited Causal Risk Diagrams

Exposure to spaceflight poses risk to human health in complex ways. To help manage this risk, the Human Systems Risk Board (HSRB) at the National Aeronautics and Space Administration (NASA) maintains a set of causal diagrams that attempt to explain how spaceflight hazards generate health risks and lead to adverse outcomes both in-mission, immediately post-mission, and over the long term. These causal risk diagrams are formulated as directed acyclic graphs (DAGs) and can function as knowledge graphs of connected risks and outcomes. These DAGs have proven useful for communication, and, through network analysis, have allowed for the identification of structurally important factors in the risk network. However, the utility these DAGs provide is directly proportional to their verisimilitude, making assessment of this trait using empirical data – whether from actual human spaceflight or various spaceflight analogue exposures and model organisms – a high priority. In this research we explore the use of machine learning algorithms to learn DAG structure from empirical data as a means of evaluating human-elicited DAG structures. To do so, we test several different graph structure-learning algorithms on data concerning changes in the bones of rats and mice after exposure to either spaceflight or a spaceflight analogue. We explore potential methods for indexing the similarity between each algorithm’s output DAG with all the others and with that of the expert-elicited DAG. We discuss next steps in this ongoing line of research and open science initiatives underway to complete them.

directed acyclic graphs↗

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins—including gene expression, chromatin occupancy and three-dimensional conformation—we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

59 BASIC BIOLOGICAL SCIENCES↗

Liana optical traits increase tropical forest albedo and reduce ecosystem productivity

Lianas are a key growth form in tropical forests. Their lack of self-supporting tissues and their vertical position on top of the canopy make them strong competitors of resources. A few pioneer studies have shown that liana optical traits differ on average from those of colocated trees. Those trait discrepancies were hypothesized to be responsible for the competitive advantage of lianas over trees. Yet, in the absence of reliable modelling tools, it is impossible to unravel their impact on the forest energy balance, light competition, and on the liana success in Neotropical forests. To bridge this gap, we performed a meta-analysis of the literature to gather all published liana leaf optical spectra, as well as all canopy spectra measured over different levels of liana infestation. We then used a Bayesian data assimilation framework applied to two radiative transfer models (RTMs) covering the leaf and canopy scales to derive tropical tree and liana trait distributions, which finally informed a full dynamic vegetation model. According to the RTMs inversion, lianas grew thinner, more horizontal leaves with lower pigment concentrations. Those traits made the lianas very efficient at light interception and significantly modified the forest energy balance and its carbon cycle. While forest albedo increased by 14% in the shortwave, light availability was reduced in the understorey (−30% of the PAR radiation) and soil temperature decreased by 0.5°C. Those liana-specific traits were also responsible for a significant reduction of tree (−19%) and ecosystem (−7%) gross primary productivity (GPP) while lianas benefited from them (their GPP increased by +27%). This study provides a novel mechanistic explanation to the increase in liana abundance, new evidence of the impact of lianas on forest functioning, and paves the way for the evaluation of the large-scale impacts of lianas on forest biogeochemical cycles.

Félicien Meunier↗

Functional succession of actively growing soil microorganisms during rewetting is shaped by precipitation history

Rewetting of seasonally dry soils induces a burst of microbial activity and carbon mineralization that changes nutrient availability and leads to succession. Yet the microbial functions that underpin this succession are not well described. Further, it’s unclear how previous precipitation frames microbial capacities after rewetting and how long these effects persist. We used isotopically-labeled water to rewet seasonally dry annual grassland soil that experienced either mean annual or reduced precipitation during the previous two years, and sampled at five subsequent time points. We used quantitative stable isotope probing (qSIP)-informed genome- resolved metagenomics to identify growing microorganisms, predict their capabilities, and analyze how these traits differed over time and between precipitation treatments. Organisms that grew after wetup were more abundant than non-growing organisms prior to the wet-up, suggesting that traits that initiate succession are pre-defined at the end of the prior plant growing season or via survival over the summer. Fast-growing organisms had fewer carbohydrate active enzyme (CAZy) genes per genome than slower-growing organisms, suggesting that although fast-growers were capable of degrading complex C, they may not specialize in this process. Differential abundance of CAZy genes in growing organisms throughout the succession implies that substrate availability varied with time. In contrast, the abundance of peptidases in growing organisms differed between precipitation treatments, but not over time following wet-up. Before wet-up, the soil organisms’ gene inventories were different between the two precipitation treatments. Surprisingly, this legacy effect waned after just one week. Thus, pre-wetup differences in microbial functional capacity converged shortly after rewetting.

Sieradzki, E T↗

CHESS 2025: Field-collected vegetation attributes and site photos

This dataset represents field observations of vegetation samples collected as part of the Colorado Headwaters Ecological Spectroscopy Study (CHESS) during June and July of 2025. Samples were collected in the field using tablet computers and digital forms, with target data differing by sample type (individual trees, individual shrubs, or 1-meter square plots of meadow and subshrub vegetation). Field samples were collected within 72 hours of airborne data collection using the National Ecological Observatory Network’s Aerial Observation Platform (NEON AOP). The NEON AOP collected waveform LiDAR (Light Detection and Ranging) and imaging spectrometer data in 426 spectral bands from the visible to shortwave infrared. Remote sensing data for the project is available on ESS-DIVE (DOI and citation to be added upon publication). Field data collected included canopy height and per-species horizontal proportional cover for meadow plots, species identity and height information for shrubs, as well as species identity, height, diameter at breast height, and health assessment information for trees. Photos of the focal site and surrounding landscape were taken for all sampling sites and are included in this archive. Green leaves or needles were collected for plant trait and foliar chemistry analysis. This data is archived separately (DOI and citation to be added upon publication). High-precision geospatial data for each sample (crown perimeter polygons for trees and shrubs, plot boundaries for meadow plots) is available here (Henderson et al., 2026). Field and remote sensing protocols largely followed those of a previous field and airborne imaging campaign performed in 2018 (described in Chadwick et al. 2020). Field data from the 2018 campaign can be found here (Chadwick et al., 2020 doi:10.15485/1618130). Because different field measurements were taken for meadow, shrub, and tree sites, data from these three sample types are archived as separate tables (chess_meadow_site_cleaned.csv, chess_shrub_site_cleaned.csv, chess_tree_site_cleaned.csv). Meadow proportional cover data is stored in a separate table (chess_meadow_cover_cleaned.csv). Taxonomy was treated identically between sample types, and the dataset shares a common set of voucher specimens (chess_voucher_IDs_cleaned.csv), as well as a single species list (chess_species_list_cleaned.csv). All taxonomic determinations were performed to the species level, and adhere to the Global Biodiversity Information Facility (GBIF) backbone taxonomy as of January 10th, 2026 (GBIF Secretariat 2023). CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgment: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Bioenergy cropping systems shape ant community composition and functional roles

R code and partial data for Haan, Helms, & Landis 2023, Bioenergy cropping systems shape ant community composition and functional roles, to be published in Frontiers in Conservation Science. The ant species list we use to generate the species matrix was uploaded previously (10.5281/zenodo.8215581) in association with another manuscript (Haan et al. 2023, Science Advances, Contrasting effects of bioenergy crops on biodiversity). Therefore here we include a small table with trait information for each species along with the R code used specifically for this manuscript.

ants↗

Origins of prokaryotes, eukaryotes, mitochondria, and chloroplasts

A computer branching model is used to analyze cellular evolution. Attention is given to certain key amino acids and nucleotide residues (ferredoxin, 5s ribosomal RNA, and c-type cytochromes) because of their commonality over a wide variety of cell types. Each amino acid or nucleotide residue is a sequence in an inherited biological trait; and the branching method is employed to align sequences so that changes reflect substitution of one residue for another. Based on the computer analysis, the symbiotic theory of cellular evolution is considered the most probable. This theory holds that organelles, e.g., mitochondria and chloroplasts invaded larger bodies, e.g., bacteria, and combined functions to form eucaryotic cells.

Schwartz, R. M.↗

Success Through Synergy (STS): Increasing Cultivation Yield and Stability with Rationally Designed Consortia

Consortia of functionally diverse microbial assemblages hold potential to enhance the production of microalgae for biofuels and bioproducts. In laboratory experiments, increases in algal yield have been achieved through the amendment of algal cultures with growth-promoting bacteria, as well as the rational design of algal-algal consortia (i.e., combining strains with complementary traits). Despite this evidence, commercial systems typically rely on monocultures because of the technical challenges that multi-taxa cultures impose. In this project, we aimed to generate algae-bacteria and algae-algae consortia in the lab and test the performance of these consortia against monoculture counterparts in the field. We also aimed to generate toolkits that could be used by the scientific community and algal farmers to advance the cultivation of consortia.

60 APPLIED LIFE SCIENCES↗

The Earliest Ion Channels in Protocellular Membranes

Cellular membranes with their hydrophobic interior are virtually impermeable to ions. Bulk of ion transport through them is enabled through ion channels. Ion channels of contemporary cells are complex protein molecules which span the membrane creating a cylindrical pore filled with water. Protocells, which are widely regarded as precursors to modern cells, had similarly impermeable membranes, but the set of proteins in their disposal was much simpler and more limited. We have been, therefore, exploring an idea that the first ion channels in protocellular membranes were formed by much smaller peptide molecules that could spontaneously selfassemble into short-lived cylindrical bundles in a membrane. Earlier studies have shown that a group of peptides known as peptaibols is capable of forming ion channels in lipid bilayers when they are exposed to an electric field. Peptaibols are small, non-genetically encoded peptides produced by some fungi as a part of their system of defense against bacteria. They are usually only 14-20 residues long, which is just enough to span the membrane. Their sequence is characterized by the presence of non-standard amino acids which, interestingly, are also expected to have existed on the early earth. In particular, the presence of 2-aminoisobutyric acid (AIB) gives peptaibols strong helix forming propensities. Association of the helices inside membranes leads to the formation of cylindrical bundles, typically containing 4 to 10 monomers. Although peptaibols are excellent candidates for models of the earliest ion channels their structures, which are stabilized only by van der Waals forces and occasional hydrogen bonds between neighboring helices, are not very stable. Although it might properly reflect protobiological reality, it is also a major obstacle in studying channel behavior. For this reason we focused on two members of the peptaibol family, trichotoxin and antiamoebin, which are characterized by a single conductance level. This indicates that their structures are unique and stable. In addition, it is also believed that the trichotoxin channel displays some selectivity between potassium and chloride ions. This makes trichotoxin and antiamoebin ideal models of the earliest ion channels that could provide insight into the origins of ion conductance and selectivity. In the absence of crystal structure of the trichotoxin and antiamoebin channels, we propose their molecular models based on experimentally determined number of monomers forming the bundles. We use molecular dynamics simulations to validate the models in terms of their conductance and selectivity. On the basis of our simulations we show that the emergence of channels built of small, alpha-helical peptides was protobiologically plausible and did not require highly specific amino acid sequences, which is a convenient evolutionary trait. Despite their simple structure, such channels could possess properties that, at the first sight, appear to require markedly larger complexity. To this end, we will discuss how the amino acid sequence and structure of primitive channels give rise to the phenomena of ionic conductance and selectivity across the earliest cell walls, which were essential functions for the emergence and early evolution of protocells. Furthermore, we will argue that even though architectures of membrane proteins are not nearly as diverse as those of water-soluble proteins, they are sufficiently flexible to adapt readily to the functional demands arising during evolution.

Mijajlovic, Milan↗

Multi-objective optimization of root phenotypes for nutrient capture using evolutionary algorithms

Root phenotypes are avenues to the development of crop cultivars with improved nutrient capture, which is an important goal for global agriculture. The fitness landscape of root phenotypes is highly complex and multidimensional. It is difficult to predict which combinations of traits (phene states) will create the best performing integrated phenotypes in various environments. Brute force methods to map the fitness landscape by simulating millions of phenotypes in multiple environments are computationally challenging. Evolutionary optimization algorithms may provide more efficient avenues to explore high dimensional domains such as the root phenotypic space. We coupled the three-dimensional functional–structural plant model, SimRoot, to the Borg Multi-Objective Evolutionary Algorithm (MOEA) and the evolutionary search over several generations facilitated the identification of optimal root phenotypes balancing trade-offs across nutrient uptake, biomass accumulation, and root carbon costs in environments varying in nutrient availability. Our results show that several combinations of root phenes generate optimal integrated phenotypes where performance in one objective comes at the cost of reduced performance in one or more of the remaining objectives, and such combinations differed for mobile and non-mobile nutrients and for maize (a monocot) and bean (a dicot). Functional–structural plant models can be used with multi-objective optimization to identify optimal root phenotypes under various environments, including future climate scenarios, which will be useful in developing the more resilient, efficient crops urgently needed in global agriculture.

59 BASIC BIOLOGICAL SCIENCES↗

The role of height–driven constraints and compensations on tree vulnerability to drought

Frequent observations of higher mortality in larger trees than in smaller ones during droughts have sparked an increasing interest in size–dependent drought–induced mortality. However, the underlying physiological mechanisms are not well understood, with height–associated hydraulic constraints often being implied as the potential mechanism driving increased drought vulnerability. Here we performed a quantitative synthesis on how key traits that drive plant water and carbon economy change with tree height within species and assessed the implications that the different constraints and compensations may have on the interacting mechanisms (hydraulic failure, carbon starvation and/or biotic–agent attacks) affecting tree vulnerability to drought. While xylem tension increases with tree height, taller trees present a range of structural and functional adjustments, including more efficient water use and transport and greater water uptake and storage capacity, that mitigate the path–length–associated drop in water potential. These adaptations allow taller trees to withstand episodic water stress. Conclusive evidence for height–dependent increased vulnerability to hydraulic failure and carbon starvation, and their coupling to defence mechanisms and pest and pathogen dynamics, is still lacking. Further research is needed, particularly at the intraspecific level, to ascertain the specific conditions and thresholds above which height hinders tree survival under drought.

54 ENVIRONMENTAL SCIENCES↗

Characterizing Natural Organic Matter Transformations by Microbial Communities in Terrestrial Subsurface Ecosystems: A Critical Review of Analytical Techniques and Challenges

Determining the mechanisms, traits, and pathways that regulate microbial transformation of natural organic matter (NOM) is critical to informing our understanding of the microbial impacts on the global carbon cycle. The capillary fringe of subsurface soils is a highly dynamic environment that remains poorly understood. Characterization of organo-mineral chemistry combined with a nuanced understanding of microbial community composition and function is necessary to understand microbial impacts on NOM speciation in the capillary fringe. We present a critical review of the popular analytical and omics techniques used for characterizing complex carbon transformation by microbial communities and focus on how complementary information obtained from the different techniques enable us to connect chemical signatures with microbial genes and pathways. This holistic approach offers a way forward for the comprehensive characterization of the formation, transformation, and mineralization of terrestrial NOM as influenced by microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Long-reads↗