SEARCH · Engineering Papers

Results for “PCA”

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2014

A statistical method for retrospective cardiac and respiratory motion gating of interventional cardiac x-ray images

Purpose: Image-guided cardiac interventions involve the use of fluoroscopic images to guide the insertion and movement of interventional devices. Cardiorespiratory gating can be useful for 3D reconstruction from multiple x-ray views and for reducing misalignments between 3D anatomical models overlaid onto fluoroscopy. Methods: The authors propose a novel and potentially clinically useful retrospective cardiorespiratory gating technique. The principal component analysis (PCA) statistical method is used in combination with other image processing operations to make our proposed masked-PCA technique suitable for cardiorespiratory gating. Unlike many previously proposed techniques, our technique is robust to varying image-content, thus it does not require specific catheters or any other optically opaque structures to be visible. Therefore, it works without any knowledge of catheter geometry. The authors demonstrate the application of our technique for the purposes of retrospective cardiorespiratory gating of normal and very low dose x-ray fluoroscopy images. Results: For normal dose x-ray images, the algorithm was validated using 28 clinical electrophysiology x-ray fluoroscopy sequences (2168 frames), from patients who underwent radiofrequency ablation (RFA) procedures for the treatment of atrial fibrillation and cardiac resynchronization therapy procedures for heart failure. The authors established end-systole, end-expiration, and end-inspiration success rates of 97.0%, 97.9%, and 97.0%, respectively. For very low dose applications, the technique was tested on ten x-ray sequences from the RFA procedures with added noise at signal to noise ratio (SNR) values of√(5)0, √(1)0, √(8), √(6), √(5), √(2), and √(1) to simulate the image quality of increasingly lower dose x-ray images. Even at the low SNR value of √(2), representing a dose reduction of more than 25 times, gating success rates of 89.1%, 88.8%, and 86.8% were established. Conclusions: The proposed technique can therefore extract useful information from interventional x-ray images while minimizing exposure to ionizing radiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

SU-E-J-01: 3D Fluoroscopic Image Estimation From Patient-Specific 4DCBCT-Based Motion Models

Purpose: 3D motion modeling derived from 4DCT images, taken days or weeks before treatment, cannot reliably represent patient anatomy on the day of treatment. We develop a method to generate motion models based on 4DCBCT acquired at the time of treatment, and apply the model to estimate 3D time-varying images (referred to as 3D fluoroscopic images). Methods: Motion models are derived through deformable registration between each 4DCBCT phase, and principal component analysis (PCA) on the resulting displacement vector fields. 3D fluoroscopic images are estimated based on cone-beam projections simulating kV treatment imaging. PCA coefficients are optimized iteratively through comparison of these cone-beam projections and projections estimated based on the motion model. Digital phantoms reproducing ten patient motion trajectories, and a physical phantom with regular and irregular motion derived from measured patient trajectories, are used to evaluate the method in terms of tumor localization, and the global voxel intensity difference compared to ground truth. Results: Experiments included: 1) assuming no anatomic or positioning changes between 4DCT and treatment time; and 2) simulating positioning and tumor baseline shifts at the time of treatment compared to 4DCT acquisition. 4DCBCT were reconstructed from the anatomy as seen at treatment time. In case 1) the tumor localization error and the intensity differences in ten patient were smaller using 4DCT-based motion model, possible due to superior image quality. In case 2) the tumor localization error and intensity differences were 2.85 and 0.15 respectively, using 4DCT-based motion models, and 1.17 and 0.10 using 4DCBCT-based models. 4DCBCT performed better due to its ability to reproduce daily anatomical changes. Conclusion: The study showed an advantage of 4DCBCT-based motion models in the context of 3D fluoroscopic images estimation. Positioning and tumor baseline shift uncertainties were mitigated by the 4DCBCT-based motion models, while they caused errors when using 4DCT-based motion models. Partially funded by Varian research grant.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

TU-C-17A-03: An Integrated Contour Evaluation Software Tool Using Supervised Pattern Recognition for Radiotherapy

Purpose: Radiotherapy (RT) contours delineated either manually or semiautomatically require verification before clinical usage. Manual evaluation is very time consuming. A new integrated software tool using supervised pattern contour recognition was thus developed to facilitate this process. Methods: The contouring tool was developed using an object-oriented programming language C# and application programming interfaces, e.g. visualization toolkit (VTK). The C# language served as the tool design basis. The Accord.Net scientific computing libraries were utilized for the required statistical data processing and pattern recognition, while the VTK was used to build and render 3-D mesh models from critical RT structures in real-time and 360° visualization. Principal component analysis (PCA) was used for system self-updating geometry variations of normal structures based on physician-approved RT contours as a training dataset. The inhouse design of supervised PCA-based contour recognition method was used for automatically evaluating contour normality/abnormality. The function for reporting the contour evaluation results was implemented by using C# and Windows Form Designer. Results: The software input was RT simulation images and RT structures from commercial clinical treatment planning systems. Several abilities were demonstrated: automatic assessment of RT contours, file loading/saving of various modality medical images and RT contours, and generation/visualization of 3-D images and anatomical models. Moreover, it supported the 360° rendering of the RT structures in a multi-slice view, which allows physicians to visually check and edit abnormally contoured structures. Conclusion: This new software integrates the supervised learning framework with image processing and graphical visualization modules for RT contour verification. This tool has great potential for facilitating treatment planning with the assistance of an automatic contour evaluation module in avoiding unnecessary manual verification for physicians/dosimetrists. In addition, its nature as a compact and stand-alone tool allows for future extensibility to include additional functions for physicians’ clinical needs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

WE-G-BRF-04: Robust Real-Time Volumetric Imaging Based On One Single Projection

Purpose: Real-time volumetric imaging is highly desirable to provide instantaneous image guidance for lung radiation therapy. This study proposes a scheme to achieve this goal using one single projection by utilizing sparse learning and a principal component analysis (PCA) based lung motion model. Methods: A patient-specific PCA-based lung motion model is first constructed by analyzing deformable vector fields (DVFs) between a reference image and 4DCT images at each phase. At the training stage, we “learn” the relationship between the DVFs and the projection using sparse learning. Specifically, we first partition the projections into patches, and then apply sparse learning to automatically identify patches that best correlate with the principal components of the DVFs. Once the relationship is established, at the application stage, we first employ a patchbased intensity correction method to overcome the problem of different intensity scale between the calculated projection in the training stage and the measured projection in the application stage. The corrected projection image is then fed to the trained model to derive a DVF, which is applied to the reference image, yielding a volumetric image corresponding to the projection. We have validated our method through a NCAT phantom simulation case and one experiment case. Results: Sparse learning can automatically select those patches containing motion information, such as those around diaphragm. For the simulation case, over 98% of the lung region pass the generalized gamma test (10HU/1mm), indicating combined accuracy in both intensity and spatial domain. For the experimental case, the average tumor localization errors projected to the imager are 0.68 mm and 0.4 mm on the axial and tangential direction, respectively. Conclusion: The proposed method is capable of accurately generating a volumetric image using one single projection. It will potentially offer real-time volumetric image guidance to facilitate lung radiotherapy. This study is supported in part by NIH (1R01CA154747-01), The Core Technology Research in Strategic Emerging Industry, Guangdong, China(2011A081402003)

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

TH-A-9A-01: Active Optical Flow Model: Predicting Voxel-Level Dose Prediction in Spine SBRT

Purpose: To predict voxel-level dose distribution and enable effective evaluation of cord dose sparing in spine SBRT. Methods: We present an active optical flow model (AOFM) to statistically describe cord dose variations and train a predictive model to represent correlations between AOFM and PTV contours. Thirty clinically accepted spine SBRT plans are evenly divided into training and testing datasets. The development of predictive model consists of 1) collecting a sequence of dose maps including PTV and OAR (spinal cord) as well as a set of associated PTV contours adjacent to OAR from the training dataset, 2) classifying data into five groups based on PTV's locations relative to OAR, two “Top”s, “Left”, “Right”, and “Bottom”, 3) randomly selecting a dose map as the reference in each group and applying rigid registration and optical flow deformation to match all other maps to the reference, 4) building AOFM by importing optical flow vectors and dose values into the principal component analysis (PCA), 5) applying another PCA to features of PTV and OAR contours to generate an active shape model (ASM), and 6) computing a linear regression model of correlations between AOFM and ASM.When predicting dose distribution of a new case in the testing dataset, the PTV is first assigned to a group based on its contour characteristics. Contour features are then transformed into ASM's principal coordinates of the selected group. Finally, voxel-level dose distribution is determined by mapping from the ASM space to the AOFM space using the predictive model. Results: The DVHs predicted by the AOFM-based model and those in clinical plans are comparable in training and testing datasets. At 2% volume the dose difference between predicted and clinical plans is 4.2±4.4% and 3.3±3.5% in the training and testing datasets, respectively. Conclusion: The AOFM is effective in predicting voxel-level dose distribution for spine SBRT. Partially supported by NIH/NCI under grant #R21CA161389 and a master research grant by Varian Medical System.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

SU-E-J-44: A Novel Approach to Quantify Patient Setup and Target Motion for Real-Time Image-Guided Radiotherapy (IGRT)

Purpose Isocenter shifts and rotations to correct patient setup errors and organ motion cannot remedy some shape changes of large targets. We are investigating new methods in quantification of target deformation for realtime IGRT of breast and chest wall cancer. Methods Ninety-five patients of breast or chest wall cancer were accrued in an IRB-approved clinical trial of IGRT using 3D surface images acquired at daily setup and beam-on time via an in-room camera. Shifts and rotations relating to the planned reference surface were determined using iterative-closest-point alignment. Local surface displacements and target deformation are measured via a ray-surface intersection and principal component analysis (PCA) of external surface, respectively. Isocenter shift, upper-abdominal displacement, and vectors of the surface projected onto the two principal components, PC1 and PC2, were evaluated for sensitivity and accuracy in detection of target deformation. Setup errors for some deformed targets were estimated by superlatively registering target volume, inner surface, or external surface in weekly CBCT or these outlines on weekly EPI. Results Setup difference according to the inner-surface, external surface, or target volume could be 1.5 cm. Video surface-guided setup agreed with EPI results to within < 0.5 cm while CBCT results were sometimes (∼20%) different from that of EPI (>0.5 cm) due to target deformation for some large breasts and some chest walls undergoing deep-breath-hold irradiation. Square root of PC1 and PC2 is very sensitive to external surface deformation and irregular breathing. Conclusion PCA of external surfaces is quick and simple way to detect target deformation in IGRT of breast and chest wall cancer. Setup corrections based on the target volume, inner surface, and external surface could be significant different. Thus, checking of target shape changes is essential for accurate image-guided patient setup and motion tracking of large deformable targets. NIH grant for the first author as cionsultant and the last author as the PI.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

TU-F-CAMPUS-T-05: Dose Escalation to Biological Tumor Volumes of Prostate Cancer Patients Using Gold Nanoparticles

Purpose: Studies have shown that radiation boosting could help reduce prostate cancer (PCa) recurrence. Biological tumor volumes (BTV) are a high priority for such radiation boosting. The purpose of this study is to investigate the potential of radiation boosting of real patient BTVs using gold nanoparticles (GNP) released from gold-loaded brachytherapy spacers (GBS) during brachytherapy. Methods: The BTVs of 12 patients having prostate adenocarcinoma identified with positron emission tomography (PET) and CT scanner using C-11 labeled tracer [11C]acetate were investigated. The initial GNP concentration and time to achieve a dose enhancement effect (DEF) of 2 was simulated using the freely downloadable software RAID APP. The investigations were carried out for low dose rate (LDR) brachytherapy sources (BTS) described in AAPM Task Group report 43: Cs-131, I-125, and Pd-103. In first case, we used 7 mg/g and 18 mg/g of GNP initial concentrations to estimate the time needed for released GNP to achieve a DEF of 2 for the different BTS, and compare with clinically relevant treatment times. In second case, we calculated the initial concentration of GNPs needed to achieve a DEF of 2 during the time the BTS would typically deliver 50%, 70% and 90% of the total dose. Results: For an initial concentration of 18 mg/g, when using Cs-131, and Pd-103, a DEF of 2 could only be achieved for BTV of 3.3 cm3 and 1 cm3 respectively. Meanwhile a DEF of 2 could be achieved for all 12 BTVs when using I-125. To achieve a DEF of 2 for all patients using Cs-131 and Pd-103, much higher initial concentrations would have to be used than have been typically employed in pre-clinical studies. Conclusion: The I-125 is the most viable BTS that can be employed with GBS to guide dose painting treatment planning for localized PCa.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

WE-G-BRD-01: A Data-Driven 4D-MRI Motion Model to Estimate Full Field-Of-View Abdominal Motion From 2D Image Navigators During MR-Linac Treatment

Purpose: To estimate full field-of-view abdominal respiratory motion from fast 2D image navigators using a 4D-MRI based motion model. This will allow for radiation dose accumulation mapping during MR-Linac treatment. Methods: Experiments were conducted on a Philips Ingenia 1.5T MRI. First, a retrospectively ordered 4D-MRI was constructed using 3D transient-bSSFP with radial in-plane sampling. Motion fields were calculated through 3D non-rigid registration. From these motion fields a PCA-based abdominal motion model was constructed and used to warp a 3D reference volume to fast 2D cine-MR image navigators that can be used for real-time tracking. To test this procedure, a time-series consisting of two interleaved orthogonal slices (sagittal and coronal), positioned on the pancreas or kidneys, were acquired for 1m38s (dynamic scan-time=0.196ms), during normal, shallow, or deep breathing. The coronal slices were used to update the optimal weights for the first two PCA components, in order to warp the 3D reference image and construct a dynamic 4D-MRI time-series. The interleaved sagittal slices served as an independent measure to test the model’s accuracy and fit. Spatial maps of the root-mean-squared error (RMSE) and histograms of the motion differences within the pancreas and kidneys were used to evaluate the method. Results: Cranio-caudal motion was accurately calculated within the pancreas using the model for normal and shallow breathing with an RMSE of 1.6mm and 1.5mm and a histogram median and standard deviation below 0.2 and 1.7mm, respectively. For deep-breathing an underestimation of the inhale amplitude was observed (RMSE=4.1mm). Respiratory-induced antero-posterior and lateral motion were correctly mapped (RMSE=0.6/0.5mm). Kidney motion demonstrated good motion estimation with RMSE-values of 0.95 and 2.4mm for the right and left kidney, respectively. Conclusion: We have demonstrated a method that can calculate dynamic 3D abdominal motion in a large volume, while acquiring real-time cine-MR images for MR-guided radiotherapy.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

WE-G-207-06: 3D Fluoroscopic Image Generation From Patient-Specific 4DCBCT-Based Motion Models Derived From Physical Phantom and Clinical Patient Images

Purpose: Respiratory-correlated cone-beam CT (4DCBCT) images acquired immediately prior to treatment have the potential to represent patient motion patterns and anatomy during treatment, including both intra- and inter-fractional changes. We develop a method to generate patient-specific motion models based on 4DCBCT images acquired with existing clinical equipment and used to generate time varying volumetric images (3D fluoroscopic images) representing motion during treatment delivery. Methods: Motion models are derived by deformably registering each 4DCBCT phase to a reference phase, and performing principal component analysis (PCA) on the resulting displacement vector fields. 3D fluoroscopic images are estimated by optimizing the resulting PCA coefficients iteratively through comparison of the cone-beam projections simulating kV treatment imaging and digitally reconstructed radiographs generated from the motion model. Patient and physical phantom datasets are used to evaluate the method in terms of tumor localization error compared to manually defined ground truth positions. Results: 4DCBCT-based motion models were derived and used to generate 3D fluoroscopic images at treatment time. For the patient datasets, the average tumor localization error and the 95th percentile were 1.57 and 3.13 respectively in subsets of four patient datasets. For the physical phantom datasets, the average tumor localization error and the 95th percentile were 1.14 and 2.78 respectively in two datasets. 4DCBCT motion models are shown to perform well in the context of generating 3D fluoroscopic images due to their ability to reproduce anatomical changes at treatment time. Conclusion: This study showed the feasibility of deriving 4DCBCT-based motion models and using them to generate 3D fluoroscopic images at treatment time in real clinical settings. 4DCBCT-based motion models were found to account for the 3D non-rigid motion of the patient anatomy during treatment and have the potential to localize tumor and other patient anatomical structures at treatment time even when inter-fractional changes occur. This project was supported, in part, through a Master Research Agreement with Varian Medical Systems, Inc., Palo Alto, CA. The project was also supported, in part, by Award Number R21CA156068 from the National Cancer Institute.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-BRC-13: Model Based Classification for Optimal Position Selection for Left-Sided Breast Radiotherapy: Free Breathing, DIBH, Or Prone

Purpose: There are clinical decision challenges to select optimal treatment positions for left-sided breast cancer patients—supine free breathing (FB), supine Deep Inspiration Breath Hold (DIBH) and prone free breathing (prone). Physicians often make the decision based on experiences and trials, which might not always result optimal OAR doses. We herein propose a mathematical model to predict the lowest OAR doses among these three positions, providing a quantitative tool for corresponding clinical decision. Methods: Patients were scanned in FB, DIBH, and prone positions under an IRB approved protocol. Tangential beam plans were generated for each position, and OAR doses were calculated. The position with least OAR doses is defined as the optimal position. The following features were extracted from each scan to build the model: heart, ipsilateral lung, breast volume, in-field heart, ipsilateral lung volume, distance between heart and target, laterality of heart, and dose to heart and ipsilateral lung. Principal Components Analysis (PCA) was applied to remove the co-linearity of the input data and also to lower the data dimensionality. Feature selection, another method to reduce dimensionality, was applied as a comparison. Support Vector Machine (SVM) was then used for classification. Thirtyseven patient data were acquired; up to now, five patient plans were available. K-fold cross validation was used to validate the accuracy of the classifier model with small training size. Results: The classification results and K-fold cross validation demonstrated the model is capable of predicting the optimal position for patients. The accuracy of K-fold cross validations has reached 80%. Compared to PCA, feature selection allows causal features of dose to be determined. This provides more clinical insights. Conclusion: The proposed classification system appeared to be feasible. We are generating plans for the rest of the 37 patient images, and more statistically significant results are to be presented.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TU-D-207B-05: Intra-Tumor Partitioning and Texture Analysis of DCE-MRI Identifies Relevant Tumor Subregions to Predict Early Pathological Response of Breast Cancer to Neoadjuvant Chemotherapy

Purpose: To predict early pathological response of breast cancer to neoadjuvant chemotherapy (NAC) based on quantitative, multi-region analysis of dynamic contrast enhancement magnetic resonance imaging (DCE-MRI). Methods: In this institution review board-approved study, 35 patients diagnosed with stage II/III breast cancer were retrospectively investigated using DCE-MR images acquired before and after the first cycle of NAC. First, principal component analysis (PCA) was used to reduce the dimensionality of the DCE-MRI data with a high-temporal resolution. We then partitioned the whole tumor into multiple subregions using k-means clustering based on the PCA-defined eigenmaps. Within each tumor subregion, we extracted four quantitative Haralick texture features based on the gray-level co-occurrence matrix (GLCM). The change in texture features in each tumor subregion between pre- and during-NAC was used to predict pathological complete response after NAC. Results: Three tumor subregions were identified through clustering, each with distinct enhancement characteristics. In univariate analysis, all imaging predictors except one extracted from the tumor subregion associated with fast wash-out were statistically significant (p< 0.05) after correcting for multiple testing, with area under the ROC curve or AUCs between 0.75 and 0.80. In multivariate analysis, the proposed imaging predictors achieved an AUC of 0.79 (p = 0.002) in leave-one-out cross validation. This improved upon conventional imaging predictors such as tumor volume (AUC=0.53) and texture features based on whole-tumor analysis (AUC=0.65). Conclusion: The heterogeneity of the tumor subregion associated with fast wash-out on DCE-MRI predicted early pathological response to neoadjuvant chemotherapy in breast cancer.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TU-FG-209-12: Treatment Site and View Recognition in X-Ray Images with Hierarchical Multiclass Recognition Models

Purpose: To investigate an approach of automatically recognizing anatomical sites and imaging views (the orientation of the image acquisition) in 2D X-ray images. Methods: A hierarchical (binary tree) multiclass recognition model was developed to recognize the treatment sites and views in x-ray images. From top to bottom of the tree, the treatment sites are grouped hierarchically from more general to more specific. Each node in the hierarchical model was designed to assign images to one of two categories of anatomical sites. The binary image classification function of each node in the hierarchical model is implemented by using a PCA transformation and a support vector machine (SVM) model. The optimal PCA transformation matrices and SVM models are obtained by learning from a set of sample images. Alternatives of the hierarchical model were developed to support three scenarios of site recognition that may happen in radiotherapy clinics, including two or one X-ray images with or without view information. The performance of the approach was tested with images of 120 patients from six treatment sites – brain, head-neck, breast, lung, abdomen and pelvis – with 20 patients per site and two views (AP and RT) per patient. Results: Given two images in known orthogonal views (AP and RT), the hierarchical model achieved a 99% average F1 score to recognize the six sites. Site specific view recognition models have 100 percent accuracy. The computation time to process a new patient case (preprocessing, site and view recognition) is 0.02 seconds. Conclusion: The proposed hierarchical model of site and view recognition is effective and computationally efficient. It could be useful to automatically and independently confirm the treatment sites and views in daily setup x-ray 2D images. It could also be applied to guide subsequent image processing tasks, e.g. site and view dependent contrast enhancement and image registration. The senior author received research grants from ViewRay Inc. and Varian Medical System.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Engineering Novosphingobium aromaticivorans to produce cis,cis -muconic acid from biomass aromatics

ABSTRACT The platform chemical cis,cis- muconic acid ( cc MA) provides facile access to a number of monomers used in the synthesis of commercial plastics. It is also a metabolic intermediate in the β-ketoadipic acid pathway of many bacteria and, therefore, a current target for microbial production from abundant renewable resources via metabolic engineering. This study investigates Novosphingobium aromaticivorans DSM12444 as a chassis for the production of cc MA from biomass aromatics. The N. aromaticivorans genome predicts that it encodes a previously uncharacterized protocatechuic acid (PCA) decarboxylase and a catechol 1,2-dioxygenase, which would be necessary for the conversion of aromatic metabolic intermediates to cc MA. This study confirmed the activity of these two enzymes in vitro and compared their activity to ones that have been previously characterized and used in cc MA production. From these results, we generated one strain that is completely derived from native genes and a second that contains genes previously used in microbial engineering synthesis of this compound. Both of these strains exhibited stoichiometric production of cc MA from PCA and produced greater than 100% yield of cc MA from the aromatic monomers that were identified in liquor derived from alkaline pretreated biomass. Our results show that a strain completely derived from native genes and one containing homologs from other hosts are both capable of stoichiometric production of cc MA from biomass aromatics. Overall, this work combines previously unknown aspects of aromatic metabolism in N. aromaticivorans and the genetic tractability of this organism to generate strains that produce cc MA from deconstructed biomass. IMPORTANCE The production of commodity chemicals from renewable resources is an important goal toward increasing the environmental and economic sustainability of industrial processes. The aromatics in plant biomass are an underutilized and abundant renewable resource for the production of valuable chemicals. However, due to the chemical composition of plant biomass, many deconstruction methods generate a heterogeneous mixture of aromatics, thus making it difficult to extract valuable chemicals using current methods. Therefore, recent efforts have focused on harnessing the pathways of microorganisms to convert a diverse set of aromatics into a single product. Novosphingobium aromaticivorans DSM12444 has the native ability to metabolize a wide range of aromatics and, thus, is a potential chassis for conversion of these abundant compounds to commodity chemicals. This study reports on new features of N. aromaticivorans that can be used to produce the commodity chemical cis,cis -muconic acid from renewable and abundant biomass aromatics.

09 BIOMASS FUELS↗

DOE OSTI2021

Redundancy in Aromatic O -Demethylation and Ring-Opening Reactions in Novosphingobium aromaticivorans and Their Impact in the Metabolism of Plant-Derived Phenolics

ABSTRACT Lignin is a plant heteropolymer composed of phenolic subunits. Because of its heterogeneity and recalcitrance, the development of efficient methods for its valorization still remains an open challenge. One approach to utilize lignin is its chemical deconstruction into mixtures of monomeric phenolic compounds, followed by biological funneling into a single product. Novosphingobium aromaticivorans DSM 12444 has been previously engineered to produce 2-pyrone-4,6-dicarboxylic acid (PDC) from depolymerized lignin by simultaneously metabolizing multiple aromatics through convergent routes involving the intermediates 3-methoxygallic acid (3-MGA) and protocatechuic acid (PCA). We investigated enzymes predicted to be responsible for O -demethylation and oxidative aromatic ring opening, two critical reactions involved in the metabolism of phenolic compounds by N. aromaticivorans . The results showed the involvement of DesA in O -demethylation of syringic and vanillic acids, LigM in O- demethylation of vanillic acid and 3-MGA, and a new O- demethylase, DmtS, in the conversion of 3-MGA into gallic acid (GA). In addition, we found that LigAB was the main aromatic ring-opening dioxygenase involved in 3-MGA, PCA, and GA metabolism and that a previously uncharacterized dioxygenase, LigAB2, had high activity with GA. Our results indicate a metabolic route not previously identified in N. aromaticivorans that involves O -demethylation of 3-MGA to GA. We predict that this pathway channels ∼15% of the carbon flow from syringic acid, with the rest following ring opening of 3-MGA. The new knowledge obtained in this study allowed for the creation of an improved engineered strain for the funneling of aromatic compounds that exhibits stoichiometric conversion of syringic acid into PDC. IMPORTANCE For lignocellulosic biorefineries to effectively contribute to reduction of fossil fuel use, they need to become efficient at producing chemicals from all major components of plant biomass. Making products from lignin will require engineering microorganisms to funnel multiple phenolic compounds to the chemicals of interest, and N. aromaticivorans is a promising chassis for this technology. The ability of N. aromaticivorans to efficiently and simultaneously degrade many phenolic compounds may be linked to having functionally redundant aromatic degradation pathways and enzymes with broad substrate specificity. A detailed knowledge of aromatic degradation pathways is thus essential to identify genetic engineering targets to maximize product yields. Furthermore, knowledge of enzyme substrate specificity is critical to redirect flow of carbon to desired pathways. This study described an uncharacterized pathway in N. aromaticivorans and the enzymes that participate in this pathway, allowing the engineering of an improved strain for production of PDC from lignin.

Biotechnology & Applied Microbiology↗

DOE OSTI2025

Engineered Methylobacterium extorquens grows well on methoxylated aromatics due to its formaldehyde metabolism and stress response

ABSTRACT Lignin is a vast yet underutilized source of renewable energy. The microbial valorization of lignin is challenging due to the toxicity of its degradation intermediates, particularly formaldehyde. In this study, we engineeredMethylobacterium extorquensPA1 to metabolize lignin-derived methoxylated aromatics, vanillate (VA) and protocatechuate (PCA), by introducing thevanandpcagene clusters. Compared toPseudomonas putida,M. extorquensPA1 exhibited better formaldehyde detoxification, enabling robust growth on VA without accumulation of formaldehyde. Genetic analyses confirmed that formaldehyde oxidation and stress response systems, rather than C 1 assimilation, were important for VA metabolism. Additionally, VA and PCA were found to disrupt membrane potential, contributing to their inherent toxicity. Our findings establishM. extorquensPA1 as a promising chassis for lignin valorization and provide a framework for engineering formaldehyde-resistant microbial platforms. IMPORTANCE In developing biotechnological solutions for a circular economy, it is critical to valorize all parts of renewable resources, such as lignocellulose from vegetative components of agricultural crops and from bioenergy feedstocks. After chemical breakdown, the aromatics arising from lignin present significant challenges for use due to their toxicity. Here, we address one component of this challenge—the methoxy groups that get released as formaldehyde—and show that existing biotechnological platform organisms with strong formaldehyde metabolism, such asMethylobacterium extorquens, can be transformed into highly capable utilizers of methoxylated aromatics.

Microbiology↗

DOE OSTI2019

A Novel Sparse Compositional Technique Reveals Microbial Perturbations

ABSTRACT The central aims of many host or environmental microbiome studies are to elucidate factors associated with microbial community compositions and to relate microbial features to outcomes. However, these aims are often complicated by difficulties stemming from high-dimensionality, non-normality, sparsity, and the compositional nature of microbiome data sets. A key tool in microbiome analysis is beta diversity, defined by the distances between microbial samples. Many different distance metrics have been proposed, all with varying discriminatory power on data with differing characteristics. Here, we propose a compositional beta diversity metric rooted in a centered log-ratio transformation and matrix completion called robust Aitchison PCA. We demonstrate the benefits of compositional transformations upstream of beta diversity calculations through simulations. Additionally, we demonstrate improved effect size, classification accuracy, and robustness to sequencing depth over the current methods on several decreased sample subsets of real microbiome data sets. Finally, we highlight the ability of this new beta diversity metric to retain the feature loadings linked to sample ordinations revealing salient intercommunity niche feature importance. IMPORTANCE By accounting for the sparse compositional nature of microbiome data sets, robust Aitchison PCA can yield high discriminatory power and salient feature ranking between microbial niches. The software to perform this analysis is available under an open-source license and can be obtained at https://github.com/biocore/DEICODE ; additionally, a QIIME 2 plugin is provided to perform this analysis at https://library.qiime2.org/plugins/q2-deicode .

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Primary Modes of Global Drop Size Distributions

Understanding drop size distribution (DSD) variability has important implications for remote sensing and numerical modeling applications. Twelve disdrometer datasets across three latitude bands are analyzed in this study, spanning a broad range of precipitation regimes: light rain, orographic, deep convective, organized midlatitude, and tropical oceanic. Principal component analysis (PCA) is used to reveal comprehensive modes of global DSD spatial and temporal variability. Although the locations contain different distributions of individual DSD parameters, all locations are found to have the same modes of variability. Based on PCA, six groups of points with unique DSD characteristics emerge. The physical processes that underpin these groups are revealed through supporting radar observations. Group 1 (group 2) is characterized by high (low) liquid water content (LWC), broad (narrow) distribution widths, and large (small) median drop diameters D 0 . Radar analysis identifies group 1 (group 2) as convective (stratiform) rainfall. Group 3 is characterized by weak, shallow radar echoes and large concentrations of small drops, indicative of warm rain showers. Group 4 identifies heavy stratiform precipitation. The low latitudes exhibit distinct bimodal distributions of the normalized intercept parameter N w , LWC, and D 0 and are found to have a clustering of points (group 5) with high rain rates, large N w , and moderate D 0 , a signature of robust warm rain processes. A distinct group associated with ice-based convection (group 6) emerges in the midlatitudes. Although all locations exhibit the same covariance of parameters associated with these groups, it is likely that the physical processes responsible for shaping the DSDs vary as a function of location.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Removing cosmic spikes using a hyperspectral upper-bound spectrum method

Cosmic ray spikes are especially problematic for hyperspectral imaging, due to the large number of spikes often present and their negative effects upon subsequent chemometric analysis. Fortunately, while the large number of spectra acquired in a hyperspectral imaging data set increases the probability and number of cosmic spikes observed, the multitude of spectra can also aid in the effective recognition and removal of the cosmic spikes. Dongmao Zhang and Dor Ben-Amotz were perhaps the first to leverage the additional spatial dimension of hyperspectral data matrices (DM). They integrated principal component analysis (PCA) into the upper bound spectrum method (UBS), resulting in a hybrid method (UBS-DM) for hyperspectral images. Here, we expand upon their use of PCA, recognizing that principal components primarily present in only a few pixels most likely correspond to cosmic spikes. Eliminating the contribution of those principal components in those pixels improves the cosmic spike removal. Both simulated and experimental hyperspectral Raman image data sets are used to test the newly developed UBS-DM-hyperspectral (UBS-DM-HS) method which extends the UBS-DM method by leveraging characteristics of hyperspectral datasets. Lastly, a comparison is provided between the performance of the UBS-DM-HS method and other methods suitable for despiking hyperspectral images, evaluating both their ability to remove cosmic ray spikes and the extent to which they introduce spectral bias.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗