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.

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.

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.

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.

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

Background: 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 after vaccination across populations and the social dimensions of prevention and treatment strategies. Methods: In the United States, 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 coronavirus disease 2019. The network comprises multidisciplinary researchers bridging gaps and fostering collaborations among immunologists, epidemiologists, virologists, clinicians and clinical laboratories, social and behavioral scientists, policymakers, 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, human immunodeficiency virus, transplant recipients, as well as otherwise healthy pregnant women, children, college students, and high-risk occupational workers (including healthcare workers and first responders). Results: 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. Conclusions: 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.

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.

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.

Unique Regulatory Mechanisms for the Human Embryonic Stem Cell Cycle.

The cell cycle in pluripotent human embryonic stem cells is governed by unique mechanisms that support unrestricted proliferation and competency for endodermal, mesodermal, and ectodermal differentiation. The abbreviated G1 period with retention of uncompromised fidelity for genetic and epigenetic mechanisms operative in control of proliferation support competency for expansion of the pluripotent cell population that is fundamental for initial stages of development. Regulatory events during the G1 period of the pluripotent cell cycle are decisive for the transition from pluripotency to lineage commitment. Recent findings indicate that a G2 cell cycle pause is present in both endodermal and mesodermal lineage cells, and is obligatory for differentiation to endoderm. J. Cell. Physiol. 232: 1254-1257, 2017. © 2016 Wiley Periodicals, Inc.

Precocious Phenotypic Transcription-Factor Expression During Early Development.

A novel role for phenotypic transcription factors in very early differentiation was recently observed and merits further study to elucidate what role this precocious expression may have in development. The RUNX1 transcription factor exhibits selective and transient upregulation during early mesenchymal differentiation. In contrast to phenotype-associated transcriptional control of gene expression to establish and sustain hematopoietic/myeloid lineage identity, precocious expression of RUNX1 is functionally linked to control of an epithelial to mesenchymal transition that is obligatory for development. This early RUNX1 expression spike provides a paradigm for precocious expression of a phenotypic transcription factor that invites detailed mechanistic study to fully understand its biological importance. J. Cell. Biochem. 118: 953-958, 2017. © 2016 Wiley Periodicals, Inc.

Transient RUNX1 Expression during Early Mesendodermal Differentiation of hESCs Promotes Epithelial to Mesenchymal Transition through TGFB2 Signaling.

The transition of human embryonic stem cells (hESCs) from pluripotency to lineage commitment is not fully understood, and a role for phenotypic transcription factors in the initial stages of hESC differentiation remains to be explored. From a screen of candidate factors, we found that RUNX1 is selectively and transiently upregulated early in hESC differentiation to mesendodermal lineages. Transcriptome profiling and functional analyses upon RUNX1 depletion established a role for RUNX1 in promoting cell motility. In parallel, we discovered a loss of repression for several epithelial genes, indicating that loss of RUNX1 impaired an epithelial to mesenchymal transition during differentiation. Cell biological and biochemical approaches revealed that RUNX1 depletion specifically compromised TGFB2 signaling. Both the decrease in motility and deregulated epithelial marker expression upon RUNX1 depletion were rescued by reintroduction of TGFB2, but not TGFB1. These findings identify roles for RUNX1-TGFB2 signaling in early events of mesendodermal lineage commitment.

MicroRNA-378-mediated suppression of Runx1 alleviates the aggressive phenotype of triple-negative MDA-MB-231 human breast cancer cells.

The Runx1 transcription factor, known for its essential role in normal hematopoiesis, was reported in limited studies to be mutated or associated with human breast tumor tissues. Runx1 increases concomitantly with disease progression in the MMTV-PyMT transgenic mouse model of breast cancer. Compelling questions relate to mechanisms that regulate Runx1 expression in breast cancer. Here, we tested the hypothesis that dysregulation of Runx1-targeting microRNAs (miRNAs) allows for pathologic increase of Runx1 during breast cancer progression. Microarray profiling of the MMTV-PyMT model revealed significant downregulation of numerous miRNAs predicted to target Runx1. One of these, miR-378, was inversely correlated with Runx1 expression during breast cancer progression in mice and in human breast cancer cell lines MCF7 and triple-negative MDA-MB-231 that represent early- and late-stage diseases, respectively. MiR-378 is nearly absent in MDA-MB-231 cells. Luciferase reporter assays revealed that miR-378 binds the Runx1 3' untranslated region (3'UTR) and inhibits Runx1 expression. Functionally, we demonstrated that ectopic expression of miR-378 in MDA-MB-231 cells inhibited Runx1 and suppressed migration and invasion, while inhibition of miR-378 in MCF7 cells increased Runx1 levels and cell migration. Depletion of Runx1 in late-stage breast cancer cells resulted in increased expression of both the miR-378 host gene PPARGC1B and pre-miR-378, suggesting a feedback loop. Taken together, our study identifies a novel and clinically relevant mechanism for regulation of Runx1 in breast cancer that is mediated by a PPARGC1B-miR-378-Runx1 regulatory pathway. Our results highlight the translational potential of miRNA replacement therapy for inhibiting Runx1 in breast cancer.

Genome-Wide Studies Reveal that H3K4me3 Modification in Bivalent Genes Is Dynamically Regulated during the Pluripotent Cell Cycle and Stabilized upon Differentiation.

Stem cell phenotypes are reflected by posttranslational histone modifications, and this chromatin-related memory must be mitotically inherited to maintain cell identity through proliferative expansion. In human embryonic stem cells (hESCs), bivalent genes with both activating (H3K4me3) and repressive (H3K27me3) histone modifications are essential to sustain pluripotency. Yet, the molecular mechanisms by which this epigenetic landscape is transferred to progeny cells remain to be established. By mapping genomic enrichment of H3K4me3/H3K27me3 in pure populations of hESCs in G2, mitotic, and G1 phases of the cell cycle, we found striking variations in the levels of H3K4me3 through the G2-M-G1 transition. Analysis of a representative set of bivalent genes revealed that chromatin modifiers involved in H3K4 methylation/demethylation are recruited to bivalent gene promoters in a cell cycle-dependent fashion. Interestingly, bivalent genes enriched with H3K4me3 exclusively during mitosis undergo the strongest upregulation after induction of differentiation. Furthermore, the histone modification signature of genes that remain bivalent in differentiated cells resolves into a cell cycle-independent pattern after lineage commitment. These results establish a new dimension of chromatin regulation important in the maintenance of pluripotency.

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.