Integrated digital pathology at scale: A solution for clinical diagnostics and cancer research at a large academic medical center.

OBJECTIVE: Broad adoption of digital pathology (DP) is still lacking, and examples for DP connecting diagnostic, research, and educational use cases are missing. We blueprint a holistic DP solution at a large academic medical center ubiquitously integrated into clinical workflows; researchapplications including molecular, genetic, and tissue databases; and educational processes. MATERIALS AND METHODS: We built a vendor-agnostic, integrated viewer for reviewing, annotating, sharing, and quality assurance of digital slides in a clinical or research context. It is the first homegrown viewer cleared by New York State provisional approval in 2020 for primary diagnosis and remote sign-out during the COVID-19 (coronavirus disease 2019) pandemic. We further introduce an interconnected Honest Broker for BioInformatics Technology (HoBBIT) to systematically compile and share large-scale DP research datasets including anonymized images, redacted pathology reports, and clinical data of patients with consent. RESULTS: The solution has been operationally used over 3 years by 926 pathologists and researchers evaluating 288 903 digital slides. A total of 51% of these were reviewed within 1 month after scanning. Seamless integration of the viewer into 4 hospital systems clearly increases the adoption of DP. HoBBIT directly impacts the translation of knowledge in pathology into effective new health measures, including artificial intelligence-driven detection models for prostate cancer, basal cell carcinoma, and breast cancer metastases, developed and validated on thousands of cases. CONCLUSIONS: We highlight major challenges and lessons learned when going digital to provide orientation for other pathologists. Building interconnected solutions will not only increase adoption of DP, but also facilitate next-generation computational pathology at scale for enhanced cancer research.

Clonal hematopoiesis in patients with dyskeratosis congenita.

Dyskeratosis congenita (DC) is a rare inherited telomeropathy most frequently caused by mutations in a number of genes all thought to be involved in telomere maintenance. The main causes of mortality in DC are bone marrow failure as well as malignancies including leukemias and solid tumors. The clinical picture including the degree of bone marrow failure is highly variable and factors that contribute to this variability are poorly understood. Based on the recent finding of frequent clonal hematopoiesis in related bone marrow failure syndromes, we hypothesized that somatic mutations may also occur in DC and may contribute at least in part to the variability in blood production. To evaluate for the presence of clonal hematopoiesis in DC, we used a combination of X-inactivation, comparative whole exome sequencing (WES) and single nucleotide polymorphism array (SNP-A) analyses. We found that clonal hematopoiesis in DC is common, as suggested by skewed X-inactivation in 8 out of 9 female patients compared to 3 out of 10 controls, and by the finding of acquired copy neutral loss-of-heterozygosity on SNP-A analysis. In addition, 3 out of 6 independent DC patients were found to have acquired somatic changes in their bone marrow by WES, including a somatic reversion in DKC1, as well as missense mutations in other protein coding genes. Our results indicate that clonal hematopoiesis is a common feature of DC, and suggest that such somatic changes, though commonly expected to indicate malignancy, may lead to improved blood cell production or stem cell survival. Am. J. Hematol. 91:1227-1233, 2016. © 2016 Wiley Periodicals, Inc.

Emergence of clonal hematopoiesis in the majority of patients with acquired aplastic anemia.

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.

Progressive increase in mtDNA 3243A>G heteroplasmy causes abrupt transcriptional reprogramming.

Variation in the intracellular percentage of normal and mutant mitochondrial DNAs (mtDNA) (heteroplasmy) can be associated with phenotypic heterogeneity in mtDNA diseases. Individuals that inherit the common disease-causing mtDNA tRNA(Leu(UUR)) 3243A>G mutation and harbor ∼10-30% 3243G mutant mtDNAs manifest diabetes and occasionally autism; individuals with ∼50-90% mutant mtDNAs manifest encephalomyopathies; and individuals with ∼90-100% mutant mtDNAs face perinatal lethality. To determine the basis of these abrupt phenotypic changes, we generated somatic cell cybrids harboring increasing levels of the 3243G mutant and analyzed the associated cellular phenotypes and nuclear DNA (nDNA) and mtDNA transcriptional profiles by RNA sequencing. Small increases in mutant mtDNAs caused relatively modest defects in oxidative capacity but resulted in sharp transitions in cellular phenotype and gene expression. Cybrids harboring 20-30% 3243G mtDNAs had reduced mtDNA mRNA levels, rounded mitochondria, and small cell size. Cybrids with 50-90% 3243G mtDNAs manifest induction of glycolytic genes, mitochondrial elongation, increased mtDNA mRNA levels, and alterations in expression of signal transduction, epigenomic regulatory, and neurodegenerative disease-associated genes. Finally, cybrids with 100% 3243G experienced reduced mtDNA transcripts, rounded mitochondria, and concomitant changes in nuclear gene expression. Thus, striking phase changes occurred in nDNA and mtDNA gene expression in response to the modest changes of the mtDNA 3243G mutant levels. Hence, a major factor in the phenotypic variation in heteroplasmic mtDNA mutations is the limited number of states that the nucleus can acquire in response to progressive changes in mitochondrial retrograde signaling.

Single nucleotide polymorphism array analysis of bone marrow failure patients reveals characteristic patterns of genetic changes.

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.

Efficient digest of high-throughput sequencing data in a reproducible report.

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.

NMNAT1 mutations cause Leber congenital amaurosis.

Leber congenital amaurosis (LCA) is an infantile-onset form of inherited retinal degeneration characterized by severe vision loss(1,2). Two-thirds of LCA cases are caused by mutations in 17 known disease-associated genes(3) (Retinal Information Network (RetNet)). Using exome sequencing we identified a homozygous missense mutation (c.25G>A, p.Val9Met) in NMNAT1 that is likely to be disease causing in two siblings of a consanguineous Pakistani kindred affected by LCA. This mutation segregated with disease in the kindred, including in three other children with LCA. NMNAT1 resides in the previously identified LCA9 locus and encodes the nuclear isoform of nicotinamide mononucleotide adenylyltransferase, a rate-limiting enzyme in nicotinamide adenine dinucleotide (NAD(+)) biosynthesis(4,5). Functional studies showed that the p.Val9Met alteration decreased NMNAT1 enzyme activity. Sequencing NMNAT1 in 284 unrelated families with LCA identified 14 rare mutations in 13 additional affected individuals. These results are the first to link an NMNAT isoform to disease in humans and indicate that NMNAT1 mutations cause LCA.

Mitochondrial genome sequence analysis: a custom bioinformatics pipeline substantially improves Affymetrix MitoChip v2.0 call rate and accuracy.

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.

Genome-wide analysis of interferon regulatory factor I binding in primary human monocytes.

IRF1 is a transcription factor that participates in interferon signaling. Previous studies of IRF1 binding have utilized in vitro assays. We used ChIP-seq in human monocytes to better define the recognition motif for IRF1. The newly identified 18bp motif (RAAASNGAAAGTGAAASY) is a refinement of the 13bp IRF1 motif commonly used. We utilized the 18bp consensus motif and identified 345 potential target genes. To compare the 18bp motif with the 13bp motif, we compared putative gene targets. Only 56 potential gene targets were defined by both consensus motifs. To compare biological effects of interferon on the 13bp and the 18bp consensus targets, we mined expression data from cells exposed to interferons or transfected with IRF1. In all cases, the 18bp consensus motif was more strongly associated with transcriptional responses than the 13bp motif. Therefore, the new 18bp consensus motif appears to have a greater association with biological activities of IRF1.

High-Resolution genomic arrays identify CNVs that phenocopy the chromosome 22q11.2 deletion syndrome.

The 22q11 Deletion Syndrome includes the overlapping phenotypes of DiGeorge/Velocardiofacial Syndromes, characterized by conotruncal heart defects, cleft palate, thymus, and parathyroid gland dysplasia. The majority (90%) of patients harbor detectable chr22q11.2 deletions, but a genetic etiology for the remainder of patients without a deletion can remain undefined despite major birth defects. We analyzed DNA from eight patients with normal 22q11 FISH studies by high-density single nucleotide polymorphism (SNP) arrays and identified potentially pathogenic copy number variants (CNVs) in four of eight patients. Two patients showed large CNVs in regions of known genomic disorders: one a deletion of distal chr22q11.2 and the other a duplication of chr5q35. A 3-Mb deletion of chr19p13.3 that includes a gene associated with conotruncal heart defects was found in a third patient. Two potentially pathogenic CNVs were found in a fourth patient: a large heterozygous deletion of chr6p24 and a smaller duplication of chr9p24. Our findings support a recent consensus statement advocating chromosomal microarray analysis as a first-line diagnostic approach for patients with multiple congenital anomalies. In patients with phenotypes suggestive of the 22q11.2 syndrome spectrum and normal FISH, microarray analysis can uncover the molecular basis of other genomic disorders whose features overlap those of 22q11.2 deletions.

Cytokine-induced monocyte characteristics in SLE.

Monocytes in SLE have been described as having aberrant behavior in a number of assays. We examined gene expression and used a genome-wide approach to study the posttranslational histone mark, H4 acetylation, to examine epigenetic changes in SLE monocytes. We compared SLE monocyte gene expression and H4 acetylation with three types of cytokine-treated monocytes to understand which cytokine effects predominated in SLE monocytes. We found that gamma-interferon and alpha-interferon both replicated a broad range of the gene expression changes seen in SLE monocytes. H4 acetylation in SLE monocytes was overall higher than in controls and there was less correlation of H4ac with cytokine-treated cells than when gene expression was compared. A set of chemokine genes had downregulated expression and H4ac. Therefore, there are significant clusters of aberrantly expressed genes in SLE which are strongly associated with altered H4ac, suggesting that these cells have experienced durable changes to their epigenome.

CNV Workshop: an integrated platform for high-throughput copy number variation discovery and clinical diagnostics.

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.

High-resolution mapping and analysis of copy number variations in the human genome: a data resource for clinical and research applications.

We present a database of copy number variations (CNVs) detected in 2026 disease-free individuals, using high-density, SNP-based oligonucleotide microarrays. This large cohort, comprised mainly of Caucasians (65.2%) and African-Americans (34.2%), was analyzed for CNVs in a single study using a uniform array platform and computational process. We have catalogued and characterized 54,462 individual CNVs, 77.8% of which were identified in multiple unrelated individuals. These nonunique CNVs mapped to 3272 distinct regions of genomic variation spanning 5.9% of the genome; 51.5% of these were previously unreported, and >85% are rare. Our annotation and analysis confirmed and extended previously reported correlations between CNVs and several genomic features such as repetitive DNA elements, segmental duplications, and genes. We demonstrate the utility of this data set in distinguishing CNVs with pathologic significance from normal variants. Together, this analysis and annotation provides a useful resource to assist with the assessment of CNVs in the contexts of human variation, disease susceptibility, and clinical molecular diagnostics.

Genomic analysis using high-density single nucleotide polymorphism-based oligonucleotide arrays and multiplex ligation-dependent probe amplification provides a comprehensive analysis of INI1/SMARCB1 in malignant rhabdoid tumors.

PURPOSE: A high-resolution genomic profiling and comprehensive targeted analysis of INI1/SMARCB1 of a large series of pediatric rhabdoid tumors was done. The aim was to identify regions of copy number change and loss of heterozygosity (LOH) that might pinpoint additional loci involved in the development or progression of rhabdoid tumors and define the spectrum of genomic alterations of INI1 in this malignancy. EXPERIMENTAL DESIGN: A multiplatform approach using Illumina single nucleotide polymorphism-based oligonucleotide arrays, multiplex ligation-dependent probe amplification, fluorescence in situ hybridization, and coding sequence analysis was used to characterize genome-wide copy number changes, LOH, and genomic alterations of INI1/SMARCB1 in a series of pediatric rhabdoid tumors. RESULTS: The biallelic alterations of INI1 that led to inactivation were elucidated in 50 of 51 tumors. INI1 inactivation was shown by a variety of mechanisms, including deletions, mutations, and LOH. The results from the array studies highlighted the complexity of rearrangements of chromosome 22 compared with the low frequency of alterations involving the other chromosomes. CONCLUSIONS: The results from the genome-wide single nucleotide polymorphism array analysis suggest that INI1 is the primary tumor suppressor gene involved in the development of rhabdoid tumors with no second locus identified. In addition, we did not identify hotspots for the breakpoints in sporadic tumors with deletions of chromosome 22q11.2. By employing a multimodality approach, the wide spectrum of alterations of INI1 can be identified in the majority of patients, which increases the clinical utility of molecular diagnostic testing.

A 3.1-Mb microdeletion of 3p21.31 associated with cortical blindness, cleft lip, CNS abnormalities, and developmental delay.

We report a 3.1-Mb de novo deletion of 3p21.31 in a 3.5-year-old female with cortical blindness, cleft lip, CNS abnormalities, and gross developmental delays. Examination of the region showed approximately 80 genes to be involved in the deletion. Functional analysis of the deleted genes suggests that several of them may be important in normal neuronal maturation and function. Thus, haploinsufficiency of one or more of these genes could potentially contribute to the observed phenotype. Our patient does not have clinical features that overlap completely with either proximal or distal 3p deletions, suggesting that the deletion seen in our patient leads to a distinct clinical phenotype not described previously.

Systematic epitope analysis of the p26 EIAV core protein.

The major core protein of equine infectious anemia virus (EIAV), p26, is one of the primary immunogenic structural proteins during a persistent infection of horses and is highly conserved among antigenically variants of viral isolates. In order to investigate its immune profile in more detail for a better diagnostic, an epitope mapping was carried out by means of two libraries of overlapping peptide fragments prepared by simultaneous and parallel SPPS on derivatized cellulose membranes (SPOT synthesis). Polyclonal equine sera from infected horses were used for the biological assay. Particularly two promising continuous epitopes (NAMRHL and MYACRD) were localized on the C-terminal extreme of p26, region 194-222. A cyclic synthetic fragment of 29 amino acid residues containing the identified epitopes was designed and studied. A significant conformational change towards a helical structure was observed when the peptide was cyclized by a bridge between Cys198 and Cys218. This observation correlated with an improvement of its ability to be recognized by specific antibodies in an EIA (Enzyme-linked Immunosorbent assay). These results suggest that the conformationally restricted synthetic antigen adequately mimics the native structure of this region of p26 core protein.

Monitoring substrate and products in a bioprocess with FTIR spectroscopy coupled to artificial neural networks enhanced with a genetic-algorithm-based method for wavelength selection.

An experiment was developed as a simple alternative to existing analytical methods for the simultaneous quantitation of glucose (substrate) and glucuronic acid (main product) in the bioprocesses Kombucha by using FTIR spectroscopy coupled to multivariate calibration (partial least-squares, PLS-1 and artificial neural networks, ANNs). Wavelength selection through a novel ranked regions genetic algorithm (RRGA) was used to enhance the predictive ability of the chemometric models. Acceptable results were obtained by using the ANNs models considering the complexity of the sample and the speediness and simplicity of the method. The accuracy on the glucuronic acid determination was calculated by analysing spiked real fermentation samples (recoveries ca. 115%).