Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PATIENTS”

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 919 records · Page 51

Prenatal Ultrasound Analysis of Umbilical‐Portal‐Systemic Venous Shunts Concurrent With Trisomy 21

Objectives A classification termed umbilical‐portal‐systemic venous shunt (UPSVS) for an abnormal umbilical vein (UV), portal vein (PV), and ductus venosus (DV) was proposed recently. According to this classification, there are 3 types of UPSVSs: types I, II, and III. Trisomy 21 associated with UV‐PV‐DV anomalies has been described, but the incidence of trisomy 21 in UPSVS cases, the relationship between UPSVS types and trisomy 21, and the pregnancy outcome are poorly documented. This study aimed to address these issues. Methods All UPSVS cases diagnosed at our department from 2016 to 2019 were retrospectively studied. The English literature describing UV‐PV‐DV anomalies and trisomy 21 from 2000 to 2019 was searched, and the retrieved cases were analyzed. Results Four of 20 UPSVS cases identified by us also had trisomy 21, with 2 type I and 2 type II UPSVSs. Ultrasound markers of Down syndrome were observed in all 4 cases that underwent termination of pregnancy (TOP). The literature search retrieved 12 reports including 279 patients, with 29 also having trisomy 21, giving a pooled trisomy 21 incidence rate of 10.4%. Of the 29 cases, 16 had type I, and 9 had type II, whereas UPSVS types in 4 were undeterminable, and 22 cases underwent TOP. Conclusions There is a high incidence of trisomy 21 in UPSVS cases. Trisomy 21 is associated with a type I or II UPSVS. Most cases with the combined defect underwent TOP. These findings may be used to direct prenatal counseling and management of the combined condition.

Dong, Xiangyi↗

Computational models of direct and indirect X‐ray breast imaging detectors for in silico trials

Abstract Background To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed. Purpose The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems. Methods We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open‐source GPU‐accelerated Monte Carlo x‐ray imaging simulation code (MC‐GPU). We build on the previous version of the MC‐GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)‐approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x‐ray photon till its first interaction within the x‐ray detector. On the other hand, the signal spread in the x‐ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron–hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth‐dependent spread and collection of optical photons and relies on a pre‐computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x‐ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data. Results We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors. Conclusions We have presented computational models of three commercially available FDA‐approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC‐GPU codes that can be used to replicate three systems are available in open source format through GitHub.

Sengupta, Aunnasha↗

Screening a knowledge‐based library of low molecular weight compounds against the proline biosynthetic enzyme 1‐pyrroline‐5‐carboxylate 1 ( PYCR1)

Abstract Δ 1 ‐pyrroline‐5‐carboxylate reductase isoform 1 (PYCR1) is the last enzyme of proline biosynthesis and catalyzes the NAD(P)H‐dependent reduction of Δ 1 ‐pyrroline‐5‐carboxylate toL‐proline. High PYCR1 gene expression is observed in many cancers and linked to poor patient outcomes and tumor aggressiveness. The knockdown of thePYCR1gene or the inhibition of PYCR1 enzyme has been shown to inhibit tumorigenesis in cancer cells and animal models of cancer, motivating inhibitor discovery. We screened a library of 71 low molecular weight compounds (average MW of 131 Da) against PYCR1 using an enzyme activity assay. Hit compounds were validated with X‐ray crystallography and kinetic assays to determine affinity parameters. The library was counter‐screened against human Δ 1 ‐pyrroline‐5‐carboxylate reductase isoform 3 and proline dehydrogenase (PRODH) to assess specificity/promiscuity. Twelve PYCR1 and one PRODH inhibitor crystal structures were determined. Three compounds inhibit PYCR1 with competitive inhibition parameter of 100 μM or lower. Among these, (S)‐tetrahydro‐2H‐pyran‐2‐carboxylic acid (70 μM) has higher affinity than the current best tool compoundN‐formyl‐l‐proline, is 30 times more specific for PYCR1 over human Δ 1 ‐pyrroline‐5‐carboxylate reductase isoform 3, and negligibly inhibits PRODH. Structure‐affinity relationships suggest that hydrogen bonding of the heteroatom of this compound is important for binding to PYCR1. The structures of PYCR1 and PRODH complexed with 1‐hydroxyethane‐1‐sulfonate demonstrate that the sulfonate group is a suitable replacement for the carboxylate anchor. This result suggests that the exploration of carboxylic acid isosteres may be a promising strategy for discovering new classes of PYCR1 and PRODH inhibitors. The structure of PYCR1 complexed withl‐pipecolate and NADH supports the hypothesis that PYCR1 has an alternative function in lysine metabolism.

Biochemistry & Molecular Biology↗

Case studies in bias reduction and inference for electronic health record data with selection bias and phenotype misclassification

Electronic health records (EHR) are not designed for population‐based research, but they provide easy and quick access to longitudinal health information for a large number of individuals. Many statistical methods have been proposed to account for selection bias, missing data, phenotyping errors, or other problems that arise in EHR data analysis. However, addressing multiple sources of bias simultaneously is challenging. We developed a methodological framework (R package, SAMBA ) for jointly handling both selection bias and phenotype misclassification in the EHR setting that leverages external data sources. These methods assume factors related to selection and misclassification are fully observed, but these factors may be poorly understood and partially observed in practice. As a follow‐up to the methodological work, we demonstrate how to apply these methods for two real‐world case studies, and we evaluate their performance. In both examples, we use individual patient‐level data collected through the University of Michigan Health System and various external population‐based data sources. In case study (a), we explore the impact of these methods on estimated associations between gender and cancer diagnosis. In case study (b), we compare corrected associations between previously identified genetic loci and age‐related macular degeneration with gold standard external summary estimates. These case studies illustrate how to utilize diverse auxiliary information to achieve less biased inference in EHR‐based research.

60 APPLIED LIFE SCIENCES↗

Use of Longitudinal Serum Analysis and Machine Learning to Develop a Classifier for Cancer Early Detection

Early detection of solid tumors through a simple screening process, such as the proteomic analysis of biofluids, has the potential to significantly alter the management and outcomes of cancers. The application of advanced targeted proteomics measurements and data analysis strategies to uniformly collected serum or plasma samples would enable longitudinal studies of cancer risk, progression, and response to therapy that have the potential to significantly reduce cancer burden in general. In this article, we describe a generalizable workflow combining robust, multiplexed targeted proteomics measurements applied to longitudinal samples from the Department of Defense Serum Repository with a Random Forest machine learning method for developing and initially evaluating the performance of candidate biomarker panels for early detection of cancers. The effectiveness of this approach was demonstrated in a cohort of 175 head and neck squamous cell carcinoma patients. The outlined protocols include methods for sample preparation, instrument analysis, and data analysis and interpretation using this workflow.

Longitudinal analysis, machine learning, cancer, e↗

Comprehensive characterization of hepatocyte-derived extracellular vesicles identifies direct miRNA-based regulation of hepatic stellate cells and DAMP-based hepatic macrophage IL-1ß and IL-17 upregulation in alcoholic hepatitis mice.

Extracellular vesicles (EVs) have been growingly recognized as biomarkers and mediators of alcoholic liver disease (ALD) in human and mice. Here we characterized hepatocyte-derived EVs (HC-EVs) and their cargo for their biological functions in a novel murine model that closely resembles liver pathology observed in patients with alcoholic hepatitis (AH), the most severe spectrum of ALD. The number of circulating EVs and HC-EVs were significantly increased by 10-fold in AH mice compared to control mice. The miRNA(miR)-seq analysis detected 20 up-regulated and 4 down-regulated miRNAs (P < 0.001-0.05) in AH-HC-EVs. Treatment of murine primary hepatic stellate cells (HSCs) with AH-HC-EVs induced ?-SMA (P < 0.05) and Col1a1 (P < 0.001). Smad7, Nr1d2, and Rora genes which regulate HSC activation were predicted targets of 17 miRs up-regulated in AH-HC-EVs. Among them were miR-27a and miR-181 which upon transfection in HSCs, repressed Nr1d2, the HSC quiescence marker. AH-HC-EVs were also enriched in organelle proteins and mitochondrial DNA (10-fold, P < 0.05) and upregulated IL-1??and IL-17 production by hepatic macrophages (HMs) from AH mice in a TLR9-dependent manner. These results demonstrate HC-EV release is intensified in AH and suggest that AH-HC-EVs orchestrate liver fibrogenesis by directly targeting HSC quiescence transcripts via a unique set of miRNAs and by amplifying HSC activation via DAMP-based induction of profibrogenic IL-1? and IL-17 by HMs.

AH, hepatocyte-derived EVs, miRNAs, DAMPs, HSC act↗

Evaluating inkjet dispenser/liquid vortex capture-mass spectrometry for single-cell metabolomics in Hep G2 steatosis caused by tamoxifen

Single-cell mass spectrometry (MS) is advancing our understanding of metabolic pathways in heterogeneous cell populations; however, many techniques are slow or require disruptive sample preparations. This study evaluated coupling a modified HP D100 single-cell inkjet dispenser with liquid vortex capture-mass spectrometry (D100/LVC-MS). The D100 is a single-cell inkjet dispenser capable of titrating solutions and isolating single cells via disposable cassettes equipped with microfluidic channels and an impedance sensor. The LVC-MS enables high-throughput capture, lysis, and ionization of analytes for mass spectrometric analysis. The D100/LVC-MS system was characterized through titration and single-cell experiments. Propranolol titration demonstrated linearity across a broad concentration range using the D100/LVC-MS system. Additionally, Hep G2 hepatocarcinoma cells and Chlamydomonas reinhardtii algae were used to showcase the D100’s high-throughput or low-buffer-volume single-cell dispensing strategies. The D100/LVC-MS system’s performance was validated by evaluating tamoxifen-induced steatosis in Hep G2 cells. Tamoxifen, associated with nonalcoholic fatty liver disease in breast cancer patients following long-term use, was tested in Hep G2 cells at 20 µM for 72 h against DMSO-treated controls. High-throughput analysis of 500 cells per condition, completed in 25 min per run, demonstrated the system’s efficiency. The D100/LVC-MS system simultaneously quantified tamoxifen and measured triglycerides and phosphatidylcholines. Triglycerides were upregulated in the tamoxifen-treated cells and the results indicated two distinct cell populations, differing in tamoxifen and phosphatidylcholines levels, suggesting heterogeneity within the treated population. In conclusion, these findings highlight the D100/LVC-MS system as a cost-effective, high-throughput platform for single-cell metabolomics and lipidomics, with significant potential for evaluating metabolic alterations.

Chromatography↗

Cohort-based pan-cancer analysis and experimental studies reveal ISG15 gene as a novel biomarker for prognosis and immunotherapy efficacy prediction

Abstract ISG15, an interferon-stimulated ubiquitin-like protein, plays a multifaceted role in tumorigenesis and immune regulation. This study comprehensively evaluates ISG15 as a prognostic biomarker and predictor of immunotherapy response through pan-cancer bioinformatics analysis and experimental validation. By integrating multiomics data from TCGA, GEO, and clinical cohorts, we found that ISG15 is significantly overexpressed in multiple cancers and generally correlates with poor prognosis. Elevated ISG15 expression is associated with increased immune checkpoint gene expression, particularly PD-L1, and immune infiltration, notably M2-like tumor-associated macrophages. Immunohistochemistry and multiplexed immunofluorescence confirmed a strong positive correlation between ISG15, PD-L1, and M2-TAM infiltration in lung and gastric cancer samples. Functional analysis at the single-cell level revealed significant associations between ISG15 and tumor proliferation, angiogenesis, and immune suppression. Immunotherapy cohort analysis demonstrated that tumors with high ISG15 expression responded favorably to PD-L1 inhibitors but exhibited resistance to CTLA-4 blockade, findings further validated in lung cancer patients receiving anti-PD-1 therapy. These results suggest that ISG15 is a promising biomarker for prognosis and immunotherapy response prediction across cancers. Its integration into clinical decision-making may enhance personalized treatment strategies, improve immunotherapy outcomes, and provide new insights into the tumor immune microenvironment, cancer progression, and potential therapeutic targets for future drug development.

Immunology↗

Continuous and discrete modeling of HIV-1 decline on therapy

Mathematical models have shed light on the dynamics of HIV- 1 infection in vivo. In this paper, we generalize continuous mathematical models of drug therapy for HIV-1 by Perelson et al. (Science 271:1582–1586, 1996) and Perelson and Nelson (SIAM Rev 41:3–44, 1999) on time scales, i.e., a nonempty closed subset of real numbers in order to derive new discrete models that predict the total concentration of plasma virus as a function of time. Here, one of our main goals is to compare discrete mathematical models with the continuous model in Perelson et al. (1996) where HIV infected patients were given protease inhibitors and sampled frequently thereafter. For the comparison, we use experimental data collected in Perelson et al. (1996) and estimate the parameters such as the virion clearance rate and the rate of loss of infected cells by fitting the total concentration of plasma virus to this data set. Our results show that discrete systems describe the best fit. In the previous models of this study, the efficacy of protease inhibitor is assumed to be perfect. Motivated by Perelson and Nelson (1999), we end the paper with a mathematical model of imperfect protease inhibitor and reverse transcriptase (RT) inhibitor combination therapy of HIV-1 infection on time scales with its stability analysis.

60 APPLIED LIFE SCIENCES↗

Domain Shift Analysis in Chest Radiographs Classification in a Veterans Healthcare Administration Population

This study aims to assess the impact of domain shift on chest X-ray classification accuracy and to analyze the influence of ground truth label quality and demographic factors such as age group, sex, and study year. We used a DenseNet121 model pre-trained MIMIC-CXR dataset for deep learning-based multi-label classification using ground truth labels from radiology reports extracted using the CheXpert and CheXbert Labeler. We compared the performance of the 14 chest X-ray labels on the MIMIC-CXR and Veterans Healthcare Administration chest X-ray dataset (VA-CXR). The validation of ground truth and the assessment of multi-label classification performance across various NLP extraction tools revealed that the VA-CXR dataset exhibited lower disagreement rates than the MIMIC-CXR datasets. Additionally, there were notable differences in AUC scores between models utilizing CheXpert and CheXbert. When evaluating multi-label classification performance across different datasets, minimal domain shift was observed in the unseen VA dataset, except for the label “Enlarged Cardiomediastinum.” The subgroup with the most significant variations in multi-label classification performance was study year. These findings underscore the importance of considering domain shift in chest X-ray classification tasks, paying particular attention to the temporality of the exam. Our study reveals the significant impact of domain shift and demographic factors on chest X-ray classification, emphasizing the need for improved transfer learning and robust model development. Addressing these challenges is crucial for advancing medical imaging research and improving patient care.

chest X-ray image classification↗

Post-Combat-Injury Opioid Prescription and Alcohol Use Disorder in the Military

Previous studies have identified combat exposure and combat traumatic experience as problematic drinking risk factors. Increasing evidence suggests that opioid use increases the risk of alcohol use disorder. This study investigated the association between opioid prescription use after injury and (1) alcohol use disorder and (2) severity of alcohol use disorder among deployed military servicemembers. Deidentified health records data of 9,029 deployed servicemembers from a retrospective cohort study were analyzed. Data were randomly selected from the Department of Defense Trauma Registry and included servicemembers with combat injuries during deployment in Iraq or Afghanistan (2002–2016). Pharmacy records and International Classification of Diseases, Ninth and Tenth Revision diagnosis codes were used. Three groups were identified (no opioid prescription use, nonpersistent opioid prescription use, and persistent opioid prescription use) and were compared on the basis of alcohol use disorder risk using Cox proportional hazard models. Data analyses were performed in 2021. Of the 9,029 servicemembers with combat injury, 2,262 developed alcohol use disorder (1,322 developed severe alcohol use disorder). Compared with no opioid prescription use, increased alcohol use disorder risk was associated with persistent opioid prescription use, with a hazard ratio of 1.13 (95% CI=1.02, 1.26). After covariate adjustment, increased risk remained statistically significant (hazards ratio=1.24; 95% CI=1.10, 1.39). There was no significant difference in alcohol use disorder risk between no opioid prescription use and nonpersistent opioid prescription use. The risk of severe alcohol use disorder did not vary by opioid use among servicemembers with alcohol use disorder diagnosis. Furthermore, the findings of the study suggest that the incidence of alcohol use disorder was higher among injured servicemembers with persistent opioid prescription use than among those without opioid use. If replicated in prospective studies, the findings highlight the need for clinicians to consider the current and history of alcohol use of patients in initiating treatment involving opioids.

60 APPLIED LIFE SCIENCES↗

2.5 Å-resolution structure of human CDK-activating kinase bound to the clinical inhibitor ICEC0942

The human CDK-activating kinase (CAK), composed of CDK7, cyclin H, and MAT1, is involved in the control of transcription initiation and the cell cycle. Because of these activities, it has been identified as a promising target for cancer chemotherapy. A number of CDK7 inhibitors have entered clinical trials, among them ICEC0942 (also known as CT7001). Structural information can aid in improving the affinity and specificity of such drugs or drug candidates, reducing side effects in patients. Here, we have determined the structure of the human CAK in complex with ICEC0942 at 2.5 Å-resolution using cryogenic electron microscopy. Our structure reveals conformational differences of ICEC0942 compared with previous X-ray crystal structures of the CDK2-bound complex, and highlights the critical ability of cryogenic electron microscopy to resolve structures of drug-bound protein complexes without the need to crystalize the protein target.

59 BASIC BIOLOGICAL SCIENCES↗

Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent Cα model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.

59 BASIC BIOLOGICAL SCIENCES↗

The AML microenvironment catalyzes a step-wise evolution to gilteritinib resistance

This study reports a stepwise evolution of gilteritinib resistance in acute myeloid leukemia (AML) relapse, and identifies therapeutic vulnerabilities that could be exploited. Our in vitro model demonstrates that extrinsic microenvironmental factors within the bone marrow nurture the growth of residual cells that underlie early gilteritinib resistance. These early resistant cells exhibit a dependency upon Aurora B (AURKB) and undergo dynamic metabolic adaptations that fostered their transition to late resistance. In line with clinical observations of gilteritinib resistance, late resistant cells displayed an expansion of NRAS mutant subclones and a unique metabolic profile that originated during early resistance. Pharmacological inhibition of AURKB resensitized early resistant cells to gilteritinib. Early resistant cells from AML patients were found to display exquisite sensitivitye to AURKB inhibition. Together, our findings support that eradication of early resistant AML cells may increase the efficacy and durability of gilteritinib.

AML, Proteomics↗

CD63-mediated cloaking of VEGF in small extracellular vesicles contributes to anti-VEGF therapy resistance

Despite wide use of anti-vascular endothelial growth factor (VEGF) therapy for many solid 30 cancers, most of them become resistant to this therapy and develop into progressive disease. 31 Therefore, new biomarkers and strategies for blocking adaptive resistance of cancer to anti-VEGF therapy are needed. As described herein, we found that cancer-derived small extracellular vesicles (EVs) package increasing quantities of VEGF and other factors in response to anti-VEGF therapy. Notably, small EV-VEGF (eVEGF) is not recognized by bevacizumab and can be transferred to endothelial cells and trigger intracrine signaling and promote angiogenesis. Also, serum eVEGF was present at higher levels in ovarian cancer mouse models with adaptive resistance to bevacizumab than in models sensitive to it. Ovarian cancer cell-derived eVEGF increased tumor growth despite treatment with bevacizumab. Notably, the eVEGF level was elevated in patient serum after bevacizumab-containing therapy. Proteomic data of small EVs from mouse models suggested that several novel biomarkers are promising for drug response prediction. Collectively, these data demonstrated a new mechanism whereby eVEGF evades recognition by therapeutic antibodies and promotes tumor angiogenesis and progression. These findings have clinical implications for biomarkers and new therapeutic strategies for ovarian cancer.

59 BASIC BIOLOGICAL SCIENCES↗

An antibody-free platform for multiplexed, sensitive quantification of protein biomarkers in complex biomatrices

Sensitive, multiplexed protein quantification remains challenging despite recent advancements in LC-MS assays for targeted protein biomarker quantification. High-sensitivity protein biomarker measurements usually require immuno-affinity enrichment of target protein; a process which is highly dependent on capture reagent and limited in capability to measure multiple analytes. Herein, we report a novel antibody-free platform, which measures multiple biomarkers from complex matrices employing a strategically optimized solid-phase extraction cleanup and orthogonal multidimensional LC-MS. Eight human protein biomarkers with different specifications were spiked into canine plasma as a model investigation system. The developed strategy achieved the desired sensitivity, robustness, and throughput via the following steps: (1) post digestion mixed-mode cation exchange-reverse phase SPE enrichment cleaned up the sample initially; (2) rapid, high-pH peptide fractionation further eliminated background components efficiently while selectively enriched signature peptides (SP) to provide sufficient sensitivity for multiple targets; and (3) trapping-micro-LC-MS analysis delivered high sensitivity comparable to a nano-LC-MS method but with much better robustness and throughput for the final analysis. Compared with a conventional LC-MS assay with direct protein digestion and limited clean-up, analysis with this antibody-free platform improved the LLOQ by 1–2 orders of magnitude for the eight protein biomarkers, reaching as low as 5 ng/mL in plasma, with feasible robustness and throughput. In conclusion, this platform was applied for the quantification of biomarkers of respiratory conditions in patients with various lung diseases, demonstrating real-world applicability.

59 BASIC BIOLOGICAL SCIENCES↗

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus’ genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

59 BASIC BIOLOGICAL SCIENCES↗