The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

A combination of HLA-DP α and ß chain polymorphisms paired with a SNP in the DPB1 3' UTR region, denoting expression levels, are associated with atopic dermatitis.

Introduction: Components of the immune response have previously been associated with the pathophysiology of atopic dermatitis (AD), specifically the Human Leukocyte Antigen (HLA) Class II region via genome-wide association studies, however the exact elements have not been identified. Methods: This study examines the genetic variation of HLA Class II genes using next generation sequencing (NGS) and evaluates the resultant amino acids, with particular attention on binding site residues, for associations with AD. The Genetics of AD cohort was used to evaluate HLA Class II allelic variation on 464 subjects with AD and 384 controls. Results: Statistically significant associations with HLA-DP α and ß alleles and specific amino acids were found, some conferring susceptibility to AD and others with a protective effect. Evaluation of polymorphic residues in DP binding pockets revealed the critical role of P1 and P6 (P1: α31M + (ß84G or ß84V) [protection]; α31Q + ß84D [susceptibility] and P6: α11A + ß11G [protection]) and were replicated with a national cohort of children consisting of 424 AD subjects. Independently, AD susceptibility-associated residues were associated with the G polymorphism of SNP rs9277534 in the 3' UTR of the HLA-DPB1 gene, denoting higher expression of these HLA-DP alleles, while protection-associated residues were associated with the A polymorphism, denoting lower expression. Discussion: These findings lay the foundation for evaluating non-self-antigens suspected to be associated with AD as they potentially interact with particular HLA Class II subcomponents, forming a complex involved in the pathophysiology of AD. It is possible that a combination of structural HLA-DP components and levels of expression of these components contribute to AD pathophysiology.

Genomic characterization of HLA class I and class II genes in ethnically diverse sub-Saharan African populations: A report on novel HLA alleles.

HLA allelic variation has been well studied and documented in many parts of the world. However, African populations have been relatively under-represented in studies of HLA variation. We have characterized HLA variation from 489 individuals belonging to 13 ethnically diverse populations from rural communities from the African countries of Botswana, Cameroon, Ethiopia, and Tanzania, known to practice traditional subsistence lifestyles using next generation sequencing (Illumina) and long-reads from Oxford Nanopore Technologies. We identified 342 distinct alleles among the 11 HLA targeted genes: HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1, with 140 of those alleles containing novel sequences that were submitted to the IPD-IMGT/HLA database. Sixteen of the 140 alleles contained novel content within the exonic regions of the genes, while 110 alleles contained novel intronic variants. Four alleles were found to be recombinants of already described HLA alleles and 10 alleles extended the sequence content of already described alleles. All 140 alleles include complete allelic sequence from the 5' UTR to the 3' UTR that are inclusive of all exons and introns. This report characterizes the HLA allelic variation from these individuals and describes the novel allelic variation present within these specific African populations.

Ancestry adjustment improves genome-wide estimates of regional intolerance.

Genomic regions subject to purifying selection are more likely to carry disease-causing mutations than regions not under selection. Cross species conservation is often used to identify such regions but with limited resolution to detect selection on short evolutionary timescales such as that occurring in only one species. In contrast, genetic intolerance looks for depletion of variation relative to expectation within a species, allowing species-specific features to be identified. When estimating the intolerance of noncoding sequence, methods strongly leverage variant frequency distributions. As the expected distributions depend on ancestry, if not properly controlled for, ancestral population source may obfuscate signals of selection. We demonstrate that properly incorporating ancestry in intolerance estimation greatly improved variant classification. We provide a genome-wide intolerance map that is conditional on ancestry and likely to be particularly valuable for variant prioritization.

Accurate long-read sequencing allows assembly of the duplicated RHD and RHCE genes harboring variants relevant to blood transfusion.

Next-generation sequencing (NGS) technologies have transformed medical genetics. However, short-read lengths pose a limitation on identification of structural variants, sequencing repetitive regions, phasing of distant nucleotide changes, and distinguishing highly homologous genomic regions. Long-read sequencing technologies may offer improvements in the characterization of genes that are currently difficult to assess. We used a combination of targeted DNA capture, long-read sequencing, and a customized bioinformatics pipeline to fully assemble the RH region, which harbors variation relevant to red cell donor-recipient mismatch, particularly among patients with sickle cell disease. RHD and RHCE are a pair of duplicated genes located within an ∼175 kb region on human chromosome 1 that have high sequence similarity and frequent structural variations. To achieve the assembly, we utilized palindrome repeats in PacBio SMRT reads to obtain consensus sequences of 2.1 to 2.9 kb average length with over 99% accuracy. We used these long consensus sequences to identify 771 assembly markers and to phase the RHD-RHCE region with high confidence. The dataset enabled direct linkage between coding and intronic variants, phasing of distant SNPs to determine RHD-RHCE haplotypes, and identification of known and novel structural variations along with the breakpoints. A limiting factor in phasing is the frequency of heterozygous assembly markers and therefore was most successful in samples from African Black individuals with increased heterogeneity at the RH locus. Overall, this approach allows RH genotyping and de novo assembly in an unbiased and comprehensive manner that is necessary to expand application of NGS technology to high-resolution RH typing.

A Stat1 bound enhancer promotes Nampt expression and function within tumor associated macrophages.

Tumor associated macrophage responses are regulated by distinct metabolic states that affect their function. However, the ability of specific signals in the local tumor microenvironment to program macrophage metabolism remains under investigation. Here, we identify NAMPT, the rate limiting enzyme in NAD salvage synthesis, as a target of STAT1 during cellular activation by interferon gamma, an important driver of macrophage polarization and antitumor responses. We demonstrate that STAT1 occupies a conserved element within the first intron of Nampt, termed Nampt-Regulatory Element-1 (NRE1). Through disruption of NRE1 or pharmacological inhibition, a subset of M1 genes is sensitive to NAMPT activity through its impact on glycolytic processes. scRNAseq is used to profile in vivo responses by NRE1-deficient, tumor-associated leukocytes in melanoma tumors through the creation of a unique mouse strain. Reduced Nampt and inflammatory gene expression are present in specific myeloid and APC populations; moreover, targeted ablation of NRE1 in macrophage lineages results in greater tumor burden. Finally, elevated NAMPT expression correlates with IFNγ responses and melanoma patient survival. This study identifies IFN and STAT1-inducible Nampt as an important factor that shapes the metabolic program and function of tumor associated macrophages.

Spatial constrains and information content of sub-genomic regions of the human genome.

Complexity metrics and machine learning (ML) models have been utilized to analyze the lengths of segmental genomic entities of DNA sequences (exonic, intronic, intergenic, repeat, unique) with the purpose to ask questions regarding the segmental organization of the human genome within the size distribution of these sequences. For this we developed an integrated methodology that is based upon the reconstructed phase space theorem, the non-extensive statistical theory of Tsallis, ML techniques, and a technical index, integrating the generated information, which we introduce and named complexity factor (COFA). Our analysis revealed that the size distribution of the genomic regions within chromosomes are not random but follow patterns with characteristic features that have been seen through its complexity character, and it is part of the dynamics of the whole genome. Finally, this picture of dynamics in DNA is recognized using ML tools for clustering, classification, and prediction with high accuracy.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes.

BACKGROUND: HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk. METHODS: One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed. RESULTS: A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1. CONCLUSIONS: The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

Genomic characterization of MICA gene using multiple next generation sequencing platforms: A validation study.

We have developed a protocol regarding the genomic characterization of the MICA gene by next generation sequencing (NGS). The amplicon includes the full length of the gene and is about 13 kb. A total of 156 samples were included in the study. Ninety-seven of these samples were previously characterized at MICA by legacy methods (Sanger or sequence specific oligonucleotide) and were used to evaluate the accuracy, precision, specificity, and sensitivity of the assay. An additional 59 DNA samples of unknown ethnicity volunteers from the United States were only genotyped by NGS. Samples were chosen to contain a diverse set of alleles. Our NGS approach included a first round of sequencing on the Illumina MiSeq platform and a second round of sequencing on the MinION platform by Oxford Nanopore Technology (ONT), on selected samples for the purpose of either characterizing new alleles or setting phase among multiple polymorphisms to resolve ambiguities or generate complete sequence for alleles that were only partially reported in the IMGT/HLA database. Complete consensus sequences were generated for every allele sequenced with ONT, extending from the 5' untranslated region (UTR) to the 3' UTR of the MICA gene. Thirty-two MICA sequences were submitted to the IMGT/HLA database including either new alleles or filling up the gaps (exonic, intronic and/or UTRs) of already reported alleles. Some of the challenges associated with the characterization of these samples are discussed.

Utilizing nanopore sequencing technology for the rapid and comprehensive characterization of eleven HLA loci; addressing the need for deceased donor expedited HLA typing.

The comprehensive characterization of human leukocyte antigen (HLA) genomic sequences remains a challenging problem. Despite the significant advantages of next-generation sequencing (NGS) in the field of Immunogenetics, there has yet to be a single solution for unambiguous, accurate, simple, cost-effective, and timely genotyping necessary for all clinical applications. This report demonstrates the benefits of nanopore sequencing introduced by Oxford Nanopore Technologies (ONT) for HLA genotyping. Samples (n = 120) previously characterized at high-resolution three-field (HR-3F) for 11 loci were assessed using ONT sequencing paired to a single-plex PCR protocol (Holotype) and to two multiplex protocols OmniType (Omixon) and NGSgo®-MX6-1 (GenDx). The results demonstrate the potential of nanopore sequencing for delivering accurate HR-3F typing with a simple, rapid, and cost-effective protocol. The protocol is applicable to time-sensitive applications, such as deceased donor typings, enabling better assessments of compatibility and epitope analysis. The technology also allows significantly shorter turnaround time for multiple samples at a lower cost. Overall, the nanopore technology appears to offer a significant advancement over current next-generation sequencing platforms as a single solution for all HLA genotyping needs.

AnthOligo: automating the design of oligonucleotides for capture/enrichment technologies.

SUMMARY: A number of methods have been devised to address the need for targeted genomic resequencing. One of these methods, region-specific extraction (RSE) is characterized by the capture of long DNA fragments (15-20 kb) by magnetic beads, after enzymatic extension of oligonucleotides hybridized to selected genomic regions. Facilitating the selection of the most appropriate capture oligos for targeting a region of interest, satisfying the properties of temperature (Tm) and entropy (ΔG), while minimizing the formation of primer-dimers in a pooled experiment, is therefore necessary. Manual design and selection of oligos becomes very challenging, complicated by factors such as length of the target region and number of targeted regions. Here we describe, AnthOligo, a web-based application developed to optimally automate the process of generation of oligo sequences used to target and capture the continuum of large and complex genomic regions. Apart from generating oligos for RSE, this program may have wider applications in the design of customizable internal oligos to be used as baits for gene panel analysis or even probes for large-scale comparative genomic hybridization array processes. AnthOligo was tested by capturing the Major Histocompatibility Complex (MHC) of a random sample.The application provides users with a simple interface to upload an input file in BED format and customize parameters for each task. The task of probe design in AnthOligo commences when a user uploads an input file and concludes with the generation of a result-set containing an optimal set of region-specific oligos. AnthOligo is currently available as a public web application with URL: http://antholigo.chop.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Resolving MiSeq-Generated Ambiguities in HLA-DPB1 Typing by Using the Oxford Nanopore Technology.

The technical limitations of current next-generation sequencing technologies, combined with an ever-increasing number of human leukocyte antigen (HLA) alleles, form the basis for the additional ambiguities encountered at an increasing rate in clinical practice. HLA-DPB1 characterization, particularly, generates a significant percentage of ambiguities (25.5%), posing a challenge for accurate and unambiguous HLA-DPB1 genotyping. Phasing of exonic heterozygous positions between exon 2 and all other downstream exons has been the major cause of ambiguities. In this study, the Oxford Nanopore MinION, a third-generation sequencing technology, was used to resolve the phasing. The accurate MiSeq sequencing data, combined with the long reads obtained from the MinION platform, allow for the resolution of the tested ambiguities.

Next-generation HLA typing of 382 International Histocompatibility Working Group reference B-lymphoblastoid cell lines: Report from the 17th International HLA and Immunogenetics Workshop.

Extended molecular characterization of HLA genes in the IHWG reference B-lymphoblastoid cell lines (B-LCLs) was one of the major goals for the 17th International HLA and Immunogenetics Workshop (IHIW). Although reference B-LCLs have been examined extensively in previous workshops complete high-resolution typing was not completed for all the classical class I and class II HLA genes. To address this, we conducted a single-blind study where select panels of B-LCL genomic DNA samples were distributed to multiple laboratories for HLA genotyping by next-generation sequencing methods. Identical cell panels comprised of 24 and 346 samples were distributed and typed by at least four laboratories in order to derive accurate consensus HLA genotypes. Overall concordance rates calculated at both 2- and 4-field allele-level resolutions ranged from 90.4% to 100%. Concordance for the class I genes ranged from 91.7 to 100%, whereas concordance for class II genes was variable; the lowest observed at HLA-DRB3 (84.2%). At the maximum allele-resolution 78 B-LCLs were defined as homozygous for all 11 loci. We identified 11 novel exon polymorphisms in the entire cell panel. A comparison of the B-LCLs NGS HLA genotypes with the HLA genotypes catalogued in the IPD-IMGT/HLA Database Cell Repository, revealed an overall allele match at 68.4%. Typing discrepancies between the two datasets were mostly due to the lower-resolution historical typing methods resulting in incomplete HLA genotypes for some samples listed in the IPD-IMGT/HLA Database Cell Repository. Our approach of multiple-laboratory NGS HLA typing of the B-LCLs has provided accurate genotyping data. The data generated by the tremendous collaborative efforts of the 17th IHIW participants is useful for updating the current cell and sequence databases and will be a valuable resource for future studies.

Transcriptome-wide profiling of multiple RNA modifications simultaneously at single-base resolution.

The breadth and importance of RNA modifications are growing rapidly as modified ribonucleotides can impact the sequence, structure, function, stability, and fate of RNAs and their interactions with other molecules. Therefore, knowing cellular RNA modifications at single-base resolution could provide important information regarding cell status and fate. A current major limitation is the lack of methods that allow the reproducible profiling of multiple modifications simultaneously, transcriptome-wide and at single-base resolution. Here we developed RBS-Seq, a modification of RNA bisulfite sequencing that enables the sensitive and simultaneous detection of m5C, Ψ, and m1A at single-base resolution transcriptome-wide. With RBS-Seq, m5C and m1A are accurately detected based on known signature base mismatches and are detected here simultaneously along with Ψ sites that show a 1-2 base deletion. Structural analyses revealed the mechanism underlying the deletion signature, which involves Ψ-monobisulfite adduction, heat-induced ribose ring opening, and Mg2+-assisted reorientation, causing base-skipping during cDNA synthesis. Detection of each of these modifications through a unique chemistry allows high-precision mapping of all three modifications within the same RNA molecule, enabling covariation studies. Application of RBS-Seq on HeLa RNA revealed almost all known m5C, m1A, and ψ sites in tRNAs and rRNAs and provided hundreds of new m5C and Ψ sites in noncoding RNAs and mRNAs. However, our results diverge greatly from earlier work, suggesting ∼10-fold fewer m5C sites in noncoding and coding RNAs and the absence of substantial m1A in mRNAs. Taken together, the approaches and refined datasets in this work will greatly enable future epitranscriptome studies.

Anti-inflammatory microRNA-146a protects mice from diet-induced metabolic disease.

Identifying regulatory mechanisms that influence inflammation in metabolic tissues is critical for developing novel metabolic disease treatments. Here, we investigated the role of microRNA-146a (miR-146a) during diet-induced obesity in mice. miR-146a is reduced in obese and type 2 diabetic patients and our results reveal that miR-146a-/- mice fed a high-fat diet (HFD) have exaggerated weight gain, increased adiposity, hepatosteatosis, and dysregulated blood glucose levels compared to wild-type controls. Pro-inflammatory genes and NF-κB activation increase in miR-146a-/- mice, indicating a role for this miRNA in regulating inflammatory pathways. RNA-sequencing of adipose tissue macrophages demonstrated a role for miR-146a in regulating both inflammation and cellular metabolism, including the mTOR pathway, during obesity. Further, we demonstrate that miR-146a regulates inflammation, cellular respiration and glycolysis in macrophages through a mechanism involving its direct target Traf6. Finally, we found that administration of rapamycin, an inhibitor of mTOR, was able to rescue the obesity phenotype in miR-146a-/- mice. Altogether, our study provides evidence that miR-146a represses inflammation and diet-induced obesity and regulates metabolic processes at the cellular and organismal levels, demonstrating how the combination of diet and miRNA genetics influences obesity and diabetic phenotypes.

Assessing the utilization of high-resolution 2-field HLA typing in solid organ transplantation.

HLA typing in solid organ transplantation (SOT) is necessary for determining HLA-matching status between donor-recipient pairs and assessing patients' anti-HLA antibody profiles. Histocompatibility has traditionally been evaluated based on serologically defined HLA antigens. The evolution of HLA typing and antibody identification technologies, however, has revealed many limitations with using serologic equivalents for assessing compatibility in SOT. The significant improvements to HLA typing introduced by next-generation sequencing (NGS) require an assessment of the impact of this technology on SOT. We have assessed the role of high-resolution 2-field HLA typing (HR-2F) in SOT by retrospectively evaluating NGS-typed pre- and post-SOT cases. HR-2F typing was highly instructive or necessary in 41% (156/385) of the cases. Several pre- and posttransplant scenarios were identified as being better served by HR-2F typing. Five different categories are presented with specific case examples. The experience of another center (Temple University Hospital) is also included, whereby 21% of the cases required HR-2F typing by Sanger sequencing, as supported by other legacy methods, to properly address posttransplant anti-HLA antibody issues.

FoxA1 and FoxA2 drive gastric differentiation and suppress squamous identity in NKX2-1-negative lung cancer.

Changes in cancer cell identity can alter malignant potential and therapeutic response. Loss of the pulmonary lineage specifier NKX2-1 augments the growth of KRAS-driven lung adenocarcinoma and causes pulmonary to gastric transdifferentiation. Here, we show that the transcription factors FoxA1 and FoxA2 are required for initiation of mucinous NKX2-1-negative lung adenocarcinomas in the mouse and for activation of their gastric differentiation program. Foxa1/2 deletion severely impairs tumor initiation and causes a proximal shift in cellular identity, yielding tumors expressing markers of the squamocolumnar junction of the gastrointestinal tract. In contrast, we observe downregulation of FoxA1/2 expression in the squamous component of both murine and human lung adenosquamous carcinoma. Using sequential in vivo recombination, we find that FoxA1/2 loss in established KRAS-driven neoplasia originating from SPC-positive alveolar cells induces keratinizing squamous cell carcinomas. Thus, NKX2-1, FoxA1 and FoxA2 coordinately regulate the growth and identity of lung cancer in a context-specific manner.

The Lineage-Defining Transcription Factors SOX2 and NKX2-1 Determine Lung Cancer Cell Fate and Shape the Tumor Immune Microenvironment.

The major types of non-small-cell lung cancer (NSCLC)-squamous cell carcinoma and adenocarcinoma-have distinct immune microenvironments. We developed a genetic model of squamous NSCLC on the basis of overexpression of the transcription factor Sox2, which specifies lung basal cell fate, and loss of the tumor suppressor Lkb1 (SL mice). SL tumors recapitulated gene-expression and immune-infiltrate features of human squamous NSCLC; such features included enrichment of tumor-associated neutrophils (TANs) and decreased expression of NKX2-1, a transcriptional regulator that specifies alveolar cell fate. In Kras-driven adenocarcinomas, mis-expression of Sox2 or loss of Nkx2-1 led to TAN recruitment. TAN recruitment involved SOX2-mediated production of the chemokine CXCL5. Deletion of Nkx2-1 in SL mice (SNL) revealed that NKX2-1 suppresses SOX2-driven squamous tumorigenesis by repressing adeno-to-squamous transdifferentiation. Depletion of TANs in SNL mice reduced squamous tumors, suggesting that TANs foster squamous cell fate. Thus, lineage-defining transcription factors determine the tumor immune microenvironment, which in turn might impact the nature of the tumor.

Characterization of 108 Genomic DNA Reference Materials for 11 Human Leukocyte Antigen Loci: A GeT-RM Collaborative Project.

The highly polymorphic human leukocyte antigen (HLA) genes, located in the human major histocompatibility complex, encode the class I and II antigen-presenting molecules, which are centrally involved in the immune response. HLA typing is used for several clinical applications, such as transplantation, pharmacogenetics, and diagnosis of autoimmune disease. HLA typing is highly complex because of the homology of HLA genes and pseudogenes and the extensive polymorphism in the population. The Centers for Disease Control and Prevention established the Genetic Testing Reference Materials Coordination Program (GeT-RM) in partnership with the genetics community to improve the availability of genomic DNA reference materials necessary for quality assurance of genetic laboratory testing. The GeT-RM together with three clinical laboratories and the Coriell Cell Repositories have characterized genomic DNA obtained from a panel of 108 cell lines for all HLA classic polymorphic loci: HLA-A, B, C, DRB1, DRB3, DRB4, DRB5, DQA1, DQB1, DPA1, and DPB1. The goal was to develop a publicly available and renewable source of well-characterized genomic DNA reference materials to support molecular HLA typing assay development, validation, and verification, quality control, and proficiency testing. These genomic DNA samples are publicly available from the National Institutes of General Medical Science Repository at the Coriell Cell Repositories.

Reversible LSD1 inhibition with HCI-2509 induces the p53 gene expression signature and disrupts the MYCN signature in high-risk neuroblastoma cells.

Lysine-Specific Demethylase 1 (LSD1) over-expression correlates with poorly differentiated neuroblastoma and predicts poor outcome despite multimodal therapy. We have studied the efficacy of reversible and specific LSD1 inhibition with HCI-2509 in neuroblastoma cell lines and particularly the effect of HCI-2509 on the transcriptomic profile in MYCN amplified NGP cells. Cell survival assays show that HCI-2509 is cytotoxic to poorly differentiated neuroblastoma cell lines in low micromole or lower doses. Transcriptional profiling of NGP cells treated with HCI-2509 shows a significant effect on p53, cell cycle, MYCN and hypoxia pathway gene sets. HCI-2509 results in increased histone methyl marks and p53 levels along with cell cycle arrest in the G2/M phase and inhibition of colony formation of NGP cells. Our findings indicate that LSD1 inhibition with HCI-2509 has a multi-target effect in neuroblastoma cell lines, mediated in part via p53. MYCN-amplified neuroblastoma cells have a targeted benefit as HCI-2509 downregulates the MYCN upregulated gene set.