Mutational Landscape of Secondary Glioblastoma Guides MET-Targeted Trial in Brain Tumor.

Low-grade gliomas almost invariably progress into secondary glioblastoma (sGBM) with limited therapeutic option and poorly understood mechanism. By studying the mutational landscape of 188 sGBMs, we find significant enrichment of TP53 mutations, somatic hypermutation, MET-exon-14-skipping (METex14), PTPRZ1-MET (ZM) fusions, and MET amplification. Strikingly, METex14 frequently co-occurs with ZM fusion and is present in ∼14% of cases with significantly worse prognosis. Subsequent studies show that METex14 promotes glioma progression by prolonging MET activity. Furthermore, we describe a MET kinase inhibitor, PLB-1001, that demonstrates remarkable potency in selectively inhibiting MET-altered tumor cells in preclinical models. Importantly, this compound also shows blood-brain barrier permeability and is subsequently applied in a phase I clinical trial that enrolls MET-altered chemo-resistant glioma patients. Encouragingly, PLB-1001 achieves partial response in at least two advanced sGBM patients with rarely significant side effects, underscoring the clinical potential for precisely treating gliomas using this therapy.

The African Society of Human Genetics successfully launches global data science workshops.

To accelerate the impact of African genomics on human health, data science skills and awareness of Africa's rich genetic diversity must be strengthened globally. We describe the first African genomics data science workshop, implemented by the African Society of Human Genetics (AfSHG) and international partners, providing a framework for future workshops.

CHARR efficiently estimates contamination from DNA sequencing data.

DNA sample contamination is a major issue in clinical and research applications of whole-genome and -exome sequencing. Even modest levels of contamination can substantially affect the overall quality of variant calls and lead to widespread genotyping errors. Currently, popular tools for estimating the contamination level use short-read data (BAM/CRAM files), which are expensive to store and manipulate and often not retained or shared widely. We propose a metric to estimate DNA sample contamination from variant-level whole-genome and -exome sequence data called CHARR, contamination from homozygous alternate reference reads, which leverages the infiltration of reference reads within homozygous alternate variant calls. CHARR uses a small proportion of variant-level genotype information and thus can be computed from single-sample gVCFs or callsets in VCF or BCF formats, as well as efficiently stored variant calls in Hail VariantDataset format. Our results demonstrate that CHARR accurately recapitulates results from existing tools with substantially reduced costs, improving the accuracy and efficiency of downstream analyses of ultra-large whole-genome and exome sequencing datasets.

Cloud-based introduction to BASH programming for biologists.

This manuscript describes the development of a resource module that is part of a learning platform named 'NIGMS Sandbox for Cloud-based Learning', https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial authored by National Institute of General Medical Sciences: NIGMS Sandbox: A Learning Platform toward Democratizing Cloud Computing for Biomedical Research at the beginning of this supplement. This module delivers learning materials introducing the utility of the BASH (Bourne Again Shell) programming language for genomic data analysis in an interactive format that uses appropriate cloud resources for data access and analyses. The next-generation sequencing revolution has generated massive amounts of novel biological data from a multitude of platforms that survey an ever-growing list of genomic modalities. These data require significant downstream computational and statistical analyses to glean meaningful biological insights. However, the skill sets required to generate these data are vastly different from the skills required to analyze these data. Bench scientists that generate next-generation data often lack the training required to perform analysis of these datasets and require support from bioinformatics specialists. Dedicated computational training is required to empower biologists in the area of genomic data analysis, however, learning to efficiently leverage a command line interface is a significant barrier in learning how to leverage common analytical tools. Cloud platforms have the potential to democratize access to the technical tools and computational resources necessary to work with modern sequencing data, providing an effective framework for bioinformatics education. This module aims to provide an interactive platform that slowly builds technical skills and knowledge needed to interact with genomics data on the command line in the Cloud. The sandbox format of this module enables users to move through the material at their own pace and test their grasp of the material with knowledge self-checks before building on that material in the next sub-module. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

Translational bioinformatics and data science for biomarker discovery in mental health: an analytical review.

Translational bioinformatics and data science play a crucial role in biomarker discovery as it enables translational research and helps to bridge the gap between the bench research and the bedside clinical applications. Thanks to newer and faster molecular profiling technologies and reducing costs, there are many opportunities for researchers to explore the molecular and physiological mechanisms of diseases. Biomarker discovery enables researchers to better characterize patients, enables early detection and intervention/prevention and predicts treatment responses. Due to increasing prevalence and rising treatment costs, mental health (MH) disorders have become an important venue for biomarker discovery with the goal of improved patient diagnostics, treatment and care. Exploration of underlying biological mechanisms is the key to the understanding of pathogenesis and pathophysiology of MH disorders. In an effort to better understand the underlying mechanisms of MH disorders, we reviewed the major accomplishments in the MH space from a bioinformatics and data science perspective, summarized existing knowledge derived from molecular and cellular data and described challenges and areas of opportunities in this space.

MM optimization: Proximal distance algorithms, path following, and trust regions.

We briefly review the majorization-minimization (MM) principle and elaborate on the closely related notion of proximal distance algorithms, a generic approach for solving constrained optimization problems via quadratic penalties. We illustrate how the MM and proximal distance principles apply to a variety of problems from statistics, finance, and nonlinear optimization. Drawing from our selected examples, we also sketch a few ideas pertinent to the acceleration of MM algorithms: a) structuring updates around efficient matrix decompositions, b) path following in proximal distance iteration, and c) cubic majorization and its connections to trust region methods. These ideas are put to the test on several numerical examples, but for the sake of brevity, we omit detailed comparisons to competing methods. The current article, which is a mix of review and current contributions, celebrates the MM principle as a powerful framework for designing optimization algorithms and reinterpreting existing ones.

How to avoid a local epidemic becoming a global pandemic.

Here, we combine international air travel passenger data with a standard epidemiological model of the initial 3 mo of the COVID-19 pandemic (January through March 2020; toward the end of which the entire world locked down). Using the information available during this initial phase of the pandemic, our model accurately describes the main features of the actual global development of the pandemic demonstrated by the high degree of coherence between the model and global data. The validated model allows for an exploration of alternative policy efficacies (reducing air travel and/or introducing different degrees of compulsory immigration quarantine upon arrival to a country) in delaying the global spread of SARS-CoV-2 and thus is suggestive of similar efficacy in anticipating the spread of future global disease outbreaks. We show that a lesson from the recent pandemic is that reducing air travel globally is more effective in reducing the global spread than adopting immigration quarantine. Reducing air travel out of a source country has the most important effect regarding the spreading of the disease to the rest of the world. Based upon our results, we propose a digital twin as a further developed tool to inform future pandemic decision-making to inform measures intended to control the spread of disease agents of potential future pandemics. We discuss the design criteria for such a digital twin model as well as the feasibility of obtaining access to the necessary online data on international air travel.

A review of genetic variant databases and machine learning tools for predicting the pathogenicity of breast cancer.

Studies continue to uncover contributing risk factors for breast cancer (BC) development including genetic variants. Advances in machine learning and big data generated from genetic sequencing can now be used for predicting BC pathogenicity. However, it is unclear which tool developed for pathogenicity prediction is most suited for predicting the impact and pathogenicity of variant effects. A significant challenge is to determine the most suitable data source for each tool since different tools can yield different prediction results with different data inputs. To this end, this work reviews genetic variant databases and tools used specifically for the prediction of BC pathogenicity. We provide a description of existing genetic variants databases and, where appropriate, the diseases for which they have been established. Through example, we illustrate how they can be used for prediction of BC pathogenicity and discuss their associated advantages and disadvantages. We conclude that the tools that are specialized by training on multiple diverse datasets from different databases for the same disease have enhanced accuracy and specificity and are thereby more helpful to the clinicians in predicting and diagnosing BC as early as possible.

Predicting relative efficiency of amide bond formation using multivariate linear regression.

Amides are ubiquitous in biologically active natural products and commercial drugs. The most common strategy for introducing this functional group is the coupling of a carboxylic acid with an amine, which requires the use of a coupling reagent to facilitate elimination of water. However, the optimal reaction conditions often appear rather arbitrary to the specific reaction. Herein, we report the development of statistical models correlating measured rates to physical organic descriptors to enable the prediction of reaction rates for untested carboxylic acid/amine pairs. The key to the success of this endeavor was the development of an end-to-end data science­based workflow to select a set of coupling partners that are appropriately distributed in chemical space to facilitate statistical model development. By using a parameterization, dimensionality reduction, and clustering protocol, a training set was identified. Reaction rates for a range of carboxylic acid and primary alkyl amine couplings utilizing carbonyldiimidazole (CDI) as the coupling reagent were measured. The collected rates span five orders of magnitude, confirming that the designed training set encompasses a wide range of chemical space necessary for effective model development. Regressing these rates with high-level density functional theory (DFT) descriptors allowed for identification of a statistical model wherein the molecular features of the carboxylic acid are primarily responsible for the observed rates. Finally, out-of-sample amide couplings are used to determine the limitations and effectiveness of the model.

Single-cell-resolution map of human retinal pigment epithelium helps discover subpopulations with differential disease sensitivity.

Regional phenotypic and functional differences in the retinal pigment epithelium (RPE) monolayer have been suggested to account for regional susceptibility in ocular diseases such as age-related macular degeneration (AMD), late-onset retinal degeneration (L-ORD), and choroideremia (CHM). However, a comprehensive description of human topographical RPE diversity is not yet available, thus limiting the understanding of regional RPE diversity and degenerative disease sensitivity in the eye. To develop a complete morphometric RPE map of the human eye, artificial intelligence­based software was trained to recognize, segment, and analyze RPE borders. Five statistically different, concentric RPE subpopulations (P1 to P5) were identified using cell area as a parameter, including a subpopulation (P4) with cell area comparable to that of macular cells in the far periphery of the eye. This work provides a complete reference map of human RPE subpopulations and their location in the eye. In addition, the analysis of cadaver non-AMD and AMD eyes and ultra-widefield fundus images of patients revealed differential vulnerability of the five RPE subpopulations to different retinal diseases.

Estimation of functional diversity and species traits from ecological monitoring data.

The twin crises of climate change and biodiversity loss define a strong need for functional diversity monitoring. While the availability of high-quality ecological monitoring data is increasing, the quantification of functional diversity so far requires the identification of species traits, for which data are harder to obtain. However, the traits that are relevant for the ecological function of a species also shape its performance in the environment and hence, should be reflected indirectly in its spatiotemporal distribution. Thus, it may be possible to reconstruct these traits from a sufficiently extensive monitoring dataset. Here, we use diffusion maps, a deterministic and de facto parameter-free analysis method, to reconstruct a proxy representation of the species' traits directly from monitoring data and use it to estimate functional diversity. We demonstrate this approach with both simulated data and real-world phytoplankton monitoring data from the Baltic Sea. We anticipate that wider application of this approach to existing data could greatly advance the analysis of changes in functional biodiversity.

A hybrid empirical and parametric approach for managing ecosystem complexity: Water quality in Lake Geneva under nonstationary futures.

Severe deterioration of water quality in lakes, characterized by overabundance of algae and declining dissolved oxygen in the deep lake (DOB), was one of the ecological crises of the 20th century. Even with large reductions in phosphorus loading, termed "reoligotrophication," DOB and chlorophyll (CHL) have often not returned to their expected pre-20th-century levels. Concurrently, management of lake health has been confounded by possible consequences of climate change, particularly since the effects of climate are not neatly separable from the effects of eutrophication. Here, using Lake Geneva as an iconic example, we demonstrate a complementary alternative to parametric models for understanding and managing lake systems. This involves establishing an empirically-driven baseline that uses supervised machine learning to capture the changing interdependencies among biogeochemical variables and then combining the empirical model with a more conventional equation-based model of lake physics to predict DOB over decadal time-scales. The hybrid model not only leads to substantially better forecasts, but also to a more actionable description of the emergent rates and processes (biogeochemical, ecological, etc.) that drive water quality. Notably, the hybrid model suggests that the impact of a moderate 3°C air temperature increase on water quality would be on the same order as the eutrophication of the previous century. The study provides a template and a practical path forward to cope with shifts in ecology to manage environmental systems for non-analogue futures.

A critique and examination of the polysocial risk score approach: predicting cognition in the Health and Retirement Study.

Polysocial risk scores were recently proposed as a strategy for improving the clinical relevance of knowledge about social determinants of health. Our objective in this study was to assess whether the polysocial risk score model improves prediction of cognition and all-cause mortality in middle-aged and older adults beyond simpler models including a smaller set of key social determinants of health. We used a sample of 13 773 individuals aged ≥50 years at baseline from the 2006-2018 waves of the Health and Retirement Study, a US population-based longitudinal cohort study. Four linear mixed models were compared: 2 simple models including a priori-selected covariates and 2 polysocial risk score models which used least absolute shrinkage and selection operator (LASSO) regularization to select covariates among 9 or 21 candidate social predictors. All models included age. Predictive accuracy was assessed via R2 and root mean-squared prediction error (RMSPE) using training/test split validation and cross-validation. For predicting cognition, the simple model including age, race, sex, and education had an R2 value of 0.31 and an RMSPE of 0.880. Compared with this, the most complex polysocial risk score selected 12 predictors (R2 = 0.35 and RMSPE = 0.858; 2.2% improvement). For all-cause mortality, the simple model including age, race, sex, and education had an area under the receiver operating characteristic curve (AUROC) of 0.747, while the most complex polysocial risk score did not demonstrate improved performance (AUROC = 0.745). Models built on a smaller set of key social determinants performed comparably to models built on a more complex set of social "risk factors."

The modern alchemy of clinical pathology: turning the output of microbiology laboratory operations into gold.

The clinical microbiology laboratory generates a huge amount of high-quality data that play a vital role in clinical care. With proper extraction, cleaning, analysis, and validation pipelines, these data can serve multiple other purposes that include supporting laboratory operations, understanding local epidemiology, informing hospital-specific policies, and public health surveillance. In this review, I use one of the core activities of the microbiology laboratory, antimicrobial susceptibility testing (AST), to illustrate several potential applications of next-generation data analytics. The first involves continuous monitoring of commercial AST systems using comparisons of minimum inhibitory concentration (MIC) distributions over time to trigger re-verification when statistically significant differences are detected. An extension of this is temporal analysis of joint MIC distributions to understand performance for multidrug-resistant organisms. More sophisticated analyses involve linking microbiologic data to clinical metadata to gain insight into the clinical validity of AST data and to inform treatment policies. The elements of a robust, validated analysis engine using routine data streams already exist, but numerous challenges must be overcome to make it a reality. Most importantly, it will require the sustained collaboration and advocacy of hospital leadership, microbiologists, clinicians, antimicrobial stewardship, data scientists, and regulatory agencies. Though no small feat, achieving this vision would provide an important resource for microbiology laboratories facing a rapidly evolving practice landscape and further cement its role as an integral part of a learning health system.

Automated analysis of the total choline resonance peak in breast proton magnetic resonance spectroscopy.

The aim of the current study was to compare the performance of fully automated software with human expert interpretation of single-voxel proton magnetic resonance spectroscopy (1H-MRS) spectra in the assessment of breast lesions. Breast magnetic resonance imaging (MRI) (including contrast-enhanced T1-weighted, T2-weighted, and diffusion-weighted imaging) and 1H-MRS images of 74 consecutive patients were acquired on a 3-T positron emission tomography-MRI scanner then automatically imported into and analyzed by SpecTec-ULR 1.1 software (LifeTec Solutions GmbH). All ensuing 117 spectra were additionally independently analyzed and interpreted by two blinded radiologists. Histopathology of at least 24 months of imaging follow-up served as the reference standard. Nonparametric Spearman's correlation coefficients for all measured parameters (signal-to-noise ratio [SNR] and integral of total choline [tCho]), Passing and Bablok regression, and receiver operating characteristic analysis, were calculated to assess test diagnostic performance, as well as to compare automated with manual reading. Based on 117 spectra of 74 patients, the area under the curve for tCho SNR and integrals ranged from 0.768 to 0.814 and from 0.721 to 0.784 to distinguish benign from malignant tissue, respectively. Neither method displayed significant differences between measurements (automated vs. human expert readers, p > 0.05), in line with the results from the univariate Spearman's rank correlation coefficients, as well as the Passing and Bablok regression analysis. It was concluded that this pilot study demonstrates that 1H-MRS data from breast MRI can be automatically exported and interpreted by SpecTec-ULR 1.1 software. The diagnostic performance of this software was not inferior to human expert readers.

Interfacing data science with cell therapy manufacturing: where we are and where we need to be.

Although several cell-based therapies have received FDA approval, and others are showing promising results, scalable, and quality-driven reproducible manufacturing of therapeutic cells at a lower cost remains challenging. Challenges include starting material and patient variability, limited understanding of manufacturing process parameter effects on quality, complex supply chain logistics, and lack of predictive, well-understood product quality attributes. These issues can manifest as increased production costs, longer production times, greater batch-to-batch variability, and lower overall yield of viable, high-quality cells. The lack of data-driven insights and decision-making in cell manufacturing and delivery is an underlying commonality behind all these problems. Data collection and analytics from discovery, preclinical and clinical research, process development, and product manufacturing have not been sufficiently utilized to develop a "systems" understanding and identify actionable controls. Experience from other industries shows that data science and analytics can drive technological innovations and manufacturing optimization, leading to improved consistency, reduced risk, and lower cost. The cell therapy manufacturing industry will benefit from implementing data science tools, such as data-driven modeling, data management and mining, AI, and machine learning. The integration of data-driven predictive capabilities into cell therapy manufacturing, such as predicting product quality and clinical outcomes based on manufacturing data, or ensuring robustness and reliability using data-driven supply-chain modeling could enable more precise and efficient production processes and lead to better patient access and outcomes. In this review, we introduce some of the relevant computational and data science tools and how they are being or can be implemented in the cell therapy manufacturing workflow. We also identify areas where innovative approaches are required to address challenges and opportunities specific to the cell therapy industry. We conclude that interfacing data science throughout a cell therapy product lifecycle, developing data-driven manufacturing workflow, designing better data collection tools and algorithms, using data analytics and AI-based methods to better understand critical quality attributes and critical-process parameters, and training the appropriate workforce will be critical for overcoming current industry and regulatory barriers and accelerating clinical translation.

Odor dilution sorting as a clinical test of olfactory function: normative values and reliability data.

Clinical assessment of an individual's sense of smell has gained prominence, but its resource-intensive nature necessitates the exploration of self-administered methods. In this study, a cohort of 68 patients with olfactory loss and 55 controls were assessed using a recently introduced olfactory test. This test involves sorting 2 odorants (eugenol and phenylethyl alcohol) in 5 dilutions according to odor intensity, with an average application time of 3.5 min. The sorting task score, calculated as the mean of Kendall's Tau between the assigned and true dilution orders and normalized to [0,1], identified a cutoff for anosmia at a score ≤ 0.7. This cutoff, which marks the 90th percentile of scores obtained with randomly ordered dilutions, had a balanced accuracy of 89% (78% to 97%) for detecting anosmia, comparable to traditional odor threshold assessments. Retest evaluations suggested a score difference of ±0.15 as a cutoff for clinically significant changes in olfactory function. In conclusion, the olfactory sorting test represents a simple, self-administered approach to the detection of anosmia or preserved olfactory function. With balanced accuracy similar to existing brief olfactory tests, this method offers a practical and user-friendly alternative for screening anosmia, addressing the need for resource-efficient assessments in clinical settings.

Diagnosed and subjectively perceived long-term effects of COVID-19 infection on olfactory function assessed by supervised machine learning.

Loss of olfactory function is a typical acute coronavirus disease 2019 (COVID-19) symptom, at least in early variants of SARS-CoV2. The time that has elapsed since the emergence of COVID-19 now allows for assessing the long-term prognosis of its olfactory impact. Participants (n = 722) of whom n = 464 reported having had COVID-19 dating back with a mode of 174 days were approached in a museum as a relatively unbiased environment. Olfactory function was diagnosed by assessing odor threshold and odor identification performance. Subjects also rated their actual olfactory function on an 11-point numerical scale [0,…10]. Neither the frequency of olfactory diagnostic categories nor olfactory test scores showed any COVID-19-related effects. Olfactory diagnostic categories (anosmia, hyposmia, or normosmia) were similarly distributed among former patients and controls (0.86%, 18.97%, and 80.17% for former patients and 1.17%, 17.51%, and 81.32% for controls). Former COVID-19 patients, however, showed differences in their subjective perception of their own olfactory function. The impact of this effect was substantial enough that supervised machine learning algorithms detected past COVID-19 infections in new subjects, based on reduced self-awareness of olfactory performance and parosmia, while the diagnosed olfactory function did not contribute any relevant information in this context. Based on diagnosed olfactory function, results suggest a positive prognosis for COVID-19-related olfactory loss in the long term. Traces of former infection are found in self-perceptions of olfaction, highlighting the importance of investigating the long-term effects of COVID-19 using reliable and validated diagnostic measures in olfactory testing.

Noninvasive Assessment of Diabetic Kidney Disease With MRI: Hype or Hope?

Owing to the increasing prevalence of diabetic mellitus, diabetic kidney disease (DKD) is presently the leading cause of chronic kidney disease and end-stage renal disease worldwide. Early identification and disease interception is of paramount clinical importance for DKD management. However, current diagnostic, disease monitoring and prognostic tools are not satisfactory, due to their low sensitivity, low specificity, or invasiveness. Magnetic resonance imaging (MRI) is noninvasive and offers a host of contrast mechanisms that are sensitive to pathophysiological changes and risk factors associated with DKD. MRI tissue characterization involves structural and functional information including renal morphology (kidney volume (TKV) and parenchyma thickness using T1- or T2-weighted MRI), renal microstructure (diffusion weighted imaging, DWI), renal tissue oxygenation (blood oxygenation level dependent MRI, BOLD), renal hemodynamics (arterial spin labeling and phase contrast MRI), fibrosis (DWI) and abdominal or perirenal fat fraction (Dixon MRI). Recent (pre)clinical studies demonstrated the feasibility and potential value of DKD evaluation with MRI. Recognizing this opportunity, this review outlines key concepts and current trends in renal MRI technology for furthering our understanding of the mechanisms underlying DKD and for supplementing clinical decision-making in DKD. Progress in preclinical MRI of DKD is surveyed, and challenges for clinical translation of renal MRI are discussed. Future directions of DKD assessment and renal tissue characterization with (multi)parametric MRI are explored. Opportunities for discovery and clinical break-through are discussed including biological validation of the MRI findings, large-scale population studies, standardization of DKD protocols, the synergistic connection with data science to advance comprehensive texture analysis, and the development of smart and automatic data analysis and data visualization tools to further the concepts of virtual biopsy and personalized DKD precision medicine. We hope that this review will convey this vision and inspire the reader to become pioneers in noninvasive assessment and management of DKD with MRI. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 2.

Applying data science methodologies with artificial intelligence variant reinterpretation to map and estimate genetic disorder prevalence utilizing clinical data.

Data science methodologies can be utilized to ascertain and analyze clinical genetic data that is often unstructured and rarely used outside of patient encounters. Genetic variants from all genetic testing resulting to a large pediatric healthcare system for a 5-year period were obtained and reinterpreted utilizing the previously validated Franklin© Artificial Intelligence (AI). Using PowerBI©, the data were further matched to patients in the electronic healthcare record to associate with demographic data to generate a variant data table and mapped by ZIP codes. Three thousand and sixty-five variants were identified and 98% were matched to patients with geographic data. Franklin© changed the interpretation for 24% of variants. One hundred and fifty-six clinically actionable variant reinterpretations were made. A total of 739 Mendelian genetic disorders were identified with disorder prevalence estimation. Mapping of variants demonstrated hot-spots for pathogenic genetic variation such as PEX6-associated Zellweger Spectrum Disorder. Seven patients were identified with Bardet-Biedl syndrome and seven patients with Rett syndrome amenable to newly FDA-approved therapeutics. Utilizing readily available software we developed a database and Exploratory Data Analysis (EDA) methodology enabling us to systematically reinterpret variants, estimate variant prevalence, identify conditions amenable to new treatments, and localize geographies enriched for pathogenic variants.