Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time sequence analysis”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Permafrost thaw with warming reduces microbial metabolic capacities in subsurface soils

Abstract Microorganisms are major constituents of the total biomass in permafrost regions, whose underlain soils are frozen for at least two consecutive years. To understand potential microbial responses to climate change, here we examined microbial community compositions and functional capacities across four soil depths in an Alaska tundra site. We showed that a 5‐year warming treatment increased soil thaw depth by 25.7% ( p = .011) within the deep organic layer (15–25 cm). Concurrently, warming reduced 37% of bacterial abundance and 64% of fungal abundances in the deep organic layer, while it did not affect microbial abundance in other soil layers (i.e., 0–5, 5–15, and 45–55 cm). Warming treatment altered fungal community composition and microbial functional structure ( p < .050), but not bacterial community composition. Using a functional gene array, we found that the relative abundances of a variety of carbon (C)‐decomposing, iron‐reducing, and sulphate‐reducing genes in the deep organic layer were decreased, which was not observed by the shotgun sequencing‐based metagenomics analysis of those samples. To explain the reduced metabolic capacities, we found that warming treatment elicited higher deterministic environmental filtering, which could be linked to water‐saturated time, soil moisture, and soil thaw duration. In contrast, plant factors showed little influence on microbial communities in subsurface soils below 15 cm, despite a 25.2% higher ( p < .05) aboveground plant biomass by warming treatment. Collectively, we demonstrate that microbial metabolic capacities in subsurface soils are reduced, probably arising from enhanced thaw by warming.

Wu, Linwei↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Simple, Cost-Effective, and Automation-Friendly Direct PCR Approach for Bacterial Community Analysis

Understanding bacterial interactions and assembly in complex microbial communities using 16S rRNA sequencing normally requires a large experimental load. However, the current DNA extraction methods, including cell disruption and genomic DNA purification, are normally biased, costly, time-consuming, labor-intensive, and not amenable to miniaturization by droplets or 1,536-well plates due to the significant DNA loss during the purification step for tiny-volume and low-cell-density samples.

16S rRNA sequencing↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2022

Optimal Control of Biomass Feedstock Processing System Under Uncertainty in Biomass Quality

Planning of biorefinery operations is complicated by the stochastic nature of physical and chemical characteristics of biomass feedstock, such as, moisture level and carbohydrate content. Biomass characteristics affect the performance of the equipment which feed the reactor and the efficiency of the conversion process in a biorefinery. We propose a stochastic optimization model to identify a blend of feedstocks, inventory levels, and operating conditions of equipment to ensure a continuous flowing of biomass to the reactor while meeting the requirements of the biochemical conversion process. We propose a sample average approximation (SAA) of the model, and develop an efficient algorithm to solve the SAA model. A feedstock preprocessing process consists of two-stage grinding and pelleting is used to develop a case study. Extensive numerical analysis are conducted which lead to a number of observations. Our main observation is that sequencing bales based on moisture level and carbohydrate content leads to robust solutions that improve processing time and processing rate of the reactor. We provide a number of managerial insights that facilitate the implementation of the model proposed. Note to Practitioners—This paper is motivated by the challenges faced in the bioenergy industry. The focus of this paper is on plants which use the biochemical conversion process to generate liquid fuels. It has been observed that variations in biomass characteristics, such as moisture content, cause variations in feeding of the system which lead to under-utilization of equipment. A requirement of biochemical conversion process is to maintain the carbohydrate content of biomass processed by the reactor, larger than a threshold. We propose a model that identifies the inventory levels and operating conditions of equipment to ensure a continuous flowing of biomass to the reactor. The goal is to improve equipment utilization while satisfying the requirements of the conversion process. The model is tested using real-life data. We found out that by sequencing bales based on moisture level and carbohydrate content, a plant can reduce variability in the system leading to improved system reliability, higher processing rates of the reactor, and higher throughput.

09 BIOMASS FUELS↗
DOE OSTI2022

Genome-wide identification and functional prediction of silicon (Si) transporters in poplar (Populus trichocarpa)

Abstract Silicon (Si) enhances plant tolerance to various biotic and abiotic stressors such as salinity, drought, and heat. In addition, Si can be biomineralized within plants to form organic carbon-containing phytoliths that can have ecosystem-level consequences by contributing to long-term carbon sequestration. Si is taken up and transported in plants via different transporter proteins such as influx transporters (e.g., Lsi1, Lsi6) and efflux transporters (e.g., Lsi2). Additionally, the imported Si can be deposited in plant leaves via silicification process using the Siliplant 1 (e.g., Slp1) protein. Functional homologs of these proteins have been reported in different food crops. Here, we performed a genome-wide analysis to identify different Si transporters and Slp1 homologs in the bioenergy crop poplar ( Populus trichocarpa Torr. and A. Gray ex W. Hook). We identified one channel-type Si influx transporter (PtLsi1; Potri.017G083300), one Si efflux transporter (PtLsi2; Potri.012G144000) and two proteins like Slp1 (PtSlp1a; Potri.004G168600 and PtSlp1b; Potri.009G129900 ) in the P. trichocarpa genome. We found a unique sequence (KPKPPVFKPPPVPI) in PtSlp1a which is repeated six times. Repeated presence of this sequence in PtSlp1a indicates that this protein might be important for silicification processes in P. trichocarpa. The mutation profiles of different Si transporters in a P. trichocarpa genome-wide association study population identified significant and impactful mutations in Potri.004G168600 and Potri.009G129900 . Using a publically accessible database ( http://bar.utoronto.ca/eplant_poplar/ ), digital expression analysis of the putative Si transporters in P. trichocarpa found low to moderate expression in the anticipated tissues, such as roots and leaves. Subcellular localization analysis found that PtLsi1/PtLsi2 are localized in the plasma membrane, whereas PtSlp1a/PtSlp1b are found in the extracellular spaces. Protein–Protein interaction analysis of PtLsi1/PtLsi2 identified Delta-1-pyrroline-5-carboxylate synthase (P5CS) as one of the main interacting partners of PtLsi2, which plays a key role in proline biosynthesis. Proline is a well-known participant in biotic and abiotic stress tolerance in plants. These findings will reinforce future efforts to modify Si accumulation for enhancing plant stress tolerance and carbon sequestration in poplar.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Legacy surface change analyses of the 1993 Rock Valley earthquake sequence for direct comparison to planned NA-22 underground conventional high-explosive experiments (Source Physics Experiment 3 (SPE3) - RV/DC Task 1.3 FY22 Final Report)

Recent work under two previous phases of the NNSA NA-22 Source Physics Experiment (SPE) have shown that underground chemical high-explosive experiments can produce detectable surface changes that differ in spatial extent and vertical magnitude depending on the geologic media at the site (Schultz-Fellenz et al., 2018; 2020; Crawford et al., 2021). Neither of these two prior phases of SPE identified natural earthquake-related surface effects in the same region that occurred at a similar depth as the explosive experiments for direct comparison. While earthquakes can also produce surface changes that are detectable using remote sensing data analyses, it is expected that the pattern and spatial extent of surface changes would vary between earthquakes and explosions. However, no direct-observed surface-change signature comparison between earthquakes and explosions has ever been performed. The SPE Phase 3 Rock Valley/Direct Comparison (RV/DC) program presents a unique opportunity to investigate and characterize co-occurrence of both earthquakes and explosions. In this project, we worked to address five tasks in a workflow, as follows: 1. Identify and obtain existing high-resolution legacy satellite and aerial imagery as close in time before and after the 1993 Rock Valley earthquake sequence to temporally constrain the analyses. 2. Transform these pre-earthquake and post-earthquake datasets into digital elevation models (DEMs) using geospatial analysis software packages (e.g., ArcGIS, Agisoft Metashape, and Google Earth Engine). 3. Perform DEM differencing analyses to assess and quantify earthquake-related changes from the 1993 sequence, and develop map products that visualize these analyses. 4. Use the analyses from (3) to: (a) assess spatial distribution and magnitude of surface changes due to the 1993 Rock Valley earthquake sequence, and (b) determine parameters of forthcoming, planned explosion-related surface change data collection from sensors mounted on unmanned aerial vehicles (UAVs) (e.g., spatial extent of collection, design and density of survey control, sensors to deploy, forward speed and line spacing of UAV flight lines, flight altitude). 5. Develop a summary report on the analyses, including how the analyses define parameters and identify focus areas for any future surface change analytical field campaigns related to the explosive experiment. Analyzing these legacy data and identifying whether they can detect any surface changes related to the earthquake sequence facilitates opportunities for direct signature comparison of surface change from explosions at one location, which has never previously been performed. Comparing the surface change signatures from a co-located and depth-equivalent earthquake and an explosion could help to advance remote sensing event discrimination techniques. This report summarizes the work completed toward this ambitious goal.

42 ENGINEERING↗
DOE OSTI2023

A mechanistic model of a PWR-based nuclear power plant in response to external hazard-induced station blackout accidents

Natural hazard-induced nuclear accidents, such as the Fukushima Daiichi Accident that occurred in Japan in 2011, have significantly increased reactor safety studies in understanding nuclear power plant (NPP) responses to external hazard events such as earthquakes and floods. Natural hazards could cause the loss of offsite power in nuclear power plants, potentially leading to a Station Blackout (SBO) accident that significantly contributes to the overall risk of nuclear power plant accidents. Despite the fact that extensive research has been conducted on the station blackout accident for nuclear power plant, further understanding of these events is needed, particularly in the context of the dynamic nature of external hazards such as external flooding. This paper estimates the progression of station blackout events for a generic pressurized water reactor (PWR) in response to external flooding events. The original RELAP5-3D model of the Westinghouse four-loop design pressurized water reactor was adopted and modified to simulate the external flood-induced station blackout accident, including the short-term and long-term station blackout scenarios. A sensitivity analysis of long-term station blackout, examining reactor operation times and analyzing key parameters over time, was also conducted in this work. The results of the analyses, especially the critical timing parameters of key event sequences, provide useful insights about the time during the external flooding event, which is important for plant operators to make timely decisions to prevent potential core damage. This paper represents significant progress toward developing an integrated risk assessment framework for further identifying and assessing the effects of the critical sources of uncertainties of nuclear power plant under external hazard-induced events.

Liu, Tao↗
DOE OSTI2022

Identification of a new circulating recombinant form of human immunodeficiency virus type 1, CRF124_cpx involving subtypes A, G, H, and CRF27_cpx in Angola

Angola, located in Central Africa, has around 320,000 (270,000–380,000) people living with human immunodeficiency virus (HIV)/AIDS, equivalent to 1% of the country’s population at the end of 2021. A previous study conducted in 2012, using Angolan samples collected between 2008 and 2010 revealed a high prevalence of HIV-1 recombinants, around 42% of sequences, with 21% showing the same UH profile in partial pol region which were grouped into a monophyletic cluster with high bootstrap support. Thus, the objective of the present work was to obtain complete genomes of those sequences and characterize them, aiming at a description of a new circulating recombinant form (CRF). Whole blood from nine HIV-1 UH pol -infected individuals had their genomic DNA extracted, and nested PCR was used to amplify seven overlapping fragments targeting the full-length HIV-1 genome. The final classification was based on maximum likelihood trees, and recombination analyses were performed using a bootscan from the Simplot program. BLAST and Los Alamos Database inspections were used to search other similar H-like pol sequences. Complete genome amplification was possible for three samples, partial genomes were obtained for the other three, and only pol was available for the remaining three sequences. Bootscan analysis of the two whole-genome and three partial genome sequences retrieved from people living with HIV/AIDS (PLHIVA) without epidemiological linkage showed the same complex recombination profile involving HIV-1 subtypes A/G/H/CRF27_cpx, with a total of six recombinant breakpoints, aiming to classify a new HIV-1 CRF124_cpx. We found no other full-length HIV-1 genomes with the same mosaic profile; however, we identified 33 partial pol sequences, mainly sampled from Angola between 2001 to 2019, with the same H-like profile. Bayesian analysis of H and H-like pol sequences indicates that CRF124_cpx probably originated in Angola at mid-1970s, indicating that this CRF has been circulating in the country for a long time. In summary, our study describes a new CRF circulating principally in Angola and highlights the importance of continuing molecular surveillance studies, especially in countries with high molecular diversity of HIV.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Automated Generation of Weather and Climate Analysis Products

Wind roses are an important part to site operations as they depict the wind speed and direction percentage over a time period. In this project, I automated the production of wind roses for SRS meteorological towers. I developed a script to sequence through the dataset by period and create a wind rose for each of those period sets. Graphics were generated for every 4-hour period of the day, for each of the 4 heights of the instruments, for every half of each month. A 10-year climatological period consisted of wind speed and directions measurements at 15-minute intervals. We used the years 2014-2024 as the wind instruments on the tower were upgraded to sonic anemometers in early 2014. I then compared these wind roses to those from a previous study done in 2003 to analyze the differences and similarities in the wind patterns. The wind roses are also used to understand environmental transport conditions at SRS. When comparing the different heights of the anemometers to the 2003 report, there is a similarity in the fact that the direction is relatively the same across the levels, with the wind speeds increasing as you get higher in the air. When comparing the season of our graphs against the 2003 report, we noticed that the winter and summer months appear to have very similar wind directions. However, in the spring we noticed that there were more southerly winds compared to the 2003 report which had more westerly winds. Finally, in the fall months, we noticed that there was a higher percentage of northeasterly winds, while the 2003 report showed more southeasterly winds. We can use these summaries to estimate the directional dependence of dose that the surrounding areas receive throughout the year according to time of day.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Image Processing Pipeline for Fluoroelastomer Crystallite Detection in Atomic Force Microscopy Images

Phase transformations in materials systems can be tracked using atomic force microscopy (AFM), enabling the examination of surface properties and macroscale morphologies. In situ measurements investigating phase transformations generate large datasets of time-lapse image sequences. The interpretation of the resulting image sequences, guided by domain-knowledge, requires manual image processing using handcrafted masks. Here this approach is time-consuming and restricts the number of images that can be processed. Her in this study, we developed an automated image processing pipeline which integrates image detection and segmentation methods. We examine five time-series AFM videos of various fluoroelastomer phase transformations. The number of image sequences per video ranges from a hundred to a thousand image sequences. The resulting image processing pipeline aims to automatically classify and analyze images to enable batch processing. Using this pipeline, the growth of each individual fluoroelastomer crystallite can be tracked through time. We incorporated statistical analysis into the pipeline to investigate trends in phase transformations between different fluoroelastomer batches. Understanding these phase transformations is crucial, as it can provide valuable insights into manufacturing processes, improve product quality, and possibly lead to the development of more advanced fluoroelastomer formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detecting Masquerade Attacks in Controller Area Networks Using Graph Machine Learning

Modern vehicles rely on a myriad of electronic control units (ECUs) interconnected via controller area networks (CANs) for critical operations. Despite their ubiquitous use and reliability, CANs are susceptible to sophisticated cyberattacks, particularly masquerade attacks, which inject false data that mimic legitimate messages at the expected frequency. These attacks pose severe risks such as unintended acceleration, brake deactivation, and rogue steering. Traditional intrusion detection systems (IDS) often struggle to detect these subtle intrusions due to their seamless integration into normal traffic. This paper introduces a novel framework for detecting masquerade attacks in the CAN bus using graph machine learning (ML). We hypothesize that the integration of shallow graph embeddings with time series features derived from CAN frames enhances the detection of masquerade attacks. We show that by representing CAN bus frames as message sequence graphs (MSGs) and enriching each node with contextual statistical attributes from time series, we can enhance detection capabilities across various attack patterns compared to using graph-based features only. Our method ensures a comprehensive and dynamic analysis of CAN frame interactions, improving robustness and efficiency. Extensive experiments on the ROAD dataset validate the effectiveness of our approach, demonstrating statistically significant improvements in the detection rates of masquerade attacks compared to a baseline that uses graph-based features only as confirmed by Mann-Whitney U and Kolmogorov-Smirnov tests (p < 0.05) .

Marfo, William [Univ. of Texas, El Paso, TX (Unite↗
DOE OSTI2021

Data for Intra- and inter-annual variability of nitrification in the rhizosphere of field-grown bioenergy sorghum

These data were collected in 2018 and 2019 at the University of Illinois Energy Farm (N 40.063607, W 88.206926). During each growing season, bulk and rhizosphere soil were collected from replicate Sorghum bicolor nitrogen use efficiency trial plots at three separate time points (approximately July 1, August 1, and September 1). We measured soil moisture, pH, soil nitrate and ammonium, potential nitrification, potential denitrification, and extracted and sequenced the V4 region of the 16S rRNA gene for microbial community analysis. All microbial sequence data is archived in the National Center for Biotechnology Information’s (NCBI) Sequence Read Archive (accession number SRP326979, project number PRJNA741261).

bioenergy↗
DOE OSTI2026

Automated Signal Timing Plan Reconstruction Using High-Resolution Event-Based Controller Data for Digital Twins

Transportation digital twins are essential tools for evaluating emerging technologies such as connected and automated vehicles, adaptive traffic signal control, and mobility optimization strategies. Realistic digital twins require accurate emulation of real-world signal controllers and detailed signal timing plans. However, signal timing plans are often unavailable or difficult to access, forcing researchers and modelers to rely on assumed fixed timings or halt their analysis. To overcome this challenge, we present a method that directly estimates signal timing plan parameters using high-resolution, event-based data from traffic signal controllers. The proposed method extracts key parameters, including cycle length, offset, phase sequence, coordinated phases, phase-specific minimum and maximum green durations, vehicle extensions, and splits under coordination. A rule-based deterministic signal timing reconstruction algorithm based on traffic signal operation rules, such as those outlined in the Signal Timing Manual, is developed and validated. We evaluate this method, which uses high-resolution controller event logs and verified signal timing plans, on 94 signalized intersections in Nashville, Tennessee, demonstrating their ability to generate accurate, simulation-ready signal timing plans for tools such as SUMO and Vissim.

Saroj, Abhilasha [ORNL] (ORCID:0000000191178063)↗
DOE OSTI2023

Effect of fiber reinforcement on the process parameters of injection molding manufacturing process

Use of fillers with the resin is common in injection molding manufacturing technique to enhance the stiffness and strength of the manufactured components. However, the addition of fillers changes the rheological and mechanical properties of the material, required adjustments in the process parameters used for the manufacturing. The present study aims to investigate this impact of fillers on the process parameters of the injection molding manufacturing process. Neat Acrylonitrile Butadiene Styrene (ABS) and ABS reinforced with short glass fibers are used as the materials to perform the comprehensive numerical analysis including complete process sequence, namely filling, packing, cooling and warpage. Parametric studies are conducted to optimize the manufacturing process for minimal warpage (or shrinkage) in the produced part. The selected process parameters for the parametric study include cooling time, packing pressure, and mold temperature, all of which directly impact the quality of the manufactured parts. The optimal values obtained for these parameters are compared to examine the effect fiber reinforcement has on the manufacturing process.

Garg, Nikhil↗
DOE OSTI2025

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗
DOE OSTI2023

classLog: Logistic regression for the classification of genetic sequences

Introduction Sequencing and phylogenetic classification have become a common task in human and animal diagnostic laboratories. It is routine to sequence pathogens to identify genetic variations of diagnostic significance and to use these data in realtime genomic contact tracing and surveillance. Under this paradigm, unprecedented volumes of data are generated that require rapid analysis to provide meaningful inference. Methods We present a machine learning logistic regression pipeline that can assign classifications to genetic sequence data. The pipeline implements an intuitive and customizable approach to developing a trained prediction model that runs in linear time complexity, generating accurate output rapidly, even with incomplete data. Our approach was benchmarked against porcine respiratory and reproductive syndrome virus (PRRSv) and swine H1 influenza A virus (IAV) datasets. Trained classifiers were tested against sequences and simulated datasets that artificially degraded sequence quality at 0, 10, 20, 30, and 40%. Results When applied to a poor-quality sequence data, the classifier achieved between >85% to 95% accuracy for the PRRSv and the swine H1 IAV HA dataset and this increased to near perfect accuracy when using the full dataset. The model also identifies amino acid positions used to determine genetic clade identity through a feature selection ranking within the model. These positions can be mapped onto a maximum-likelihood phylogenetic tree, allowing for the inference of clade defining mutations. Discussion Our approach is implemented as a python package with code available at https://github.com/flu-crew/classLog .

Zeller, Michael A.↗
DOE OSTI2022

A high-throughput skim-sequencing approach for genotyping, dosage estimation and identifying translocations

The development of next-generation sequencing (NGS) enabled a shift from array-based genotyping to directly sequencing genomic libraries for high-throughput genotyping. Even though whole-genome sequencing was initially too costly for routine analysis in large populations such as breeding or genetic studies, continued advancements in genome sequencing and bioinformatics have provided the opportunity to capitalize on whole-genome information. As new sequencing platforms can routinely provide high-quality sequencing data for sufficient genome coverage to genotype various breeding populations, a limitation comes in the time and cost of library construction when multiplexing a large number of samples. Here we describe a high-throughput whole-genome skim-sequencing (skim-seq) approach that can be utilized for a broad range of genotyping and genomic characterization. Using optimized low-volume Illumina Nextera chemistry, we developed a skim-seq method and combined up to 960 samples in one multiplex library using dual index barcoding. With the dual-index barcoding, the number of samples for multiplexing can be adjusted depending on the amount of data required, and could be extended to 3,072 samples or more. Panels of doubled haploid wheat lines ( Triticum aestivum , CDC Stanley x CDC Landmark), wheat-barley ( T . aestivum x Hordeum vulgare ) and wheat-wheatgrass ( Triticum durum x Thinopyrum intermedium ) introgression lines as well as known monosomic wheat stocks were genotyped using the skim-seq approach. Bioinformatics pipelines were developed for various applications where sequencing coverage ranged from 1 × down to 0.01 × per sample. Using reference genomes, we detected chromosome dosage, identified aneuploidy, and karyotyped introgression lines from the skim-seq data. Leveraging the recent advancements in genome sequencing, skim-seq provides an effective and low-cost tool for routine genotyping and genetic analysis, which can track and identify introgressions and genomic regions of interest in genetics research and applied breeding programs.

60 APPLIED LIFE SCIENCES↗