Interpreting single-cell messages in normal and aberrant hematopoiesis with the Cell Marker Accordion

Single-cell technologies offer a unique opportunity to explore cellular heterogeneity in hematopoiesis, reveal malignant hematopoietic cells with clinically significant features and measure gene signatures linked to pathological pathways. However, reliable identification of cell types is a crucial bottleneck in single-cell analysis. Available databases contain dissimilar nomenclature and non-concurrent marker sets, leading to inconsistent annotations and poor interpretability. Furthermore, current tools focus mostly on physiological cell types, lacking extensive applicability in disease. We developed the Cell Marker Accordion, a user-friendly platform for the automatic annotation and biological interpretation of single-cell populations based on consistency weighted markers. We validated our approach on peripheral blood and bone marrow single-cell datasets, using surface markers and expert-based annotation as the ground truth. In all cases, we significantly improved the accuracy in identifying cell types with respect to any single source database. Moreover, the Cell Marker Accordion can identify disease-critical cells and pathological processes, extracting potential biomarkers in a wide variety of contexts in human and murine single-cell datasets. It characterizes leukemia stem cell subtypes, including therapy-resistant cells in acute myeloid leukemia patients; it identifies malignant plasma cells in multiple myeloma samples; it dissects cell type alterations in splicing factor-mutant cells from myelodysplastic syndrome patients; it discovers activation of innate immunity pathways in bone marrow from mice treated with METTL3 inhibitors. The breadth of these applications elevates the Cell Marker Accordion as a flexible, faithful and standardized tool to annotate and interpret hematopoietic populations in single-cell datasets focused on the study of hematopoietic development and disease.

Impact of Memory T Cells on SARS-COV-2 Vaccine Response in Hematopoietic Stem Cell Transplant

During the COVID-19 pandemic, hematopoietic stem cell transplant (HSCT) recipients faced an elevated mortality rate from SARS-CoV-2 infection, ranging between 10-40%. The SARS-CoV-2 mRNA vaccines are important tools in preventing severe disease, yet their efficacy in the post-transplant setting remains unclear, especially in patients subjected to myeloablative chemotherapy and immunosuppression. We evaluated the humoral and adaptive immune responses to the SARS-CoV-2 mRNA vaccination series in 42 HSCT recipients and 5 healthy controls. Peripheral blood mononuclear nuclear cells and serum were prospectively collected before and after each dose of the SARS-CoV-2 vaccine. Post-vaccination responses were assessed by measuring anti-spike IgG and nucleocapsid titers, and antigen specific T cell activity, before and after vaccination. In order to examine mechanisms behind a lack of response, pre-and post-vaccine samples were selected based on humoral and cellular responses for single-cell RNA sequencing with TCR and BCR sequencing. Our observations revealed that while all participants eventually mounted a humoral response, transplant recipients had defects in memory T cell populations that were associated with an absence of T cell response, some of which could be detected pre-vaccination.

Microfluidic Immuno-Serolomic Assay Reveals Systems Level Association with COVID-19 Pathology and Vaccine Protection

How to develop highly informative serology assays to evaluate the quality of immune protection against coronavirus disease-19 (COVID-19) has been a global pursuit over the past years. Here, a microfluidic high-plex immuno-serolomic assay is developed to simultaneously measure50 plasma or serum samples for50 soluble markers including 35proteins, 11 anti-spike/receptor binding domian (RBD) IgG antibodies spanningmajor variants, and controls. This assay demonstrates the quintuplicate test in a single run with high throughput, low sample volume, high reproducibilityand accuracy. It is applied to the measurement of 1012 blood samples including in-depth analysis of sera from 127 patients and 21 healthy donors over multiple time points, either with acute COVID infection or vaccination. The protein analysis reveals distinct immune mediator modules that exhibit a reduced degree of diversity in protein-protein cooperation in patients with hematologic malignancies or receiving B cell depletion therapy. Serological analysis identifies that COVID-infected patients with hematologic malignancies display impaired anti-RBD antibody response despite high level of anti-spike IgG, which can be associated with limited clonotype diversity and functional deficiency in B cells. These findings underscore the importance to individualize immunization strategies for these high-risk patients and provide an informative tool to monitor their responses at the systems level.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Abstract 798: SeroNet Pooling Project of immunocompromised populations

Introduction: COVID-19 vaccination substantially reduces morbidity and mortality associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and severe illness. However, despite effective COVID-19 vaccines many questions remain about the efficacy of vaccines and the durability and robustness of immune responses, especially in immunocompromised persons. The NCI-funded Serological Sciences Network (SeroNet) is a coordinated effort including 11 sites to advance research on the immune response to SARS-CoV-2 infection and COVID-19 vaccination among diverse and vulnerable populations. The goals of the Pooling Project are: (1) to conduct real-world data (RWD) analyses using electronic medical records (EMR) data from four health care systems (Kaiser Permanente Northern California, Northwell Health, Veterans Affairs-Case Western, and Cedars-Sinai) to determine vaccine effectiveness in (a) cancer patients; (b) autoimmune diseases and (c) solid organ transplant recipients (SOTR); (2) to conduct meta-analyses of prospective cohort studies from eight SeroNet institutions (Cedars-Sinai, Johns Hopkins, Northwell Health, Emory University, University of Minnesota, Mount Sinai, Yale University) to determine post-vaccine immune responses in (a) lung cancer patients; (b) hematologic cancers/hematopoietic stem cell transplant (HSCT) recipients; (c) SOTR; (d) lupus.

Methods: For our RWD analyses, data is extracted from EMR using standardized algorithms using ICD-10 codes to identify immunocompromised persons (hematologic and solid organ malignancy; SOTR; autoimmune disease, including inflammatory bowel disease, rheumatoid arthritis, and SLE). We use common case definitions to extract data on demographic, laboratory values, clinical co-morbidity, COVID-19 vaccination, SARS-CoV-2 infection and severe COVID-19, and disease-specific variables. In addition, we pool individual-level data from prospective cohorts enrolling patients with cancer and other immunosuppressed conditions from across network. Surveys and biospecimens from serology and immune profiling are collected at pre-specified timepoints across longitudinal cohorts.

Results: Currently, we have EMR data extracted from 4 health systems including >715,000 cancer patients, >9,500 SOTR and >180,000 with autoimmune conditions. Prospective cohorts across the network have longitudinal data on >450 patients with lung cancer, >1,200 patients with hematologic malignancies, >400 SOTR and >400 patients with lupus. We will report results examining vaccine effectiveness for prevention of SARS-CoV-2 infection, severe COVID-19 and post-acute sequelae of COVID-19 (PAS-C or long COVID) in cancer patients compared to other immunocompromised conditions.

Conclusion: Our goal is to inform public health guidelines on COVID-19 vaccine and boosters to reduce SARS-CoV-2 infection and severe illness in immunocompromised populations.

Citation Format: Elham Kazemain, Jane Figueiredo, Jacek Skarbinski, Russell McBride, Viviana Simon, Amy B. Karger, F. Eun-Hyung Lee, Fred R. Hirsch, Andrea Cox, Sabra Klein, Rong Fan, Stephanie Halene, David A. Zidar, James M. Crawford, Bharat Thyagarajan, Charles Gleason, Alex Mathson, Komal Srivastava, Puleng Moshele, Toby Amoss, Martin Runnstrom, Susanne Linderman, Ananda M. Rodilla, Philip C. Mack, Yu Shyr, Anna Yin, Patrick Shea, Jennifer VanOudenhove, Hinnah Siddiqui, Brigid M. Wilson, Eric P. Elkin, Crystal A. Hsiao, Yonah Ziemba, Cheryl B. Schleicher, Sharon Fox, Lawrence H. Kushi, Karen Reckamp, Akil Merchant, Noah Merin. SeroNet Pooling Project of immunocompromised populations [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 798.

Is it the time to integrate novel sequencing technologies into clinical practice?

PURPOSE OF REVIEW

The aim of this study was to provide insight into how novel next-generation sequencing (NGS) techniques are set to revolutionize clinical practice.

RECENT FINDINGS

Advances in sequencing technologies have focused on improved capture of mutations and reads and cellular resolution. Both short and long read DNA sequencing technology are being refined and combined in novel ways with other multiomic approaches to gain unprecedented biological insight into disease. Single-cell (sc)DNA-seq and integrated scDNA-seq with immunophenotyping provide granular information on disease composition such as clonal hierarchy, co-mutation status, zygosity, clonal diversity and genotype phenotype correlations. These and other techniques can identify rare cell populations providing the opportunity for increased sensitivity in measurable residual disease monitoring and precise characterization of residual clones permitting distinction of leukemic from pre/nonmalignant clones.

SUMMARY

Increasing genetics-based mechanistic insights and classification of myeloid diseases along with a decrease in the cost of high-throughput NGS mean novel sequencing technologies are closer to being a reality in standard clinical practice. These technologies are poised to improve diagnostics, our ability to monitor treatment response and minimal residual disease and allow the study of premalignant conditions such as clonal haematopoiesis.

Microfluidic immuno-serology assay revealed a limited diversity of protection against COVID-19 in patients with altered immunity

The immune response to SARS-CoV-2 for patients with altered immunity such as hematologic malignancies and autoimmune disease may differ substantially from that in general population. These patients remain at high risk despite wide-spread adoption of vaccination. It is critical to examine the differences at the systems level between the general population and the patients with altered immunity in terms of immunologic and serological responses to COVID-19 infection and vaccination. Here, we developed a novel microfluidic chip for high-plex immuno-serological assay to simultaneously measure up to 50 plasma or serum samples for up to 50 soluble markers including 35 plasma proteins, 11 anti-spike/RBD IgG antibodies spanning all major variants, and controls. Our assay demonstrated the quintuplicate test in a single run with high throughput, low sample volume input, high reproducibility and high accuracy. It was applied to the measurement of 1,012 blood samples including in-depth analysis of sera from 127 patients and 21 healthy donors over multiple time points, either with acute COVID infection or vaccination. The protein association matrix analysis revealed distinct immune mediator protein modules that exhibited a reduced degree of diversity in protein-protein cooperation in patients with hematologic malignancies and patients with autoimmune disorders receiving B cell depletion therapy. Serological analysis identified that COVID infected patients with hematologic malignancies display impaired anti-RBD antibody response despite high level of anti-spike IgG, which could be associated with limited clonotype diversity and functional deficiency in B cells and was further confirmed by single-cell BCR and transcriptome sequencing. These findings underscore the importance to individualize immunization strategy for these high-risk patients and provide an informative tool to monitor their responses at the systems level.

Zebrafish models of alx-linked frontonasal dysplasia reveal a role for Alx1 and Alx3 in the anterior segment and vasculature of the developing eye

The cellular and genetic mechanisms that coordinate formation of facial sensory structures with surrounding skeletal and soft tissue elements remain poorly understood. Alx1, a homeobox transcription factor, is a key regulator of midfacial morphogenesis. ALX1 mutations in humans are linked to severe congenital anomalies of the facial skeleton (frontonasal dysplasia, FND) with malformation or absence of eyes and orbital contents (micro- and anophthalmia). Zebrafish with loss-of-function alx1 mutations develop with craniofacial and ocular defects of variable penetrance, likely due to compensatory upregulation in expression of a paralogous gene, alx3. Here we show that zebrafish alx1;alx3 mutants develop with highly penetrant cranial and ocular defects that resemble human ALX1-linked FND. alx1 and alx3 are expressed in anterior cranial neural crest (aCNC), which gives rise to the anterior neurocranium (ANC), anterior segment structures of the eye and vascular pericytes. Consistent with a functional requirement for alx genes in aCNC, alx1; alx3 mutants develop with nearly absent ANC and grossly aberrant hyaloid vasculature and ocular anterior segment, but normal retina. In vivo lineage labeling identified a requirement for alx1 and alx3 during aCNC migration, and transcriptomic analysis suggested oxidative stress response as a key target mechanism of this function. Oxidative stress is a hallmark of fetal alcohol toxicity, and we found increased penetrance of facial and ocular malformations in alx1 mutants exposed to ethanol, consistent with a protective role for alx1 against ethanol toxicity. Collectively, these data demonstrate a conserved role for zebrafish alx genes in controlling ocular and facial development, and a novel role in protecting these key midfacial structures from ethanol toxicity during embryogenesis. These data also reveal novel roles for alx genes in ocular anterior segment formation and vascular development and suggest that retinal deficits in alx mutants may be secondary to aberrant ocular vascularization and anterior segment defects. This study establishes robust zebrafish models for interrogating conserved genetic mechanisms that coordinate facial and ocular development, and for exploring gene--environment interactions relevant to fetal alcohol syndrome.

Mission, Organization and Future Direction of the Serological Sciences Network for COVID-19 (SeroNet) Epidemiologic Cohort Studies

Global efforts are needed to elucidate the epidemiology of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the underlying cause of coronavirus disease 2019 (COVID-19) including seroprevalence, risk factors and long-term sequelae, as well as immune responses following vaccination across populations and the social dimensions of prevention and treatment strategies. In the U.S., the National Cancer Institute in partnership with the National Institute of Allergy and Infectious Diseases, established the SARS-CoV-2 Serological Sciences Network (SeroNet) as the nation’s largest coordinated effort to study COVID-19. The network is comprised of multidisciplinary researchers bridging gaps and fostering collaborations between immunologists, epidemiologists, virologists, clinicians and clinical laboratories, social and behavioral scientists, policy makers, data scientists, and community members. In total, 49 institutions form the SeroNet consortium to study individuals with cancer, autoimmune disease, inflammatory bowel diseases, cardiovascular diseases, HIV, transplant recipients, as well as otherwise healthy pregnant women, children, college students, and high-risk occupational workers (including health care workers and first responders). Several studies focus on underrepresented populations, including ethnic minorities and rural communities. To support integrative data analyses across SeroNet studies, efforts are underway to define common data elements for standardized serology measurements, cellular and molecular assays, self-reported data, treatment, and clinical outcomes. In this paper, we discuss the overarching framework for SeroNet epidemiology studies, critical research questions under investigation, and data accessibility for the worldwide scientific community. Lessons learned will help inform preparedness and responsiveness to future emerging diseases.

Sarcomere function activates a p53-dependent DNA damage response that promotes polyploidization and limits in vivo cell engraftment

Human cardiac regeneration is limited by low cardiomyocyte replicative rates and progressive polyploidization by unclear mechanisms. To study this process, we engineer a human cardiomyocyte model to track replication and polyploidization using fluorescently tagged cyclin B1 and cardiac troponin T. Using time-lapse imaging, in vitro cardiomyocyte replication patterns recapitulate the progressive mononuclear polyploidization and replicative arrest observed in vivo. Single-cell transcriptomics and chromatin state analyses reveal that polyploidization is preceded by sarcomere assembly, enhanced oxidative metabolism, a DNA damage response, and p53 activation. CRISPR knockout screening reveals p53 as a driver of cell-cycle arrest and polyploidization. Inhibiting sarcomere function, or scavenging ROS, inhibits cell-cycle arrest and polyploidization. Finally, we show that cardiomyocyte engraftment in infarcted rat hearts is enhanced 4-fold by the increased proliferation of troponin-knockout cardiomyocytes. Thus, the sarcomere inhibits cell division through a DNA damage response that can be targeted to improve cardiomyocyte replacement strategies.

Epigenomic and Transcriptomic Dynamics During Human Heart Organogenesis

RATIONALE

There is growing evidence that common variants and rare sequence alterations in regulatory sequences can result in birth defects or predisposition to disease. Congenital heart defects are the most common birth defect and have a clear genetic component, yet only a third of cases can be attributed to structural variation in the genome or a mutation in a gene. The remaining unknown cases could be caused by alterations in regulatory sequences.

OBJECTIVE

Identify regulatory sequences and gene expression networks that are active during organogenesis of the human heart. Determine whether these sites and networks are enriched for disease-relevant genes and associated genetic variation.

METHODS AND RESULTS

We characterized ChromHMM (chromatin state) and gene expression dynamics during human heart organogenesis. We profiled 7 histone modifications in embryonic hearts from each of 9 distinct Carnegie stages (13-14, 16-21, and 23), annotated chromatin states, and compared these maps to over 100 human tissues and cell types. We also generated RNA-sequencing data, performed differential expression, and constructed weighted gene coexpression networks. We identified 177 412 heart enhancers; 12 395 had not been previously annotated as strong enhancers. We identified 92% of all functionally validated heart-positive enhancers (n=281; 7.5× enrichment; P<2.2×10-16). Integration of these data demonstrated novel heart enhancers are enriched near genes expressed more strongly in cardiac tissue and are enriched for variants associated with ECG measures and atrial fibrillation. Our gene expression network analysis identified gene modules strongly enriched for heart-related functions, regulatory control by heart-specific enhancers, and putative disease genes.

CONCLUSIONS

Well-connected hub genes with heart-specific expression targeted by embryonic heart-specific enhancers are likely disease candidates. Our functional annotations will allow for better interpretation of whole genome sequencing data in the large number of patients affected by congenital heart defects.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

PSA regulates androgen receptor expression in prostate cancer cells

BACKGROUND

Prostate-specific antigen (PSA) is a pivotal downstream target gene of the androgen receptor (AR), and a serum biomarker to monitor prostate cancer (PrCa) progression. It has been reported that PSA transactivates AR, but the mechanistic requirements of this response have not been investigated.

METHODS

We studied the localization of PSA, AR, and Src in intracellular compartments of synthetic androgen (R1881)-stimulated LNCaP and C4-2B PrCa cells, using immunofluorescence and subcellular fractionation approaches. We also investigated the effect of downregulation of PSA on AR expression by immunoblotting and real-time PCR using short hairpin RNA (shRNA) and small interfering RNA (siRNA). Src activity was analyzed by immunoblotting.

RESULTS

R1881 stimulation induced nuclear localization of both PSA and AR in LNCaP and C4-2B PrCa cells as well as increased phosphorylation of Src. Stable shRNA or transient siRNA knockdown of PSA resulted in reduced AR protein levels as well as AR mRNA levels in C4-2B cells. Similar to C4-2B cells, ablation of AR levels upon silencing of PSA was also confirmed in VCaP cells, another androgen-independent cell line. Silencing of PSA did not cause significant changes in Src activation; besides, Src regulation by integrins did not appear to affect AR transcriptional activity.

CONCLUSIONS

PSA localizes to nuclei of androgen-stimulated PrCa cells, and controls AR mRNA and protein levels. This regulatory loop is specific for PSA, does not involve known AR activators such as Src and AKT, and may contribute to AR signaling under conditions of increasing PSA levels in patients.

Analysis of PKR structure by small-angle scattering

Protein kinase R (PKR) is a key component of the interferon antiviral defense pathway. Upon binding double-stranded RNA, PKR undergoes autophosphorylation reactions that activate the kinase. PKR contains an N-terminal double-stranded RNA binding domain, which consists of two tandem double-stranded RNA binding motifs, and a C-terminal kinase domain. We have used small-angle X-ray scattering and small-angle neutron scattering to define the conformation of latent PKR in solution. Guinier analysis indicates a radius of gyration of about 35 A. The p(r) distance distribution function exhibits a peak near 30 A, with a broad shoulder extending to longer distances. Good fits to the scattering data require models that incorporate multiple compact and extended conformations of the two interdomain linker regions. Thus, PKR belongs to the growing family of proteins that contain intrinsically unstructured regions. We propose that the flexible linkers may allow PKR to productively dimerize upon interaction with RNA activators that have diverse structures.