ABSTRACT
Sequencing the tyrosine phosphoproteome using MS-based proteomics is challenging due to the low abundance of tyrosine phosphorylation in cells, a challenge compounded in scarce samples like primary cells or clinical samples. The broad-spectrum optimisation of selective triggering (BOOST) method was recently developed to increase phosphotyrosine sequencing in low protein input samples by leveraging tandem mass tags (TMT), phosphotyrosine enrichment, and a phosphotyrosine-loaded carrier channel. Here, we demonstrate the viability of BOOST in T cell receptor (TCR)-stimulated primary murine T cells by benchmarking the accuracy and precision of the BOOST method and discerning significant alterations in the phosphoproteome associated with receptor stimulation. Using 1 mg of protein input (about 20 million cells) and BOOST, we identify and precisely quantify more than 2000 unique pY sites compared to about 300 unique pY sites in non-BOOST control samples. We show that although replicate variation increases when using the BOOST method, BOOST does not jeopardise quantitative precision or the ability to determine statistical significance for peptides measured in triplicate. Many pY previously uncharacterised sites on important T cell signalling proteins are quantified using BOOST, and we identify new TCR responsive pY sites observable only with BOOST. Finally, we determine that the phase-spectrum deconvolution method on Orbitrap instruments can impair pY quantitation in BOOST experiments.
ABSTRACT
DNA methylation (DNAm)-based deconvolution estimates contain relative data, forming a composition, that standard methods (testing directly on cell proportions) are ill-suited to handle. In this study we examined the performance of an alternative method, analysis of compositions of microbiomes (ANCOM), for the analysis of DNAm-based deconvolution estimates. We performed two different simulation studies comparing ANCOM to a standard approach (two sample t-test performed directly on cell proportions) and analyzed a real-world data from the Women's Health Initiative to evaluate the applicability of ANCOM to DNAm-based deconvolution estimates. Our findings indicate that ANCOM can effectively account for the compositional nature of DNAm-based deconvolution estimates. ANCOM adequately controls the false discovery rate while maintaining statistical power comparable to that of standard methods.
DNA methylation (DNAm)-based deconvolution provides highly accurate estimates of the proportion of each cell type in a mixed-cell type biological sample (e.g., whole-blood). These estimates can be used for examining the association between cell type proportions and biological or clinical end points; for example, comparing the estimated neutrophil proportion in whole blood between smokers and non-smokers. Cell proportion data has unique features which present challenges for traditional and widely used statistical methods. In response to this issue, our work presents two simulation studies and a real-world analysis that benchmark the performance of current standard statistical methods against an alternative method called analysis composition of microbes (ANCOM), which was originally developed for the analysis of microbiome data. In our real-world analysis we used DNAm data collected from Women's Health Initiative Long Life Study I and compared the results of each method against a gold-standard that is typically not available for these analyses. In each of our simulation studies, ANCOM was able to detect true differences in cell proportions between the groups being compared but had a much lower rate of false discovery compared with the standard statistical methods. Our real-world analysis demonstrated similar findings. Overall, our study highlights the potential of ANCOM as a powerful and robust method for analyzing DNAm-derived deconvolution estimates when the interest is comparisons of cell type proportions and biological or clinical end points. ANCOM's ability to minimize false discovery while maintaining robust statistical power positions it as a valuable addition to the epigenomic analysis toolkit.
ABSTRACT
BACKGROUND: Although the clinical signs of inflammatory breast cancer (IBC) resemble acute inflammation, the role played by infiltrating immune and stromal cells in this aggressive disease is uncharted. The tumor microenvironment (TME) presents molecular alterations, such as epimutations, prior to morphological abnormalities. These changes affect the distribution and the intricate communication between the TME components related to cancer prognosis and therapy response. Herein, we explored the global DNA methylation profile of IBC and surrounding tissues to estimate the microenvironment cellular composition and identify epigenetically dysregulated markers. METHODS: We used the HiTIMED algorithm to deconvolve the bulk DNA methylation data of 24 IBC and six surrounding non-tumoral tissues (SNT) (GSE238092) and determine their cellular composition. The prognostic relevance of cell types infiltrating IBC and their relationship with clinicopathological variables were investigated. CD34 (endothelial cell marker) and CD68 (macrophage marker) immunofluorescence staining was evaluated in an independent set of 17 IBC and 16 non-IBC samples. RESULTS: We found lower infiltration of endothelial, stromal, memory B, dendritic, and natural killer cells in IBC than in SNT samples. Higher endothelial cell (EC) and stromal cell content were related to better overall survival. EC proportions positively correlated with memory B and memory CD8+ T infiltration in IBC. Immune and EC markers exhibited distinct DNA methylation profiles between IBC and SNT samples, revealing hypermethylated regions mapped to six genes (CD40, CD34, EMCN, HLA-G, PDPN, and TEK). We identified significantly higher CD34 and CD68 protein expression in IBC compared to non-IBC. CONCLUSIONS: Our findings underscored cell subsets that distinguished patients with better survival and dysregulated markers potentially actionable through combinations of immunotherapy and epigenetic drugs.
Subject(s)
DNA Methylation , Inflammatory Breast Neoplasms , Tumor Microenvironment , Humans , DNA Methylation/genetics , Tumor Microenvironment/genetics , Female , Inflammatory Breast Neoplasms/genetics , Inflammatory Breast Neoplasms/pathology , Inflammatory Breast Neoplasms/metabolism , Treatment Outcome , Middle Aged , Prognosis , Molecular Targeted Therapy , Gene Expression Regulation, NeoplasticABSTRACT
PURPOSE: The measurement of slice sensitivity profile (SSP) in non-helical CT is conventionally performed by repeated scans with moving a micro-coin phantom little by little in the longitudinal direction at a small interval, which is reliable but laborious and time-consuming. The purpose of this study was to propose a simple method for measuring the SSP in non-helical CT based on a previous method that measured the slice thickness using a tilted metal wire. METHODS: In the proposed method, a CT image was obtained by scanning a wire tilted at an angle θ=30° to the scan plane. By deconvolving the image with the point spread function (PSF) measured at the scanner, we obtained an image that was not affected by the PSF blurring. The CT value profile along the wire was acquired on the obtained image. The SSP was determined by multiplying the profile by tan θ. In addition, the SSP was measured by the conventional method using a micro-coin phantom and compared with the SSP obtained by the proposed method. RESULTS: The SSP measured by the proposed method agreed well with that measured by the conventional method. The full-width at half-maximum values of these SSPs were the same. CONCLUSION: The proposed method was demonstrated to easily and accurately measure the SSP in non-helical CT.
ABSTRACT
Recent technological advancements have enabled spatially resolved transcriptomic profiling but at a multicellular resolution that is more cost-effective. The task of cell type deconvolution has been introduced to disentangle discrete cell types from such multicellular spots. However, existing benchmark datasets for cell type deconvolution are either generated from simulation or limited in scale, predominantly encompassing data on mice and are not designed for human immuno-oncology. To overcome these limitations and promote comprehensive investigation of cell type deconvolution for human immuno-oncology, we introduce a large-scale spatial transcriptomic deconvolution benchmark dataset named SpatialCTD, encompassing 1.8 million cells and 12,900 pseudo spots from the human tumor microenvironment across the lung, kidney, and liver. In addition, SpatialCTD provides more realistic reference than those generated from single-cell RNA sequencing (scRNA-seq) data for most reference-based deconvolution methods. To utilize the location-aware SpatialCTD reference, we propose a graph neural network-based deconvolution method (i.e., GNNDeconvolver). Extensive experiments show that GNNDeconvolver often outperforms existing state-of-the-art methods by a substantial margin, without requiring scRNA-seq data. To enable comprehensive evaluations of spatial transcriptomics data from flexible protocols, we provide an online tool capable of converting spatial transcriptomic data from various platforms (e.g., 10× Visium, MERFISH, and sci-Space) into pseudo spots, featuring adjustable spot size. The SpatialCTD dataset and GNNDeconvolver implementation are available at https://github.com/OmicsML/SpatialCTD, and the online converter tool can be accessed at https://omicsml.github.io/SpatialCTD/.
ABSTRACT
Immunologic treatment options are uncommon in low-grade gliomas, although such therapies might be beneficial for inoperable and aggressive cases. Knowledge of the immune and stromal cells in low-grade gliomas is highly relevant for such approaches but still needs to be improved. Published gene-expression data from 400 low-grade gliomas and 193 high-grade gliomas were gathered to quantify 10 microenvironment cell populations with a deconvolution method designed explicitly for brain tumors. First, we investigated general differences in the microenvironment of low- and high-grade gliomas. Lower-grade and high-grade tumors cluster together, respectively, and show a general similarity within and distinct differences between these groups, the main difference being a higher infiltration of fibroblasts and T cells in high-grade gliomas. Among the analyzed entities, gangliogliomas and pleomorphic xanthoastrocytomas presented the highest overall immune cell infiltration. Further analyses of the low-grade gliomas presented three distinct microenvironmental signatures of immune cell infiltration, which can be divided into T-cell/dendritic/natural killer cell-, neutrophilic/B lineage/natural killer cell-, and monocytic/vascular/stromal-cell-dominated immune clusters. These clusters correlated with tumor location, age, and histological diagnosis but not with sex or progression-free survival. A survival analysis showed that the prognosis can be predicted from gene expression, clinical data, and a combination of both with a support vector machine and revealed the negative prognostic relevance of vascular markers. Overall, our work shows that low- and high-grade gliomas can be characterized and differentiated by their immune cell infiltration. Low-grade gliomas cluster into three distinct immunologic tumor microenvironments, which may be of further interest for upcoming immunotherapeutic research.
Subject(s)
Brain Neoplasms , Glioma , Tumor Microenvironment , Humans , Tumor Microenvironment/immunology , Tumor Microenvironment/genetics , Glioma/genetics , Glioma/immunology , Glioma/pathology , Brain Neoplasms/genetics , Brain Neoplasms/pathology , Brain Neoplasms/immunology , Child , Female , Male , Neoplasm Grading , Gene Expression Profiling , Transcriptome , Child, Preschool , Adolescent , Stromal Cells/pathology , Stromal Cells/metabolism , Stromal Cells/immunologyABSTRACT
Spatial transcriptomics technology has been an essential and powerful method for delineating tissue architecture at the molecular level. However, due to the limitations of the current spatial techniques, the cellular information cannot be directly measured but instead spatial spots typically varying from a diameter of 0.2 to 100 µm are characterized. Therefore, it is vital to apply computational strategies for inferring the cellular composition within each spatial spot. The main objective of this review is to summarize the most recent progresses to estimate the exact cellular proportions for each spatial spot, and to prospect the future directions of this field.
ABSTRACT
Spatially resolved transcriptomics integrates high-throughput transcriptome measurements with preserved spatial cellular organization information. However, many technologies cannot reach single-cell resolution. We present STdGCN, a graph model leveraging single-cell RNA sequencing (scRNA-seq) as reference for cell-type deconvolution in spatial transcriptomic (ST) data. STdGCN incorporates expression profiles from scRNA-seq and spatial localization from ST data for deconvolution. Extensive benchmarking on multiple datasets demonstrates that STdGCN outperforms 17 state-of-the-art models. In a human breast cancer Visium dataset, STdGCN delineates stroma, lymphocytes, and cancer cells, aiding tumor microenvironment analysis. In human heart ST data, STdGCN identifies changes in endothelial-cardiomyocyte communications during tissue development.
Subject(s)
Single-Cell Analysis , Transcriptome , Humans , Single-Cell Analysis/methods , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Tumor Microenvironment , Gene Expression Profiling/methods , RNA-Seq/methods , Sequence Analysis, RNA/methodsABSTRACT
Objective Time-of-flight (TOF) scatter rejection requires a total timing jitter, including the detector timing jitter and the X-ray source's pulses width, of 50 ps or less to mitigate most of the effects of scattered photons in radiography and CT imaging. However, since the total contribution of the source and detector to the timing jitter can be retrieved during an acquisition with nothing between the source and detector, it can be demonstrated that this contribution may be partially removed to improve the image quality. Approach A scatter correction method using iterative deconvolution of the measured time point-spread function estimates the number of scattered photons detected in each pixel. To evaluate the quality of the estimation, GATE was used to simulate the radiography of a water cylinder with bone inserts, and a head and torso in a system with total timing jitters from 100 ps up to 500 ps full-width-at-half-maximum. Main results With a total timing jitter of 200 ps, 89% of the contrast degradation caused by scattered photons was recovered in a head and torso radiography, compared to 28% with a simple time threshold method. Corrected images using the estimation have a percent root-mean square error between 2 and 14% in both phantoms with timing jitters from 100 to 500 ps which is lower than the error achieved with scatter rejection alone at 100 ps. Significance TOF X-ray imaging has the potential to mitigate the effects of the scattering contribution and offers an alternative to anti-scatter grids that avoids loss of primary photons. Compare to simple TOF scatter rejection using only a threshold, the deconvolution estimation approach has lower requirements on both the source and detector. These requirements are now within reach of state-of-the-art systems.
ABSTRACT
Background: This study introduces a novel prognostic tool, the Disulfidoptosis-Related lncRNA Index (DRLI), integrating the molecular signatures of disulfidoptosis and long non-coding RNAs (lncRNAs) with the cellular heterogeneity of the tumor microenvironment, to predict clinical outcomes in patients with clear cell renal cell carcinoma (ccRCC). Methods: We analyzed 530 tumor and 72 normal samples from The Cancer Genome Atlas (TCGA), employing k-means clustering based on disulfidoptosis-associated gene expression to stratify ccRCC samples into prognostic groups. lncRNAs correlated with disulfidoptosis were identified and used to construct the DRLI, which was validated by Kaplan-Meier and receiver operating characteristic curves. We utilized single-cell deconvolution analysis to estimate the proportion of immune cell types within the tumor microenvironment, while the ESTIMATE and TIDE algorithms were employed to assess immune infiltration and potential response to immunotherapy. Results: The Disulfidoptosis-Related lncRNA Index (DRLI) effectively stratified ccRCC patients into high and low-risk groups, significantly impacting survival outcomes (P < 0.001). High-risk patients, marked by a unique lncRNA profile associated with disulfidoptosis, faced worse prognoses. Single-cell analysis revealed marked tumor microenvironment heterogeneity, especially in immune cell makeup, correlating with patient risk levels. In prognostic predictions, DRLI outperformed traditional clinical indicators, achieving AUC values of 0.779, 0.757, and 0.779 for 1-year, 3-year, and 5-year survival in the training set, and 0.746, 0.734, and 0.750 in the validation set. Notably, while the constructed nomogram showed exceptional predictive capability for short-term prognosis (AUC = 0.877), the DRLI displayed remarkable long-term predictive accuracy, with its AUC value reaching 0.823 for 10-year survival, closely approaching the nomogram's performance. Conclusions: The study introduces the DRLI as a groundbreaking molecular stratification tool for ccRCC, enhancing prognostic precision and potentially guiding personalized treatment strategies. This advancement is particularly significant in the context of long-term survival predictions. Our findings also elucidate the complex interplay between disulfidoptosis, lncRNAs, and the immune microenvironment in ccRCC, offering a comprehensive perspective on its pathogenesis and progression. The DRLI and the nomogram together represent significant strides in ccRCC research, highlighting the importance of molecular-based assessments in predicting patient outcomes.
ABSTRACT
Inferring cell type proportions from bulk transcriptome data is crucial in immunology and oncology. Here, we introduce guided LDA deconvolution (GLDADec), a bulk deconvolution method that guides topics using cell type-specific marker gene names to estimate topic distributions for each sample. Through benchmarking using blood-derived datasets, we demonstrate its high estimation performance and robustness. Moreover, we apply GLDADec to heterogeneous tissue bulk data and perform comprehensive cell type analysis in a data-driven manner. We show that GLDADec outperforms existing methods in estimation performance and evaluate its biological interpretability by examining enrichment of biological processes for topics. Finally, we apply GLDADec to The Cancer Genome Atlas tumor samples, enabling subtype stratification and survival analysis based on estimated cell type proportions, thus proving its practical utility in clinical settings. This approach, utilizing marker gene names as partial prior information, can be applied to various scenarios for bulk data deconvolution. GLDADec is available as an open-source Python package at https://github.com/mizuno-group/GLDADec.
Subject(s)
Software , Humans , Gene Expression Profiling/methods , Algorithms , Transcriptome , Computational Biology/methods , Neoplasms/genetics , Biomarkers, Tumor/genetics , Genetic MarkersABSTRACT
Peptidoglycan is a major and essential component of the bacterial cell envelope that confers cell shape and provides protection against internal osmotic pressure. This complex macromolecule is made of glycan strands cross-linked by short peptides, and its structure is continually modified throughout growth via a process referred to as "remodeling." Peptidoglycan remodeling allows cells to grow, adapt to their environment, and release fragments that can act as signaling molecules during host-pathogen interactions. Preparing peptidoglycan samples for structural analysis first requires purification of the peptidoglycan sacculus, followed by its enzymatic digestion into disaccharide peptides (muropeptides). These muropeptides can then be characterized by liquid chromatography coupled mass spectrometry (LC-MS) and used to infer the structure of intact peptidoglycan sacculi. Due to the presence of unusual crosslinks, noncanonical amino acids, and amino sugars, the analysis of peptidoglycan LC-MS datasets cannot be handled by traditional proteomics software. In this chapter, we describe a protocol to perform the analysis of peptidoglycan LC-MS datasets using the open-source software PGFinder. We provide a step-by-step strategy to deconvolute data from various mass spectrometry instruments, generate muropeptide databases, perform a PGFinder search, and process the data output.
Subject(s)
Peptidoglycan , Software , Peptidoglycan/chemistry , Peptidoglycan/metabolism , Peptidoglycan/analysis , Chromatography, Liquid/methods , Mass Spectrometry/methods , Glycomics/methods , Proteomics/methods , Bacteria/metabolism , Bacteria/chemistry , Liquid Chromatography-Mass SpectrometryABSTRACT
Monitoring of volatile organic compounds (VOCs) in air is crucial for understanding their atmospheric impacts and advancing their emission reduction plans. This study presents an innovative integrated methodology suitable for achieving semireal-time high spatiotemporal resolution three-dimensional measurements of VOCs from ground to hundreds of meters above ground. The methodology integrates an active AirCore sampler, custom-designed for deployment from unmanned aerial vehicles (UAV), a proton-transfer-reaction mass spectrometry (PTR-MS) for sample analysis, and a data deconvolution algorithm for improved time resolution for measurements of multiple VOCs in air. The application of the deconvolution technique significantly improves the signal strength of data from PTR-MS analysis of AirCore samples and enhances their temporal resolution by 4 to 8 times to 4-11 s. A case study demonstrates that the methodology can achieve sample collection and analysis of VOCs within 45 min, resulting in >120-360 spatially resolved data points for each VOC measured and achieving a horizontal resolution of 20-55 m at a UAV flight speed of 5 m/s and a vertical resolution of 5 m. This methodology presents new possibilities for acquiring 3-dimensional spatial distributions of VOC concentrations, effectively tackling the longstanding challenge of characterizing three-dimensional VOC distributions in the lowest portion of the atmospheric boundary layer.
Subject(s)
Air Pollutants , Environmental Monitoring , Volatile Organic Compounds , Volatile Organic Compounds/analysis , Environmental Monitoring/methods , Air Pollutants/analysis , Mass Spectrometry/methods , Algorithms , AircraftABSTRACT
Purpose: The study aimed to (1) assess the feasibility constrained spherical deconvolution (CSD) tractography to reconstruct crossing fiber bundles with unsedated neonatal diffusion MRI (dMRI), and (2) demonstrate the impact of spatial and angular resolution and processing settings on tractography and derived quantitative measures. Methods: For the purpose of this study, the term-equivalent dMRIs (single-shell b800, and b2000, both 5 b0, and 45 gradient directions) of two moderate-late preterm infants (with and without motion artifacts) from a local cohort [Brain Imaging in Moderate-late Preterm infants (BIMP) study; Calgary, Canada] and one infant from the developing human connectome project with high-quality dMRI (using the b2600 shell, comprising 20 b0 and 128 gradient directions, from the multi-shell dataset) were selected. Diffusion tensor imaging (DTI) and CSD tractography were compared on b800 and b2000 dMRI. Varying image resolution modifications, (pre-)processing and tractography settings were tested to assess their impact on tractography. Each experiment involved visualizing local modeling and tractography for the corpus callosum and corticospinal tracts, and assessment of morphological and diffusion measures. Results: Contrary to DTI, CSD enabled reconstruction of crossing fibers. Tractography was susceptible to image resolution, (pre-) processing and tractography settings. In addition to visual variations, settings were found to affect streamline count, length, and diffusion measures (fractional anisotropy and mean diffusivity). Diffusion measures exhibited variations of up to 23%. Conclusion: Reconstruction of crossing fiber bundles using CSD tractography with unsedated neonatal dMRI data is feasible. Tractography settings affected streamline reconstruction, warranting careful documentation of methods for reproducibility and comparison of cohorts.
ABSTRACT
BACKGROUND: Computational cell type deconvolution enables the estimation of cell type abundance from bulk tissues and is important for understanding tissue microenviroment, especially in tumor tissues. With rapid development of deconvolution methods, many benchmarking studies have been published aiming for a comprehensive evaluation for these methods. Benchmarking studies rely on cell-type resolved single-cell RNA-seq data to create simulated pseudobulk datasets by adding individual cells-types in controlled proportions. RESULTS: In our work, we show that the standard application of this approach, which uses randomly selected single cells, regardless of the intrinsic difference between them, generates synthetic bulk expression values that lack appropriate biological variance. We demonstrate why and how the current bulk simulation pipeline with random cells is unrealistic and propose a heterogeneous simulation strategy as a solution. The heterogeneously simulated bulk samples match up with the variance observed in real bulk datasets and therefore provide concrete benefits for benchmarking in several ways. We demonstrate that conceptual classes of deconvolution methods differ dramatically in their robustness to heterogeneity with reference-free methods performing particularly poorly. For regression-based methods, the heterogeneous simulation provides an explicit framework to disentangle the contributions of reference construction and regression methods to performance. Finally, we perform an extensive benchmark of diverse methods across eight different datasets and find BayesPrism and a hybrid MuSiC/CIBERSORTx approach to be the top performers. CONCLUSIONS: Our heterogeneous bulk simulation method and the entire benchmarking framework is implemented in a user friendly package https://github.com/humengying0907/deconvBenchmarking and https://doi.org/10.5281/zenodo.8206516 , enabling further developments in deconvolution methods.
Subject(s)
Benchmarking , Single-Cell Analysis , Single-Cell Analysis/methods , Humans , Computer Simulation , RNA-Seq/methods , Computational Biology/methodsABSTRACT
Mass spectrometry-based omics technologies are increasingly used in perturbation studies to map drug effects to biological pathways by identifying significant molecular events. Significance is influenced by fold change and variation of each molecular parameter, but also by multiple testing corrections. While the fold change is largely determined by the biological system, the variation is determined by experimental workflows. Here, it is shown that memory effects of prior subculture can influence the variation of perturbation profiles using the two colon carcinoma cell lines SW480 and HCT116. These memory effects are largely driven by differences in growth states that persist into the perturbation experiment. In SW480 cells, memory effects combined with moderate treatment effects amplify the variation in multiple omics levels, including eicosadomics, proteomics, and phosphoproteomics. With stronger treatment effects, the memory effect was less pronounced, as demonstrated in HCT116 cells. Subculture homogeneity was controlled by real-time monitoring of cell growth. Controlled homogeneous subculture resulted in a perturbation network of 321 causal conjectures based on combined proteomic and phosphoproteomic data, compared to only 58 causal conjectures without controlling subculture homogeneity in SW480 cells. Some cellular responses and regulatory events were identified that extend the mode of action of arsenic trioxide (ATO) only when accounting for these memory effects. Controlled prior subculture led to the finding of a synergistic combination treatment of ATO with the thioredoxin reductase 1 inhibitor auranofin, which may prove useful in the management of NRF2-mediated resistance mechanisms.
Subject(s)
Proteomics , Humans , Proteomics/methods , Cell Line, Tumor , HCT116 Cells , Cell Culture Techniques/methods , Colonic Neoplasms/metabolism , Colonic Neoplasms/drug therapy , Colonic Neoplasms/pathology , Arsenic Trioxide/pharmacology , Auranofin/pharmacology , Cell Proliferation/drug effects , Mass Spectrometry/methodsABSTRACT
During the SARS-CoV-2 pandemic, genome-based wastewater surveillance sequencing has been a powerful tool for public health to monitor circulating and emerging viral variants. As a medium, wastewater is very complex because of its mixed matrix nature, which makes the deconvolution of wastewater samples more difficult. Here we introduce a gold standard dataset constructed from synthetic viral control mixtures of known composition, spiked into a wastewater RNA matrix and sequenced on the Oxford Nanopore Technologies platform. We compare the performance of eight of the most commonly used deconvolution tools in identifying SARS-CoV-2 variants present in these mixtures. The software evaluated was primarily chosen for its relevance to the CDC wastewater surveillance reporting protocol, which until recently employed a pipeline that incorporates results from four deconvolution methods: Freyja, kallisto, Kraken 2/Bracken, and LCS. We also tested Lollipop, a deconvolution method used by the Swiss SARS-CoV-2 Sequencing Consortium, and three additional methods not used in the C-WAP pipeline: lineagespot, Alcov, and VaQuERo. We found that the commonly used software Freyja outperformed the other CDC pipeline tools in correct identification of lineages present in the control mixtures, and that the VaQuERo method was similarly accurate, with minor differences in the ability of the two methods to avoid false negatives and suppress false positives. Our results also provide insight into the effect of the tiling primer scheme and wastewater RNA extract matrix on viral sequencing and data deconvolution outcomes.
Subject(s)
COVID-19 , SARS-CoV-2 , Wastewater , SoftwareABSTRACT
BACKGROUND: Changing cell-type proportions can confound studies of differential gene expression or DNA methylation (DNAm) from peripheral blood mononuclear cells (PBMCs). We examined how cell-type proportions derived from the transcriptome versus the methylome (DNAm) influence estimates of differentially expressed genes (DEGs) and differentially methylated positions (DMPs). METHODS: Transcriptome and DNAm data were obtained from PBMC RNA and DNA of Kenyan children (n = 8) before, during, and 6 weeks following uncomplicated malaria. DEGs and DMPs between time points were detected using cell-type adjusted modeling with Cibersortx or IDOL, respectively. RESULTS: Most major cell types and principal components had moderate to high correlation between the two deconvolution methods (r = 0.60-0.96). Estimates of cell-type proportions and DEGs or DMPs were largely unaffected by the method, with the greatest discrepancy in the estimation of neutrophils. CONCLUSION: Variation in cell-type proportions is captured similarly by both transcriptomic and methylome deconvolution methods for most major cell types.
ABSTRACT
Background: The recently launched DNA methylation profiling platform, Illumina MethylationEPIC BeadChip Infinium microarray v2.0 (EPICv2), is highly correlated with measurements obtained from its predecessor MethylationEPIC BeadChip Infinium microarray v1.0 (EPICv1). However, the concordance between the two versions in the context of DNA methylation-based tools, including cell type deconvolution algorithms, epigenetic clocks, and inflammation and lifestyle biomarkers has not yet been investigated. Findings: We profiled DNA methylation on both EPIC versions using matched venous blood samples from individuals spanning early to late adulthood across three cohorts. On combining the DNA methylomes of the cohorts, we observed that samples primarily clustered by the EPIC version they were measured on. Within each cohort, when we calculated cell type proportions, epigenetic age acceleration (EAA), rate of aging estimates, and biomarker scores for the matched samples on each version, we noted significant differences between EPICv1 and EPICv2 in the majority of these estimates. These differences were not significant, however, when estimates were adjusted for EPIC version or when EAAs were calculated separately for each EPIC version. Conclusions: Our findings indicate that EPIC version differences predominantly explain DNA methylation variation and influence estimates of DNA methylation-based tools, and therefore we recommend caution when combining cohorts run on different versions. We demonstrate the importance of calculating DNA methylation-based estimates separately for each EPIC version or accounting for EPIC version either as a covariate in statistical models or by using version correction algorithms.
ABSTRACT
A variety of enabling formulations has been developed to address poor oral drug absorption caused by insufficient dissolution in the gastrointestinal tract. As the in vivo performance of these formulations is a result of a complex interplay between dissolution, digestion and permeation, development of suitable in vitro assays that captures these phenomena are called for. The enabling-absorption (ENA) device, consisting of a donor and receiver chamber separated by a semipermeable membrane, has successfully been used to study the performance of lipid-based formulations. In this work, the ENA device was prepared with two different setups (a Caco-2 cell monolayer and an artificial lipid membrane) to study the performance of a lipid-based formulation (LBF), an amorphous solid dispersion (ASD) and the potential benefit of combining the two formulation strategies. An in vivo pharmacokinetic study in rats was performed to evaluate the in vitro-in vivo correlation. In the ENA, high drug concentrations in the donor chamber did not translate to a high mass transfer, which was particularly evident for the ASD as compared to the LBF. The solubility of the polymer used in the ASD was strongly affected by pH-shifts in vitro, and the ph_dependence resulted in poor in vivo performance of the formulation. The dissolution was however increased in vitro when the ASD was combined with a blank lipid-based formulation. This beneficial effect was also observed in vivo, where the drug exposure of the ASD increased significantly when the ASD was co-administered with the blank LBF. To conclude, the in vitro model managed to capture solubility limitations and strategies to overcome these for one of the formulations studied. The correlation between the in vivo exposure of the drug exposure and AUC in the ENA was good for the non pH-sensitive formulations. The deconvoluted pharmacokinetic data indicated that the receiver chamber was a better predictor for the in vivo performance of the drug, however both chambers provided valuable insights to the observed outcome in vivo. This shows that the advanced in vitro setting used herein successfully could explain absorption differences of highly complex formulations.