Search

Internal tandem duplication mutations in fms-like tyrosine kinase 3 (FLT3-ITD) are recurrent in acute myeloid leukemia (AML) and increase the risk of relapse. Clinical responses to FLT3 inhibitors (FLT3i) include myeloid differentiation of the FLT3-ITD clone in nearly half of patients through an unknown mechanism. We identified enhancer of zeste homolog 2 (EZH2), a component of polycomb repressive complex 2 (PRC2), as a mediator of this effect using a proteomic-based screen. FLT3i downregulated EZH2 protein expression and PRC2 activity on H3K27me3. FLT3-ITD and loss-of-function mutations in EZH2 are mutually exclusive in human AML. We demonstrated that FLT3i increase myeloid maturation with reduced stem/progenitor cell populations in murine Flt3-ITD AML. Combining EZH1/2 inhibitors with FLT3i increased terminal maturation of leukemic cells and reduced leukemic burden. Our data suggest that reduced EZH2 activity following FLT3 inhibition promotes myeloid differentiation of FLT3-ITD leukemic cells, providing a mechanistic explanation for the clinical observations. These results demonstrate that in addition to its known cell survival and proliferation signaling, FLT3-ITD has a second, previously undefined function to maintain a myeloid stem/progenitor cell state through modulation of PRC2 activity. Our findings support exploring EZH1/2 inhibitors as therapy for FLT3-ITD AML., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community., Competing Interests: C.S.G.’s spouse was an employee of Alex’s Lemonade Stand Foundation, which was a sponsor of this research. J.A.S., C.L.S., C.J.B., S.J.S., and J.N.T. are or were employees of Alex’s Lemonade Stand Foundation, a sponsor of this research. A.J.W. is a member of the Scientific Advisory boards for Alexion and DayOne Biopharmaceuticals., (© 2023 The Author(s).)

Background and Aims: A key histopathological feature of inflammatory bowel disease is damage to the mucosa, including breakdown of the epithelial barrier. Human enteroids and colonoids are a critical bench-to-bedside tool for studying the epithelium in inflammatory bowel disease. The goal of the current study was to define transcriptional differences in healthy versus diseased subjects that are sustained in enteroids and colonoids, including from disease-spared tissue., Methods: Biopsies and matching enteroid or colonoid cultures from pediatric patients with ileal Crohn disease (N = 6) and control subjects (N = 17) were subjected to RNA sequencing followed by bioinformatic and machine learning analyses. Late passage enteroids were exposed to cytokines to assess durable transcriptional differences., Results: We observed substantial overlap of pathways upregulated in Crohn disease in enteroids and ileal biopsies, as well as colonoids and rectal biopsies. KEGG pathways for cytokine-cytokine receptor interaction, chemokine signaling, protein export, and Toll-like receptor signaling were upregulated in both ileal and rectal biopsies, as well as enteroids and colonoids. In vitro cytokine exposure reactivated genes previously increased in biopsies. Machine learning predicted biopsy location (100% accuracy) and donor disease status (83% accuracy). A random forest classifier generated using ileal enteroids identified rectal colonoids from ileal Crohn disease subjects with 80% accuracy., Conclusion: We confirmed transcriptional profiles of Crohn disease biopsies are expressed in enteroids and colonoids. Furthermore, transcriptomic data from disease-spared rectal tissue can identify patients with ileal Crohn disease. Our data support the use of patient enteroids and colonoids as critical translational tools for the study of inflammatory bowel disease., Competing Interests: Conflicts of Interest: The authors disclose no conflicts.

Heterozygous pathogenic variants in DNM1 cause developmental and epileptic encephalopathy (DEE) as a result of a dominant-negative mechanism impeding vesicular fission. Thus far, pathogenic variants in DNM1 have been studied with a canonical transcript that includes the alternatively spliced exon 10b. However, after performing RNA sequencing in 39 pediatric brain samples, we find the primary transcript expressed in the brain includes the downstream exon 10a instead. Using this information, we evaluated genotype-phenotype correlations of variants affecting exon 10a and identified a cohort of eleven previously unreported individuals. Eight individuals harbor a recurrent de novo splice site variant, c.1197-8G>A (GenBank: NM_001288739.1), which affects exon 10a and leads to DEE consistent with the classical DNM1 phenotype. We find this splice site variant leads to disease through an unexpected dominant-negative mechanism. Functional testing reveals an in-frame upstream splice acceptor causing insertion of two amino acids predicted to impair oligomerization-dependent activity. This is supported by neuropathological samples showing accumulation of enlarged synaptic vesicles adherent to the plasma membrane consistent with impaired vesicular fission. Two additional individuals with missense variants affecting exon 10a, p.Arg399Trp and p.Gly401Asp, had a similar DEE phenotype. In contrast, one individual with a missense variant affecting exon 10b, p.Pro405Leu, which is less expressed in the brain, had a correspondingly less severe presentation. Thus, we implicate variants affecting exon 10a as causing the severe DEE typically associated with DNM1-related disorders. We highlight the importance of considering relevant isoforms for disease-causing variants as well as the possibility of splice site variants acting through a dominant-negative mechanism., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Acquired aplastic anemia (AA) is caused by autoreactive T cell-mediated destruction of early hematopoietic cells. Somatic loss of human leukocyte antigen (HLA) class I alleles was identified as a mechanism of immune escape in surviving hematopoietic cells of some patients with AA. However, pathogenicity, structural characteristics, and clinical impact of specific HLA alleles in AA remain poorly understood. Here, we evaluated somatic HLA loss in 505 patients with AA from 2 multi-institutional cohorts. Using a combination of HLA mutation frequencies, peptide-binding structures, and association with AA in an independent cohort of 6,323 patients from the National Marrow Donor Program, we identified 19 AA risk alleles and 12 non-risk alleles and established a potentially novel AA HLA pathogenicity stratification. Our results define pathogenicity for the majority of common HLA-A/B alleles across diverse populations. Our study demonstrates that HLA alleles confer different risks of developing AA, but once AA develops, specific alleles are not associated with response to immunosuppression or transplant outcomes. However, higher pathogenicity alleles, particularly HLA-B*14:02, are associated with higher rates of clonal evolution in adult patients with AA. Our study provides insights into the immune pathogenesis of AA, opening the door to future autoantigen identification and improved understanding of clonal evolution in AA.

Angiosarcoma is a rare, devastating malignancy with few curative options for disseminated disease. We analyzed a recently published genomic data set of 48 angiosarcomas and noticed recurrent amplifications of HOXA -cluster genes in 33% of patients. HOXA genes are master regulators of embryonic vascular development and adult neovascularization, which provides a molecular rationale to suspect that amplified HOXA genes act as oncogenes in angiosarcoma. HOXA amplifications typically affected multiple pro-angiogenic HOXA genes and co-occurred with amplifications of CD36 and KDR, whereas the overall mutation rate in these tumors was relatively low. HOXA amplifications were found most commonly in angiosarcomas located in the breast and were rare in angiosarcomas arising in sun-exposed areas on the head, neck, face and scalp. Our data suggest that HOXA -amplified angiosarcoma is a distinct molecular subgroup. Efforts to develop therapies targeting oncogenic HOX gene expression in AML and other sarcomas may have relevance for HOXA -amplified angiosarcoma.

Bone marrow (BM) niche-derived signals are critical for facilitating engraftment after hematopoietic stem cell (HSC) transplantation (HSCT). HSCT is required for restoration of hematopoiesis in patients with inherited BM failure syndromes (iBMFSs). Shwachman-Diamond syndrome (SDS) is a rare iBMFS associated with mutations in SBDS. Previous studies have demonstrated that SBDS deficiency in osteolineage niche cells causes BM dysfunction that promotes leukemia development. However, it is unknown whether BM niche defects caused by SBDS deficiency also impair efficient engraftment of healthy donor HSC after HSCT, a hypothesis that could explain morbidity noted after clinical HSCT for patients with SDS. Here, we report a mouse model with inducible Sbds deletion in hematopoietic and osteolineage cells. Primary and secondary BM transplantation (BMT) studies demonstrated that SBDS deficiency within BM niches caused poor donor hematopoietic recovery and specifically poor HSC engraftment after myeloablative BMT. We have also identified multiple molecular and cellular defects within niche populations that are driven by SBDS deficiency and are accentuated by or develop specifically after myeloablative conditioning. These abnormalities include altered frequencies of multiple niche cell subsets, including mesenchymal lineage cells, macrophages, and endothelial cells; disruption of growth factor signaling, chemokine pathway activation, and adhesion molecule expression; and p53 pathway activation and signals involved in cell cycle arrest. Taken together, this study demonstrates that SBDS deficiency profoundly impacts recipient hematopoietic niche function in the setting of HSCT, suggesting that novel therapeutic strategies targeting host niches could improve clinical HSCT outcomes for patients with SDS., (© 2022 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

Conotruncal defects with normally related great vessels (CTD-NRGVs) occur in both patients with and without 22q11.2 deletion syndrome (22q11.2DS), but it is unclear to what extent the genetically complex etiologies of these heart defects may overlap across these two groups, potentially involving variation within and/or outside of the 22q11.2 region. To explore this potential overlap, we conducted genome-wide SNP-level, gene-level, and gene set analyses using common variants, separately in each of five cohorts, including two with 22q11.2DS ( N = 1472 total cases) and three without 22q11.2DS ( N = 935 total cases). Results from the SNP-level analyses were combined in meta-analyses, and summary statistics from these analyses were also used in gene and gene set analyses. Across all these analyses, no association was significant after correction for multiple comparisons. However, several SNPs, genes, and gene sets with suggestive evidence of association were identified. For common inherited variants, we did not identify strong evidence for shared genomic mechanisms for CTD-NRGVs across individuals with and without 22q11.2 deletions. Nevertheless, several of our top gene-level and gene set results have been linked to cardiogenesis and may represent candidates for future work.

Meningioma-1 (MN1) overexpression in AML is associated with poor prognosis, and forced expression of MN1 induces leukemia in mice. We sought to determine how MN1 causes AML. We found that overexpression of MN1 can be induced by translocations that result in hijacking of a downstream enhancer. Structure predictions revealed that the entire MN1 coding frame is disordered. We identified the myeloid progenitor-specific BAF complex as the key interaction partner of MN1. MN1 over-stabilizes BAF on enhancer chromatin, a function directly linked to the presence of a long polyQ-stretch within MN1. BAF over-stabilization at binding sites of transcription factors regulating a hematopoietic stem/progenitor program prevents the developmentally appropriate decommissioning of these enhancers and results in impaired myeloid differentiation and leukemia. Beyond AML, our data detail how the overexpression of a polyQ protein, in the absence of any coding sequence mutation, can be sufficient to cause malignant transformation., Competing Interests: Declaration of interests K.M.B. holds a patent on the use of DOT1L inhibitors for MN1 leukemia. K.M.B. has received research funding from Syndax and has previously consulted for Agios. I.H. received research funding from Bristol-Myers Squibb (BMS), Celgene, Merck Sharpe & Dohme (MSD), Astellas, Otsuka, Ono, Kyowa Kirin, Sanofi, Shionogi, Zenyaku, Daiichi Sankyo, Taiho, Takeda, Chugai, Eli Lilly, Nihon Shinyaku, Novartis, Pfizer, Fujimoto, Tanabe-Mitsubishi, Fukuyu Hospital, and Yamada Yohojo and received honoraria from or holds membership on an entity’s board of directors, speaker’s bureau, or its advisory committees for Celgene, Janssen, Takeda, Ono, BMS, Novartis, Daiichi Sankyo, Kyowa Kirin, Eisai, Nihon-Shinyaku, Pfizer, AbbVie, Otsuka, Shionogi, Mundi, CSL, and MSD. All other authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Translocations of Meningioma-1 (MN1) occur in a subset of acute myeloid leukemias (AML) and result in high expression of MN1, either as a full-length protein, or as a fusion protein that includes most of the N-terminus of MN1. High levels of MN1 correlate with poor prognosis. When overexpressed in murine hematopoietic progenitors, MN1 causes an aggressive AML characterized by an aberrant myeloid precursor-like gene expression program that shares features of KMT2A-rearranged (KMT2A-r) leukemia, including high levels of Hoxa and Meis1 gene expression. Compounds that target a critical KMT2A-Menin interaction have proven effective in KMT2A-r leukemia. Here, we demonstrate that Menin (Men1) is also critical for the self-renewal of MN1-driven AML through the maintenance of a distinct gene expression program. Genetic inactivation of Men1 led to a decrease in the number of functional leukemia-initiating cells. Pharmacologic inhibition of the KMT2A-Menin interaction decreased colony-forming activity, induced differentiation programs in MN1-driven murine leukemia and decreased leukemic burden in a human AML xenograft carrying an MN1-ETV6 translocation. Collectively, these results nominate Menin inhibition as a promising therapeutic strategy in MN1-driven leukemia.

Understanding mechanisms of cooperation between oncogenes is critical for the development of novel therapies and rational combinations. Acute myeloid leukemia (AML) cells with KMT2A-fusions and KMT2A partial tandem duplications (KMT2APTD) are known to depend on the histone methyltransferase DOT1L, which methylates histone 3 lysine 79 (H3K79). About 30% of KMT2APTD AMLs carry mutations in IDH1/2 (mIDH1/2). Previous studies showed that 2-hydroxyglutarate produced by mIDH1/2 increases H3K79 methylation, and mIDH1/2 patient samples are sensitive to DOT1L inhibition. Together, these findings suggested that stabilization or increases in H3K79 methylation associated with IDH mutations support the proliferation of leukemias dependent on this mark. However, we found that mIDH1/2 and KMT2A alterations failed to cooperate in an experimental model. Instead, mIDH1/2 and 2-hydroxyglutarate exert toxic effects, specifically on KMT2A-rearranged AML cells (fusions/partial tandem duplications). Mechanistically, we uncover an epigenetic barrier to efficient cooperation; mIDH1/2 expression is associated with high global histone 3 lysine 79 dimethylation (H3K79me2) levels, whereas global H3K79me2 is obligate low in KMT2A-rearranged AML. Increasing H3K79me2 levels, specifically in KMT2A-rearrangement leukemias, resulted in transcriptional downregulation of KMT2A target genes and impaired leukemia cell growth. Our study details a complex genetic and epigenetic interaction of 2 classes of oncogenes, IDH1/2 mutations and KMT2A rearrangements, that is unexpected based on the high percentage of IDH mutations in KMT2APTD AML. KMT2A rearrangements are associated with a trend toward lower response rates to mIDH1/2 inhibitors. The substantial adaptation that has to occur for 2 initially counteracting mutations to be tolerated within the same leukemic cell may provide at least a partial explanation for this observation., (© 2020 by The American Society of Hematology.)

Over 50% of patients with 22q11.2 deletion syndrome (DS) have a conotruncal or related cardiac defect (CTRD). We hypothesized that similar genetic variants, developmental pathways and biological functions, contribute to disease risk for CTRD in patients without a 22q11.2 deletion (ND-CTRD) and with a 22q11.2 deletion (DS-CTRD). To test this hypothesis, we performed rare CNV (rCNV)-based analyses on 630 ND-CTRD cases and 602 DS-CTRD cases with comparable cardiac lesions separately and jointly. First, we detected a collection of heart development related pathways from Gene Ontology and Mammalian Phenotype Ontology analysis. We then constructed gene regulation networks using unique genes collected from the rCNVs found in the ND-CTRD and DS-CTRD cohorts. These gene networks were clustered and their predicted functions were examined. We further investigated expression patterns of those unique genes using publicly available mouse embryo microarray expression data from single-cell embryos to fully developed hearts. By these bioinformatics approaches, we identified a commonly shared gene expression pattern in both the ND-CTRD and DS-CTRD cohorts. Computational analysis of gene functions characterized with this expression pattern revealed a collection of significantly enriched terms related to cardiovascular development. By our combined analysis of rCNVs in the ND-CTRD and DS-CTRD cohorts, a group of statistically significant shared pathways, biological functions, and gene expression patterns were identified that can be tested in future studies for their biological relevance., (© 2019 Wiley Periodicals, Inc.)

Acquired aplastic anemia (aAA) is an acquired deficiency of early hematopoietic cells, characterized by inadequate blood production, and a predisposition to myelodysplastic syndrome (MDS) and leukemia. Although its exact pathogenesis is unknown, aAA is thought to be driven by Human Leukocyte Antigen (HLA)-restricted T cell immunity, with earlier studies favoring HLA class II-mediated pathways. Using whole exome sequencing (WES), we recently identified two aAA patients with somatic mutations in HLA class I genes. We hypothesized that HLA class I mutations are pathognomonic for autoimmunity in aAA, but were previously underappreciated because the Major Histocompatibility Complex (MHC) region is notoriously difficult to analyze by WES. Using a combination of targeted deep sequencing of HLA class I genes and single nucleotide polymorphism array (SNP-A) genotyping we screened 66 aAA patients for somatic HLA class I loss. We found somatic HLA loss in eleven patients (17%), with thirteen loss-of-function mutations in HLA-A *33:03, HLA-A *68:01, HLA-B *14:02 and HLA-B *40:02 alleles. Three patients had more than one mutation targeting the same HLA allele. Interestingly, HLA-B *14:02 and HLA-B *40:02 were significantly overrepresented in aAA patients, compared to ethnicity-matched controls. Patients who inherited the targeted HLA alleles, regardless of HLA mutation status, had a more severe disease course with more frequent clonal complications as assessed by WES, SNP-A, and metaphase cytogenetics, and more frequent secondary MDS. The finding of recurrent HLA class I mutations provides compelling evidence for a predominant HLA class I-driven autoimmunity in aAA, and establishes a novel link between aAA patients' immunogenetics and clonal evolution., Competing Interests: Conflict-of-interest disclosure: The authors declare no competing financial interests.

Purpose: Despite substantial progress in sequencing, current strategies can genetically solve only approximately 55-60% of inherited retinal degeneration (IRD) cases. This can be partially attributed to elusive mutations in the known IRD genes, which are not easily identified by the targeted next-generation sequencing (NGS) or Sanger sequencing approaches. We hypothesized that copy-number variations (CNVs) are a major contributor to the elusive genetic causality of IRDs., Methods: Twenty-eight cases previously unsolved with a targeted NGS were investigated with whole-genome single-nucleotide polymorphism (SNP) and comparative genomic hybridization (CGH) arrays., Results: Deletions in the IRD genes were detected in 5 of 28 families, including a de novo deletion. We suggest that the de novo deletion occurred through nonallelic homologous recombination (NAHR) and we constructed a genomic map of NAHR-prone regions with overlapping IRD genes. In this article, we also report an unusual case of recessive retinitis pigmentosa due to compound heterozygous mutations in SNRNP200, a gene that is typically associated with the dominant form of this disease., Conclusions: CNV mapping substantially increased the genetic diagnostic rate of IRDs, detecting genetic causality in 18% of previously unsolved cases. Extending the search to other structural variations will probably demonstrate an even higher contribution to genetic causality of IRDs.Genet Med advance online publication 13 October 2016.

Background: Previous studies using different cardiac phenotypes, technologies and designs suggest a burden of large, rare or de novo copy number variants (CNVs) in subjects with congenital heart defects. We sought to identify disease-related CNVs, candidate genes, and functional pathways in a large number of cases with conotruncal and related defects that carried no known genetic syndrome., Methods: Cases and control samples were divided into two cohorts and genotyped to assess each subject's CNV content. Analyses were performed to ascertain differences in overall CNV prevalence and to identify enrichment of specific genes and functional pathways in conotruncal cases relative to healthy controls., Results: Only findings present in both cohorts are presented. From 973 total conotruncal cases, a burden of rare CNVs was detected in both cohorts. Candidate genes from rare CNVs found in both cohorts were identified based on their association with cardiac development or disease, and/or their reported disruption in published studies. Functional and pathway analyses revealed significant enrichment of terms involved in either heart or early embryonic development., Conclusion: Our study tested one of the largest cohorts specifically with cardiac conotruncal and related defects. These results confirm and extend previous findings that CNVs contribute to disease risk for congenital heart defects in general and conotruncal defects in particular. As disease heterogeneity renders identification of single recurrent genes or loci difficult, functional pathway and gene regulation network analyses appear to be more informative. Birth Defects Research 109:271-295, 2017. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)

Hepatitis-associated aplastic anemia (HAA) is a variant of acquired aplastic anemia (AA) in which immune-mediated bone marrow failure (BMF) develops following an acute episode of seronegative hepatitis. Dyskeratosis congenita (DC) is an inherited BMF syndrome characterized by the presence of short telomeres, mucocutaneous abnormalities, and cancer predisposition. While both conditions may cause BMF and hepatic impairment, therapeutic approaches are distinct, making it imperative to establish the correct diagnosis. In clinical practice, lymphocyte telomere lengths (TL) are used as a first-line screen to rule out inherited telomeropathies before initiating treatment for AA. To evaluate the reliability of TL in the HAA population, we performed a retrospective analysis of TL in 10 consecutively enrolled HAA patients compared to 19 patients with idiopathic AA (IAA). HAA patients had significantly shorter telomeres than IAA patients (P = 0.009), including four patients with TL at or below the 1st percentile for age-matched controls. HAA patients had no clinical features of DC and did not carry disease-causing mutations in known genes associated with inherited telomere disorders. Instead, short TLs were significantly correlated with severe lymphopenia and skewed lymphocyte subsets, features characteristic of HAA. Our results indicate the importance of caution in the interpretation of TL measurements in HAA, because, in this patient population, short telomeres have limited specificity., (© 2015 The Authors. American Journal of Hematology Published by Wiley Periodicals, Inc.)

Acquired aplastic anemia (aAA) is a nonmalignant disease caused by autoimmune destruction of early hematopoietic cells. Clonal hematopoiesis is a late complication, seen in 20-25% of older patients. We hypothesized that clonal hematopoiesis in aAA is a more general phenomenon, which can arise early in disease, even in younger patients. To evaluate clonal hematopoiesis in aAA, we used comparative whole exome sequencing of paired bone marrow and skin samples in 22 patients. We found somatic mutations in 16 patients (72.7%) with a median disease duration of 1 year; of these, 12 (66.7%) were patients with pediatric-onset aAA. Fifty-eight mutations in 51 unique genes were found primarily in pathways of immunity and transcriptional regulation. Most frequently mutated was PIGA, with seven mutations. Only two mutations were in genes recurrently mutated in myelodysplastic syndrome. Two patients had oligoclonal loss of the HLA alleles, linking immune escape to clone emergence. Two patients had activating mutations in key signaling pathways (STAT5B (p.N642H) and CAMK2G (p.T306M)). Our results suggest that clonal hematopoiesis in aAA is common, with two mechanisms emerging-immune escape and increased proliferation. Our findings expand conceptual understanding of this nonneoplastic blood disorder. Future prospective studies of clonal hematopoiesis in aAA will be critical for understanding outcomes and for designing personalized treatment strategies., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Background: We sought to characterize the landscape of structural variation associated with the subset of congenital cardiac defects characterized by left-sided obstruction., Methods: Cases with left-sided cardiac defects (LSCD) and pediatric controls were uniformly genotyped and assessed for copy number variant (CNV) calls. Significance testing was performed to ascertain differences in overall CNV incidence, and for CNV enrichment of specific genes and gene functions in LSCD cases relative to controls., Results: A total of 257 cases of European descent and 962 ethnically matched, disease-free pediatric controls were included. Although there was no difference in CNV rate between cases and controls, a significant enrichment in rare LSCD CNVs was detected overall (p=7.30 × 10(-3) , case/control ratio=1.26) and when restricted either to deletions (p=7.58 × 10(-3) , case/control ratio=1.20) or duplications (3.02 × 10(-3) , case/control ratio=1.43). Neither gene-based, functional nor knowledge-based analyses identified genes, loci or pathways that were significantly enriched in cases as compared to controls when appropriate corrections for multiple tests were applied. However, several genes of interest were identified by virtue of their association with cardiac development, known human conditions, or reported disruption by CNVs in other patient cohorts., Conclusion: This study examines the largest cohort to date with LSCD for structural variation. These data suggest that CNVs play a role in disease risk and identify numerous genes disrupted by CNVs of potential disease relevance. These findings further highlight the genetic heterogeneity and complexity of these disorders., (© 2014 Wiley Periodicals, Inc.)

Setting: Left-sided cardiac lesions have a birth prevalence of approximately 1 in 1000 and have been shown to be heritable in pedigree studies. A large microdeletion at 16p12.1 is associated with childhood developmental delay, and initial studies describing this deletion identified left-sided lesions as an enriched phenotype compared with a control population., Objective: The aim of this study is to determine whether patients with left-sided cardiac lesions have an increased frequency of 16p12.1 microdeletions as compared with control populations., Design: A cohort of 262 probands with left-sided lesions, including 53 with isolated aortic stenosis/bicuspid aortic valve, 83 with coarctation of the aorta with or without aortic stenosis/bicuspid aortic valve, and 126 with hypoplastic left heart syndrome were assessed for copy number variation at 16p12.1. The control cohort included 595 patients with conotruncal defects as a cardiac control and 971 healthy children., Results: We detected one patient in the left-sided lesion cohort with a large duplication partially overlapping the reported 16p12.1 microdeletion, along with one patient each in the conotruncal and control cohorts with a deletion in the same region. None of these patients had dysmorphic features, extracardiac malformations, or developmental delay., Conclusion: In our cohort, structural variation at 16p12.1 was not identified with increased frequency in patients with left-sided lesions as compared with controls., (© 2013 Wiley Periodicals, Inc.)

The bone marrow failure syndromes (BMFS) are a heterogeneous group of rare blood disorders characterized by inadequate haematopoiesis, clonal evolution, and increased risk of leukaemia. Single nucleotide polymorphism arrays (SNP-A) have been proposed as a tool for surveillance of clonal evolution in BMFS. To better understand the natural history of BMFS and to assess the clinical utility of SNP-A in these disorders, we analysed 124 SNP-A from a comprehensively characterized cohort of 91 patients at our BMFS centre. SNP-A were correlated with medical histories, haematopathology, cytogenetic and molecular data. To assess clonal evolution, longitudinal analysis of SNP-A was performed in 25 patients. We found that acquired copy number-neutral loss of heterozygosity (CN-LOH) was significantly more frequent in acquired aplastic anaemia (aAA) than in other BMFS (odds ratio 12·2, P < 0·01). Homozygosity by descent was most common in congenital BMFS, frequently unmasking autosomal recessive mutations. Copy number variants (CNVs) were frequently polymorphic, and we identified CNVs enriched in neutropenia and aAA. Our results suggest that acquired CN-LOH is a general phenomenon in aAA that is probably mechanistically and prognostically distinct from typical CN-LOH of myeloid malignancies. Our analysis of clinical utility of SNP-A shows the highest yield of detecting new clonal haematopoiesis at diagnosis and at relapse., (© 2013 John Wiley & Sons Ltd.)

Background: Mitochondrial DNA (mtDNA) diseases are rare disorders whose prevalence is estimated around 1 in 5000. Patients are usually tested only for deletions and for common mutations of mtDNA which account for 5-40% of cases, depending on the study. However, the prevalence of rare mtDNA mutations is not known., Methods: We analysed the whole mtDNA in a cohort of 743 patients suspected of manifesting a mitochondrial disease, after excluding deletions and common mutations. Both heteroplasmic and homoplasmic variants were identified using two complementary strategies (Surveyor and MitoChip). Multiple correspondence analyses followed by hierarchical ascendant cluster process were used to explore relationships between clinical spectrum, age at onset and localisation of mutations., Results: 7.4% of deleterious mutations and 22.4% of novel putative mutations were identified. Pathogenic heteroplasmic mutations were more frequent than homoplasmic mutations (4.6% vs 2.8%). Patients carrying deleterious mutations showed symptoms before 16 years of age in 67% of cases. Early onset disease (<1 year) was significantly associated with mutations in protein coding genes (mainly in complex I) while late onset disorders (>16 years) were associated with mutations in tRNA genes. MTND5 and MTND6 genes were identified as 'hotspots' of mutations, with Leigh syndrome accounting for the large majority of associated phenotypes., Conclusions: Rare mitochondrial DNA mutations probably account for more than 7.4% of patients with respiratory chain deficiency. This study shows that a comprehensive analysis of mtDNA is essential, and should include young children, for an accurate diagnosis that is now accessible with the development of next generation sequencing technology.

Background: High-throughput sequencing (HTS) technologies are spearheading the accelerated development of biomedical research. Processing and summarizing the large amount of data generated by HTS presents a non-trivial challenge to bioinformatics. A commonly adopted standard is to store sequencing reads aligned to a reference genome in SAM (Sequence Alignment/Map) or BAM (Binary Alignment/Map) files. Quality control of SAM/BAM files is a critical checkpoint before downstream analysis. The goal of the current project is to facilitate and standardize this process., Results: We developed bamchop, a robust program to efficiently summarize key statistical metrics of HTS data stored in BAM files, and to visually present the results in a formatted report. The report documents information about various aspects of HTS data, such as sequencing quality, mapping to a reference genome, sequencing coverage, and base frequency. Bamchop uses the R language and Bioconductor packages to calculate statistical matrices and the Sweave utility and associated LaTeX markup for documentation. Bamchop's efficiency and robustness were tested on BAM files generated by local sequencing facilities and the 1000 Genomes Project. Source code, instruction and example reports of bamchop are freely available from https://github.com/CBMi-BiG/bamchop., Conclusions: Bamchop enables biomedical researchers to quickly and rigorously evaluate HTS data by providing a convenient synopsis and user-friendly reports.

Background: Mitochondrial genome sequence analysis is critical to the diagnostic evaluation of mitochondrial disease. Existing methodologies differ widely in throughput, complexity, cost efficiency, and sensitivity of heteroplasmy detection. Affymetrix MitoChip v2.0, which uses a sequencing-by-genotyping technology, allows potentially accurate and high-throughput sequencing of the entire human mitochondrial genome to be completed in a cost-effective fashion. However, the relatively low call rate achieved using existing software tools has limited the wide adoption of this platform for either clinical or research applications. Here, we report the design and development of a custom bioinformatics software pipeline that achieves a much improved call rate and accuracy for the Affymetrix MitoChip v2.0 platform. We used this custom pipeline to analyze MitoChip v2.0 data from 24 DNA samples representing a broad range of tissue types (18 whole blood, 3 skeletal muscle, 3 cell lines), mutations (a 5.8 kilobase pair deletion and 6 known heteroplasmic mutations), and haplogroup origins. All results were compared to those obtained by at least one other mitochondrial DNA sequence analysis method, including Sanger sequencing, denaturing HPLC-based heteroduplex analysis, and/or the Illumina Genome Analyzer II next generation sequencing platform., Results: An average call rate of 99.75% was achieved across all samples with our custom pipeline. Comparison of calls for 15 samples characterized previously by Sanger sequencing revealed a total of 29 discordant calls, which translates to an estimated 0.012% for the base call error rate. We successfully identified 4 known heteroplasmic mutations and 24 other potential heteroplasmic mutations across 20 samples that passed quality control., Conclusions: Affymetrix MitoChip v2.0 analysis using our optimized MitoChip Filtering Protocol (MFP) bioinformatics pipeline now offers the high sensitivity and accuracy needed for reliable, high-throughput and cost-efficient whole mitochondrial genome sequencing. This approach provides a viable alternative of potential utility for both clinical diagnostic and research applications to traditional Sanger and other emerging sequencing technologies for whole mitochondrial genome analysis.

Background: Recent studies have shown that copy number variations (CNVs) are frequent in higher eukaryotes and associated with a substantial portion of inherited and acquired risk for various human diseases. The increasing availability of high-resolution genome surveillance platforms provides opportunity for rapidly assessing research and clinical samples for CNV content, as well as for determining the potential pathogenicity of identified variants. However, few informatics tools for accurate and efficient CNV detection and assessment currently exist., Results: We developed a suite of software tools and resources (CNV Workshop) for automated, genome-wide CNV detection from a variety of SNP array platforms. CNV Workshop includes three major components: detection, annotation, and presentation of structural variants from genome array data. CNV detection utilizes a robust and genotype-specific extension of the Circular Binary Segmentation algorithm, and the use of additional detection algorithms is supported. Predicted CNVs are captured in a MySQL database that supports cohort-based projects and incorporates a secure user authentication layer and user/admin roles. To assist with determination of pathogenicity, detected CNVs are also annotated automatically for gene content, known disease loci, and gene-based literature references. Results are easily queried, sorted, filtered, and visualized via a web-based presentation layer that includes a GBrowse-based graphical representation of CNV content and relevant public data, integration with the UCSC Genome Browser, and tabular displays of genomic attributes for each CNV., Conclusions: To our knowledge, CNV Workshop represents the first cohesive and convenient platform for detection, annotation, and assessment of the biological and clinical significance of structural variants. CNV Workshop has been successfully utilized for assessment of genomic variation in healthy individuals and disease cohorts and is an ideal platform for coordinating multiple associated projects., Availability and Implementation: Available on the web at: http://sourceforge.net/projects/cnv.