DNA Methylation signatures underpinning blood neutrophil to lymphocyte ratio during first week of human life.

Understanding of newborn immune ontogeny in the first week of life will enable age-appropriate strategies for safeguarding vulnerable newborns against infectious diseases. Here we conducted an observational study exploring the immunological profile of infants longitudinally throughout their first week of life. Our Expanded Program on Immunization - Human Immunology Project Consortium (EPIC-HIPC) studies the epigenetic regulation of systemic immunity using small volumes of peripheral blood samples collected from West African neonates on days of life (DOL) 0, 1, 3, and 7. Genome-wide DNA methylation and single nucleotide polymorphism markers are examined alongside matched transcriptomic and flow cytometric data. Integrative analysis reveals that a core network of transcription factors mediates dynamic shifts in neutrophil-to-lymphocyte ratios (NLR), which are underpinned by cell-type specific methylation patterns in the two cell types. Genetic variants are associated with lower NLRs at birth, and healthy newborns with lower NLRs at birth are more likely to subsequently develop sepsis. These findings provide valuable insights into the early-life determinants of immune system development.

Breastfeeding and Neonatal Age Influence Neutrophil-Driven Ontogeny of Blood Cell Populations in the First Week of Human Life.

The first few days of life are characterized by rapid external and internal changes that require substantial immune system adaptations. Despite growing evidence of the impact of this period on lifelong immune health, this period remains largely uncharted. To identify factors that may impact the trajectory of immune development, we conducted stringently standardized, high-throughput phenotyping of peripheral white blood cell (WBC) populations from 796 newborns across two distinct cohorts (The Gambia, West Africa; Papua New Guinea, Melanesia) in the framework of a Human Immunology Project Consortium (HIPC) study. Samples were collected twice from each newborn during the first week of life, first at Day of Life 0 (at birth) and then subsequently at Day of Life 1, 3, or 7 depending on the randomization group the newborn belongs to. The subsequent analysis was conducted at an unprecedented level of detail using flow cytometry and an unbiased automated gating algorithm. The results showed that WBC composition in peripheral blood changes along patterns highly conserved across populations and environments. Changes across days of life were most pronounced in the innate myeloid compartment. Breastfeeding, and at a smaller scale neonatal vaccination, were associated with changes in peripheral blood neutrophil and monocyte cell counts. Our results suggest a common trajectory of immune development in newborns and possible association with timing of breastfeeding initiation, which may contribute to immune-mediated protection from infection in early life. These data begin to outline a specific window of opportunity for interventions that could deliberately direct WBC composition, and with that, immune trajectory and thus ontogeny in early life. This trial is registered with NCT03246230.

Recommendations for using artificial intelligence in clinical flow cytometry.

Flow cytometry is a key clinical tool in the diagnosis of many hematologic malignancies and traditionally requires close inspection of digital data by hematopathologists with expert domain knowledge. Advances in artificial intelligence (AI) are transferable to flow cytometry and have the potential to improve efficiency and prioritization of cases, reduce errors, and highlight fundamental, previously unrecognized associations with underlying biological processes. As a multidisciplinary group of stakeholders, we review a range of critical considerations for appropriately applying AI to clinical flow cytometry, including use case identification, low and high risk use cases, validation, revalidation, computational considerations, and the present regulatory frameworks surrounding AI in clinical medicine. In particular, we provide practical guidance for the development, implementation, and suggestions for potential regulation of AI-based methods in the clinical flow cytometry laboratory. We expect these recommendations to be a helpful initial framework of reference, which will also require additional updates as the field matures.

MAGIC-DR: An interpretable machine-learning guided approach for acute myeloid leukemia measurable residual disease analysis.

Multiparameter flow cytometry is widely used for acute myeloid leukemia minimal residual disease testing (AML MRD) but is time consuming and demands substantial expertise. Machine learning offers potential advancements in accuracy and efficiency, but has yet to be widely adopted for this application. To explore this, we trained single cell XGBoost classifiers from 98 diagnostic AML cell populations and 30 MRD negative samples. Performance was assessed by cross-validation. Predictions were integrated with UMAP as a heatmap parameter for an augmented/interactive AML MRD analysis framework, which was benchmarked against traditional MRD analysis for 25 test cases. The results showed that XGBoost achieved a median AUC of 0.97, effectively distinguishing diverse AML cell populations from normal cells. When integrated with UMAP, the classifiers highlighted MRD populations against the background of normal events. Our pipeline, MAGIC-DR, incorporated classifier predictions and UMAP into flow cytometry standard (FCS) files. This enabled a human-in-the-loop machine learning guided MRD workflow. Validation against conventional analysis for 25 MRD samples showed 100% concordance in myeloid blast detection, with MAGIC-DR also identifying several immature monocytic populations not readily found by conventional analysis. In conclusion, Integrating a supervised classifier with unsupervised dimension reduction offers a robust method for AML MRD analysis that can be seamlessly integrated into conventional workflows. Our approach can support and augment human analysis by highlighting abnormal populations that can be gated on for quantification and further assessment. This has the potential to speed up MRD analysis, and potentially improve detection sensitivity for certain AML immunophenotypes.

flowSim: Near duplicate detection for flow cytometry data.

The analysis of large amounts of data is important for the development of machine learning (ML) models. flowSim is the first algorithm designed to visualize, detect and remove highly redundant information in flow cytometry (FCM) training sets to decrease the computational time for training and increase the performance of ML algorithms by reducing overfitting. flowSim performs near duplicate image detection by combining community detection algorithms with the density analysis of the marker expression values. flowSim clustering compared to consensus manual clustering on a dataset composed of 160 images of bivariate FCM data had a mean Adjusted Rand Index of 0.90, demonstrating its efficiency in identifying similar patterns. flowSim selectively discarded near duplicate files in datasets constructed with known redundancy, and removed 92.6% of FCM images in a dataset of over 500,000 drawn from public repositories.

flowCut: An R package for automated removal of outlier events and flagging of files based on time versus fluorescence analysis.

Technical artifacts such as clogging that occur during the data acquisition process of flow cytometry data can cause spurious events and fluorescence intensity shifting that impact the quality of the data and its analysis results. These events should be identified and potentially removed before being passed to the next stage of analysis. flowCut, an R package, automatically detects anomaly events in flow cytometry experiments and flags files for potential review. Its results are on par with manual analysis and it outperforms existing automated approaches.

Single-cell profiling reveals a memory B cell-like subtype of follicular lymphoma with increased transformation risk.

Follicular lymphoma (FL) is an indolent cancer of mature B-cells but with ongoing risk of transformation to more aggressive histology over time. Recurrent mutations associated with transformation have been identified; however, prognostic features that can be discerned at diagnosis could be clinically useful. We present here comprehensive profiling of both tumor and immune compartments in 155 diagnostic FL biopsies at single-cell resolution by mass cytometry. This revealed a diversity of phenotypes but included two recurrent patterns, one which closely resembles germinal center B-cells (GCB) and another which appears more related to memory B-cells (MB). GCB-type tumors are enriched for EZH2, TNFRSF14, and MEF2B mutations, while MB-type tumors contain increased follicular helper T-cells. MB-type and intratumoral phenotypic diversity are independently associated with increased risk of transformation, supporting biological relevance of these features. Notably, a reduced 26-marker panel retains sufficient information to allow phenotypic profiling of future cohorts by conventional flow cytometry.

Automated identification of maximal differential cell populations in flow cytometry data.

We introduce a new cell population score called SpecEnr (specific enrichment) and describe a method that discovers robust and accurate candidate biomarkers from flow cytometry data. Our approach identifies a new class of candidate biomarkers we define as driver cell populations, whose abundance is associated with a sample class (e.g., disease), but not as a result of a change in a related population. We show that the driver cell populations we find are also easily interpretable using a lattice-based visualization tool. Our method is implemented in the R package flowGraph, freely available on GitHub (github.com/aya49/flowGraph) and on BioConductor.

ISAC Probe Tag Dictionary: Standardized Nomenclature for Detection and Visualization Labels Used in Cytometry and Microscopy Imaging.

Since the advent of microscopy imaging and flow cytometry, there has been an explosion in the number of probes, consisting of a component binding to an analyte and a detectable tag, to mark areas of interest in or on cells and tissue. Probe tags have been created to detect and/or visualize probes. Over time, these probe tags have increased in number. The expansion has resulted in arbitrarily created synonyms of probe tags used in publications and software. The synonyms are problematic for readability of publications, accuracy of text/data mining, and bridging data from multiple platforms, protocols, and databases for Big Data analysis. Development and implementation of a universal language for probe tags will ensure equivalent quality and level of data being reported or extracted for clinical/scientific evaluation as well as help connect data from many platforms. The International Society for Advancement of Cytometry Data Standards Task Force composed of academic scientists and industry hardware/software/reagent manufactures have developed recommendations for a standardized nomenclature for probe tags used in cytometry and microscopy imaging. These recommendations are shared in this technical note in the form of a Probe Tag Dictionary. © 2020 International Society for Advancement of Cytometry.

Data File Standard for Flow Cytometry, Version FCS 3.2.

FCS 3.2 is a revision of the flow cytometry data standard based on a decade of suggested improvements from the community as well as industry needs to capture instrument conditions and measurement features more precisely. The unchanged goal of the standard is to provide a uniform file format that allows files created by one type of acquisition hardware and software to be analyzed by any other type. The standard retains the overall FCS file structure and most features of previous versions, but also contains a few changes that were required to support new types of data and use cases efficiently. These changes are incompatible with existing FCS file readers. Notably, FCS 3.2 supports mixed data types to, for example, allow FCS measurements that are intrinsically integers (e.g., indices or class assignments) or measurements that are commonly captured as integers (e.g., time ticks) to be more represented as integer values, while capturing other measurements as floating-point values in the same FCS data set. In addition, keywords explicitly specifying dyes, detectors, and analytes were added to avoid having to extract those heuristically and unreliably from measurement names. Types of measurements were formalized, several keywords added, others removed, or deprecated, and various aspects of the specification were clarified. A reference implementation of the cyclic redundancy check (CRC) calculation is provided in two programming languages since a correct CRC implementation was problematic for many vendors. © 2020 International Society for Advancement of Cytometry.

Corrigendum: Clinical Protocol for a Longitudinal Cohort Study Employing Systems Biology to Identify Markers of Vaccine Immunogenicity in Newborn Infants in The Gambia and Papua New Guinea.

[This corrects the article DOI: 10.3389/fped.2020.00197.].

Multi-Omic Data Integration Allows Baseline Immune Signatures to Predict Hepatitis B Vaccine Response in a Small Cohort.

Background: Vaccination remains one of the most effective means of reducing the burden of infectious diseases globally. Improving our understanding of the molecular basis for effective vaccine response is of paramount importance if we are to ensure the success of future vaccine development efforts. Methods: We applied cutting edge multi-omics approaches to extensively characterize temporal molecular responses following vaccination with hepatitis B virus (HBV) vaccine. Data were integrated across cellular, epigenomic, transcriptomic, proteomic, and fecal microbiome profiles, and correlated to final HBV antibody titres. Results: Using both an unsupervised molecular-interaction network integration method (NetworkAnalyst) and a data-driven integration approach (DIABLO), we uncovered baseline molecular patterns and pathways associated with more effective vaccine responses to HBV. Biological associations were unravelled, with signalling pathways such as JAK-STAT and interleukin signalling, Toll-like receptor cascades, interferon signalling, and Th17 cell differentiation emerging as important pre-vaccination modulators of response. Conclusion: This study provides further evidence that baseline cellular and molecular characteristics of an individual's immune system influence vaccine responses, and highlights the utility of integrating information across many parallel molecular datasets.

"Age Related Differences in the Biology of Chronic Graft-Versus-Host Disease After Hematopoietic Stem Cell Transplantation".

It is established that pediatric hematopoietic stem cell transplant (HSCT) recipients have a lower rate of chronic graft-versus-host disease (cGvHD) compared to adults. Our group has previously published immune profiles changes associated with cGvHD of clinically well-defined adult and pediatric HSCT cohorts. Since all analyses were performed by the same research group and analyzed using identical methodology, we first compared our previous immune profile analyses between adults and children. We then performed additional analyses comparing the T cell populations across age groups, and a sub-analysis of the impact of the estimated pubertal status at time of HSCT in our pediatric cohort. In all analyses, we corrected for clinical covariates including total body irradiation and time of onset of cGvHD. Three consistent findings were seen in both children and adults, including elevations of ST2 and naive helper T (Th) cells and depression of NKreg cells. However, significant differences exist between children and adults in certain cytokines, B cell, and Treg populations. In children, we saw a broad suppression of newly formed B (NF-B) cells, whereas adults exhibited an increase in T1-CD21lo B cells and a decrease in T1-CD24hiCD38hi B cells. Prepubertal children had elevations of aminopeptidase N (sCD13) and ICAM-1. Treg abnormalities in children appeared to be primarily in memory Treg cells, whereas in adults the abnormalities were in naïve Treg cells. In adults, the loss of PD1 expression in naïve Treg and naïve Th cells was associated with cGvHD. We discuss the possible mechanisms for these age-related differences, and how they might theoretically impact on different therapeutic approaches to cGvHD between children and adults.

Systems Biology Methods Applied to Blood and Tissue for a Comprehensive Analysis of Immune Response to Hepatitis B Vaccine in Adults.

Conventional vaccine design has been based on trial-and-error approaches, which have been generally successful. However, there have been some major failures in vaccine development and we still do not have highly effective licensed vaccines for tuberculosis, HIV, respiratory syncytial virus, and other major infections of global significance. Approaches at rational vaccine design have been limited by our understanding of the immune response to vaccination at the molecular level. Tools now exist to undertake in-depth analysis using systems biology approaches, but to be fully realized, studies are required in humans with intensive blood and tissue sampling. Methods that support this intensive sampling need to be developed and validated as feasible. To this end, we describe here a detailed approach that was applied in a study of 15 healthy adults, who were immunized with hepatitis B vaccine. Sampling included ~350 mL of blood, 12 microbiome samples, and lymph node fine needle aspirates obtained over a ~7-month period, enabling comprehensive analysis of the immune response at the molecular level, including single cell and tissue sample analysis. Samples were collected for analysis of immune phenotyping, whole blood and single cell gene expression, proteomics, lipidomics, epigenetics, whole blood response to key immune stimuli, cytokine responses, in vitro T cell responses, antibody repertoire analysis and the microbiome. Data integration was undertaken using different approaches-NetworkAnalyst and DIABLO. Our results demonstrate that such intensive sampling studies are feasible in healthy adults, and data integration tools exist to analyze the vast amount of data generated from a multi-omics systems biology approach. This will provide the basis for a better understanding of vaccine-induced immunity and accelerate future rational vaccine design.

Clinical Protocol for a Longitudinal Cohort Study Employing Systems Biology to Identify Markers of Vaccine Immunogenicity in Newborn Infants in The Gambia and Papua New Guinea.

Background: Infection contributes to significant morbidity and mortality particularly in the very young and in low- and middle-income countries. While vaccines are a highly cost-effective tool against infectious disease little is known regarding the cellular and molecular pathways by which vaccines induce protection at an early age. Immunity is distinct in early life and greater precision is required in our understanding of mechanisms of early life protection to inform development of new pediatric vaccines. Methods and Analysis: We will apply transcriptomic, proteomic, metabolomic, multiplex cytokine/chemokine, adenosine deaminase, and flow cytometry immune cell phenotyping to delineate early cellular and molecular signatures that correspond to vaccine immunogenicity. This approach will be applied to a neonatal cohort in The Gambia (N ~ 720) receiving at birth: (1) Hepatitis B (HepB) vaccine alone, (2) Bacille Calmette Guerin (BCG) vaccine alone, or (3) HepB and BCG vaccines, (4) HepB and BCG vaccines delayed till day 10 at the latest. Each study participant will have a baseline peripheral blood sample drawn at DOL0 and a second blood sample at DOL1,-3, or-7 as well as late timepoints to assess HepB vaccine immunogenicity. Blood will be fractionated via a "small sample big data" standard operating procedure that enables multiple downstream systems biology assays. We will apply both univariate and multivariate frameworks and multi-OMIC data integration to identify features associated with anti-Hepatitis B (anti-HB) titer, an established correlate of protection. Cord blood sample collection from a subset of participants will enable human in vitro modeling to test mechanistic hypotheses identified in silico regarding vaccine action. Maternal anti-HB titer and the infant microbiome will also be correlated with our findings which will be validated in a smaller cohort in Papua New Guinea (N ~ 80). Ethics and Dissemination: The study has been approved by The Gambia Government/MRCG Joint Ethics Committee and The Boston Children's Hospital Institutional Review Board. Ethics review is ongoing with the Papua New Guinea Medical Research Advisory Committee. All de-identified data will be uploaded to public repositories following submission of study output for publication. Feedback meetings will be organized to disseminate output to the study communities. Clinical Trial Registration: Clinicaltrials.gov Registration Number: NCT03246230.

The miR-185/PAK6 axis predicts therapy response and regulates survival of drug-resistant leukemic stem cells in CML.

Overcoming drug resistance and targeting cancer stem cells remain challenges for curative cancer treatment. To investigate the role of microRNAs (miRNAs) in regulating drug resistance and leukemic stem cell (LSC) fate, we performed global transcriptome profiling in treatment-naive chronic myeloid leukemia (CML) stem/progenitor cells and identified that miR-185 levels anticipate their response to ABL tyrosine kinase inhibitors (TKIs). miR-185 functions as a tumor suppressor: its restored expression impaired survival of drug-resistant cells, sensitized them to TKIs in vitro, and markedly eliminated long-term repopulating LSCs and infiltrating blast cells, conferring a survival advantage in preclinical xenotransplantation models. Integrative analysis with mRNA profiles uncovered PAK6 as a crucial target of miR-185, and pharmacological inhibition of PAK6 perturbed the RAS/MAPK pathway and mitochondrial activity, sensitizing therapy-resistant cells to TKIs. Thus, miR-185 presents as a potential predictive biomarker, and dual targeting of miR-185-mediated PAK6 activity and BCR-ABL1 may provide a valuable strategy for overcoming drug resistance in patients.

Single Cell Phenotypic Profiling of 27 DLBCL Cases Reveals Marked Intertumoral and Intratumoral Heterogeneity.

Diffuse large B-cell lymphoma (DLBCL) is the most common histologic subtype of non-Hodgkin lymphoma and is notorious for its clinical heterogeneity. Patient outcomes can be predicted by cell-of-origin (COO) classification, demonstrating that the underlying transcriptional signature of malignant B-cells informs biological behavior in the context of standard combination chemotherapy regimens. In the current study, we used mass cytometry (CyTOF) to examine tumor phenotypes at the protein level with single cell resolution in a collection of 27 diagnostic DLBCL biopsy specimens from treatment naïve patients. We found that malignant B-cells from each patient occupied unique regions in 37-dimensional phenotypic space with no apparent clustering of samples into discrete subtypes. Interestingly, variable MHC class II expression was found to be the greatest contributor to phenotypic diversity. Within individual tumors, a subset of cases showed multiple phenotypic subpopulations, and in one case, we were able to demonstrate direct correspondence between protein-level phenotypic subsets and DNA mutation-defined subclones. In summary, CyTOF analysis can resolve both intertumoral and intratumoral heterogeneity among primary samples and reveals that each case of DLBCL is unique and may be comprised of multiple, genetically distinct subclones. © 2019 International Society for Advancement of Cytometry.

High-throughput phenotyping reveals expansive genetic and structural underpinnings of immune variation.

By developing a high-density murine immunophenotyping platform compatible with high-throughput genetic screening, we have established profound contributions of genetics and structure to immune variation (http://www.immunophenotype.org). Specifically, high-throughput phenotyping of 530 unique mouse gene knockouts identified 140 monogenic 'hits', of which most had no previous immunologic association. Furthermore, hits were collectively enriched in genes for which humans show poor tolerance to loss of function. The immunophenotyping platform also exposed dense correlation networks linking immune parameters with each other and with specific physiologic traits. Such linkages limit freedom of movement for individual immune parameters, thereby imposing genetically regulated 'immunologic structures', the integrity of which was associated with immunocompetence. Hence, we provide an expanded genetic resource and structural perspective for understanding and monitoring immune variation in health and disease.