Background and Aims: Metabolic dysfunction–associated steatotic liver disease affects about a third of adults worldwide and is projected soon to be the leading cause of liver cirrhosis. It occurs when fat accumulates in hepatocytes and can progress to metabolic dysfunction–associated steatohepatitis, liver cirrhosis, and HCC. Metabolic dysfunction–associated steatotic liver disease pathogenesis is believed to involve a combination of genetic and environmental risk factors. Single nucleotide polymorphisms have been implicated, but non-syndromic monogenic causes are lacking. Approach and Results: We identified a novel genetic variant in a familial case of metabolic dysfunction–associated steatohepatitis and performed deep variant functional analysis, including protein modeling, dynamics, and cell-based assays to assess molecular mechanisms of dysfunction and altered cellular signaling. We analyzed exome sequencing data of 3904 individuals with steatotic liver disease (SLD) to identify additional cases and establish the link between specific gene variants and SLD diagnosis. We discovered and functionally validated the NM_000245.4:c.3505A>T; p.(Ile1169Phe) variant in the MET (mesenchymal-epithelial transition) kinase domain as a monogenic cause of SLD. Subsequently, we detected additional ultra-rare, previously uninterpreted, and likely deleterious variants in MET from screening sequencing data. Among individuals with confirmed SLD based on electronic record review, 1.1% (45/3904) had rare predicted deleterious MET variants. Eight of 45 (17.7%) individuals had predicted deleterious variants in the MET kinase domain confirmed to be functionally like the familial case variant. Conclusions: We report the first germline nonmalignant rare MET -driven disease, a monogenic form of SLD.
Epigenetic disease markers in primary sclerosing cholangitis and primary biliary cholangitis—methylomics of cholestatic liver disease
Background: The epigenome, the set of modifications to DNA and associated molecules that control gene expression, cellular identity, and function, plays a major role in mediating cellular responses to outside factors. Thus, evaluation of the epigenetic state can provide insights into cellular adaptions occurring over the course of disease. Methods: We performed epigenome-wide association studies of primary sclerosing cholangitis (PSC) and primary biliary cholangitis (PBC) using the Illumina MethylationEPIC Bead Chip. Results: We found evidence of increased epigenetic age acceleration and differences in predicted immune cell composition in patients with PSC and PBC. Epigenetic profiles demonstrated differences in predicted protein levels including increased levels of tumor necrosis factor receptor superfamily member 1B in patients with cirrhotic compared to noncirrhotic PSC and PBC. Epigenome-wide association studies of PSC discovered strongly associated 5′-C-phosphate-G-3′ sites in genes including vacuole membrane protein 1 and SOCS3, and epigenome-wide association studies of PBC found strong 5′-C-phosphate-G-3′ associations in genes including NOD-like receptor family CARD domain containing 5, human leukocyte antigen-E, and PSMB8. Analyses identified disease-associated canonical pathways and upstream regulators involved with immune signaling and activation of macrophages and T-cells. A comparison of PSC and PBC data found relatively little overlap at the 5′-C-phosphate-G-3′ and gene levels with slightly more overlap at the level of pathways and upstream regulators. Conclusions: This study provides insights into methylation profiles of patients that support current concepts of disease mechanisms and provide novel data to inspire future research. Studies to corroborate our findings and expand into other -omics layers will be invaluable to further our understanding of these rare diseases with the goal to improve and individualize prognosis and treatment.
Exome Sequencing Identifies Carriers of the Autosomal Dominant Cancer Predisposition Disorders Beyond Current Practice Guideline Recommendations
PURPOSE The autosomal dominant cancer predisposition disorders hereditary breast and ovarian cancer (HBOC) and Lynch syndrome (LS) are genetic conditions for which early identification and intervention have a positive effect on the individual and public health. The goals of this study were to determine whether germline genetic screening using exome sequencing could be used to efficiently identify carriers of HBOC and LS. METHODS Participants were recruited from three geographically and racially diverse sites in the United States (Rochester, MN; Phoenix, AZ; Jacksonville, FL). Participants underwent Exome+ sequencing (Helix Inc, San Mateo, CA) and return of results for specific genetic findings: HBOC ( BRCA1 and BRCA1) and LS ( MLH1, MSH2, MSH6, PMS2, and EPCAM). Chart review was performed to collect demographics and personal and family cancer history. RESULTS To date, 44,306 participants have enrolled in Tapestry. Annotation and interpretation of all variants in genes for HBOC and LS resulted in the 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). More than half of these participants (52.1%) were newly diagnosed carriers with HBOC and LS. In all, 39.2% of HBOC/LS carriers did not satisfy National Comprehensive Cancer Network (NCCN) criteria for genetic evaluation. NCCN criteria were less commonly met in underrepresented minority populations versus self-reported White race (51.5% v 37.5%, P = .028). CONCLUSION Our results emphasize the need for wider utilization of germline genetic sequencing for enhanced screening and detection of individuals who have LS and HBOC cancer predisposition syndromes.
IL-17 signaling in primary sclerosing cholangitis patient-derived organoids
Background: The pathogenesis of primary sclerosing cholangitis (PSC) is unclear, although studies implicate IL-17A as an inflammatory mediator in this disease. However, a direct assessment of IL-17 signaling in PSC cholangiocytes is lacking. In this study, we aimed to investigate and characterize the response of PSC extrahepatic cholangiocyte organoids (ECO) to IL-17A stimulation. Methods: Cholangiocytes obtained from patients with PSC and without PSC by endoscopic retrograde cholangiography were cultured as ECO. The ECO were treated with vehicle or IL-17A and assessed by transcriptomics, secretome analysis, and genome sequencing. Results: Unsupervised clustering of all integrated single-cell RNA sequencing data identified 8 cholangiocyte clusters that did not differ between PSC and non-PSC ECO. However, PSC ECO cells demonstrated a robust response to IL-17 treatment, as noted by an increased number of differentially expressed genes by transcriptomics and more abundant chemokine and cytokine expression and secretion. After rigorous filtering, genome sequencing identified candidate somatic variants shared among PSC ECO from unrelated individuals. However, no candidate rare variants in genes regulating the IL-17 pathway were identified, but rare variants regulating the MAPK signaling pathway were present in all PSC ECO. Conclusions: PSC and non-PSC patient-derived ECO respond differently to IL-17 stimulation, implicating this pathway in the pathogenesis of PSC.
Clinomics Implementation in the Mayo Clinic Nephrology Practice
AbstractBackgroundHow variations predicted by pharmacogenomic testing to alter drug metabolism and therapeutic response affect outcomes for patients with disorders of gut‐ brain interaction is unclear.AimsTo assess the prevalence of pharmacogenomics‐predicted drug‐gene interactions and symptom outcomes for patients with disorders of gut‐brain interaction.MethodsPatients who were treated in our clinical practice for functional dyspepsia/bowel disorder underwent pharmacogenomic testing. The change in symptoms from baseline to 6 months was compared for patients with variations in CYP2D6 and CYP2C19, which metabolize neuromodulators, and SLC6A4, which encodes the sodium‐ dependent serotonin transporter.ResultsAt baseline, 79 of 94 participants (84%) had at least one predicted major drug‐ gene interaction, and all 94 (100%) had at least one predicted moderate interaction. For the 44 participants who completed a survey of their symptoms at 6 months, the mean (SD) irritable bowel syndrome‐symptom severity score decreased from 284 (71) at baseline to 231 (95) at 6 months (p < 0.001). Among patients taking selective serotonin reuptake inhibitors, the decrease in symptom severity (p = 0.03) and pain (p = 0.002) scores from baseline to 6 months was greater for patients with a homozygous SLC6A4 long/long genotype (n = 30) (ie, increased serotonin transporter activity) than for patients with homozygous short/short or heterozygous long/short genotypes (n = 64). Symptom outcomes were not affected by CYP2D6 or CYP2C19 variations.ConclusionsThe homozygous SLC6A4 long/long genotype confers better symptom resolution for patients with disorders of gut‐brain interaction who take selective serotonin reuptake inhibitors than do the homozygous short/short or heterozygous long/short genotypes.
Distribution of Chromosome-Arm Specific Telomere Length in Patients with Telomere Biology Disorders
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 <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 <5 kb and <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>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 <5kb and <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.
Functional Testing of Variants of Uncertain Significance in TERC, TERT,and RTEL1 from Adult Patients with Telomere Biology Disorders
INTRODUCTION Telomere biology disorders (TBDs) are multisystem diseases caused by accelerated telomere shortening that can present with bone marrow failure and pulmonary fibrosis, among other abnormalities. TBD diagnosis is considered when average blood telomere length (TL) measured by clinical flowFISH is <1 st percentile and/or a pathogenic variant in telomere maintenance genes (most commonly TERT, TERC, RTEL1) is identified. Pathogenic variants in genes coding for the telomerase protein ( TERT) and RNA component ( TERC) decrease telomerase activity. Similarly, variants in RTEL1 result in defective unwinding of the t-loop at telomere ends, increasing the levels of free t-circles in cells. Variants in these genes are often classified as variants of uncertain significance (VUSs by ACMG criteria) with additional efforts needed for variant curation. To address this, we describe results of functional testing performed in TBD patients with VUSs involving TERT, TERC, and RTEL1. METHODS TBD suspected patients were identified through the Pre-Myeloid Cancer and Bone Marrow Failure Clinic (Mayo Clinic) where they underwent clinical flowFISH and genetic testing for TBD related genes. Informed consent was obtained with approval from the Mayo Clinic Institutional Review Board. Viable peripheral blood and bone marrow samples were collected along with age matched healthy volunteers. Telomerase activity was measured in samples with TERT and TERC variants through telomeric repeat amplification protocol (TRAP) assays using the TRAPeze XL Telomerase Detection Kit (MilliporeSigma) following manufacturer instructions. HCT116 cells were used as positive controls while heat inactivated samples served as negative controls. Fluorescence was measured in a SpectraMax plate reader. Percentage of total product generated (TPG) relative to healthy age matched controls was calculated. T-circle detection assays were used on RTEL1 samples based on the protocol by Zhang et al. (2017). Final southern blotting was completed using the TeloTAGGG Telomere Length Assay Kit (Roche), and t-circle bands were visualized in a c600 GEL Imaging System (Azure Biosystems). RESULTS Ten patients (mean age = 53 years (27-71), 50% female) were included in our study. Patients 1 to 5 presented VUS in TERT (n=4) or TERC (n=1) and TL <1 st percentile in lymphocytes and/or granulocytes. Patient 6 carried a known pathogenic variant in RTEL1 and was included as a positive control. Patients 7 to 9 with RTEL1 variants presented TL close to the 10 th or between the 1 st-10 th percentile. No TL results were available for patient 10 but was tested due to the presence of a VUS in RTEL1 (Table 1). Similarly, while patient 7 carried a likely pathogenic variant, it was tested due to presenting lymphocyte TL between the 1 st-10 th percentile. TRAP results indicated decreased telomerase activity compared to age matched controls in all samples with TERT and TERC variants suggesting a pathogenic effect (52.%, 10.7%, 48.7%, 4.4% and 5.2% of TPG compared to age matched controls, respectively) (Table 1). On the other hand, t-circles were only seen in patient 6 (Figure 2). DISCUSSION We used TRAP and t-circle assays to explore the telomerase activity and presence of t-circles in patient samples with TERT/TERC and RTEL1 variants. In the case of TERT/TERC, we observed a reduction of telomerase activity in all samples. Together with blood TL below the 1 st percentile, this suggests that these variants are indeed deleterious (Table 1). However, since patient 4 presented two TERT VUSs, additional testing is needed to conclude how each variant impacts telomerase activity independently. In regard to RTEL1 variants, we observed t-circles in patient 6 as expected. No t-circles could be seen on all remaining samples, including patient 7 who carried a variant classified as likely pathogenic (Figure 1). This in-frame deletion has been described once in a TBD family (Cogan et al. 2015) but was not functionally tested. This result, together with the borderline TL seen in this patient, may suggest that some compensatory mechanism exists on this sample, or that the deleterious effect of the variant is not strong enough to be detected by our assay. Further studies are currently ongoing by our group to validate these results using a cell-based transfection disease model. In summary, we show that functional testing can be a helpful tool to determine pathogenicity of VUSs in TERT/TERC and RTEL1.
IL-17 Signaling in Primary Sclerosing Cholangitis Patient-Derived Organoids
Abstract The pathogenesis of primary sclerosing cholangitis (PSC) is unclear, although studies implicate IL-17A as an inflammatory mediator in this disease. However, a direct assessment of IL-17 signaling in PSC cholangiocytes is lacking. In this study we aimed to investigate the response of PSC extrahepatic cholangiocyte organoids (ECO) to IL-17A stimulation. Cholangiocytes obtained from PSC and non-PSC patients by endoscopic retrograde cholangiography (ERC) were cultured as ECO. The ECO were treated with vehicle or IL-17A and assessed by transcriptomics, secretome analysis, and genome sequencing (GS). Unsupervised clustering of all integrated scRNA-seq data identified 8 cholangiocyte clusters which did not differ between PSC and non-PSC ECO. However, PSC ECO cells demonstrated a robust response to IL-17 treatment, noted by an increased number of differentially expressed genes (DEG) by transcriptomics, and more abundant chemokine and cytokine expression and secretion. After rigorous filtering, GS identified candidate somatic variants shared among PSC ECO from unrelated individuals. However, no candidate rare variants in genes regulating the IL-17 pathway were identified, but rare variants regulating the MAPK signaling pathway were present in all PSC ECO. In conclusion, PSC and non-PSC patient derived ECO respond differently to IL-17 stimulation implicating this pathway in the pathogenesis of PSC.