Removing Batch Effects from Longitudinal Gene Expression - Quantile Normalization Plus ComBat as Best Approach for Microarray Transcriptome Data

Review on radiomic analysis in 18F-fluorodeoxyglucose positron emission tomography for prediction of melanoma outcomes

Over the past decade, several strategies have revolutionized the clinical management of patients with cutaneous melanoma (CM), including immunotherapy and targeted tyrosine kinase inhibitor (TKI)-based therapies. Indeed, immune checkpoint inhibitors (ICIs), alone or in combination, represent the standard of care for patients with advanced disease without an actionable mutation. Notably BRAF combined with MEK inhibitors represent the therapeutic standard for disease disclosing BRAF mutation. At the same time, FDG PET/CT has become part of the routine staging and evaluation of patients with cutaneous melanoma. There is growing interest in using FDG PET/CT measurements to predict response to ICI therapy and/or target therapy. While semiquantitative values such as standardized uptake value (SUV) are limited for predicting outcome, new measures including tumor metabolic volume, total lesion glycolysis and radiomics seem promising as potential imaging biomarkers for nuclear medicine. The aim of this review, prepared by an interdisciplinary group of experts, is to take stock of the current literature on radiomics approaches that could improve outcomes in CM.

Quantification Quality Control Emerges as a Crucial Factor to Enhance Single-Cell Proteomics Data Analysis

Mass spectrometry (MS)-based single-cell proteomics (SCP) provides us the opportunity to unbiasedly explore biological variability within cells without the limitation of antibody availability. This field is rapidly developed with the main focuses on instrument advancement, sample preparation refinement, and signal boosting methods; however, the optimal data processing and analysis are rarely investigated which holds an arduous challenge because of the high proportion of missing values and batch effect. Here, we introduced a quantification quality control to intensify the identification of differentially expressed proteins (DEPs) by considering both within and across SCP data. Combining quantification quality control with isobaric matching between runs (IMBR) and PSM-level normalization, an additional 12% and 19% of proteins and peptides, with more than 90% of proteins/peptides containing valid values, were quantified. Clearly, quantification quality control was able to reduce quantification variations and q-values with the more apparent cell type separations. In addition, we found that PSM-level normalization performed similar to other protein-level normalizations but kept the original data profiles without the additional requirement of data manipulation. In proof of concept of our refined pipeline, six uniquely identified DEPs exhibiting varied fold-changes and playing critical roles for melanoma and monocyte functionalities were selected for validation using immunoblotting. Five out of six validated DEPs showed an identical trend with the SCP dataset, emphasizing the feasibility of combining the IMBR, cell quality control, and PSM-level normalization in SCP analysis, which is beneficial for future SCP studies.

Comparison of the effectiveness of different normalization methods for metagenomic cross-study phenotype prediction under heterogeneity

The human microbiome, comprising microorganisms residing within and on the human body, plays a crucial role in various physiological processes and has been linked to numerous diseases. To analyze microbiome data, it is essential to account for inherent heterogeneity and variability across samples. Normalization methods have been proposed to mitigate these variations and enhance comparability. However, the performance of these methods in predicting binary phenotypes remains understudied. This study systematically evaluates different normalization methods in microbiome data analysis and their impact on disease prediction. Our findings highlight the strengths and limitations of scaling, compositional data analysis, transformation, and batch correction methods. Scaling methods like TMM show consistent performance, while compositional data analysis methods exhibit mixed results. Transformation methods, such as Blom and NPN, demonstrate promise in capturing complex associations. Batch correction methods, including BMC and Limma, consistently outperform other approaches. However, the influence of normalization methods is constrained by population effects, disease effects, and batch effects. These results provide insights for selecting appropriate normalization approaches in microbiome research, improving predictive models, and advancing personalized medicine. Future research should explore larger and more diverse datasets and develop tailored normalization strategies for microbiome data analysis.

Label-aware distance mitigates temporal and spatial variability for clustering and visualization of single-cell gene expression data

Clustering and visualization are essential parts of single-cell gene expression data analysis. The Euclidean distance used in most distance-based methods is not optimal. The batch effect, i.e., the variability among samples gathered from different times, tissues, and patients, introduces large between-group distance and obscures the true identities of cells. To solve this problem, we introduce Label-Aware Distance (LAD), a metric using temporal/spatial locality of the batch effect to control for such factors. We validate LAD on simulated data as well as apply it to a mouse retina development dataset and a lung dataset. We also found the utility of our approach in understanding the progression of the Coronavirus Disease 2019 (COVID-19). LAD provides better cell embedding than state-of-the-art batch correction methods on longitudinal datasets. It can be used in distance-based clustering and visualization methods to combine the power of multiple samples to help make biological findings.

Short-term predictor for COVID-19 severity from a longitudinal multi-omics study for practical application in intensive care units

BACKGROUND

The COVID-19 pandemic challenged the management of technical and human resources in intensive care units (ICU) across the world. Several long-term predictors for COVID-19 disease progression have been discovered. However, predictors to support short-term planning of resources and medication that can be translated to future pandemics are still missing. A workflow was established to identify a predictor for short-term COVID-19 disease progression in the acute phase of intensive care patients to support clinical decision-making.

METHODS

Thirty-two patients with SARS-CoV-2 infection were recruited on admission to the ICU and clinical data collected. During their hospitalization, plasma samples were acquired from each patient on multiple occasions, excepting one patient for which only one time point was possible, and the proteome (Inflammation, Immune Response and Organ Damage panels from Olink® Target 96), metabolome and lipidome (flow injection analysis and liquid chromatography-mass spectrometry) analyzed for each sample. Patient visits were grouped according to changes in disease severity based on their respiratory and organ function, and evaluated using a combination of statistical analysis and machine learning. The resulting short-term predictor from this multi-omics approach was compared to the human assessment of disease progression. Furthermore, the potential markers were compared to the baseline levels of 50 healthy subjects with no known SARS-CoV-2 or other viral infections.

RESULTS

A total of 124 clinical parameters, 271 proteins and 782 unique metabolites and lipids were assessed. The dimensionality of the dataset was reduced, selecting 47 from the 1177 parameters available following down-selection, to build the machine learning model. Subsequently, two proteins (C-C motif chemokine 7 (CCL7) and carbonic anhydrase 14 (CA14)) and one lipid (hexosylceramide 18:2; O2/20:0) were linked to disease progression in the studied SARS-CoV-2 infections. Thus, a predictor delivering the prognosis of an upcoming worsening of the patient's condition up to five days in advance with a reasonable accuracy (79 % three days prior to event, 84 % four to five days prior to event) was found. Interestingly, the predictor's performance was complementary to the clinicians' capabilities to foresee a worsening of a patient.

CONCLUSION

This study presents a workflow to identify omics-based biomarkers to support clinical decision-making and resource management in the ICU. This was successfully applied to develop a short-term predictor for aggravation of COVID-19 symptoms. The applied methods can be adapted for future small cohort studies.

Disulfidptosis status influences prognosis and therapeutic response in clear cell renal cell carcinoma

Disulfidptosis is a recently identified type of programmed cell death. It is characterized by aberrant accumulation of intracellular disulfides. The clinical implications of disulfidptosis in clear cell renal cell carcinoma (ccRCC) remain unclear. A series of bioinformatics approaches were employed to analyze ten disulfidptosis-related molecules. Firstly, the expression patterns of the disulfidptosis-related molecules were different between normal and ccRCC tissues. A comprehensive cohort of patients with ccRCC was then assembled from three public databases and subjected to cluster analysis based on disulfidptosis-related molecules. Consensus cluster analysis revealed three distinct disulfidptosis clusters. We then conducted weighted gene co-expression network analysis (WGCNA) to identify highly correlated genes. 267 hub genes were screened out through WGCNA, and three gene clusters were then determined. Finally, we identified 87 genes with prognostic value and then used them to develop a disulfidptosis scoring (DSscore) system, which was proven to independently predict survival in ccRCC. Patients in the high-DSscore group exhibited a significant survival advantage and better immunotherapeutic responses compared with those in the low-DSscore group. However, the patients in the low-DSscore group exhibited a greater degree of chemotherapeutic response. In addition, the expression of disulfidptosis-related molecules was validated by qRT-PCR, and the potential of disulfidptosis-related molecules to indicate distinct cell subtypes were validated by single-cell RNA-sequencing. In conclusion, DSscore is a promising index for predicting the prognosis and efficacy of immunotherapy in patients with ccRCC and may provide a basis for novel strategies for future studies.

Unraveling Epigenetic Interplay between Inflammation, Thrombosis, and Immune-Related Disorders through a Network Meta-analysis

Inflammation and thrombosis are two distinct yet interdependent physiological processes. The inflammation results in the activation of the coagulation system that directs the immune system and its activation, resulting in the initiation of the pathophysiology of thrombosis, a process termed immune-thrombosis. Still, the shared underlying molecular mechanism related to the immune system and coagulation has not yet been explored extensively. Inspired to answer this, we carried out a comprehensive gene expression meta-analysis using publicly available datasets of four diseases, including venous thrombosis, systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease. A total of 609 differentially expressed genes (DEGs) shared by all four datasets were identified based on the combined effect size approach. The pathway enrichment analysis of the DEGs showed enrichment of various epigenetic pathways such as histone-modifying enzymes, posttranslational protein modification, chromatin organization, chromatin-modifying enzymes, HATs acetylate proteins. Network-based protein-protein interaction analysis showed epigenetic enzyme coding genes dominating among the top hub genes. The miRNA-interacting partner of the top 10 hub genes was determined. The predomination of epitranscriptomics regulation opens a layout for the meta-analysis of miRNA datasets of the same four diseases. We identified 30 DEmiRs shared by these diseases. There were 9 common DEmiRs selected from the list of miRNA-interacting partners of top 10 hub genes and shared significant DEmiRs from microRNAs dataset acquisition. These common DEmiRs were found to regulate genes involved in epigenetic modulation and indicate a promising epigenetic aspect that needs to be explored for future molecular studies in the context of immunothrombosis and inflammatory disease.

Investigation of normalization procedures for transcriptome profiles of compounds oriented toward practical study design

The transcriptome profile is a representative phenotype-based descriptor of compounds, widely acknowledged for its ability to effectively capture compound effects. However, the presence of batch differences is inevitable. Despite the existence of sophisticated statistical methods, many of them presume a substantial sample size. How should we design a transcriptome analysis to obtain robust compound profiles, particularly in the context of small datasets frequently encountered in practical scenarios? This study addresses this question by investigating the normalization procedures for transcriptome profiles, focusing on the baseline distribution employed in deriving biological responses as profiles. Firstly, we investigated two large GeneChip datasets, comparing the impact of different normalization procedures. Through an evaluation of the similarity between response profiles of biological replicates within each dataset and the similarity between response profiles of the same compound across datasets, we revealed that the baseline distribution defined by all samples within each batch under batch-corrected condition is a good choice for large datasets. Subsequently, we conducted a simulation to explore the influence of the number of control samples on the robustness of response profiles across datasets. The results offer insights into determining the suitable quantity of control samples for diminutive datasets. It is crucial to acknowledge that these conclusions stem from constrained datasets. Nevertheless, we believe that this study enhances our understanding of how to effectively leverage transcriptome profiles of compounds and promotes the accumulation of essential knowledge for the practical application of such profiles.

Radiation-sensitive genetic prognostic model identifies individuals at risk for radiation resistance in head and neck squamous cell carcinoma

BACKGROUND

The advantages of radiotherapy for head and neck squamous cell carcinoma (HNSCC) depend on the radiation sensitivity of the patient. Here, we established and verified radiological factor-related gene signature and built a prognostic risk model to predict whether radiotherapy would be beneficial.

METHODS

Data from The Cancer Genome Atlas, Gene Expression Omnibus, and RadAtlas databases were subjected to LASSO regression, univariate COX regression, and multivariate COX regression analyses to integrate genomic and clinical information from patients with HNSCC. HNSCC radiation-related prognostic genes were identified, and patients classified into high- and low-risk groups, based on risk scores. Variations in radiation sensitivity according to immunological microenvironment, functional pathways, and immunotherapy response were investigated. Finally, the expression of HNSCC radiation-related genes was verified by qRT-PCR.

RESULTS

We built a clinical risk prediction model comprising a 15-gene signature and used it to divide patients into two groups based on their susceptibility to radiation: radiation-sensitive and radiation-resistant. Overall survival was significantly greater in the radiation-sensitive than the radiation-resistant group. Further, our model was an independent predictor of radiotherapy response, outperforming other clinical parameters, and could be combined with tumor mutational burden, to identify the target population with good predictive value for prognosis at 1, 2, and 3 years. Additionally, the radiation-resistant group was more vulnerable to low levels of immune infiltration, which are significantly associated with DNA damage repair, hypoxia, and cell cycle regulation. Tumor Immune Dysfunction and Exclusion scores also suggested that the resistant group would respond less favorably to immunotherapy.

CONCLUSIONS

Our prognostic model based on a radiation-related gene signature has potential for application as a tool for risk stratification of radiation therapy for patients with HNSCC, helping to identify candidates for radiation therapy and overcome radiation resistance.

Correcting batch effects in large-scale multiomics studies using a reference-material-based ratio method

BACKGROUND

Batch effects are notoriously common technical variations in multiomics data and may result in misleading outcomes if uncorrected or over-corrected. A plethora of batch-effect correction algorithms are proposed to facilitate data integration. However, their respective advantages and limitations are not adequately assessed in terms of omics types, the performance metrics, and the application scenarios.

RESULTS

As part of the Quartet Project for quality control and data integration of multiomics profiling, we comprehensively assess the performance of seven batch effect correction algorithms based on different performance metrics of clinical relevance, i.e., the accuracy of identifying differentially expressed features, the robustness of predictive models, and the ability of accurately clustering cross-batch samples into their own donors. The ratio-based method, i.e., by scaling absolute feature values of study samples relative to those of concurrently profiled reference material(s), is found to be much more effective and broadly applicable than others, especially when batch effects are completely confounded with biological factors of study interests. We further provide practical guidelines for implementing the ratio based approach in increasingly large-scale multiomics studies.

CONCLUSIONS

Multiomics measurements are prone to batch effects, which can be effectively corrected using ratio-based scaling of the multiomics data. Our study lays the foundation for eliminating batch effects at a ratio scale.

Exploring the relationship between lactate metabolism and immunological function in colorectal cancer through genes identification and analysis

Introduction: Metabolic dysregulation is a widely acknowledged contributor for the development and tumorigenesis of colorectal cancer (CRC), highlighting the need for reliable prognostic biomarkers in this malignancy. Methods: Herein, we identified key genes relevant to CRC metabolism through a comprehensive analysis of lactate metabolism-related genes from GSEA MsigDB, employing univariate Cox regression analysis and random forest algorithms. Clinical prognostic analysis was performed following identification of three key genes, and consistent clustering enabled the classification of public datasets into three patterns with significant prognostic differences. The molecular pathways and tumor microenvironment (TME) of these patterns were then investigated through correlation analyses. Quantitative PCR was employed to quantify the mRNA expression levels of the three pivotal genes in CRC tissue. Single-cell RNA sequencing data and fluorescent multiplex immunohistochemistry were utilized to analyze relevant T cells and validate the correlation between key genes and CD4+ T cells. Results: Our analysis revealed that MPC1, COQ2, and ADAMTS13 significantly stratify the cohort into three patterns with distinct prognoses. Additionally, the immune infiltration and molecular pathways were significantly different for each pattern. Among the key genes, MPC1 and COQ2 were positively associated with good prognosis, whereas ADAMTS13 was negatively associated with good prognosis. Single-cell RNA sequencing (scRNA-seq) data illustrated that the relationship between three key genes and T cells, which was further confirmed by the results of fluorescent multiplex immunohistochemistry demonstrating a positive correlation between MPC1 and COQ2 with CD4+ T cells and a negative correlation between ADAMTS13 and CD4+ T cells. Discussion: These findings suggest that the three key lactate metabolism genes, MPC1, COQ2, and ADAMTS13, may serve as effective prognostic biomarkers and support the link between lactate metabolism and the immune microenvironment in CRC.

malbacR: A Package for Standardized Implementation of Batch Correction Methods for Omics Data

Mass spectrometry is a powerful tool for identifying and analyzing biomolecules such as metabolites and lipids in complex biological samples. Liquid chromatography and gas chromatography mass spectrometry studies quite commonly involve large numbers of samples, which can require significant time for sample preparation and analyses. To accommodate such studies, the samples are commonly split into batches. Inevitably, variations in sample handling, temperature fluctuation, imprecise timing, column degradation, and other factors result in systematic errors or biases of the measured abundances between the batches. Numerous methods are available via R packages to assist with batch correction for omics data; however, since these methods were developed by different research teams, the algorithms are available in separate R packages, each with different data input and output formats. We introduce the malbacR package, which consolidates 11 common batch effect correction methods for omics data into one place so users can easily implement and compare the following: pareto scaling, power scaling, range scaling, ComBat, EigenMS, NOMIS, RUV-random, QC-RLSC, WaveICA2.0, TIGER, and SERRF. The malbacR package standardizes data input and output formats across these batch correction methods. The package works in conjunction with the pmartR package, allowing users to seamlessly include the batch effect correction in a pmartR workflow without needing any additional data manipulation.

Making human immune systems more interpretable through systems immunology

The human immune system is a distributed system of specialized cell populations with unique functions that collectively give rise to immune responses to infections and during immune-mediated diseases. Cell composition, plasma proteins, and functional responses vary among individuals, making the system difficult to interpret, but this variation is nonrandom. With careful analyses using novel experimental and computational tools, human immune system composition and function carry interpretable information. Here, we propose that systems-level analyses offer an opportunity to make human immune responses more interpretable in the future, and we discuss herein important considerations and lessons learned to this end. Predictable human immunology holds implications for better diagnostic and curative precision in patients with infectious and immune-associated diseases.

Molecular signature of neutrophil extracellular trap mediating disease module in idiopathic inflammatory myopathy

The rarity and heterogeneity of idiopathic inflammatory myopathy (IIM) pose challenges for researching IIM in affected individuals. We analyzed integrated transcriptomic datasets obtained using muscle tissues from patients with five distinct IIM subtypes to investigate the shared and distinctive cellular and molecular characteristics. A transcriptomic dataset of muscle tissues from normal controls (n = 105) and patients with dermatomyositis (n = 89), polymyositis (n = 33), inclusion body myositis (n = 121), immune-mediated necrotizing myositis (n = 75), and anti-synthetase syndrome (n = 18) was used for differential gene-expression analysis, functional-enrichment analysis, gene set-enrichment analysis, disease-module identification, and kernel-based diffusion scoring. Damage-associated molecular pattern-associated pathways and neutrophil-mediated immunity were significantly enriched across different IIM subtypes, although their activities varied. Interferons-signaling pathways were differentially activated across all five IIM subtypes. In particular, neutrophil extracellular trap (NET) formation was significantly activated and correlated with Fcγ R-mediated signaling pathways. NET formation-associated genes were key for establishing disease modules, and FCGRs, C1QA, and SERPINE1 markedly perturbed the disease modules. Integrated transcriptomic analysis of muscle tissues identified NETs as key components of neutrophil-mediated immunity involved in the pathogenesis of IIM subtypes and, thus, has therapeutically targetable value.

Methodology to identify a gene expression signature by merging microarray datasets

A vast number of microarray datasets have been produced as a way to identify differentially expressed genes and gene expression signatures. A better understanding of these biological processes can help in the diagnosis and prognosis of diseases, as well as in the therapeutic response to drugs. However, most of the available datasets are composed of a reduced number of samples, leading to low statistical, predictive and generalization power. One way to overcome this problem is by merging several microarray datasets into a single dataset, which is typically a challenging task. Statistical methods or supervised machine learning algorithms are usually used to determine gene expression signatures. Nevertheless, statistical methods require an arbitrary threshold to be defined, and supervised machine learning methods can be ineffective when applied to high-dimensional datasets like microarrays. We propose a methodology to identify gene expression signatures by merging microarray datasets. This methodology uses statistical methods to obtain several sets of differentially expressed genes and uses supervised machine learning algorithms to select the gene expression signature. This methodology was validated using two distinct research applications: one using heart failure and the other using autism spectrum disorder microarray datasets. For the first, we obtained a gene expression signature composed of 117 genes, with a classification accuracy of approximately 98%. For the second use case, we obtained a gene expression signature composed of 79 genes, with a classification accuracy of approximately 82%. This methodology was implemented in R language and is available, under the MIT licence, at https://github.com/bioinformatics-ua/MicroGES.

Crosstalk between Gut Microbiota and Epigenetic Markers in Obesity Development: Relationship between Ruminococcus, BMI, and MACROD2/SEL1L2 Methylation

Changes in gut microbiota composition and in epigenetic mechanisms have been proposed to play important roles in energy homeostasis, and the onset and development of obesity. However, the crosstalk between epigenetic markers and the gut microbiome in obesity remains unclear. The main objective of this study was to establish a link between the gut microbiota and DNA methylation patterns in subjects with obesity by identifying differentially methylated DNA regions (DMRs) that could be potentially regulated by the gut microbiota. DNA methylation and bacterial DNA sequencing analysis were performed on 342 subjects with a BMI between 18 and 40 kg/m². DNA methylation analyses identified a total of 2648 DMRs associated with BMI, while ten bacterial genera were associated with BMI. Interestingly, only the abundance of Ruminococcus was associated with one BMI-related DMR, which is located between the MACROD2/SEL1L2 genes. The Ruminococcus abundance negatively correlated with BMI, while the hypermethylated DMR was associated with reduced MACROD2 protein levels in serum. Additionally, the mediation test showed that 19% of the effect of Ruminococcus abundance on BMI is mediated by the methylation of the MACROD2/SEL1L2 DMR. These findings support the hypothesis that a crosstalk between gut microbiota and epigenetic markers may be contributing to obesity development.

Present and future challenges for the investigation of transgenerational epigenetic inheritance

Epigenetic pathways are essential in different biological processes and in phenotype-environment interactions in response to different stressors and they can induce phenotypic plasticity. They encompass several processes that are mitotically and, in some cases, meiotically heritable, so they can be transferred to subsequent generations via the germline. Transgenerational Epigenetic Inheritance (TEI) describes the phenomenon that phenotypic traits, such as changes in fertility, metabolic function, or behavior, induced by environmental factors (e.g., parental care, pathogens, pollutants, climate change), can be transferred to offspring generations via epigenetic mechanisms. Investigations on TEI contribute to deciphering the role of epigenetic mechanisms in adaptation, adversity, and evolution. However, molecular mechanisms underlying the transmission of epigenetic changes between generations, and the downstream chain of events leading to persistent phenotypic changes, remain unclear. Therefore, inter-, (transmission of information between parental and offspring generation via direct exposure) and transgenerational (transmission of information through several generations with disappearance of the triggering factor) consequences of epigenetic modifications remain major issues in the field of modern biology. In this article, we review and describe the major gaps and issues still encountered in the TEI field: the general challenges faced in epigenetic research; deciphering the key epigenetic mechanisms in inheritance processes; identifying the relevant drivers for TEI and implement a collaborative and multi-disciplinary approach to study TEI. Finally, we provide suggestions on how to overcome these challenges and ultimately be able to identify the specific contribution of epigenetics in transgenerational inheritance and use the correct tools for environmental science investigation and biomarkers identification.

Antibody Profiling and In Silico Functional Analysis of Differentially Reactive Antibody Signatures of Glioblastomas and Meningiomas

Studies on tumor-associated antigens in brain tumors are sparse. There is scope for enhancing our understanding of molecular pathology, in order to improve on existing forms, and discover new forms, of treatment, which could be particularly relevant to immuno-oncological strategies. To elucidate immunological differences, and to provide another level of biological information, we performed antibody profiling, based on a high-density protein array (containing 8173 human transcripts), using IgG isolated from the sera of n = 12 preoperative and n = 16 postoperative glioblastomas, n = 26 preoperative and n = 29 postoperative meningiomas, and n = 27 healthy, cancer-free controls. Differentially reactive antigens were compared to gene expression data from an alternate public GBM data set from OncoDB, and were analyzed using the Reactome pathway browser. Protein array analysis identified approximately 350-800 differentially reactive antigens, and revealed different antigen profiles in the glioblastomas and meningiomas, with approximately 20-30%-similar and 10-15%-similar antigens in preoperative and postoperative sera, respectively. Seroreactivity did not correlate with OncoDB-derived gene expression. Antigens in the preoperative glioblastoma sera were enriched for signaling pathways, such as signaling by Rho-GTPases, COPI-mediated anterograde transport and vesicle-mediated transport, while the infectious disease, SRP-dependent membrane targeting cotranslational proteins were enriched in the meningiomas. The pre-vs. postoperative seroreactivity in the glioblastomas was enriched for antigens, e.g., platelet degranulation and metabolism of lipid pathways; in the meningiomas, the antigens were enriched in infectious diseases, metabolism of amino acids and derivatives, and cell cycle. Antibody profiling in both tumor entities elucidated several hundred antigens and characteristic signaling pathways that may provide new insights into molecular pathology and may be of interest for the development of new treatment strategies.

Normalization of Large-Scale Transcriptome Data Using Heuristic Methods

In this study, we introduce an artificial intelligent method for addressing the batch effect of a transcriptome data. The method has several clear advantages in comparison with the alternative methods presently in use. Batch effect refers to the discrepancy in gene expression data series, measured under different conditions. While the data from the same batch (measurements performed under the same conditions) are compatible, combining various batches into 1 data set is problematic because of incompatible measurements. Therefore, it is necessary to perform correction of the combined data (normalization), before performing biological analysis. There are numerous methods attempting to correct data set for batch effect. These methods rely on various assumptions regarding the distribution of the measurements. Forcing the data elements into pre-supposed distribution can severely distort biological signals, thus leading to incorrect results and conclusions. As the discrepancy between the assumptions regarding the data distribution and the actual distribution is wider, the biases introduced by such “correction methods” are greater. We introduce a heuristic method to reduce batch effect. The method does not rely on any assumptions regarding the distribution and the behavior of data elements. Hence, it does not introduce any new biases in the process of correcting the batch effect. It strictly maintains the integrity of measurements within the original batches.

Batch correction methods for nontarget chemical analysis data: application to a municipal wastewater collection system

Nontarget chemical analysis using high-resolution mass spectrometry has increasingly been used to discern spatial patterns and temporal trends in anthropogenic chemical abundance in natural and engineered systems. A critical experimental design consideration in such applications, especially those monitoring complex matrices over long time periods, is a choice between analyzing samples in multiple batches as they are collected, or in one batch after all samples have been processed. While datasets acquired in multiple analytical batches can include the effects of instrumental variability over time, datasets acquired in a single batch risk compound degradation during sample storage. To assess the influence of batch effects on the analysis and interpretation of nontarget data, this study examined a set of 56 samples collected from a municipal wastewater system over 7 months. Each month's samples included 6 from sites within the collection system, one combined influent, and one treated effluent sample. Samples were analyzed using liquid chromatography high-resolution mass spectrometry in positive electrospray ionization mode in multiple batches as the samples were collected and in a single batch at the conclusion of the study. Data were aligned and normalized using internal standard scaling and ComBat, an empirical Bayes method developed for estimating and removing batch effects in microarrays. As judged by multiple lines of evidence, including comparing principal variance component analysis between single and multi-batch datasets and through patterns in principal components and hierarchical clustering analyses, ComBat appeared to significantly reduce the influence of batch effects. For this reason, we recommend the use of more, small batches with an appropriate batch correction step rather than acquisition in one large batch.

Whole-blood methylation signatures are associated with and accurately classify multiple sclerosis disease severity

BACKGROUND

The variation in multiple sclerosis (MS) disease severity is incompletely explained by genetics, suggesting genetic and environmental interactions are involved. Moreover, the lack of prognostic biomarkers makes it difficult for clinicians to optimise care. DNA methylation is one epigenetic mechanism by which gene-environment interactions can be assessed. Here, we aimed to identify DNA methylation patterns associated with mild and severe relapse-onset MS (RMS) and to test the utility of methylation as a predictive biomarker.

METHODS

We conducted an epigenome-wide association study between 235 females with mild (n = 119) or severe (n = 116) with RMS. Methylation was measured with the Illumina methylationEPIC array and analysed using logistic regression. To generate hypotheses about the functional consequence of differential methylation, we conducted gene set enrichment analysis using ToppGene. We compared the accuracy of three machine learning models in classifying disease severity: (1) clinical data available at baseline (age at onset and first symptoms) built using elastic net (EN) regression, (2) methylation data using EN regression and (3) a weighted methylation risk score of differentially methylated positions (DMPs) from the main analysis using logistic regression. We used a conservative 70:30 test:train split for classification modelling. A false discovery rate threshold of 0.05 was used to assess statistical significance.

RESULTS

Females with mild or severe RMS had 1472 DMPs in whole blood (839 hypermethylated, 633 hypomethylated in the severe group). Differential methylation was enriched in genes related to neuronal cellular compartments and processes, and B-cell receptor signalling. Whole-blood methylation levels at 1708 correlated CpG sites classified disease severity more accurately (machine learning model 2, AUC = 0.91) than clinical data (model 1, AUC = 0.74) or the wMRS (model 3, AUC = 0.77). Of the 1708 selected CpGs, 100 overlapped with DMPs from the main analysis at the gene level. These overlapping genes were enriched in neuron projection and dendrite extension, lending support to our finding that neuronal processes, rather than immune processes, are implicated in disease severity.

CONCLUSION

RMS disease severity is associated with whole-blood methylation at genes related to neuronal structure and function. Moreover, correlated whole-blood methylation patterns can assign disease severity in females with RMS more accurately than clinical data available at diagnosis.

Diagnostic signature, subtype classification, and immune infiltration of key m6A regulators in osteomyelitis patients

Background: As a recurrent inflammatory bone disease, the treatment of osteomyelitis is always a tricky problem in orthopaedics. N6-methyladenosine (m6A) regulators play significant roles in immune and inflammatory responses. Nevertheless, the function of m6A modification in osteomyelitis remains unclear. Methods: Based on the key m6A regulators selected by the GSE16129 dataset, a nomogram model was established to predict the incidence of osteomyelitis by using the random forest (RF) method. Through unsupervised clustering, osteomyelitis patients were divided into two m6A subtypes, and the immune infiltration of these subtypes was further evaluated. Validating the accuracy of the diagnostic model for osteomyelitis and the consistency of clustering based on the GSE30119 dataset. Results: 3 writers of Methyltransferase-like 3 (METTL3), RNA-binding motif protein 15B (RBM15B) and Casitas B-lineage proto-oncogene like 1 (CBLL1) and three readers of YT521-B homology domain-containing protein 1 (YTHDC1), YT521-B homology domain-containing family 3 (YTHDF2) and Leucine-rich PPR motif-containing protein (LRPPRC) were identified by difference analysis, and their Mean Decrease Gini (MDG) scores were all greater than 10. Based on these 6 significant m6A regulators, a nomogram model was developed to predict the incidence of osteomyelitis, and the fitting curve indicated a high degree of fit in both the test and validation groups. Two m6A subtypes (cluster A and cluster B) were identified by the unsupervised clustering method, and there were significant differences in m6A scores and the abundance of immune infiltration between the two m6A subtypes. Among them, two m6A regulators (METTL3 and LRPPRC) were closely related to immune infiltration in patients with osteomyelitis. Conclusion: m6A regulators play key roles in the molecular subtypes and immune response of osteomyelitis, which may provide assistance for personalized immunotherapy in patients with osteomyelitis.

Integrating bioinformatic strategies in spatial life science research

As space exploration programs progress, manned space missions will become more frequent and farther away from Earth, putting a greater emphasis on astronaut health. Through the collaborative efforts of researchers from various countries, the effect of the space environment factors on living systems is gradually being uncovered. Although a large number of interconnected research findings have been produced, their connection seems to be confused, and many unknown effects are left to be discovered. Simultaneously, several valuable data resources have emerged, accumulating data measuring biological effects in space that can be used to further investigate the unknown biological adaptations. In this review, the previous findings and their correlations are sorted out to facilitate the understanding of biological adaptations to space and the design of countermeasures. The biological effect measurement methods/data types are also organized to provide references for experimental design and data analysis. To aid deeper exploration of the data resources, we summarized common characteristics of the data generated from longitudinal experiments, outlined challenges or caveats in data analysis and provided corresponding solutions by recommending bioinformatics strategies and available models/tools.

Machine learning on syngeneic mouse tumor profiles to model clinical immunotherapy response

Most patients with cancer are refractory to immune checkpoint blockade (ICB) therapy, and proper patient stratification remains an open question. Primary patient data suffer from high heterogeneity, low accessibility, and lack of proper controls. In contrast, syngeneic mouse tumor models enable controlled experiments with ICB treatments. Using transcriptomic and experimental variables from >700 ICB-treated/control syngeneic mouse tumors, we developed a machine learning framework to model tumor immunity and identify factors influencing ICB response. Projected on human immunotherapy trial data, we found that the model can predict clinical ICB response. We further applied the model to predicting ICB-responsive/resistant cancer types in The Cancer Genome Atlas, which agreed well with existing clinical reports. Last, feature analysis implicated factors associated with ICB response. In summary, our computational framework based on mouse tumor data reliably stratified patients regarding ICB response, informed resistance mechanisms, and has the potential for wide applications in disease treatment studies.

reComBat: batch-effect removal in large-scale multi-source gene-expression data integration

Motivation

With the steadily increasing abundance of omics data produced all over the world under vastly different experimental conditions residing in public databases, a crucial step in many data-driven bioinformatics applications is that of data integration. The challenge of batch-effect removal for entire databases lies in the large number of batches and biological variation, which can result in design matrix singularity. This problem can currently not be solved satisfactorily by any common batch-correction algorithm.

Results

We present reComBat, a regularized version of the empirical Bayes method to overcome this limitation and benchmark it against popular approaches for the harmonization of public gene-expression data (both microarray and bulkRNAsq) of the human opportunistic pathogen Pseudomonas aeruginosa. Batch-effects are successfully mitigated while biologically meaningful gene-expression variation is retained. reComBat fills the gap in batch-correction approaches applicable to large-scale, public omics databases and opens up new avenues for data-driven analysis of complex biological processes beyond the scope of a single study.

Availability and implementation

The code is available at https://github.com/BorgwardtLab/reComBat, all data and evaluation code can be found at https://github.com/BorgwardtLab/batchCorrectionPublicData.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Dynamics of Methylation of CpG Sites Associated With Systemic Lupus Erythematosus Subtypes in a Longitudinal Cohort

OBJECTIVE

Findings from cross-sectional studies have revealed associations between DNA methylation and systemic lupus erythematosus (SLE) outcomes. This study was undertaken to investigate the dynamics of DNA methylation by examining participants from an SLE longitudinal cohort using samples collected at 2 time points.

METHODS

A total of 101 participants from the California Lupus Epidemiology Study were included in our analysis. DNA was extracted from blood samples collected at the time of enrolment in the cohort and samples collected after 2 years and was analyzed using Illumina EPIC BeadChip kit. Paired t-tests were used to identify genome-wide changes which included 256 CpG sites previously found to be associated with SLE subtypes. Linear mixed models were developed to understand the relationship between DNA methylation and disease activity, medication use, and sample cell-type proportions, adjusted for age, sex, and genetic principal components.

RESULTS

The majority of CpGs that were previously determined to be associated with SLE subtypes remained stable over 2 years (185 CpGs [72.3%]; t-test false discovery rate >0.05). Compared to background genome-wide methylation, there was an enrichment of SLE subtype-associated CpGs that changed over time (27.7% versus 0.34%). Changes in cell-type proportions were associated with changes at 67 CpGs (P < 2.70 × 10^-5 ), and 15 CpGs had at least 1 significant association with immunosuppressant use.

CONCLUSION

In this longitudinal SLE cohort, we identified a subset of SLE subtype-associated CpGs that remained stable over time and may be useful as biomarkers of disease subtypes. Another subset of SLE subtype-associated CpGs changed at a higher proportion compared to the genome-wide methylome. Additional studies are needed to understand the etiology and impact of these changes on methylation of SLE-associated CpGs.

Microarray meta-analysis reveals IL6 and p38β/MAPK11 as potential targets of hsa-miR-124 in endothelial progenitor cells: Implications for stent re-endothelization in diabetic patients

Circulating endothelial progenitor cells (EPCs) play an important role in the repair processes of damaged vessels, favoring re-endothelization of stented vessels to minimize restenosis. EPCs number and function is diminished in patients with type 2 diabetes, a known risk factor for restenosis. Considering the impact of EPCs in vascular injury repair, we conducted a meta-analysis of microarray to assess the transcriptomic profile and determine target genes during the differentiation process of EPCs into mature ECs. Five microarray datasets, including 13 EPC and 12 EC samples were analyzed, using the online tool ExpressAnalyst. Differentially expressed genes (DEGs) analysis was done by Limma method, with an | log2FC| &gt; 1 and FDR &lt; 0.05. Combined p-value by Fisher exact method was computed for the intersection of datasets. There were 3,267 DEGs, 1,539 up-regulated and 1,728 down-regulated in EPCs, with 407 common DEGs in at least four datasets. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed enrichment for terms related to “AGE-RAGE signaling pathway in diabetic complications.” Intersection of common DEGs, KEGG pathways genes and genes in protein-protein interaction network (PPI) identified four key genes, two up-regulated (IL1B and STAT5A) and two down-regulated (IL6 and MAPK11). MicroRNA enrichment analysis of common DEGs depicted five hub microRNA targeting 175 DEGs, including STAT5A, IL6 and MAPK11, with hsa-miR-124 as common regulator. This group of genes and microRNAs could serve as biomarkers of EPCs differentiation during coronary stenting as well as potential therapeutic targets to improve stent re-endothelization, especially in diabetic patients.

Using Machine Learning to Identify Biomarkers Affecting Fat Deposition in Pigs by Integrating Multisource Transcriptome Information

Fat deposition in pigs is not only closely related to pig production efficiency and pork quality but also an ideal model for human obesity. Transcriptome sequencing is widely used to study fat deposition. However, due to small sample sizes, high false positive rates, and poor consistency of results from different studies, new strategies are urgently needed. Machine learning, a new analysis method, can effectively fit complex data and accurately identify samples and genes. In this study, 36 samples of adipose tissue, muscle tissue, and liver tissue were collected from Songliao black pigs and Landrace pigs, and the mRNA of all the samples was sequenced. In addition, we collected transcriptome data for 64 samples in the GEO database from four different sources. After standardization and imputation of missing values in the data set comprising 100 samples, traditional differential expression analysis was carried out, and different numbers of expressed genes were selected as features for the training model of eight machine learning methods. In the 1000 replications of fourfold cross validation with 100 samples, AdaBoost performed best, with an average prediction accuracy greater than 93% and the highest mean area under the curve in predicting the high- and low-fat content groups among the eight ML methods. According to their performance-based ranks inferred by AdaBoost, 12 genes related to fat deposition were identified; among them, FASN and APOD were specifically expressed in adipose tissue, and APOA1 was specifically expressed in the liver, which could be important candidate biomarkers affecting fat deposition.

Evaluation of the Stellae-123 prognostic gene expression signature in acute myeloid leukemia

Risk stratification in acute myeloid leukemia (AML) has been extensively improved thanks to the incorporation of recurrent cytogenomic alterations into risk stratification guidelines. However, mortality rates among fit patients assigned to low or intermediate risk groups are still high. Therefore, significant room exists for the improvement of AML prognostication. In a previous work, we presented the Stellae-123 gene expression signature, which achieved a high accuracy in the prognostication of adult patients with AML. Stellae-123 was particularly accurate to restratify patients bearing high-risk mutations, such as ASXL1, RUNX1 and TP53. The intention of the present work was to evaluate the prognostic performance of Stellae-123 in external cohorts using RNAseq technology. For this, we evaluated the signature in 3 different AML cohorts (2 adult and 1 pediatric). Our results indicate that the prognostic performance of the Stellae-123 signature is reproducible in the 3 cohorts of patients. Additionally, we evidenced that the signature was superior to the European LeukemiaNet 2017 and the pediatric clinical risk scores in the prediction of survival at most of the evaluated time points. Furthermore, integration with age substantially enhanced the accuracy of the model. In conclusion, Stellae-123 is a reproducible machine learning algorithm based on a gene expression signature with promising utility in the field of AML.

Spatial intra-tumor heterogeneity is associated with survival of lung adenocarcinoma patients

Intra-tumor heterogeneity (ITH) of human tumors is important for tumor progression, treatment response, and drug resistance. However, the spatial distribution of ITH remains incompletely understood. Here, we present spatial analysis of ITH in lung adenocarcinomas from 147 patients using multi-region mass spectrometry of >5,000 regions, single-cell copy number sequencing of ~2,000 single cells, and cyclic immunofluorescence of >10 million cells. We identified two distinct spatial patterns among tumors, termed clustered and random geographic diversification (GD). These patterns were observed in the same samples using both proteomic and genomic data. The random proteomic GD pattern, which is characterized by decreased cell adhesion and lower levels of tumor-interacting endothelial cells, was significantly associated with increased risk of recurrence or death in two independent patient cohorts. Our study presents comprehensive spatial mapping of ITH in lung adenocarcinoma and provides insights into the mechanisms and clinical consequences of GD.

Longitudinal phenotypic aging metrics in the Baltimore Longitudinal Study of Aging

To define metrics of phenotypic aging, it is essential to identify biological and environmental factors that influence the pace of aging. Previous attempts to develop aging metrics were hampered by cross-sectional designs and/or focused on younger populations. In the Baltimore Longitudinal Study of Aging (BLSA), we collected longitudinally across the adult age range a comprehensive list of phenotypes within four domains (body composition, energetics, homeostatic mechanisms and neurodegeneration/neuroplasticity) and functional outcomes. We integrated individual deviations from population trajectories into a global longitudinal phenotypic metric of aging and demonstrate that accelerated longitudinal phenotypic aging is associated with faster physical and cognitive decline, faster accumulation of multimorbidity and shorter survival. These associations are more robust compared with the use of phenotypic and epigenetic measurements at a single time point. Estimation of these metrics required repeated measures of multiple phenotypes over time but may uniquely facilitate the identification of mechanisms driving phenotypic aging and subsequent age-related functional decline.

HarmonizR enables data harmonization across independent proteomic datasets with appropriate handling of missing values

Dataset integration is common practice to overcome limitations in statistically underpowered omics datasets. Proteome datasets display high technical variability and frequent missing values. Sophisticated strategies for batch effect reduction are lacking or rely on error-prone data imputation. Here we introduce HarmonizR, a data harmonization tool with appropriate missing value handling. The method exploits the structure of available data and matrix dissection for minimal data loss, without data imputation. This strategy implements two common batch effect reduction methods—ComBat and limma (removeBatchEffect()). The HarmonizR strategy, evaluated on four exemplarily analyzed datasets with up to 23 batches, demonstrated successful data harmonization for different tissue preservation techniques, LC-MS/MS instrumentation setups, and quantification approaches. Compared to data imputation methods, HarmonizR was more efficient and performed superior regarding the detection of significant proteins. HarmonizR is an efficient tool for missing data tolerant experimental variance reduction and is easily adjustable for individual dataset properties and user preferences.

Dataset integration is common practice to overcome limitations in statistically underpowered omics datasets. Here the authors present “HarmonizR”, a tool for missing data tolerant experimental variance reduction in large, integrated but independently generated datasets without data imputation, adjustable for individual dataset modalities, correction algorithm, and user preferences.

Data harmonisation for information fusion in digital healthcare: A state-of-the-art systematic review, meta-analysis and future research directions

Removing the bias and variance of multicentre data has always been a challenge in large scale digital healthcare studies, which requires the ability to integrate clinical features extracted from data acquired by different scanners and protocols to improve stability and robustness. Previous studies have described various computational approaches to fuse single modality multicentre datasets. However, these surveys rarely focused on evaluation metrics and lacked a checklist for computational data harmonisation studies. In this systematic review, we summarise the computational data harmonisation approaches for multi-modality data in the digital healthcare field, including harmonisation strategies and evaluation metrics based on different theories. In addition, a comprehensive checklist that summarises common practices for data harmonisation studies is proposed to guide researchers to report their research findings more effectively. Last but not least, flowcharts presenting possible ways for methodology and metric selection are proposed and the limitations of different methods have been surveyed for future research.

Decision Theory versus Conventional Statistics for Personalized Therapy of Breast Cancer

Estrogen and progesterone receptors being present or not represents one of the most important biomarkers for therapy selection in breast cancer patients. Conventional measurement by immunohistochemistry (IHC) involves errors, and numerous attempts have been made to increase precision by additional information from gene expression. This raises the question of how to fuse information, in particular, if there is disagreement. It is the primary domain of Dempster–Shafer decision theory (DST) to deal with contradicting evidence on the same item (here: receptor status), obtained through different techniques. DST is widely used in technical settings, such as self-driving cars and aviation, and is also promising to deliver significant advantages in medicine. Using data from breast cancer patients already presented in previous work, we focus on comparing DST with classical statistics in this work, to pave the way for its application in medicine. First, we explain how DST not only considers probabilities (a single number per sample), but also incorporates uncertainty in a concept of ‘evidence’ (two numbers per sample). This allows for very powerful displays of patient data in so-called ternary plots, a novel and crucial advantage for medical interpretation. Results are obtained according to conventional statistics (ODDS) and, in parallel, according to DST. Agreement and differences are evaluated, and the particular merits of DST discussed. The presented application demonstrates how decision theory introduces new levels of confidence in diagnoses derived from medical data.

Large-Scale Meta-Longitudinal Microbiome Data with a Known Batch Factor

Data contamination in meta-approaches where multiple biological samples are combined considerably affects the results of subsequent downstream analyses, such as differential abundance tests comparing multiple groups at a fixed time point. Little has been thoroughly investigated regarding the impact of the lurking variable of various batch sources, such as different days or different laboratories, in more complicated time series experimental designs, for instance, repeatedly measured longitudinal data and metadata. We highlight that the influence of batch factors is significant on subsequent downstream analyses, including longitudinal differential abundance tests, by performing a case study of microbiome time course data with two treatment groups and a simulation study of mimic microbiome longitudinal counts.

Modern Beijing sublineage of Mycobacterium tuberculosis shift macrophage into a hyperinflammatory status

The high prevalence of the modern Beijing sublineage of Mycobacterium tuberculosis may be related to increased virulence, although the responsible mechanisms remain poorly understood. We previously described enhanced triacylglycerol accumulation in modern Beijing strains. Here we show that modern Beijing strains grow faster in vitro and trigger a vigorous immune response and pronounced macrophage infiltration. Transcriptomic analysis of bone marrow derived macrophages infected with modern Beijing lineage strains revealed a significant enrichment of infection, cholesterol homeostasis and amino acid metabolic pathways. The upregulation of proinflammatory / bactericidal cytokines was confirmed by RT–PCR analysis, which is also in consistent with the reduced bacterial burden in modern strains infected macrophages. These results suggest that modern Beijing strains elicit a hyperinflammatory response which might indicate a stronger virulence and contribute to their extensive global prevalence.

Personally Tailored Survival Prediction of Patients With Follicular Lymphoma Using Machine Learning Transcriptome-Based Models

Follicular Lymphoma (FL) has a 10-year mortality rate of 20%, and this is mostly related to lymphoma progression and transformation to higher grades. In the era of personalized medicine it has become increasingly important to provide patients with an optimal prediction about their expected outcomes. The objective of this work was to apply machine learning (ML) tools on gene expression data in order to create individualized predictions about survival in patients with FL. Using data from two different studies, we were able to create a model which achieved good prediction accuracies in both cohorts (c-indexes of 0.793 and 0.662 in the training and test sets). Integration of this model with m7-FLIPI and age rendered high prediction accuracies in the test set (cox c-index 0.79), and a simplified approach identified 4 groups with remarkably different outcomes in terms of survival. Importantly, one of the groups comprised 27.35% of patients and had a median survival of 4.64 years. In summary, we have created a gene expression-based individualized predictor of overall survival in FL that can improve the predictions of the m7-FLIPI score.

pwrBRIDGE: a user-friendly web application for power and sample size estimation in batch-confounded microarray studies with dependent samples

Batch effect Reduction of mIcroarray data with Dependent samples usinG Empirical Bayes (BRIDGE) is a recently developed statistical method to address the issue of batch effect correction in batch-confounded microarray studies with dependent samples. The key component of the BRIDGE methodology is the use of samples run as technical replicates in two or more batches, "bridging samples", to inform batch effect correction/attenuation. While previously published results indicate a relationship between the number of bridging samples, M, and the statistical power of downstream statistical testing on the batch-corrected data, there is of yet no formal statistical framework or user-friendly software, for estimating M to achieve a specific statistical power for hypothesis tests conducted on the batch-corrected data. To fill this gap, we developed pwrBRIDGE, a simulation-based approach to estimate the bridging sample size, M, in batch-confounded longitudinal microarray studies. To illustrate the use of pwrBRIDGE, we consider a hypothetical, longitudinal batch-confounded study whose goal is to identify Alzheimer's disease (AD) progression-associated genes from amnestic mild cognitive impairment (aMCI) to AD in human blood after a 5-year follow-up. pwrBRIDGE helps researchers design and plan batch-confounded microarray studies with dependent samples to avoid over- or under-powered studies.

Batch effect reduction of microarray data with dependent samples using an empirical Bayes approach (BRIDGE)

Batch-effects present challenges in the analysis of high-throughput molecular data and are particularly problematic in longitudinal studies when interest lies in identifying genes/features whose expression changes over time, but time is confounded with batch. While many methods to correct for batch-effects exist, most assume independence across samples; an assumption that is unlikely to hold in longitudinal microarray studies. We propose Batch effect Reduction of mIcroarray data with Dependent samples usinGEmpirical Bayes (BRIDGE), a three-step parametric empirical Bayes approach that leverages technical replicate samples profiled at multiple timepoints/batches, so-called "bridge samples", to inform batch-effect reduction/attenuation in longitudinal microarray studies. Extensive simulation studies and an analysis of a real biological data set were conducted to benchmark the performance of BRIDGE against both ComBat and longitudinalComBat. Our results demonstrate that while all methods perform well in facilitating accurate estimates of time effects, BRIDGE outperforms both ComBat and longitudinal ComBat in the removal of batch-effects in data sets with bridging samples, and perhaps as a result, was observed to have improved statistical power for detecting genes with a time effect. BRIDGE demonstrated competitive performance in batch effect reduction of confounded longitudinal microarray studies, both in simulated and a real data sets, and may serve as a useful preprocessing method for researchers conducting longitudinal microarray studies that include bridging samples.

Integrative analysis of lung molecular signatures reveals key drivers of idiopathic pulmonary fibrosis

Background

Idiopathic pulmonary fibrosis (IPF) is a devastating disease with a high clinical burden. The molecular signatures of IPF were analyzed to distinguish molecular subgroups and identify key driver genes and therapeutic targets.

Methods

Thirteen datasets of lung tissue transcriptomics including 585 IPF patients and 362 normal controls were obtained from the databases and subjected to filtration of differentially expressed genes (DEGs). A functional enrichment analysis, agglomerative hierarchical clustering, network-based key driver analysis, and diffusion scoring were performed, and the association of enriched pathways and clinical parameters was evaluated.

Results

A total of 2,967 upregulated DEGs was filtered during the comparison of gene expression profiles of lung tissues between IPF patients and healthy controls. The core molecular network of IPF featured p53 signaling pathway and cellular senescence. IPF patients were classified into two molecular subgroups (C1, C2) via unsupervised clustering. C1 was more enriched in the p53 signaling pathway and ciliated cells and presented a worse prognostic score, while C2 was more enriched for cellular senescence, profibrosing pathways, and alveolar epithelial cells. The p53 signaling pathway was closely correlated with a decline in forced vital capacity and carbon monoxide diffusion capacity and with the activation of cellular senescence. CDK1/2, CKDNA1A, CSNK1A1, HDAC1/2, FN1, VCAM1, and ITGA4 were the key regulators as evidence by high diffusion scores in the disease module. Currently available and investigational drugs showed differential diffusion scores in terms of their target molecules.

Conclusions

An integrative molecular analysis of IPF lungs identified two molecular subgroups with distinct pathobiological characteristics and clinical prognostic scores. Inhibition against CDKs or HDACs showed great promise for controlling lung fibrosis. This approach provided molecular insights to support the prediction of clinical outcomes and the selection of therapeutic targets in IPF patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12890-021-01749-3.

Epigenome-wide association studies: current knowledge, strategies and recommendations

The aetiology and pathophysiology of complex diseases are driven by the interaction between genetic and environmental factors. The variability in risk and outcomes in these diseases are incompletely explained by genetics or environmental risk factors individually. Therefore, researchers are now exploring the epigenome, a biological interface at which genetics and the environment can interact. There is a growing body of evidence supporting the role of epigenetic mechanisms in complex disease pathophysiology. Epigenome-wide association studies (EWASes) investigate the association between a phenotype and epigenetic variants, most commonly DNA methylation. The decreasing cost of measuring epigenome-wide methylation and the increasing accessibility of bioinformatic pipelines have contributed to the rise in EWASes published in recent years. Here, we review the current literature on these EWASes and provide further recommendations and strategies for successfully conducting them. We have constrained our review to studies using methylation data as this is the most studied epigenetic mechanism; microarray-based data as whole-genome bisulphite sequencing remains prohibitively expensive for most laboratories; and blood-based studies due to the non-invasiveness of peripheral blood collection and availability of archived DNA, as well as the accessibility of publicly available blood-cell-based methylation data. Further, we address multiple novel areas of EWAS analysis that have not been covered in previous reviews: (1) longitudinal study designs, (2) the chip analysis methylation pipeline (ChAMP), (3) differentially methylated region (DMR) identification paradigms, (4) methylation quantitative trait loci (methQTL) analysis, (5) methylation age analysis and (6) identifying cell-specific differential methylation from mixed cell data using statistical deconvolution.

Childhood adversity correlates with stable changes in DNA methylation trajectories in children and converges with epigenetic signatures of prenatal stress

Childhood maltreatment (CM) is an established major risk factor for a number of negative health outcomes later in life. While epigenetic mechanisms, such as DNA methylation (DNAm), have been proposed as a means of embedding this environmental risk factor, little is known about its timing and trajectory, especially in very young children. It is also not clear whether additional environmental adversities, often experienced by these children, converge on similar DNAm changes. Here, we calculated a cumulative adversity score, which additionally to CM includes socioeconomic status (SES), other life events, parental psychopathology and epigenetic biomarkers of prenatal smoking and alcohol consumption. We investigated the effects of CM alone as well as the adversity score on longitudinal DNAm trajectories in the Berlin Longitudinal Child Study. This is a cohort of 173 children aged 3-5 years at baseline of whom 86 were exposed to CM. These children were followed-up for 2 years with extensive psychometric and biological assessments as well as saliva collection at 5 time points providing genome-wide DNAm levels. Overall, only a few DNAm patterns were stable over this timeframe, but less than 10 DNAm regions showed significant changes. At baseline, neither CM nor the adversity score associated with DNAm changes. However, in 6 differentially methylated regions (DMRs), CM and the adversity score significantly moderated DNAm trajectories over time. A number of these DMRs have previously been associated with adverse prenatal exposures. In our study, children exposed to CM also presented with epigenetic signatures indicative of increased prenatal exposure to tobacco and alcohol, as compared to non-CM exposed children. These epigenetic signatures of prenatal exposure strongly correlate with DNAm regions associated with CM and the adversity score. Finally, weighted correlation network analysis revealed a module of CpGs exclusively associated with CM. While our study identifies DNAm loci specifically associated with CM, especially within long non-coding RNAs, the majority of associations were found with the adversity score with convergent association with indicators of adverse prenatal exposures. This study highlights the importance of mapping not only of the epigenome but also the exposome and extending the observational timeframe to well before birth.

The pediatric buccal epigenetic clock identifies significant ageing acceleration in children with internalizing disorder and maltreatment exposure

Background

Studies reporting accelerated ageing in children with affective disorders or maltreatment exposure have relied on algorithms for estimating epigenetic age derived from adult samples. These algorithms have limited validity for epigenetic age estimation during early development. We here use a pediatric buccal epigenetic (PedBE) clock to predict DNA methylation-based ageing deviation in children with and without internalizing disorder and assess the moderating effect of maltreatment exposure. We further conduct a gene set enrichment analysis to assess the contribution of glucocorticoid signaling to PedBE clock-based results.

Method

DNA was isolated from saliva of 158 children [73 girls, 85 boys; mean age (SD) = 4.25 (0.8) years] including children with internalizing disorder and maltreatment exposure. Epigenetic age was estimated based on DNA methylation across 94 CpGs of the PedBE clock. Residuals of epigenetic age regressed against chronological age were contrasted between children with and without internalizing disorder. Maltreatment was coded in 3 severity levels and entered in a moderation model. Genome-wide dexamethasone-responsive CpGs were derived from an independent sample and enrichment of these CpGs within the PedBE clock was identified.

Results

Children with internalizing disorder exhibited significant acceleration of epigenetic ageing as compared to children without internalizing disorder (F_1,147 = 6.67, p = .011). This association was significantly moderated by maltreatment severity (b = 0.49, 95% CI [0.073, 0.909], t = 2.322, p = .022). Children with internalizing disorder who had experienced maltreatment exhibited ageing acceleration relative to children with no internalizing disorder (1–2 categories: b = 0.50, 95% CI [0.170, 0.821], t = 3.008, p = .003; 3 or more categories: b = 0.99, 95% CI [0.380, 1.593], t = 3.215, p = .002). Children with internalizing disorder who were not exposed to maltreatment did not show epigenetic ageing acceleration. There was significant enrichment of dexamethasone-responsive CpGs within the PedBE clock (OR = 4.36, p = 1.65*10–6). Among the 94 CpGs of the PedBE clock, 18 (19%) were responsive to dexamethasone.

Conclusion

Using the novel PedBE clock, we show that internalizing disorder is associated with accelerated epigenetic ageing in early childhood. This association is moderated by maltreatment severity and may, in part, be driven by glucocorticoids. Identifying developmental drivers of accelerated epigenetic ageing after maltreatment will be critical to devise early targeted interventions.

DNA methylation-based biomarkers of aging were slowed down in a two-year diet and physical activity intervention trial: the DAMA study

Several biomarkers of healthy aging have been proposed in recent years, including the epigenetic clocks, based on DNA methylation (DNAm) measures, which are getting increasingly accurate in predicting the individual biological age. The recently developed "next-generation clock" DNAmGrimAge outperforms "first-generation clocks" in predicting longevity and the onset of many age-related pathological conditions and diseases. Additionally, the total number of stochastic epigenetic mutations (SEMs), also known as the epigenetic mutation load (EML), has been proposed as a complementary DNAm-based biomarker of healthy aging. A fundamental biological property of epigenetic, and in particular DNAm modifications, is the potential reversibility of the effect, raising questions about the possible slowdown of epigenetic aging by modifying one's lifestyle. Here, we investigated whether improved dietary habits and increased physical activity have favorable effects on aging biomarkers in healthy postmenopausal women. The study sample consists of 219 women from the "Diet, Physical Activity, and Mammography" (DAMA) study: a 24-month randomized factorial intervention trial with DNAm measured twice, at baseline and the end of the trial. Women who participated in the dietary intervention had a significant slowing of the DNAmGrimAge clock, whereas increasing physical activity led to a significant reduction of SEMs in crucial cancer-related pathways. Our study provides strong evidence of a causal association between lifestyle modification and slowing down of DNAm aging biomarkers. This randomized trial elucidates the causal relationship between lifestyle and healthy aging-related epigenetic mechanisms.

TISMO: syngeneic mouse tumor database to model tumor immunity and immunotherapy response

Syngeneic mouse models are tumors derived from murine cancer cells engrafted on genetically identical mouse strains. They are widely used tools for studying tumor immunity and immunotherapy response in the context of a fully functional murine immune system. Large volumes of syngeneic mouse tumor expression profiles under different immunotherapy treatments have been generated, although a lack of systematic collection and analysis makes data reuse challenging. We present Tumor Immune Syngeneic MOuse (TISMO), a database with an extensive collection of syngeneic mouse model profiles with interactive visualization features. TISMO contains 605 in vitro RNA-seq samples from 49 syngeneic cancer cell lines across 23 cancer types, of which 195 underwent cytokine treatment. TISMO also includes 1518 in vivo RNA-seq samples from 68 syngeneic mouse tumor models across 19 cancer types, of which 832 were from immune checkpoint blockade (ICB) studies. We manually annotated the sample metadata, such as cell line, mouse strain, transplantation site, treatment, and response status, and uniformly processed and quality-controlled the RNA-seq data. Besides data download, TISMO provides interactive web interfaces to investigate whether specific gene expression, pathway enrichment, or immune infiltration level is associated with differential immunotherapy response. TISMO is available at http://tismo.cistrome.org.

Insights Into Transcriptome Profiles Associated With Wooden Breast Myopathy in Broilers Slaughtered at the Age of 6 or 7 Weeks

Transcriptomes associated with wooden breast (WB) were characterized in broilers at two different market ages. Breasts (Pectoralis major) were collected, 20-min postmortem, from male Ross 308 broilers slaughtered at 6 and 7 weeks of age. The breasts were classified as "non-WB" or "WB" based on palpation hardness scoring (non-WB = no abnormal hardness, WB = consistently hardened). Total RNA was isolated from 16 samples (n = 3 for 6 week non-WB, n = 3 for 6 week WB; n = 5 for 7 week non-WB, n = 5 for 7 week WB). Transcriptome was profiled using a chicken gene expression microarray with one-color hybridization technique, and compared between non-WB and WB samples of the same age. Among 6 week broilers, 910 transcripts were differentially expressed (DE) (false discovery rate, FDR < 0.05). Pathway analysis underlined metabolisms of glucose and lipids along with gap junctions, tight junction, and focal adhesion (FA) signaling as the top enriched pathways. For the 7 week broilers, 1,195 transcripts were identified (FDR < 0.05) with regulation of actin cytoskeleton, mitogen-activated protein kinase (MAPK) signaling, protein processing in endoplasmic reticulum and FA signaling highlighted as the enriched affected pathways. Absolute transcript levels of eight genes (actinin-1 - ACTN1, integrin-linked kinase - ILK, integrin subunit alpha 8 - ITGA8, integrin subunit beta 5 - ITGB5, protein tyrosine kinase 2 - PTK2, paxillin - PXN, talin 1 - TLN1, and vinculin - VCL) of FA signaling pathway were further elucidated using a droplet digital polymerase chain reaction. The results indicated that, in 6 week broilers, ITGA8 abundance in WB was greater than that of non-WB samples (p < 0.05). Concerning 7 week broilers, greater absolute levels of ACTN1, ILK, ITGA8, and TLN1, accompanied with a reduced ITGB5 were found in WB compared with non-WB (p < 0.05). Transcriptional modification of FA signaling underlined the potential of disrupted cell-cell communication that may incite aberrant molecular events in association with development of WB myopathy.

A prognostic model based on seven immune-related genes predicts the overall survival of patients with hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a disease with a high incidence and a poor prognosis. Growing amounts of evidence have shown that the immune system plays a critical role in the biological processes of HCC such as progression, recurrence, and metastasis, and some have discussed using it as a weapon against a variety of cancers. However, the impact of immune-related genes (IRGs) on the prognosis of HCC remains unclear.

METHODS

Based on The Cancer Gene Atlas (TCGA) and Immunology Database and Analysis Portal (ImmPort) datasets, we integrated the ribonucleic acid (RNA) sequencing profiles of 424 HCC patients with IRGs to calculate immune-related differentially expressed genes (DEGs). Survival analysis was used to establish a prognostic model of survival- and immune-related DEGs. Based on genomic and clinicopathological data, we constructed a nomogram to predict the prognosis of HCC patients. Gene set enrichment analysis further clarified the signalling pathways of the high-risk and low-risk groups constructed based on the IRGs in HCC. Next, we evaluated the correlation between the risk score and the infiltration of immune cells, and finally, we validated the prognostic performance of this model in the GSE14520 dataset.

RESULTS

A total of 100 immune-related DEGs were significantly associated with the clinical outcomes of patients with HCC. We performed univariate and multivariate least absolute shrinkage and selection operator (Lasso) regression analyses on these genes to construct a prognostic model of seven IRGs (Fatty Acid Binding Protein 6 (FABP6), Microtubule-Associated Protein Tau (MAPT), Baculoviral IAP Repeat Containing 5 (BIRC5), Plexin-A1 (PLXNA1), Secreted Phosphoprotein 1 (SPP1), Stanniocalcin 2 (STC2) and Chondroitin Sulfate Proteoglycan 5 (CSPG5)), which showed better prognostic performance than the tumour/node/metastasis (TNM) staging system. Moreover, we constructed a regulatory network related to transcription factors (TFs) that further unravelled the regulatory mechanisms of these genes. According to the median value of the risk score, the entire TCGA cohort was divided into high-risk and low-risk groups, and the low-risk group had a better overall survival (OS) rate. To predict the OS rate of HCC, we established a gene- and clinical factor-related nomogram. The receiver operating characteristic (ROC) curve, concordance index (C-index) and calibration curve showed that this model had moderate accuracy. The correlation analysis between the risk score and the infiltration of six common types of immune cells showed that the model could reflect the state of the immune microenvironment in HCC tumours.

CONCLUSION

Our IRG prognostic model was shown to have value in the monitoring, treatment, and prognostic assessment of HCC patients and could be used as a survival prediction tool in the near future.

A Meta-Analysis of the Effects of High-LET Ionizing Radiations in Human Gene Expression

The use of high linear energy transfer (LET) ionizing radiation (IR) is progressively being incorporated in radiation therapy due to its precise dose localization and high relative biological effectiveness. At the same time, these benefits of particle radiation become a high risk for astronauts in the case of inevitable cosmic radiation exposure. Nonetheless, DNA Damage Response (DDR) activated via complex DNA damage in healthy tissue, occurring from such types of radiation, may be instrumental in the induction of various chronic and late effects. An approach to elucidating the possible underlying mechanisms is studying alterations in gene expression. To this end, we identified differentially expressed genes (DEGs) in high Z and high energy (HZE) particle-, γ-ray- and X-ray-exposed healthy human tissues, utilizing microarray data available in public repositories. Differential gene expression analysis (DGEA) was conducted using the R programming language. Consequently, four separate meta-analyses were conducted, after DEG lists were grouped depending on radiation type, radiation dose and time of collection post-irradiation. To highlight the biological background of each meta-analysis group, functional enrichment analysis and biological network construction were conducted. For HZE particle exposure at 8–24 h post-irradiation, the most interesting finding is the variety of DNA repair mechanisms that were downregulated, a fact that is probably correlated with complex DNA damage formation. Simultaneously, after X-ray exposure during the same hours after irradiation, DNA repair mechanisms continue to take place. Finally, in a further comparison of low- and high-LET radiation effects, the most prominent result is that autophagy mechanisms seem to persist and that adaptive immune induction seems to be present. Such bioinformatics approaches may aid in obtaining an overview of the cellular response to high-LET particles. Understanding these response mechanisms can consequently aid in the development of countermeasures for future space missions and ameliorate heavy ion treatments.

Differentially methylated regions (DMRs) in PON3 gene between responders and non-responders to a weight loss dietary intervention: a new tool for precision management of obesity

Differentially methylated regions (DMR) are genomic regions with different methylation status. The aim of this research was to identify DMRs in subjects with obesity that predict the response to a weight-loss dietary intervention and its association with metabolic variables. Based on the change in body mass index (BMI), 201 subjects with overweight and obesity were categorized in tertiles according to their response to a hypocaloric diet: Responders (R; n = 64) and Non-Responders (NR; n = 63). The R group lost 4.55 ± 0.91 BMI units (kg/m²) and the NR group lost 1.95 ± 0.73 kg/m² (p < 0.001). DNA methylation was analysed in buffy coat through a methylation array at baseline. DMRs were analysed using a function of ChAMP (Chip Analysis Methylation Pipeline) in R software. Baseline DNA methylation analysis between R and NR exhibited a DMR located at paraoxonase 3 gene (PON3) consisting of 13 CpG sites, eleven of them significantly hypermethylated in R. To analyse the implication of these 11 CpGs on weight loss, a z-score was performed as a measure of DMR methylation. This analysis showed a correlation between PON3 DNA methylation and BMI loss. This z-score negatively correlated with PON3 protein serum levels. Total paraoxonase activity in serum was not different between groups, but PON enzymatic activity positively correlated with oxidized LDL levels. The present study identified a DMR within PON3 gene that is related to PON3 protein levels in serum, and that could be used as a potential biomarker to predict the response to weight-loss dietary interventions.