Automated machine learning for genome wide association studies.

MOTIVATION: Genome-wide association studies (GWAS) present several computational and statistical challenges for their data analysis, including knowledge discovery, interpretability, and translation to clinical practice. RESULTS: We develop, apply, and comparatively evaluate an automated machine learning (AutoML) approach, customized for genomic data that delivers reliable predictive and diagnostic models, the set of genetic variants that are important for predictions (called a biosignature), and an estimate of the out-of-sample predictive power. This AutoML approach discovers variants with higher predictive performance compared to standard GWAS methods, computes an individual risk prediction score, generalizes to new, unseen data, is shown to better differentiate causal variants from other highly correlated variants, and enhances knowledge discovery and interpretability by reporting multiple equivalent biosignatures. AVAILABILITY AND IMPLEMENTATION: Code for this study is available at: https://github.com/mensxmachina/autoML-GWAS. JADBio offers a free version at: https://jadbio.com/sign-up/. SNP data can be downloaded from the EGA repository (https://ega-archive.org/). PRS data are found at: https://www.aicrowd.com/challenges/opensnp-height-prediction. Simulation data to study population structure can be found at: https://easygwas.ethz.ch/data/public/dataset/view/1/.

A robust clustering strategy for stratification unveils unique patient subgroups in acutely decompensated cirrhosis.

BACKGROUND: Patient heterogeneity poses significant challenges for managing individuals and designing clinical trials, especially in complex diseases. Existing classifications rely on outcome-predicting scores, potentially overlooking crucial elements contributing to heterogeneity without necessarily impacting prognosis. METHODS: To address patient heterogeneity, we developed ClustALL, a computational pipeline that simultaneously faces diverse clinical data challenges like mixed types, missing values, and collinearity. ClustALL enables the unsupervised identification of patient stratifications while filtering for stratifications that are robust against minor variations in the population (population-based) and against limited adjustments in the algorithm's parameters (parameter-based). RESULTS: Applied to a European cohort of patients with acutely decompensated cirrhosis (n = 766), ClustALL identified five robust stratifications, using only data at hospital admission. All stratifications included markers of impaired liver function and number of organ dysfunction or failure, and most included precipitating events. When focusing on one of these stratifications, patients were categorized into three clusters characterized by typical clinical features; notably, the 3-cluster stratification showed a prognostic value. Re-assessment of patient stratification during follow-up delineated patients' outcomes, with further improvement of the prognostic value of the stratification. We validated these findings in an independent prospective multicentre cohort of patients from Latin America (n = 580). CONCLUSIONS: By applying ClustALL to patients with acutely decompensated cirrhosis, we identified three patient clusters. Following these clusters over time offers insights that could guide future clinical trial design. ClustALL is a novel and robust stratification method capable of addressing the multiple challenges of patient stratification in most complex diseases.

Thymic function and survival at advance ages in nursing home residents from Southern Italy.

BACKGROUND: Immunosenescence is a complex process characterized by an age-related remodelling of immune system. The prominent effects of the immunosenescence process is the thymic involution and, consequently, the decreased numbers and functions of T cells. Since thymic involution results in a collapse of the T-cell receptor (TCR) repertoire, a reliable biomarker of its activity is represented by the quantification of signal joint T-cell receptor rearrangement excision circles (sjTRECs) levels. Although it is reasonable to think that thymic function could play a crucial role on elderly survival, only a few studies investigated the relationship between an accurate measurement of human thymic function and survival at old ages. METHODS AND FINDINGS: By quantifying the amount sjTRECs by real-time polymerase chain reaction (PCR), the decrease in thymic output in 241 nursing home residents from Calabria (Southern Italy) was evaluated to investigate the relationship between thymic function and survival at old ages. We found that low sjTREC levels were associated with a significant increased risk of mortality at older ages. Nursing home residents with lower sjTREC exhibit a near 2-fold increase in mortality risk compared to those with sjTREC levels in a normal range. CONCLUSION: Thymic function failure is an independent predictor of mortality among elderly nursing home residents. sjTREC represents a biomarker of effective ageing as its blood levels could anticipate individuals at high risk of negative health outcomes. The identification of these subjects is crucial to manage older people's immune function and resilience, such as, for instance, to plan more efficient vaccinal campaigns in older populations.

An ELOVL2-Based Epigenetic Clock for Forensic Age Prediction: A Systematic Review.

The prediction of chronological age from methylation-based biomarkers represents one of the most promising applications in the field of forensic sciences. Age-prediction models developed so far are not easily applicable for forensic caseworkers. Among the several attempts to pursue this objective, the formulation of single-locus models might represent a good strategy. The present work aimed to develop an accurate single-locus model for age prediction exploiting ELOVL2, a gene for which epigenetic alterations are most highly correlated with age. We carried out a systematic review of different published pyrosequencing datasets in which methylation of the ELOVL2 promoter was analysed to formulate age prediction models. Nine of these, with available datasets involving 2298 participants, were selected. We found that irrespective of which model was adopted, a very strong relationship between ELOVL2 methylation levels and age exists. In particular, the model giving the best age-prediction accuracy was the gradient boosting regressor with a prediction error of about 5.5 years. The findings reported here strongly support the use of ELOVL2 for the formulation of a single-locus epigenetic model, but the inclusion of additional, non-redundant markers is a fundamental requirement to apply a molecular model to forensic applications with more robust results.

Feedforward regulation of Myc coordinates lineage-specific with housekeeping gene expression during B cell progenitor cell differentiation.

The differentiation of self-renewing progenitor cells requires not only the regulation of lineage- and developmental stage-specific genes but also the coordinated adaptation of housekeeping functions from a metabolically active, proliferative state toward quiescence. How metabolic and cell-cycle states are coordinated with the regulation of cell type-specific genes is an important question, because dissociation between differentiation, cell cycle, and metabolic states is a hallmark of cancer. Here, we use a model system to systematically identify key transcriptional regulators of Ikaros-dependent B cell-progenitor differentiation. We find that the coordinated regulation of housekeeping functions and tissue-specific gene expression requires a feedforward circuit whereby Ikaros down-regulates the expression of Myc. Our findings show how coordination between differentiation and housekeeping states can be achieved by interconnected regulators. Similar principles likely coordinate differentiation and housekeeping functions during progenitor cell differentiation in other cell lineages.

Age-related changes in medial septal cholinergic and GABAergic projection neurons and hippocampal neurotransmitter receptors: relationship with memory impairment.

The hippocampus, which provides cognitive functions, has been shown to become highly vulnerable during aging. One important modulator of the hippocampal neural network is the medial septum (MS). The present study attempts to determine how age-related mnemonic dysfunction is associated with neurochemical changes in the septohippocampal (SH) system, using behavioral and immunochemical experiments performed on young-adult, middle-aged and aged rats. According to these behavioral results, the aged and around 52.8% of middle-aged rats (within the "middle-aged-impaired" sub-group) showed both impaired spatial reference memory in the Morris water maze and habituation in the open field. Immunohistochemical studies revealed a significant decrease in the number of MS choline acetyltransferase immunoreactive cells in the aged and all middle-aged rats, in comparison to the young; however the number of gamma-aminobutyric acid-ergic (GABAergic) parvalbumin immunoreactive cells was higher in middle-aged-impaired and older rats compared to young and middle-aged-unimpaired rats. Western Blot analysis moreover showed a decrease in the level of expression of cholinergic, GABAergic and glutamatergic receptors in the hippocampus of middle-aged-impaired and aged rats in contrast to middle-aged-unimpaired and young rats. The present results demonstrate for the first time that a decrease in the expression level of hippocampal receptors in naturally aged rats with impaired cognitive abilities occurs in parallel with an increase in the number of GABAergic neurons in the MS, and it highlights the particular importance of inhibitory signaling in the SH network for memory function.

Myo-Inositol Limits Kainic Acid-Induced Epileptogenesis in Rats.

Epilepsy is a severe neurological disease characterized by spontaneous recurrent seizures (SRS). A complex pathophysiological process referred to as epileptogenesis transforms a normal brain into an epileptic one. Prevention of epileptogenesis is a subject of intensive research. Currently, there are no clinically approved drugs that can act as preventive medication. Our previous studies have revealed highly promising antiepileptogenic properties of a compound-myo-inositol (MI) and the present research broadens previous results and demonstrates the long-term disease-modifying effect of this drug, as well as the amelioration of cognitive comorbidities. For the first time, we show that long-term treatment with MI: (i) decreases the frequency and duration of electrographic SRS in the hippocampus; (ii) has an ameliorating effect on spatial learning and memory deficit associated with epileptogenesis, and (iii) attenuates cell loss in the hippocampus. MI treatment also alters the expression of the glial fibrillary acidic protein, LRRC8A subunit of volume-regulated anion channels, and protein tyrosine phosphatase receptor type R, all expected to counteract the epileptogenesis. All these effects are still present even 4 weeks after MI treatment ceased. This suggests that MI may exert multiple actions on various epileptogenesis-associated changes in the brain and, therefore, could be considered as a candidate target for prevention of epileptogenesis.

Short and Long Time Bloodstains Age Determination by Colorimetric Analysis: A Pilot Study.

Bloodstains found at crime scenes represent a crucial source of information for investigative purposes. However, in forensic practice, no technique is currently used to estimate the time from deposition of bloodstains. This preliminary study focuses on the age estimation of bloodstains by exploiting the color variations over time due to the oxidation of the blood. For this purpose, we used a colorimetric methodology in order to easily obtain objective, univocal and reproducible results. We developed two bloodstain age prediction algorithms: a short-term and a long-term useful model for the first 24h and 60 days, respectively. Both models showed high levels of classification accuracy, particularly for the long-term model. Although a small-scale study, these results improve the potential application of colorimetric analysis in the time-line reconstruction of violent criminal events.

Challenges in the Multivariate Analysis of Mass Cytometry Data: The Effect of Randomization.

Cytometry by time-of-flight (CyTOF) has emerged as a high-throughput single cell technology able to provide large samples of protein readouts. Already, there exists a large pool of advanced high-dimensional analysis algorithms that explore the observed heterogeneous distributions making intriguing biological inferences. A fact largely overlooked by these methods, however, is the effect of the established data preprocessing pipeline to the distributions of the measured quantities. In this article, we focus on randomization, a transformation used for improving data visualization, which can negatively affect multivariate data analysis methods such as dimensionality reduction, clustering, and network reconstruction algorithms. Our results indicate that randomization should be used only for visualization purposes, but not in conjunction with high-dimensional analytical tools. © 2019 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

SCENERY: a web application for (causal) network reconstruction from cytometry data.

Flow and mass cytometry technologies can probe proteins as biological markers in thousands of individual cells simultaneously, providing unprecedented opportunities for reconstructing networks of protein interactions through machine learning algorithms. The network reconstruction (NR) problem has been well-studied by the machine learning community. However, the potentials of available methods remain largely unknown to the cytometry community, mainly due to their intrinsic complexity and the lack of comprehensive, powerful and easy-to-use NR software implementations specific for cytometry data. To bridge this gap, we present Single CEll NEtwork Reconstruction sYstem (SCENERY), a web server featuring several standard and advanced cytometry data analysis methods coupled with NR algorithms in a user-friendly, on-line environment. In SCENERY, users may upload their data and set their own study design. The server offers several data analysis options categorized into three classes of methods: data (pre)processing, statistical analysis and NR. The server also provides interactive visualization and download of results as ready-to-publish images or multimedia reports. Its core is modular and based on the widely-used and robust R platform allowing power users to extend its functionalities by submitting their own NR methods. SCENERY is available at scenery.csd.uoc.gr or http://mensxmachina.org/en/software/.

Feature selection for high-dimensional temporal data.

BACKGROUND: Feature selection is commonly employed for identifying collectively-predictive biomarkers and biosignatures; it facilitates the construction of small statistical models that are easier to verify, visualize, and comprehend while providing insight to the human expert. In this work we extend established constrained-based, feature-selection methods to high-dimensional "omics" temporal data, where the number of measurements is orders of magnitude larger than the sample size. The extension required the development of conditional independence tests for temporal and/or static variables conditioned on a set of temporal variables. RESULTS: The algorithm is able to return multiple, equivalent solution subsets of variables, scale to tens of thousands of features, and outperform or be on par with existing methods depending on the analysis task specifics. CONCLUSIONS: The use of this algorithm is suggested for variable selection with high-dimensional temporal data.

A comparative evaluation of data-merging and meta-analysis methods for reconstructing gene-gene interactions.

BACKGROUND: We address the problem of integratively analyzing multiple gene expression, microarray datasets in order to reconstruct gene-gene interaction networks. Integrating multiple datasets is generally believed to provide increased statistical power and to lead to a better characterization of the system under study. However, the presence of systematic variation across different studies makes network reverse-engineering tasks particularly challenging. We contrast two approaches that have been frequently used in the literature for addressing systematic biases: meta-analysis methods, which first calculate opportune statistics on single datasets and successively summarize them, and data-merging methods, which directly analyze the pooled data after removing eventual biases. This comparative evaluation is performed on both synthetic and real data, the latter consisting of two manually curated microarray compendia comprising several E. coli and Yeast studies, respectively. Furthermore, the reconstruction of the regulatory network of the transcription factor Ikaros in human Peripheral Blood Mononuclear Cells (PBMCs) is presented as a case-study. RESULTS: The meta-analysis and data-merging methods included in our experimentations provided comparable performances on both synthetic and real data. Furthermore, both approaches outperformed (a) the naïve solution of merging data together ignoring possible biases, and (b) the results that are expected when only one dataset out of the available ones is analyzed in isolation. Using correlation statistics proved to be more effective than using p-values for correctly ranking candidate interactions. The results from the PBMC case-study indicate that the findings of the present study generalize to different types of network reconstruction algorithms. CONCLUSIONS: Ignoring the systematic variations that differentiate heterogeneous studies can produce results that are statistically indistinguishable from random guessing. Meta-analysis and data merging methods have proved equally effective in addressing this issue, and thus researchers may safely select the approach that best suit their specific application.

Long-term effects of myo-inositol on traumatic brain injury: Epigenomic and transcriptomic studies.

Background and purpose: Traumatic brain injury (TBI) and its consequences remain great challenges for neurology. Consequences of TBI are associated with various alterations in the brain but little is known about long-term changes of epigenetic DNA methylation patterns. Moreover, nothing is known about potential treatments that can alter these epigenetic changes in beneficial ways. Therefore, we have examined myo-inositol (MI), which has positive effects on several pathological conditions. Methods: TBI was induced in mice by controlled cortical impact (CCI). One group of CCI animals received saline injections for two months (TBI+SAL), another CCI group received MI treatment (TBI+MI) for the same period and one group served as a sham-operated control. Mice were sacrificed 4 months after CCI and changes in DNA methylome and transcriptomes were examined. Results: For the first time we: (i) provide comprehensive map of long-term DNA methylation changes after CCI in the hippocampus; (ii) identify differences by methylation sites between the groups; (iii) characterize transcriptome changes; (iv) provide association between DNA methylation sites and gene expression. MI treatment is linked with upregulation of genes covering 33 biological processes, involved in immune response and inflammation. In support of these findings, we have shown that expression of BATF2, a transcription factor involved in immune-regulatory networks, is upregulated in the hippocampus of the TBI+MI group where the BATF2 gene is demethylated. Conclusion: TBI is followed by long-term epigenetic and transcriptomic changes in hippocampus. MI treatment has a significant effect on these processes by modulation of immune response and biological pathways of inflammation.

Src-NADH dehydrogenase subunit 2 complex and recognition memory of imprinting in domestic chicks.

Src is a non-receptor tyrosine kinase participating in a range of neuronal processes, including synaptic plasticity. We have recently shown that the amounts of total Src and its two phosphorylated forms, at tyrosine-416 (activated) and tyrosine-527 (inhibited), undergoes time-dependent, region-specific learning-related changes in the domestic chick forebrain after visual imprinting. These changes occur in the intermediate medial mesopallium (IMM), a site of memory formation for visual imprinting, but not the posterior pole of the nidopallium (PPN), a control brain region not involved in imprinting. Src interacts with mitochondrial genome-coded NADH dehydrogenase subunit 2 (NADH2), a component of mitochondrial respiratory complex I. This interaction occurs at brain excitatory synapses bearing NMDA glutamate receptors. The involvement of Src-NADH2 complexes in learning and memory is not yet explored. We show for the first time that, independently of changes in total Src or total NADH2, NADH2 bound to Src immunoprecipitated from the P2 plasma membrane-mitochondrial fraction: (i) is increased in a learning-related manner in the left IMM 1 h after the end of training; (ii), is decreased in the right IMM in a learning-related way 24 h after training. These changes occurred in the IMM but not the PPN. They are attributable to learning occurring during training rather than a predisposition to learn. Learning-related changes in Src-bound NADH2 are thus time- and region-dependent.

Improving Lung Cancer Screening Selection: The HUNT Lung Cancer Risk Model for Ever-Smokers Versus the NELSON and 2021 United States Preventive Services Task Force Criteria in the Cohort of Norway: A Population-Based Prospective Study.

Background: Improving the method for selecting participants for lung cancer (LC) screening is an urgent need. Here, we compared the performance of the Helseundersøkelsen i Nord-Trøndelag (HUNT) Lung Cancer Model (HUNT LCM) versus the Dutch-Belgian lung cancer screening trial (Nederlands-Leuvens Longkanker Screenings Onderzoek (NELSON)) and 2021 United States Preventive Services Task Force (USPSTF) criteria regarding LC risk prediction and efficiency. Methods: We used linked data from 10 Norwegian prospective population-based cohorts, Cohort of Norway. The study included 44,831 ever-smokers, of which 686 (1.5%) patients developed LC; the median follow-up time was 11.6 years (0.01-20.8 years). Results: Within 6 years, 222 (0.5%) individuals developed LC. The NELSON and 2021 USPSTF criteria predicted 37.4% and 59.5% of the LC cases, respectively. By considering the same number of individuals as the NELSON and 2021 USPSTF criteria selected, the HUNT LCM increased the LC prediction rate by 41.0% and 12.1%, respectively. The HUNT LCM significantly increased sensitivity (p < 0.001 and p = 0.028), and reduced the number needed to predict one LC case (29 versus 40, p < 0.001 and 36 versus 40, p = 0.02), respectively. Applying the HUNT LCM 6-year 0.98% risk score as a cutoff (14.0% of ever-smokers) predicted 70.7% of all LC, increasing LC prediction rate with 89.2% and 18.9% versus the NELSON and 2021 USPSTF, respectively (both p < 0.001). Conclusions: The HUNT LCM was significantly more efficient than the NELSON and 2021 USPSTF criteria, improving the prediction of LC diagnosis, and may be used as a validated clinical tool for screening selection.

Promising microRNAs in pre-diagnostic serum associated with lung cancer up to eight years before diagnosis: a HUNT study.

INTRODUCTION: Blood biomarkers for early detection of lung cancer (LC) are in demand. There are few studies of the full microRNome in serum of asymptomatic subjects that later develop LC. Here we searched for novel microRNA biomarkers in blood from non-cancer, ever-smokers populations up to eight years before diagnosis. METHODS: Serum samples from 98,737 subjects from two prospective population studies, HUNT2 and HUNT3, were considered initially. Inclusion criteria for cases were: ever-smokers; no known cancer at study entrance; 0-8 years from blood sampling to LC diagnosis. Each future LC case had one control matched to sex, age at study entrance, pack-years, smoking cessation time, and similar HUNT Lung Cancer Model risk score. A total of 240 and 72 serum samples were included in the discovery (HUNT2) and validation (HUNT3) datasets, respectively, and analysed by next-generation sequencing. The validated serum microRNAs were also tested in two pre-diagnostic plasma datasets from the prospective population studies NOWAC (n = 266) and NSHDS (n = 258). A new model adding clinical variables was also developed and validated. RESULTS: Fifteen unique microRNAs were discovered and validated in the pre-diagnostic serum datasets when all cases were contrasted against all controls, all with AUC > 0.60. In combination as a 15-microRNAs signature, the AUC reached 0.708 (discovery) and 0.703 (validation). A non-small cell lung cancer signature of six microRNAs showed AUC 0.777 (discovery) and 0.806 (validation). Combined with clinical variables of the HUNT Lung Cancer Model (age, gender, pack-years, daily cough parts of the year, hours of indoor smoke exposure, quit time in years, number of cigarettes daily, body mass index (BMI)) the AUC reached 0.790 (discovery) and 0.833 (validation). These results could not be validated in the plasma samples. CONCLUSION: There were a few significantly differential expressed microRNAs in serum up to eight years before diagnosis. These promising microRNAs alone, in concert, or combined with clinical variables have the potential to serve as early diagnostic LC biomarkers. Plasma is not suitable for this analysis. Further validation in larger prospective serum datasets is needed.

The HUNT lung-SNP model: genetic variants plus clinical variables improve lung cancer risk assessment over clinical models.

PURPOSE: The HUNT Lung Cancer Model (HUNT LCM) predicts individualized 6-year lung cancer (LC) risk among individuals who ever smoked cigarettes with high precision based on eight clinical variables. Can the performance be improved by adding genetic information? METHODS: A polygenic model was developed in the prospective Norwegian HUNT2 study with clinical and genotype data of individuals who ever smoked cigarettes (n = 30749, median follow up 15.26 years) where 160 LC were diagnosed within six years. It included the variables of the original HUNT LCM plus 22 single nucleotide polymorphisms (SNPs) highly associated with LC. External validation was performed in the prospective Norwegian Tromsø Study (n = 2663). RESULTS: The novel HUNT Lung-SNP model significantly improved risk ranking of individuals over the HUNT LCM in both HUNT2 (p < 0.001) and Tromsø (p < 0.05) cohorts. Furthermore, detection rate (number of participants selected to detect one LC case) was significantly better for the HUNT Lung-SNP vs. HUNT LCM in both cohorts (42 vs. 48, p = 0.003 and 11 vs. 14, p = 0.025, respectively) as well as versus the NLST, NELSON and 2021 USPSTF criteria. The area under the receiver operating characteristic curve (AUC) was higher for the HUNT Lung-SNP in both cohorts, but significant only in HUNT2 (AUC 0.875 vs. 0.844, p < 0.001). However, the integrated discrimination improvement index (IDI) indicates a significant improvement of LC risk stratification by the HUNT Lung-SNP in both cohorts (IDI 0.019, p < 0.001 (HUNT2) and 0.013, p < 0.001 (Tromsø)). CONCLUSION: The HUNT Lung-SNP model could have a clinical impact on LC screening and has the potential to replace the HUNT LCM as well as the NLST, NELSON and 2021 USPSTF criteria in a screening setting. However, the model should be further validated in other populations and evaluated in a prospective trial setting.

Physical performance strongly predicts all-cause mortality risk in a real-world population of older diabetic patients: machine learning approach for mortality risk stratification.

Background: Prognostic risk stratification in older adults with type 2 diabetes (T2D) is important for guiding decisions concerning advance care planning. Materials and methods: A retrospective longitudinal study was conducted in a real-world sample of older diabetic patients afferent to the outpatient facilities of the Diabetology Unit of the IRCCS INRCA Hospital of Ancona (Italy). A total of 1,001 T2D patients aged more than 70 years were consecutively evaluated by a multidimensional geriatric assessment, including physical performance evaluated using the Short Physical Performance Battery (SPPB). The mortality was assessed during a 5-year follow-up. We used the automatic machine-learning (AutoML) JADBio platform to identify parsimonious mathematical models for risk stratification. Results: Of 977 subjects included in the T2D cohort, the mean age was 76.5 (SD: 4.5) years and 454 (46.5%) were men. The mean follow-up time was 53.3 (SD:15.8) months, and 209 (21.4%) patients died by the end of the follow-up. The JADBio AutoML final model included age, sex, SPPB, chronic kidney disease, myocardial ischemia, peripheral artery disease, neuropathy, and myocardial infarction. The bootstrap-corrected concordance index (c-index) for the final model was 0.726 (95% CI: 0.687-0.763) with SPPB ranked as the most important predictor. Based on the penalized Cox regression model, the risk of death per unit of time for a subject with an SPPB score lower than five points was 3.35 times that for a subject with a score higher than eight points (P-value <0.001). Conclusion: Assessment of physical performance needs to be implemented in clinical practice for risk stratification of T2D older patients.

Learning biologically-interpretable latent representations for gene expression data: Pathway Activity Score Learning Algorithm.

Molecular gene-expression datasets consist of samples with tens of thousands of measured quantities (i.e., high dimensional data). However, lower-dimensional representations that retain the useful biological information do exist. We present a novel algorithm for such dimensionality reduction called Pathway Activity Score Learning (PASL). The major novelty of PASL is that the constructed features directly correspond to known molecular pathways (genesets in general) and can be interpreted as pathway activity scores. Hence, unlike PCA and similar methods, PASL's latent space has a fairly straightforward biological interpretation. PASL is shown to outperform in predictive performance the state-of-the-art method (PLIER) on two collections of breast cancer and leukemia gene expression datasets. PASL is also trained on a large corpus of 50000 gene expression samples to construct a universal dictionary of features across different tissues and pathologies. The dictionary validated on 35643 held-out samples for reconstruction error. It is then applied on 165 held-out datasets spanning a diverse range of diseases. The AutoML tool JADBio is employed to show that the predictive information in the PASL-created feature space is retained after the transformation. The code is available at https://github.com/mensxmachina/PASL.

scAEGAN: Unification of single-cell genomics data by adversarial learning of latent space correspondences.

Recent progress in Single-Cell Genomics has produced different library protocols and techniques for molecular profiling. We formulate a unifying, data-driven, integrative, and predictive methodology for different libraries, samples, and paired-unpaired data modalities. Our design of scAEGAN includes an autoencoder (AE) network integrated with adversarial learning by a cycleGAN (cGAN) network. The AE learns a low-dimensional embedding of each condition, whereas the cGAN learns a non-linear mapping between the AE representations. We evaluate scAEGAN using simulated data and real scRNA-seq datasets, different library preparations (Fluidigm C1, CelSeq, CelSeq2, SmartSeq), and several data modalities as paired scRNA-seq and scATAC-seq. The scAEGAN outperforms Seurat3 in library integration, is more robust against data sparsity, and beats Seurat 4 in integrating paired data from the same cell. Furthermore, in predicting one data modality from another, scAEGAN outperforms Babel. We conclude that scAEGAN surpasses current state-of-the-art methods and unifies integration and prediction challenges.