Klee, Eric W.

Introduction: Common genetic variation has been linked to an increased risk of developing dilated cardiomyopathy (DCM). However, the phenotypic profile, outcomes, and biological mechanisms by which such genetic background may predispose individuals to DCM are still unknown. Hypothesis: The primary objectives of this study were to develop a DCM polygenic risk score (PRS) based on prior case control-DCM GWAS and to apply this PRS to a well clinical characterized general population. Methods: The DCM-PRS was derived from 13 DCM risk alleles that were assayed in an unselected cohort (n=1897) of unrelated Olmsted County residents age >45 years who did not have a history of heart failure (HF) or CM and underwent physical examination, echocardiography, and blood analyses. Results: The entire population was equally divided into quartiles based on the DCM-PRS. No significant differences were present in terms of age (median 62 years) and sex (48% male). However, a significant correlation was detected between the DCM-PRS and left ventricle ejection fraction (LVEF) (Beta: -2.68, p<0.001), LV end-systolic diameter (Beta: 0.22, p<0.001), and LV end-diastolic diameter (Beta: 0.14, p=0.01). During a follow-up of 16 years (IQR: 15-17), 47 subjects developed DCM. Using a 3:1 non-DCM vs DCM matched cohort, a logistic regression analysis, found a significant correlation between median DCM-PRS and DCM development [OR=2.27 (1.06, 4.85), p=0.03]. However, there was no association between the DCM-PRS and HF biomarkers (aldosterone, natriuretic peptides, ST2, Galectin-3). Conclusions: DCM-PRS correlates with decreased LVEF, increased LV dimensions, and DCM development in a general population without significant differences in baseline clinical features. The absence of correlation between the DCM-PRS and available HF biomarkers suggests that this genetic background affects LV traits and DCM development independent of these biological mechanisms.

Introduction: Telomere biology disorders (TBD) are multisystem diseases often associated with bone marrow failure (BMF) and pulmonary fibrosis. TBD are clinically diagnosed when mean telomere lengths (TL) in blood cells measured by clinical flowFISH assays fall &lt;1 st percentile, in comparison to age-matched reference controls. This information has limitations, e.g. TL in blood cells may not represent TL in other organs (e.g. lung, liver) and moreover, individuals with high clinical suspicion of TBD can show borderline average TL which complicates their diagnosis. The likely reason for this is that cell senescence can be triggered by shortening of a subset, or even a single chromosome arm telomere, a situation that may be missed when measuring average TL. Previous approaches to measure individual chromosome-arm telomere length (CA-TL) are only available for a small number of chromosome arms or can only estimate the number of shortened telomeres without matching them to specific chromosome arms. Taking advantage of recent developments in long-read next generation sequencing (NGS) and the availability of the recently published telomere-to-telomere (T2T) reference human genome, our group developed Telogator, (Stephens Z. et al, 2022) a bioinformatic tool that estimates CA-TL for most chromosome arms using long-read NGS data. Here we use Telogator to describe for the first time the CA-TL pattern in TBD patients in comparison with a control non-TBD cohort. Methods: TBD patients were identified through the Pre-Myeloid and Bone Marrow Failure Specialty Clinic at Mayo Clinic. All patients underwent flowFISH and genetic testing. Age-matched non-TBD patients identified at the Center for Individualized Medicine (Mayo Clinic) were used as reference controls. Informed consent was obtained from all individuals with approval from the Mayo Clinic Institutional Review Board. DNA from peripheral blood mononuclear cells was sequenced using PacBio Sequel II following the manufacturer indications for HiFi sequencing. BAM files were aligned to the T2T reference and CA-TL was calculated using Telogator as described in Stephens Z. et al. 2022. Average TL, CA-TL, and percent of short telomeres (both defined as &lt;5 kb and &lt;3.2 kb) were calculated. Comparisons of average TL and percent of short telomeres were performed using t-test, while Wilcoxon rank-sum test was used to compare CA-TL. Results: Four TBD patients [median age = 58.5 (52-60), 50% female] and 6 non-TBD controls [Median age = 56 (50-63), 33% female] were included in the analysis. All TBD patients presented average TL below the 1 st centile in both lymphocytes and granulocytes, and TBD-related phenotype (100% with pulmonary fibrosis, 50% with BMF, 25% with liver fibrosis). Three had a pathogenic variant in a TBD related genes [2 in TERT (NM_198253.3: c.2030G&gt;A, p.(Gly677Asp) and c.2515_2533del, p.(Thr839Alafs*29)) and 1 in NHP2 (NM_017838.3: c.459_460delAT; p.(X154Argext29))]. The average TL calculated by Telogator showed high correlation with flowFISH values (r 2 = 0.9883, p = 0.0059) and was lower in TBD patients than controls (p = 0.0012). Next, when comparing individual CA-TL between groups, only half of CA-TL was statistically different in the TBD group compared to controls (1p, 3q, 4p, 5p, 5q, 6q, 7p, 7q, 8q, 9p, 9q, 10q, 11p, 12q, 12q, 16p, 16q, 18q, 19p, 19q, 20p, 20q, 21q - Figure 1). The percentage of telomeres &lt;5kb and &lt;3 kb was significantly increased in TBD patients compared to controls (p= 0.014 and 0.0333, respectively). Discussion: Using a novel informatic tool, we replicate previous results obtained using classic approaches and were able to expand on these studies by measuring specific CA-TL in TBD patients. We identified statistically significant differences at the CA-TL level in these patients when compared to age-matched controls, while other CA-TL were not statistically different. Future studies and the inclusion of additional TBD patients will help define the clinical relevance of the differential CA-TL shortening and its prevalence among TBD patients. While limited by the small number of samples in the study, we show the feasibility of using this approach to study TL with unprecedented resolution and show that Telogator offers important advantages compared to classic methods to study telomere biology that can be used in other relevant fields like aging and/or cancer. Our results pave the way to advance the mechanistic understanding of TBD.

Introduction: Despite the improved clinical outcomes with novel therapies in patients with chronic lymphocytic leukemia (CLL), Richter transformation (RT) of CLL remains a major clinical challenge. Our group showed a selective response (~40%) with the PD-1 blocking antibody pembrolizumab in patients with RT, particularly after prior exposure to ibrutinib (Ding et al, Blood, 2017). It has been established that PD-1 inhibits cytotoxic T-cell functions via the interaction with its ligands, PD-L1 and PD-L2. However, the role of PD-1 in tumor B-cells is understudied. Here we show the expression and functional implications of the PD-1 signaling axis in B-cell pathobiology. Methods: CLL-involved lymph node (LN) (n=26) and RT (with diffuse large B-cell lymphoma [DLBCL] histology) -involved LN (n=20) samples were tested for PD-1 expression by immunohistochemistry (IHC) (mouse clone NAT105, Abcam). Then, PD-1 expression was examined in 10 lymphoma cell lines and 1 CLL-like cell line (MEC-1) by both flow cytometry (FC) and Western blot (WB) analysis. PD-1 overexpression using pLEX-lentiviral packaging system (Thermo Scientific) was evaluated for its impact on in vitro cell cycle regulation and in vivo tumor growth. CRISPR/CAS9 mediated PD-1 knockdown evaluated on cell growth in MEC-1 and Mino cells. Role of PD-1 was further examined on the regulation of cell signaling pathways using human phospho kinase array kit (R&D) and WB analysis. Additionally, gene expression signatures were evaluated by RNA-seq using pLEX-lentiviral-PD-1 transfected DLBCL cell lines OCI-LY7 and OCI-LY19. Results: PD-1 expression was significantly increased in RT-LN compared to CLL-LN (mean ± SEM in RT vs. CLL, 30.6% ± 4.7% vs. 11.5% ± 2.8%, p&lt;0.001) as determined by IHC. PD-1 expression was highest in patients with RT where the last prior CLL therapy was ibrutinib. To test the role of PD-1 in tumor B-cells, its expression was assessed in lymphoma and CLL cell lines. The expression of PD-1 was found to be variable by both FC and WB analysis, but at very low-levels in DLBCL cell lines OCI-LY7 and OCI-LY19. Constitutive lentiviral (pLEX-PD-1)-mediated overexpression of PD-1 in OCI-LY7 and OCI-LY19 cells led to increased cell growth (1.4 and 1.9 fold, respectively, compared to original lines after 4 days of culturing of equal number of cells), which was further confirmed by cell cycle analysis which showed an increase in S phase by 20.4% (p&lt;0.05) and 24.5% (p&lt;0.05) in OCI-LY7 and OCI-LY19 cells, respectively. When the luciferase+ PD-1 overexpressing OCI-LY7 cells (1X10 6) were injected intravenously into NSG mice, increased tumor growth (p&lt;0.01) was observed by bioluminescence imaging after 22 days of follow-up. PD-1 knockdown in MEC-1 and Mino cells (high basal PD-1 expression) led to significant inhibition of cell growth (85% and 79.5%, respectively, p&lt;0.05). Using phospho-kinase array, an overall decrease of phosphorylation was detected on multiple sites of p53 (S392, S15, S46) and Chk-2 (T68) in PD-1 overexpressing OCI-LY7 and OCI-LY19 cell lines. This finding was further confirmed by WB analyses. In addition, decreased total p53 and increased phosphorylation on both SHP-1 and SHP-2 were found in PD-1 overexpressing OCI-LY7 and OCI-LY19 cell lines. SHP-1 (TPI-1) and SHP-1+SHP-2 (TPI-1+TNO155) inhibitor treatment promoted OCI-LY19 cell death and led to recovery of phosphorylation on ATM, Chk-2 and p53 in these cells. RNA sequencing data in PD-1 overexpressing OCI-LY7 and OCI-LY19 cells revealed increased expression of genes that promote cell cycle progression and DNA replication. Conclusion: Our data showed robust PD-1 expression in patients with Richter transformation from CLL to DLBCL. PD-1 overexpression in DLBCL lymphoma cell lines enhanced cell proliferation in vitro and in vivo. Further investigation identified that PD-1 modified the phosphorylation/function of SHP phosphatases and thereby regulated p53 pathways ( Figure 1). These findings indicate that RT with DLBCL cells have the ability to “hijack” the PD-1 pathway which results in downregulation of p53 mediated DNA repair. This novel finding provides for new strategies to further evaluate the interaction of PD-1 signals with DNA repair pathways and suggest potential novel targets for RT.

ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of virus evolution and spread and to inform policy decisions. We sequenced viral genomes from over 22,000 patient samples tested at Mayo Clinic Laboratories between 2020 and 2022 and used Bayesian phylodynamics to describe county and regional spread in Minnesota. The earliest calculated introduction into Minnesota was to Hennepin County from a domestic source around 22 January 2020; 6 weeks before the first confirmed case in the state. This led to the virus spreading to Northern Minnesota and, eventually, the rest of the state. International introductions were most abundant in Hennepin (home to the Minneapolis/St. Paul International Airport) totaling 45 (out of 107) over the 2-year period. Southern Minnesota counties were most common for domestic introductions, with 19 (out of 64), potentially driven by bordering states such as Iowa and Wisconsin as well as Illinois, which is nearby. Hennepin also was, by far, the most dominant source of in-state transmissions to other Minnesota locations ( n = 772) over the 2-year period. We also analyzed the diversity of the location source of SARS-CoV-2 viruses in each county and noted the timing of state-wide policies as well as trends in clinical cases. Neither the number of clinical cases nor the major policy decisions, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, appeared to have an impact on virus diversity across each individual county. IMPORTANCE We analyzed over 22,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes of patient samples tested at Mayo Clinic Laboratories during a 2-year period in the COVID-19 pandemic, which included Alpha, Delta, and Omicron variants of concern to examine the roles and relationships of Minnesota virus transmission. We found that Hennepin County, the most populous county, drove the transmission of SARS-CoV-2 viruses in the state after including the formation of earlier clades including 20A, 20C, and 20G, as well as variants of concern Alpha and Delta. We also found that Hennepin County was the source for most of the county-to-county introductions after an initial predicted introduction with the virus in early 2020 from an international source, while other counties acted as transmission “sinks.” In addition, major policies, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, did not appear to have an impact on virus diversity across individual counties.

ABSTRACT Cytomegalovirus (CMV) is a significant cause of morbidity and mortality among immunocompromised hosts, including transplant recipients. Antiviral prophylaxis or treatment is used to reduce the incidence of CMV disease in this patient population; however, there is concern about increasing antiviral resistance. Detection of antiviral resistance in CMV was traditionally accomplished using Sanger sequencing of UL54 and UL97 genes, in which specific mutations may result in reduced antiviral activity. In this study, a novel next-generation sequencing (NGS) method was developed and validated to detect mutations in UL54 / UL97 associated with antiviral resistance. Plasma samples ( n = 27) submitted for antiviral resistance testing by Sanger sequencing were also analyzed using the NGS method. When compared to Sanger sequencing, the NGS assay demonstrated 100% (27/27) overall agreement for determining antiviral resistance/susceptibility and 88% (22/25) agreement at the level of resistance-associated mutations. The limit of detection of the NGS method was determined to be 500 IU/mL, and the lower threshold for detecting mutations associated with resistance was established at 15%. The NGS assay represents a novel laboratory tool that assists healthcare providers in treating patients who are infected with CMV harboring resistance-associated mutations and who may benefit from tailored antiviral therapy.

AbstractBackgroundThe PI3K/AKT pathway, extensively studied in cancer, is vital for regulating cell metabolism, differentiation, and proliferation. Pathogenic variants in the PIK3R1 gene, which encodes three regulatory units of class IA PI3Ks, have been found in affected tissue of individuals with vascular lesions. These variants predominantly occur in the iSH2 domain, disrupting inhibitory contacts with the catalytic unit and leading to PI3K activation. Germline variants in this gene are also linked to an immunological condition called Activated PI3K delta syndrome type 2 (APDS2).MethodsThis is a case report and literature review. Clinical data were retrieved from medical records.ResultsA male patient presented with extensive vascular malformation covering over 90% of his body, along with complete 2–3 toe syndactyly, suggesting a vascular malformation syndrome called PROS. Low levels of IgA and IgG were detected. The patient achieved his developmental milestones and had above‐average weight, height, and head circumference. Exome sequencing of skin and blood DNA revealed a de novo variant in PIK3R1 (c.1746‐2A>G, p.?) in 9% of the patient's blood cells and 25% of cultured fibroblasts. Initially, classified as a variant of uncertain significance, this variant was later confirmed to be the cause.ConclusionsThis is the first intronic SNV in a canonical splice site within iSH2 described, highlighting the importance of iSH2 in the regulation of the PI3K/AKT pathway and its involvement in the development of vascular overgrowth and antibody deficiency.

Background Pathogenic variants in the zinc finger protein coding genes are rare causes of intellectual disability and congenital malformations. Mutations in the ZNF148 gene causing GDACCF syndrome (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies; MIM #617260) have been reported in five individuals so far. Methods As a result of an international collaboration using GeneMatcher Phenome Central Repository and personal communications, here we describe the clinical and molecular genetic characteristics of 22 previously unreported individuals. Results The core clinical phenotype is characterised by developmental delay particularly in the domain of speech development, postnatal growth retardation, microcephaly and facial dysmorphism. Corpus callosum abnormalities appear less frequently than suggested by previous observations. The identified mutations concerned nonsense or frameshift variants that were mainly located in the last exon of the ZNF148 gene. Heterozygous deletion including the entire ZNF148 gene was found in only one case. Most mutations occurred de novo, but were inherited from an affected parent in two families. Conclusion The GDACCF syndrome is clinically diverse, and a genotype-first approach, that is, exome sequencing is recommended for establishing a genetic diagnosis rather than a phenotype-first approach. However, the syndrome may be suspected based on some recurrent, recognisable features. Corpus callosum anomalies were not as constant as previously suggested, we therefore recommend to replace the term ‘GDACCF syndrome’ with ‘ZNF148-related neurodevelopmental disorder’.

Abstract Introduction: Most rare disease patients (75 − 50%) undergoing genomic sequencing remain unsolved, often due to lack of information about variants identified. Data review over time can leverage novel information regarding disease-causing variants and genes, increasing this diagnostic yield. However, time and resource constraints have limited reanalysis of genetic data in clinical laboratories setting. We developed RENEW, (REannotation of NEgative WES/WGS) an automated reannotation procedure that uses relevant new information in on-line genetic databases to enable rapid review of genetic findings. Methods We tested RENEW in an unselected cohort of 1,066 undiagnosed cases from the Mayo Clinic Center for Individualized Medicine using new information in ClinVar, HGMD and OMIM between the date of previous analysis/testing and April of 2022. Results 5,741 variants prioritized by RENEW were rapidly reviewed by variant interpretation specialists. Mean analysis time was approximately 20 seconds per variant (32 hours total time). Reviewed cases were classified as: 879 (93.0%) undiagnosed, 63 (6.6%) putatively diagnosed, and 4 (0.4%) definitively diagnosed. Discussion New strategies are needed to enable efficient review of genomic findings in unsolved cases. We report on a fast and practical approach to address this need and improve overall diagnostic success in patient testing through a recurrent reannotation process.

Abstract Two inherited autosomal dominant cancer predisposition conditions - BRCA related hereditary breast and ovarian cancer (HBOC) and Lynch syndrome (LS) - are termed Centers for Disease Control and Prevention Tier 1 (CDCT1) genetic conditions, for which early identification and intervention have a positive impact on public health. Selection of individuals for genetic testing for these conditions is based on personal and family history. The goal of this study was to evaluate whether screening in a tertiary medical center using exome sequencing could efficiently identify carriers of HBOC and LS and to determine the frequency of incremental carriers identified outside of traditional clinical practice guidelines. Participants from three geographically diverse Mayo Clinic practices in the USA (Rochester MN, Phoenix AZ, Jacksonville FL) consented to clinical “Exome+” sequencing (Helix, San Mateo, CA; Mayo GeneGuide, Rochester, MN) in the TAPESTRY study (NCT05212428), which links sequencing data with electronic health records and the return of CDCT1 genetic findings. For this study we focused on two inherited cancer predisposition conditions: BRCA-related HBOC (BRCA1 and BRCA2) and LS (MLH1, MSH2, MSH6, PMS2 and EPCAM). Detailed chart review to collect demographic information, personal and family history, and assessment of clinical practice guidelines for genetic evaluation (National Cancer Control Network 2021). To date, 44,306 patients have enrolled and sequenced in TAPESTRY. Annotation and interpretation of all variants in 7 genes for HBOC and LS resulted in identification of 550 carriers (prevalence 1.24%) which included 387 with HBOC (27.2% BRCA1, 42.8% BRCA2) and 163 with LS (12.3% MSH6, 8.8% PMS2, 4.5% MLH1, 3.8% MSH2 and 0.2% EPCAM). Demographics of the cohort included: 62.7% female, mean age 55.2 years, non-white race 9.6% and 3.8% Hispanic/Latino ethnicity. A personal history of cancer was present in 46.4% of test positive patients, including 22.5% of HBOC patients with breast or ovarian cancer and 20.9% of LS patients with colorectal or uterine cancer. More than half of the patients (52.1%) were newly diagnosed with HBOC and LS based on the results of this study. Overall, 39.2% of HBOC/LS carriers identified through Exome+ sequencing did not satisfy NCCN criteria for genetic evaluation, this was higher for LS (56.2%) compared to HBOC (32%). Amongst those newly diagnosed with HBOC/LS (n=286), NCCN criteria were not satisfied in 60% of cases (78% for LS and 51% for HBOC). Of the pathogenic germline variant carriers who met NCCN guidelines for testing, 34.2% were not aware of their diagnosis prior to participation in this study. Our results emphasize the need for wide genomic screening for CDCT1 cancer predisposition syndromes. Such screening could identify at-risk carriers of HBOC and LS, who would not otherwise be identified through clinical practice guidelines. Citation Format: Emily Gay, Niloy Jewel Samadder, Michelle L. Bublitz, Melanie M. Peterson, Tammy A. Wilson, Lorelei A. Bandel, Sebastian M. Armasu, Robert A. Vierkant, Matthew J. Ferber, Eric W. Klee, Nicholas B. Larson, Teresa M. Kruisselbrink, Timothy B. Curry, Jan B. Egan, Jennifer L. Kemppainen, Jessa S. Bidwell, Jennifer L. Anderson, Tammy M. McAllister, T'Nita S. Walker, Katie L. Kunze, Vanda Lindpere, Michael A. Golafshar, Margaret Klint, Richard J. Presutti, William V. Bobo, Aleksander Sekulic, Jolene M. Summer Bolster, Cheryl L. Willman, Konstantinos N. Lazaridis. Genetic screening in a tertiary medical center identifies carriers of cancer predisposition diseases that would be missed by clinical guidelines [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5768.

Abstract Background Large β-globin gene cluster deletions (hereditary persistence of fetal hemoglobin [Hb] or β-, δβ-, γδβ-, and ϵγδβ-thalassemia), are associated with widely disparate phenotypes, including variable degrees of microcytic anemia and Hb F levels. When present, increased Hb A2 is used as a surrogate marker for β-thalassemia. Notably, ϵγδβ-thalassemias lack the essential regulatory locus control region (LCR) and cause severe transient perinatal anemia but normal newborn screen (NBS) results and Hb A2 levels. Herein, we report a novel deletion of the ϵ, Aγ, Gγ, and ψβ loci with intact LCR, δ-, and β-regions in 2 women and newborn twins. Methods Capillary electrophoresis (CE), high-performance liquid chromatography (HPLC), DNA sequencing, multiplex ligation-dependent probe amplification (MLPA), gap-polymerase chain reaction (gap-PCR), and long-read sequencing (LRS) were performed. Results NBS showed an Hb A &gt; Hb F pattern for both twins. At 20 months, Hb A2 was increased similarly to that in the mother and an unrelated woman. Unexplained microcytosis was absent and the twins lacked severe neonatal anemia. MLPA, LRS, and gap-PCR confirmed a 32 599 base pair deletion of ϵ (HBE1) through ψβ (HBBP1) loci. Conclusions This deletion represents a hemoglobinopathy category with a distinct phenotype that has not been previously described, an ϵγ-thalassemia. Both the NBS Hb A &gt; F pattern and the subsequent increased Hb A2 without microcytosis are unusual. A similar deletion should be considered when this pattern is encountered and appropriate test methods selected for detection. Knowledge of the clinical impact of this new category will improve genetic counselling, with distinction from the severe transient anemia associated with ϵγδβ-thalassemia.