Faculty Bibliography

Myeloid and erythroid precursor vacuolation is a common dysplastic finding associated with myeloid malignancies, toxins, drug, and nutritional deficiencies. It has been described as a core morphologic feature in VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome. We sought to determine the number of cases attributable to VEXAS syndrome in bone marrow biopsies and aspirates (BAMB) reporting myeloid precursor vacuolation. We reviewed 1318 individual BAMB reports from January 2020 to July 2021 where "vacuole(s)," "vacuolation," or "vacuolated" was reported. Bone marrow biopsies with vacuolation confined to blasts or those completed as routine workup prior to stem cell transplant or post induction chemotherapy for AML (acute myeloid leukemia) were excluded. Myeloid and erythroid precursor vacuolation was noted in 219 reports representing 210 patients. The most common etiology was myelodysplastic syndrome (MDS) (38.6%), AML (16.7%), lymphoproliferative disorders and multiple myeloma (7.6%), drug or toxin exposure (5.2%) myeloproliferative neoplasm (MPN) or MPN/MDS overlap syndrome (4.3%). VEXAS syndrome was determined to be the etiology in 2.9% of patients. Two additional cases of VEXAS syndrome with bone marrow biopsies reported in the specified time frame did not explicitly report myeloid or erythroid precursor vacuolation but were identified based on clinical suspicion and repeat BAMB review. Myeloid and erythroid precursor vacuolation is a dysplastic feature attributable to VEXAS syndrome in at least 2.9% of cases. Standardized reporting of vacuolization, triaging of molecular sequencing and optimal treatment of this disorder are critical issues facing those seeing patients with suspected VEXAS syndrome.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

Detecting mutations from single DNA molecules is crucial in many fields but challenging. Next-generation sequencing (NGS) affords tremendous throughput but cannot directly sequence double-stranded DNA molecules ('single duplexes') to discern the true mutations on both strands. Here we present Concatenating Original Duplex for Error Correction (CODEC), which confers single duplex resolution to NGS. CODEC affords 1,000-fold higher accuracy than NGS, using up to 100-fold fewer reads than duplex sequencing. CODEC revealed mutation frequencies of 2.72 × 10^-8 in sperm of a 39-year-old individual, and somatic mutations acquired with age in blood cells. CODEC detected genome-wide, clonal hematopoiesis mutations from single DNA molecules, single mutated duplexes from tumor genomes and liquid biopsies, microsatellite instability with 10-fold greater sensitivity and mutational signatures, and specific tumor mutations with up to 100-fold fewer reads. CODEC enables more precise genetic testing and reveals biologically significant mutations, which are commonly obscured by NGS errors.

Congenital myasthenic syndrome (CMS) is a heterogeneous condition associated with 34 different genes, including SLC5A7, which encodes the high affinity choline transporter 1 (CHT1). CHT1 is expressed in presynaptic neurons of the neuromuscular junction where it uses the inward sodium gradient to re-uptake choline. Bi-allelic CHT1 mutations often lead to neonatal lethality, and less commonly to non-lethal motor weakness and developmental delays. Here, we report detailed biochemical characterization of two novel mutations in CHT1, p.I294T and p.D349N, that we identified in an 11 year-old patient with a history of neonatal respiratory distress, and subsequent hypotonia and global developmental delay. Heterologous expression of each CHT1 mutant in human embryonic kidney cells showed two different mechanisms of reduced protein function. The p.I294T CHT1 mutant transporter function was detectable, but its abundance and half-life were significantly reduced. In contrast, the p.D349N CHT1 mutant was abundantly expressed at the cell membrane, but transporter function was absent. The residual function of the p.I294T CHT1 mutant may explain the non-lethal form of CMS in this patient, and the divergent mechanisms of reduced CHT1 function that we identified may guide future functional studies of the CHT1 myasthenic syndrome. Based on these in vitro studies that provided a diagnosis, treatment with cholinesterase inhibitor together with physical and occupational therapy significantly improved the patient's strength and quality of life.

IMPORTANCE:VEXAS (vacuoles, E1-ubiquitin-activating enzyme, X-linked, autoinflammatory, somatic) syndrome is a disease with rheumatologic and hematologic features caused by somatic variants in UBA1. Pathogenic variants are associated with a broad spectrum of clinical manifestations. Knowledge of prevalence, penetrance, and clinical characteristics of this disease have been limited by ascertainment biases based on known phenotypes. OBJECTIVE:To determine the prevalence of pathogenic variants in UBA1 and associated clinical manifestations in an unselected population using a genomic ascertainment approach. DESIGN, SETTING, AND PARTICIPANTS:This retrospective observational study evaluated UBA1 variants in exome data from 163 096 participants within the Geisinger MyCode Community Health Initiative. Clinical phenotypes were determined from Geisinger electronic health record data from January 1, 1996, to January 1, 2022. EXPOSURES:Exome sequencing was performed. MAIN OUTCOMES AND MEASURES:Outcome measures included prevalence of somatic UBA1 variation; presence of rheumatologic, hematologic, pulmonary, dermatologic, and other findings in individuals with somatic UBA1 variation on review of the electronic health record; review of laboratory data; bone marrow biopsy pathology analysis; and in vitro enzymatic assays. RESULTS:In 163 096 participants (mean age, 52.8 years; 94% White; 61% women), 11 individuals harbored likely somatic variants at known pathogenic UBA1 positions, with 11 of 11 (100%) having clinical manifestations consistent with VEXAS syndrome (9 male, 2 female). A total of 5 of 11 individuals (45%) did not meet criteria for rheumatologic and/or hematologic diagnoses previously associated with VEXAS syndrome; however, all individuals had anemia (hemoglobin: mean, 7.8 g/dL; median, 7.5 g/dL), which was mostly macrocytic (10/11 [91%]) with concomitant thrombocytopenia (10/11 [91%]). Among the 11 patients identified, there was a pathogenic variant in 1 male participant prior to onset of VEXAS-related signs or symptoms and 2 female participants had disease with heterozygous variants. A previously unreported UBA1 variant (c.1861A>T; p.Ser621Cys) was found in a symptomatic patient, with in vitro data supporting a catalytic defect and pathogenicity. Together, disease-causing UBA1 variants were found in 1 in 13 591 unrelated individuals (95% CI, 1:7775-1:23 758), 1 in 4269 men older than 50 years (95% CI, 1:2319-1:7859), and 1 in 26 238 women older than 50 years (95% CI, 1:7196-1:147 669). CONCLUSIONS AND RELEVANCE:This study provides an estimate of the prevalence and a description of the clinical manifestations of UBA1 variants associated with VEXAS syndrome within a single regional health system in the US. Additional studies are needed in unselected and genetically diverse populations to better define general population prevalence and phenotypic spectrum.

OBJECTIVE/UNASSIGNED:, the gene encoding the anti-inflammatory IL-1 receptor antagonist (IL-1Ra), on the cytokine release syndrome and mortality. METHODS/UNASSIGNED:gene were assessed for association with laboratory markers of the cytokine release syndrome (CRS) and mortality. RESULTS/UNASSIGNED:rs419598 CC SNV exhibited lower inflammatory biomarker levels, and was associated with reduced mortality compared to the CT/TT genotype in men (OR 0.49 (0.23 - 1.00); 0.052), with the most pronounced effect observed between the ages of 55-74 [5.5% vs. 18.4%, p<0.001]. CONCLUSION/UNASSIGNED:modulates the COVID-19 cytokine release syndrome via endogenous " anti-inflammatory" mechanisms. SIGNIFICANCE STATEMENT/UNASSIGNED:merits further evaluation in severe SARS-CoV-2 infection.

Chromatin accessibility assays are central to the genome-wide identification of gene regulatory elements associated with transcriptional regulation. However, the data have highly variable quality arising from several biological and technical factors. To surmount this problem, we developed a sequence-based machine learning method to evaluate and refine chromatin accessibility data. Our framework, gapped k-mer SVM quality check (gkmQC), provides the quality metrics for a sample based on the prediction accuracy of the trained models. We tested 886 DNase-seq samples from the ENCODE/Roadmap projects to demonstrate that gkmQC can effectively identify "high-quality" (HQ) samples with low conventional quality scores owing to marginal read depths. Peaks identified in HQ samples are more accurately aligned at functional regulatory elements, show greater enrichment of regulatory elements harboring functional variants, and explain greater heritability of phenotypes from their relevant tissues. Moreover, gkmQC can optimize the peak-calling threshold to identify additional peaks, especially for rare cell types in single-cell chromatin accessibility data.

Numerous applications in molecular biology and genomics require characterization of mutant DNA molecules present at low levels within a larger sample of non-mutant DNA. This is often achieved either by selectively amplifying mutant DNA, or by sequencing all the DNA followed by computational identification of the mutant DNA. However, selective amplification is challenging for insertions and deletions (indels). Additionally, sequencing all the DNA in a sample may not be cost effective when only the presence of a mutation needs to be ascertained rather than its allelic fraction. The MutS protein evolved to detect DNA heteroduplexes in which the two DNA strands are mismatched. Prior methods have utilized MutS to enrich mutant DNA by hybridizing mutant to non-mutant DNA to create heteroduplexes. However, the purity of heteroduplex DNA these methods achieve is limited because they can only feasibly perform one or two enrichment cycles. We developed a MutS-magnetic bead system that enables rapid serial enrichment cycles. With six cycles, we achieve complete purification of heteroduplex indel DNA originally present at a 5% fraction and over 40-fold enrichment of heteroduplex DNA originally present at a 1% fraction. This system may enable novel approaches for enriching mutant DNA for targeted sequencing. This article is protected by copyright. All rights reserved.