Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease.

BACKGROUND: X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique. RESULTS: In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis. CONCLUSIONS: This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Genetic Testing Goes Beyond Imaging and Histological Evaluation in Pleuroparenchymal Fibroelastosis.

BACKGROUND: Lung biopsy remains the gold standard in the diagnosis of fibrotic interstitial lung disease (F-ILD), but there is a growing appreciation of the role of pathogenic gene variants in telomere and surfactant protein genes, especially in familial pulmonary fibrosis (FPF). Pleuroparenchymal fibroelastosis (PPFE) is a rare disease that can coexist with different patterns of F-ILD, including FPF. It can be progressive and often leads to respiratory failure and death. This study tested the hypothesis that genetic testing goes beyond radiological and histological findings in PPFE and other F-ILD further informing clinical decision-making for patients and affected family members by identifying pathological gene variants in telomere and surfactant protein genes. METHODS: This is a retrospective review of 70 patients with F-ILD in the setting of FPF or premature lung fibrosis. Six out of 70 patients were diagnosed with PPFE based on radiological or histological characteristics. All patients underwent telomere length evaluation in peripheral blood by Flow-FISH or genetic testing using a customized exome-based panel that included telomere and surfactant protein genes associated with lung fibrosis. RESULTS: Herein, we identified six individuals where radiographic or histopathological analyses of PPFE were linked with telomere biology disorders (TBD) or variants in surfactant protein genes. Each case involved individuals with either personal early-onset lung fibrosis or a family history of the disease. Assessments of telomere length and genetic testing offered insights beyond traditional radiological and histopathological evaluations. CONCLUSION: Detecting anomalies in TBD-related or surfactant protein genes can significantly refine the diagnosis and treatment strategies for individuals with PPFE and other F-ILD.

Semiautomated approach focused on new genomic information results in time and effort-efficient reannotation of negative exome data.

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 genomic databases to enable rapid review of genomic findings. We tested RENEW in an unselected cohort of 1066 undiagnosed cases with a broad spectrum of phenotypes 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. 5741 variants prioritized by RENEW were rapidly reviewed by variant interpretation specialists. Mean analysis time was approximately 20 s per variant (32 h total time). Reviewed cases were classified as: 879 (93.0%) undiagnosed, 63 (6.6%) putatively diagnosed, and 4 (0.4%) definitively diagnosed. 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.

Prognostic relevance of clonal hematopoiesis in myeloid neoplastic transformation in patients with follicular lymphoma treated with radioimmunotherapy.

While novel radioisotope therapies continue to advance cancer care, reports of therapy-related myeloid neoplasms (t-MN) have generated concern. The prevalence and role of clonal hematopoiesis (CH) in this process remain to be defined. We hypothesized that: (i) CH is prevalent in relapsed follicular lymphoma and is associated with t-MN transformation, and (ii) radiation in the form of radioimmunotherapy (RIT) plays a role in clonal progression. In this retrospective cohort study, we evaluated the prevalence and prognostic impact of CH on clinical outcomes in 58 heavily pre-treated follicular lymphoma patients who received RIT. Patients had been given a median of four lines of therapy before RIT. The prevalence of CH prior to RIT was 46%, while it was 67% (P=0.15) during the course of RIT and subsequent therapies in the paired samples. Fourteen (24%) patients developed t-MN. Patients with t-MN had a higher variant allele fraction (38% vs. 15%; P=0.02) and clonal complexity (P=0.03) than those without. The spectrum of CH differed from that in age-related CH, with a high prevalence of DNA damage repair and response pathway mutations, absence of spliceosome mutations, and a paucity of signaling mutations. While there were no clear clinical associations between RIT and t-MN, or overall survival, patients with t-MN had a higher mutant clonal burden, along with extensive chromosomal abnormalities (median survival, afer t-MN diagnosis, 0.9 months). The baseline prevalence of CH was high, with an increase in prevalence on exposure to RIT and subsequent therapies. The high rates of t-MN with marked clonal complexities and extensive chromosomal damage underscore the importance of better identifying and studying genotoxic stressors accentuated by therapeutic modalities.

Experimental and Computational Approaches to Measure Telomere Length: Recent Advances and Future Directions.

PURPOSE OF REVIEW: The length of telomeres, protective structures at the chromosome ends, is a well-established biomarker for pathological conditions including multisystemic syndromes called telomere biology disorders. Approaches to measure telomere length (TL) differ on whether they estimate average, distribution, or chromosome-specific TL, and each presents their own advantages and limitations. RECENT FINDINGS: The development of long-read sequencing and publication of the telomere-to-telomere human genome reference has allowed for scalable and high-resolution TL estimation in pre-existing sequencing datasets but is still impractical as a dedicated TL test. As sequencing costs continue to fall and strategies for selectively enriching telomere regions prior to sequencing improve, these approaches may become a promising alternative to classic methods. Measurement methods rely on probe hybridization, qPCR or more recently, computational methods using sequencing data. Refinements of existing techniques and new approaches have been recently developed but a test that is accurate, simple, and scalable is still lacking.

Patients with telomere biology disorders show context specific somatic mosaic states with high frequency of U2AF1 variants.

Somatic mosaic states in telomere biology disorders are characterized by somatic variants in the spliceosome and DNA damage response and repair pathways. A likely maladaptive response to short telomeres that may lead to increased hematological cancer.

Clinical and molecular correlates of somatic and germline DDX41 variants in patients and families with myeloid neoplasms.

The diagnosis of germline predisposition to myeloid neoplasms (MN) secondary to DDX41 variants is currently hindered by the long latency period, variable family histories and the frequent occurrence of DDX41 variants of uncertain significance (VUS). We reviewed 4,524 consecutive patients who underwent targeted sequencing for suspected or known MN and analyzed the clinical impact and relevance of DDX41VUS in comparison to DDX41path variants. Among 107 patients (44 [0.9%] DDX41path and 63 DDX41VUS [1.4%; 11 patients with both DDX41path and DDX41VUS]), we identified 17 unique DDX41path and 45 DDX41VUS variants: 24 (23%) and 77 (72%) patients had proven and presumed germline DDX41 variants, respectively. The median age was similar between DDX41path and DDX41VUS (66 vs. 62 years; P=0.41). The median variant allele frequency (VAF) (47% vs. 48%; P=0.62), frequency of somatic myeloid co-mutations (34% vs 25%; P= 0.28), cytogenetic abnormalities (16% vs. 12%; P=>0.99) and family history of hematological malignancies (20% vs. 33%; P=0.59) were comparable between the two groups. Time to treatment in months (1.53 vs. 0.3; P=0.16) and proportion of patients progressing to acute myeloid leukemia (14% vs. 11%; P=0.68), were similar. The median overall survival in patients with high-risk myelodysplastic syndrome/acute myloid leukemia was 63.4 and 55.7 months in the context of DDX41path and DDX41VUS, respectively (P=0.93). Comparable molecular profiles and clinical outcomes among DDX41path and DDX41VUS patients highlights the need for a comprehensive DDX41 variant interrogation/classification system, to improve surveillance and management strategies in patients and families with germline DDX41 predisposition syndromes.

TELOMERE LENGTH OF PERIPHERAL BLOOD MONONUCLEAR CELLS IS ASSOCIATED WITH DISCHARGE DISPOSITION IN OLDER TRAUMA PATIENTS.

Introduction: Little is known regarding peripheral blood mononuclear cell telomere length (PBMC-TL) and response to traumatic injury. The objective of this study was to characterize the role of PBMC-TL in coagulation and clinical outcomes after injury. Methods: Plasma and buffy coats were prospectively collected from trauma patients and healthy volunteers. DNA was purified and PBMC-TL quantified by quantitative polymerase chain reaction. Thrombin generation kinetics were expressed as lag time (in minutes), peak height (in nanometers), time to peak (in minutes), and endogenous thrombin potential (in nM × min). Results are in median and quartiles [Q1, Q3]. P < 0.05 was considered significant (Wilcoxon rank sum testing). Results: Forty-two younger patients (21 [20, 22] years, 69% were male) and 39 older patients (62 [61, 64] years, 79% were male) were included. There was no significant difference in Clinical Frailty Scores between groups. Younger patients had longer total PBMC-TL (0.40 Mb [0.30, 0.49] vs. 0.29 Mb [0.23, 0.33], P < 0.001) and longer average PBMC-TL per chromosome (4.3 kb [3.3, 5.3] vs. 3.2 kb [2.5, 3.7], P < 0.001). When older patients were stratified by 50th percentile of PBMC-TL, there were no differences in thrombin generation; however, those with shorter telomeres were less likely to be discharged home (29% vs. 77%, P = 0.004). Older patients in the bottom quartile of PBMC-TL had shorter lag time (2.78 min [2.33, 3.00] vs. 3.33 min [3.24, 3.89], P = 0.030) and were less likely to be discharged home (22% vs. 90%, P = 0.006) than those in the top quartile of PBMC-TL. Multivariable logistic regression models revealed both increased age and shorter PBMC-TL to be independent predictors of discharge disposition other than home. Conclusion: In older trauma patients, shorter PBMC-TL is associated with accelerated initiation of thrombin generation and lower likelihood of being discharged to home.

Early expression of mature αß TCR in CD4-CD8- T cell progenitors enables MHC to drive development of T-ALL bearing NOTCH mutations.

During normal T cell development in mouse and human, a low-frequency population of immature CD4-CD8- double-negative (DN) thymocytes expresses early, mature αß T cell antigen receptor (TCR). We report that these early αß TCR+ DN (EADN) cells are DN3b-DN4 stage and require CD3δ but not major histocompatibility complex (MHC) for their generation/detection. When MHC - is present, however, EADN cells can respond to it, displaying a degree of coreceptor-independent MHC reactivity not typical of mature, conventional αß T cells. We found these data to be connected with observations that EADN cells were susceptible to T cell acute lymphoblastic leukemia (T-ALL) transformation in both humans and mice. Using the OT-1 TCR transgenic system to model EADN-stage αß TCR expression, we found that EADN leukemogenesis required MHC to induce development of T-ALL bearing NOTCH1 mutations. This leukemia-driving MHC requirement could be lost, however, upon passaging the tumors in vivo, even when matching MHC was continuously present in recipient animals and on the tumor cells themselves. These data demonstrate that MHC:TCR signaling can be required to initiate a cancer phenotype from an understudied developmental state that appears to be represented in the mouse and human disease spectrum.

Heterozygous variants in PRPF8 are associated with neurodevelopmental disorders.

The pre-mRNA-processing factor 8, encoded by PRPF8, is a scaffolding component of a spliceosome complex involved in the removal of introns from mRNA precursors. Previously, heterozygous pathogenic variants in PRPF8 have been associated with autosomal dominant retinitis pigmentosa. More recently, PRPF8 was suggested as a candidate gene for autism spectrum disorder due to the enrichment of sequence variants in this gene in individuals with neurodevelopmental disorders. We report 14 individuals with various forms of neurodevelopmental conditions, found to have heterozygous, predominantly de novo, missense, and loss-of-function variants in PRPF8. These individuals have clinical features that may represent a new neurodevelopmental syndrome.

Clonal Hematopoiesis and Myeloid Neoplasms in the Context of Telomere Biology Disorders.

PURPOSE OF REVIEW: Telomere biology disorders (TBDs) are cancer-predisposing multisystemic diseases that portend a higher risk of transforming into myeloid neoplasms (MNs). Due to the rarity and high variability of clinical presentations, TBD-specific characteristics of MN and the mechanisms behind this predisposition are not well defined. Herein, we review recent studies on TBD patient cohorts describing myeloid transformation events and summarize efforts to develop screening and treatment guidelines for these patients. RECENT FINDINGS: Preliminary studies have indicated that TBD patients have a higher prevalence of somatic genetic alterations in hematopoietic cells, an age-related phenomenon, also known as clonal hematopoiesis; increasing predisposition to MN. The CH mutational landscape in TBD differs from that observed in non-TBD patients and preliminary data suggest a higher frequency of somatic mutations in the DNA repair mechanism pathway. Although initial studies did not observe specific features of MN in TBD patients, certain events are common in TBD, such as hypocellular bone marrows. The mechanisms of MN development need further elucidation. Current management options for MN-TBD patients need to be individualized and tailored as per the clinical context. Because of the high sensitivity to alkylator chemotherapy and radiation conferred by short telomeres, non-cytotoxic targeted therapies and immunotherapy are ideal therapeutic options, but these therapies are still being tested in clinical trials. Defining the mechanisms of CH evolution in TBD and identifying risk factors leading to MN evolution will allow for the development of screening and treatment guidelines for these patients.

Telogator: a method for reporting chromosome-specific telomere lengths from long reads.

MOTIVATION: Telomeres are the repetitive sequences found at the ends of eukaryotic chromosomes and are often thought of as a 'biological clock,' with their average length shortening during division in most cells. In addition to their association with senescence, abnormal telomere lengths are well known to be associated with multiple cancers, short telomere syndromes and as risk factors for a broad range of diseases. While a majority of methods for measuring telomere length will report average lengths across all chromosomes, it is known that aberrations in specific chromosome arms are biomarkers for certain diseases. Due to their repetitive nature, characterizing telomeres at this resolution is prohibitive for short read sequencing approaches, and is challenging still even with longer reads. RESULTS: We present Telogator: a method for reporting chromosome-specific telomere length from long read sequencing data. We demonstrate Telogator's sensitivity in detecting chromosome-specific telomere length in simulated data across a range of read lengths and error rates. Telogator is then applied to 10 germline samples, yielding a high correlation with short read methods in reporting average telomere length. In addition, we investigate common subtelomere rearrangements and identify the minimum read length required to anchor telomere/subtelomere boundaries in samples with these haplotypes. AVAILABILITY AND IMPLEMENTATION: Telogator is written in Python3 and is available at github.com/zstephens/telogator. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Inverse Association of Telomere Length With Liver Disease and Mortality in the US Population.

Physiologic aging leads to attrition of telomeres and replicative senescence. An acceleration of this process has been hypothesized in the progression of chronic liver disease. We sought to examine the association of telomere length (TL) with liver disease and its impact on mortality risk. A cohort of 7,072 adults with leukocyte TL measurements from the National Health and Nutrition Examination Survey 1999-2002 with mortality follow-up through 2015 was analyzed. Liver disease was defined by aminotransferase levels and classified into etiology-based and advanced fibrosis categories. Multivariable-adjusted linear regression models estimated effect sizes, with 95% confidence intervals (CIs), of the presence of liver disease on TL. Cox regression models evaluated associations between TL and all-cause mortality risk using adjusted hazard ratios (HRs). The cohort was representative of the US population with mean age 46.1 years and mean TL 5.79 kilobase pairs. No overall association between TL and liver disease was found; however, there was a significant negative association of TL and advanced liver fibrosis in individuals aged 65 and above. The liver disease cohort (HR 1.22, 95% CI 0.99-1.51) and those with metabolic syndrome (HR 1.26, 95% CI 0.96-1.67) had increased mortality risk with shorter TL. The relationship between TL and all-cause mortality was stronger in women (HR 1.51, 95% CI 1.02-2.23) and in non-Hispanic Whites (HR 1.37, 95% CI 1.02-1.84). Conclusion: Shortened leukocyte TL is independently associated with advanced liver disease at older ages, and with a higher risk of all-cause mortality in those with liver disease. These associations reaffirm the need to better understand the role of telomeres in the progression of liver disease.

Clinicoradiographic and genetic features of cerebral small vessel disease indicate variability in mode of inheritance for monoallelic HTRA1 variants.

BACKGROUND: Biallelic pathogenic variants in HTRA1 cause CARASIL. More recently, monoallelic variants have been associated with the autosomal dominant disorder CADASIL2 but not all carriers develop disease manifestations. We describe the clinicoradiologic and mutation spectrum of four new CADASIL2 individuals. METHODS: Medical records at Mayo Clinic between 2013 and 2020 were retrospectively reviewed to identify patients with cerebral small vessel disease related to monoallelic HTRA1 variants. RESULTS: Four patients met the study inclusion criteria for cerebral small vessel disease related to HTRA1 monoallelic variants. The mean age at onset of first clinical stroke was 51.25 years (range 41-64 years). The mean disease duration was 6.5 years (range 4-12). All individuals had recurrent strokes within the duration of follow-up with a mean number of strokes per patient being 5.5 (range 2-12). Three individuals had leukoencephalopathy with brain stem involvement. Microhemorrhages were seen on brain MRI in three patients. HTRA1 monoallelic variants identified in our cohort were missense variants in three patients and a novel frameshift variation in one patient. Interestingly, two of these missense variants were previously reported in an autosomal recessive pattern of inheritance and here are associated with a dominant effect. CONCLUSIONS: Clinicoradiologic characteristics of heterozygous HTRA1-related CSVD may overlap with sporadic CSVD. Heterozygous HTRA1 variants can contribute to dominant or recessive disease mechanisms.

Spectrum of hematological malignancies, clonal evolution and outcomes in 144 Mayo Clinic patients with germline predisposition syndromes.

Germline predisposition syndromes (GPS) result from constitutional aberrations in tumor suppressive and homeostatic genes, increasing risk for neoplasia in affected kindred. In this study, we present clinical and genomic data on 144 Mayo Clinic patients with GPS; 59 evaluated prospectively using an algorithm-based diagnostic approach in the setting of a dedicated GPS/ inherited bone marrow failure syndrome (IBMFS) clinic. Seventy-two (50%) patients had IBMFS (telomere biology disorders-32,Fanconi anemia-18, Diamond Blackfan Anemia - 11, congenital neutropenia-5, Schwachman-Diamond Syndrome-5 and Bloom Syndrome-1), 27 (19%) had GPS with antecedent thrombocytopenia (RUNX1-FPD-15, ANKRD26-6, ETV6-2, GATA1-1, MPL-3), 28 (19%) had GPS without antecedent thrombocytopenia (GATA2 haploinsufficiency-16, DDX41-10, CBL-1 and CEBPA-1) and 17 (12%) had general cancer predisposition syndromes (ataxia telangiectasia-7, heterozygous ATM variants-3, CHEK2-2, TP53-2, CDK2NA-1, NF1-1 and Nijmegen Breakage Syndrome-1). Homozygous and heterozygous ATM pathogenic variants were exclusively associated with lymphoproliferative disorders (LPD), while DDX41 GPS was associated with LPD and myeloid neoplasms. The use of somatic NGS-testing identified clonal evolution in GPS patients, with ASXL1, RAS pathway genes, SRSF2 and TET2 being most frequently mutated. Fifty-two (91%) of 59 prospectively identified GPS patients had a change in their management approach, including additional GPS-related screening in 42 (71%), referral for allogenic HSCT workup and screening of related donors in 16 (27%), medication initiation and selection of specific conditioning regimens in 14 (24%), and genetic counseling with specific intent of fertility preservation and preconceptual counseling in 10 (17%) patients; highlighting the importance of dedicated GPS screening, detection and management programs for patients with hematological neoplasms.

Pathology, Radiology, and Genetics of Interstitial Lung Disease in Patients With Shortened Telomeres.

Interstitial lung diseases (ILDs) in patients with shortened telomeres have not been well characterized. We describe demographic, radiologic, histopathologic, and molecular features, and p16 expression in patients with telomeres ≤10th percentile (shortened telomeres) and compare them to patients with telomere length >10th percentile. Lung explants, wedge biopsies, and autopsy specimens of patients with telomere testing were reviewed independently by 3 pathologists using defined parameters. High-resolution computed tomography scans were reviewed by 3 radiologists. p16-positive fibroblast foci were quantified. A multidisciplinary diagnosis was recorded. Patients with shortened telomeres (N=26) were morphologically diagnosed as usual interstitial pneumonia (UIP) (N=11, 42.3%), chronic hypersensitivity pneumonitis (N=6, 23.1%), pleuroparenchymal fibroelastosis, fibrotic nonspecific interstitial pneumonia, desquamative interstitial pneumonia (N=1, 3.8%, each), and fibrotic interstitial lung disease (fILD), not otherwise specified (N=6, 23.1%). Patients with telomeres >10th percentile (N=18) showed morphologic features of UIP (N=9, 50%), chronic hypersensitivity pneumonitis (N=3, 16.7%), fibrotic nonspecific interstitial pneumonia (N=2, 11.1%), or fILD, not otherwise specified (N=4, 22.2%). Patients with shortened telomeres had more p16-positive foci (P=0.04). The number of p16-positive foci correlated with outcome (P=0.0067). Thirty-nine percent of patients with shortened telomeres harbored telomere-related gene variants. Among 17 patients with shortened telomeres and high-resolution computed tomography features consistent with or probable UIP, 8 (47.1%) patients showed morphologic features compatible with UIP; multidisciplinary diagnosis most commonly was idiopathic pulmonary fibrosis (N=7, 41.2%) and familial pulmonary fibrosis (N=5, 29%) in these patients. In conclusion, patients with shortened telomeres have a spectrum of fILDs. They often demonstrate atypical and discordant features on pathology and radiology leading to diagnostic challenges.

TSPEAR variants are primarily associated with ectodermal dysplasia and tooth agenesis but not hearing loss: A novel cohort study.

Biallelic loss-of-function variants in the thrombospondin-type laminin G domain and epilepsy-associated repeats (TSPEAR) gene have recently been associated with ectodermal dysplasia and hearing loss. The first reports describing a TSPEAR disease association identified this gene is a cause of nonsyndromic hearing loss, but subsequent reports involving additional affected families have questioned this evidence and suggested a stronger association with ectodermal dysplasia. To clarify genotype-phenotype associations for TSPEAR variants, we characterized 13 individuals with biallelic TSPEAR variants. Individuals underwent either exome sequencing or panel-based genetic testing. Nearly all of these newly reported individuals (11/13) have phenotypes that include tooth agenesis or ectodermal dysplasia, while three newly reported individuals have hearing loss. Of the individuals displaying hearing loss, all have additional variants in other hearing-loss-associated genes, specifically TMPRSS3, GJB2, and GJB6, that present competing candidates for their hearing loss phenotype. When presented alongside previous reports, the overall evidence supports the association of TSPEAR variants with ectodermal dysplasia and tooth agenesis features but creates significant doubt as to whether TSPEAR variants are a monogenic cause of hearing loss. Further functional evidence is needed to evaluate this phenotypic association.