Impact of genome build on RNA-seq interpretation and diagnostics.

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

The Undiagnosed Diseases Network: Characteristics of solvable applicants and diagnostic suggestions for non-accepted ones.

PURPOSE: Can certain characteristics identify as solvable some undiagnosed patients who seek extensive evaluation and thorough record review, like by the Undiagnosed Diseases Network (UDN)? METHODS: The UDN is a national research resource to solve medical mysteries through team science. Applicants provide informed consent to access to their medical records. After review, expert panels assess if applicants meet inclusion and exclusion criteria to select participants. When not accepting applicants, UDN experts may offer suggestions for diagnostic efforts. Using minimal information from initial applications, we compare features in applicants not accepted with those accepted and either solved or still not solved by the UDN. The diagnostic suggestions offered to non-accepted applicants and their clinicians were tallied. RESULTS: Non-accepted applicants were more often female, older at first symptoms and application, and longer in review than accepted applicants. The accepted and successfully diagnosed applicants were younger in ages, shorter in review time, more often non-white, of Hispanic ethnicity, and presenting with nervous system features. Half of non-accepted applicants were given suggestions for further local diagnostic evaluation. A few seemed to have two major diagnoses or a provocative environmental exposure history. CONCLUSION: Comprehensive UDN record review generates possibly helpful advice.

Loss of function of FAM177A1, a Golgi complex localized protein, causes a novel neurodevelopmental disorder.

PURPOSE: The function of FAM177A1 and its relationship to human disease is largely unknown. Recent studies have demonstrated FAM177A1 to be a critical immune-associated gene. One previous case study has linked FAM177A1 to a neurodevelopmental disorder in 4 siblings. METHODS: We identified 5 individuals from 3 unrelated families with biallelic variants in FAM177A1. The physiological function of FAM177A1 was studied in a zebrafish model organism and human cell lines with loss-of-function variants similar to the affected cohort. RESULTS: These individuals share a characteristic phenotype defined by macrocephaly, global developmental delay, intellectual disability, seizures, behavioral abnormalities, hypotonia, and gait disturbance. We show that FAM177A1 localizes to the Golgi complex in mammalian and zebrafish cells. Intersection of the RNA sequencing and metabolomic data sets from FAM177A1-deficient human fibroblasts and whole zebrafish larvae demonstrated dysregulation of pathways associated with apoptosis, inflammation, and negative regulation of cell proliferation. CONCLUSION: Our data shed light on the emerging function of FAM177A1 and defines FAM177A1-related neurodevelopmental disorder as a new clinical entity.

Digitized Electrocardiography Measurements Support the Biological Plausibility of the Pathological Significance of ST Segments in Athletes.

OBJECTIVE: ST segment deviations around the isoelectric line are common findings in manifest cardiovascular disease. In athletes, ST elevation is common, while ST depression is considered rare. However, clinical studies in athletes have associated ST depression with myocardial fibrosis and fatty infiltration and ST elevation with pericarditis and myocarditis. This study aims to explore the association between resting ST segment deviations and resting heart rate, an indicator of training and autonomic tone and electrocardiography (ECG) markers of exercise training effect and cardiovascular health R and T wave amplitude. DESIGN: Retrospective analysis of digitized ECG data. SETTING: Institutional setting. PARTICIPANTS: Seven thousand eight hundred thirty-six (male athletes = 4592, female athletes = 3244) healthy asymptomatic athletes (14-35 years). MAIN OUTCOME MEASURES: A series of correlations and regressions were conducted between ST depression (<0.0 µV) and ST elevation (>0.0 µV), on R and T wave amplitudes, and heart rate in leads V2, V5, and aVF. RESULTS: Positive correlations between ST elevation and R and T wave (S wave in V2) amplitudes and leads V5, V2, and aVF in male and female athletes (range of r = 0.1-0.54). In addition, there was a negative correlation between ST elevation and HR for male and female athletes. Finally, there was a negative correlation between ST depression and R wave and HR for male and female athletes in V5 ( P < 0.01). CONCLUSIONS: In athletes, ST segment elevation is correlated with R and T wave amplitudes and negatively correlated with HR. In addition, ST segment elevation is correlated with low heart rate, consistent with its higher prevalence in athletes. ST segment depression is not influenced by HR but is negatively associated with R and T wave amplitudes.

Quantifying assumptions underlying peak oxygen consumption equations across the body mass spectrum.

The goal of this study is to quantify the assumptions associated with the Wasserman-Hansen (WH) and Fitness Registry and the Importance of Exercise: A National Database (FRIEND) predictive peak oxygen consumption (pVO2) equations across body mass index (BMI). Assumptions in pVO2 for both equations were first determined using a simulation and then evaluated using exercise data from the Stanford Exercise Testing registry. We calculated percent-predicted VO2 (ppVO2) values for both equations and compared them using the Bland-Altman method. Assumptions associated with pVO2 across BMI categories were quantified by comparing the slopes of age-adjusted VO2 ratios (pVO2/pre-exercise VO2) and ppVO2 values for different BMI categories. The simulation revealed lower predicted fitness among adults with obesity using the FRIEND equation compared to the WH equations. In the clinical cohort, we evaluated 2471 patients (56.9% male, 22% with BMI >30 kg/m2, pVO2 26.8 mlO2/kg/min). The Bland-Altman plot revealed an average relative difference of -1.7% (95% CI: -2.1 to -1.2%) between WH and FRIEND ppVO2 values with greater differences among those with obesity. Analysis of the VO2 ratio to ppVO2 slopes across the BMI spectrum confirmed the assumption of lower fitness in those with obesity, and this trend was more pronounced using the FRIEND equation. Peak VO2 estimations between the WH and FRIEND equations differed significantly among individuals with obesity. The FRIEND equation resulted in a greater attributable reduction in pVO2 associated with obesity relative to the WH equations. The outlined relationships between BMI and predicted VO2 may better inform the clinical interpretation of ppVO2 values during cardiopulmonary exercise test evaluations.

Knowledge and attitudes on implementing cardiovascular pharmacogenomic testing.

Pharmacogenomics has the potential to inform drug dosing and selection, reduce adverse events, and improve medication efficacy; however, provider knowledge of pharmacogenomic testing varies across provider types and specialties. Given that many actionable pharmacogenomic genes are implicated in cardiovascular medication response variability, this study aimed to evaluate cardiology providers' knowledge and attitudes on implementing clinical pharmacogenomic testing. Sixty-one providers responded to an online survey, including pharmacists (46%), physicians (31%), genetic counselors (15%), and nurses (8%). Most respondents (94%) reported previous genetics education; however, only 52% felt their genetics education prepared them to order a clinical pharmacogenomic test. In addition, most respondents (66%) were familiar with pharmacogenomics, with genetic counselors being most likely to be familiar (p < 0.001). Only 15% of respondents had previously ordered a clinical pharmacogenomic test and a total of 36% indicated they are likely to order a pharmacogenomic test in the future; however, the vast majority of respondents (89%) were interested in pharmacogenomic testing being incorporated into diagnostic cardiovascular genetic tests. Moreover, 84% of providers preferred pharmacogenomic panel testing compared to 16% who preferred single gene testing. Half of the providers reported being comfortable discussing pharmacogenomic results with their patients, but the majority (60%) expressed discomfort with the logistics of test ordering. Reported barriers to implementation included uncertainty about the clinical utility and difficulty choosing an appropriate test. Taken together, cardiology providers have moderate familiarity with pharmacogenomics and limited experience with test ordering; however, they are interested in incorporating pharmacogenomics into diagnostic genetic tests and ordering pharmacogenomic panels.

Regional Variation in Cardiovascular Genes Enables a Tractable Genome Editing Strategy.

BACKGROUND: To realize the potential of genome engineering therapeutics, tractable strategies must be identified that balance personalized therapy with the need for off-the-shelf availability. We hypothesized that regional clustering of pathogenic variants can inform the design of rational prime editing therapeutics to treat the majority of genetic cardiovascular diseases with a limited number of reagents. METHODS: We collated 2435 high-confidence pathogenic/likely pathogenic (P/LP) variants in 82 cardiovascular disease genes from ClinVar. We assessed the regional density of these variants by defining a regional clustering index. We then combined a highly active base editor with prime editing to demonstrate the feasibility of a P/LP hotspot-directed genome engineering therapeutic strategy in vitro. RESULTS: P/LP variants in cardiovascular disease genes display higher regional density than rare variants found in the general population. P/LP missense variants displayed higher average regional density than P/LP truncating variants. Following hypermutagenesis at a pathogenic hotspot, mean prime editing efficiency across introduced variants was 57±27%. CONCLUSIONS: Designing therapeutics that target pathogenic hotspots will not only address known missense P/LP variants but also novel P/LP variants identified in these hotspots as well. Moreover, the clustering of P/LP missense rather than truncating variants in these hotspots suggests that prime editing technology is particularly valuable for dominant negative disease. Although prime editing technology in relation to cardiac health continues to improve, this study presents an approach to targeting the most impactful regions of the genome for inherited cardiovascular disease.

Impact of genome build on RNA-seq interpretation and diagnostics.

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network (UDN) and Genomics Research to Elucidate the Genetics of Rare Disease (GREGoR) Consortium. We identified 2,800 genes with build-dependent quantification across six routinely-collected biospecimens, including 1,391 protein-coding genes and 341 known rare disease genes. We further observed multiple genes that only have detectable expression in a subset of genome builds. Finally, we characterized how genome build impacts the detection of outlier transcriptomic events. Combined, we provide a database of genes impacted by build choice, and recommend that transcriptomics-guided analyses and diagnoses are cross-referenced with these data for robustness.

Integration of transcriptomics and long-read genomics prioritizes structural variants in rare disease.

Rare structural variants (SVs) - insertions, deletions, and complex rearrangements - can cause Mendelian disease, yet they remain difficult to accurately detect and interpret. We sequenced and analyzed Oxford Nanopore long-read genomes of 68 individuals from the Undiagnosed Disease Network (UDN) with no previously identified diagnostic mutations from short-read sequencing. Using our optimized SV detection pipelines and 571 control long-read genomes, we detected 716 long-read rare (MAF < 0.01) SV alleles per genome on average, achieving a 2.4x increase from short-reads. To characterize the functional effects of rare SVs, we assessed their relationship with gene expression from blood or fibroblasts from the same individuals, and found that rare SVs overlapping enhancers were enriched (LOR = 0.46) near expression outliers. We also evaluated tandem repeat expansions (TREs) and found 14 rare TREs per genome; notably these TREs were also enriched near overexpression outliers. To prioritize candidate functional SVs, we developed Watershed-SV, a probabilistic model that integrates expression data with SV-specific genomic annotations, which significantly outperforms baseline models that don't incorporate expression data. Watershed-SV identified a median of eight high-confidence functional SVs per UDN genome. Notably, this included compound heterozygous deletions in FAM177A1 shared by two siblings, which were likely causal for a rare neurodevelopmental disorder. Our observations demonstrate the promise of integrating long-read sequencing with gene expression towards improving the prioritization of functional SVs and TREs in rare disease patients.

Mavacamten Treatment for Symptomatic Obstructive Hypertrophic Cardiomyopathy: Interim Results From the MAVA-LTE Study, EXPLORER-LTE Cohort.

BACKGROUND: Data assessing the long-term safety and efficacy of mavacamten treatment for symptomatic obstructive hypertrophic cardiomyopathy are needed. OBJECTIVES: The authors sought to evaluate interim results from the EXPLORER-Long Term Extension (LTE) cohort of MAVA-LTE (A Long-Term Safety Extension Study of Mavacamten in Adults Who Have Completed EXPLORER-HCM; NCT03723655). METHODS: After mavacamten or placebo withdrawal at the end of the parent EXPLORER-HCM (Clinical Study to Evaluate Mavacamten [MYK-461] in Adults With Symptomatic Obstructive Hypertrophic Cardiomyopathy; NCT03470545), patients could enroll in MAVA-LTE. Patients received mavacamten 5 mg once daily; adjustments were made based on site-read echocardiograms. RESULTS: Between April 9, 2019, and March 5, 2021, 231 of 244 eligible patients (94.7%) enrolled in MAVA-LTE (mean age: 60 years; 39% female). At data cutoff (August 31, 2021) 217 (93.9%) remained on treatment (median time in study: 62.3 weeks; range: 0.3-123.9 weeks). At 48 weeks, patients showed improvements in left ventricular outflow tract (LVOT) gradients (mean change ± SD from baseline: resting: -35.6 ± 32.6 mm Hg; Valsalva: -45.3 ± 35.9 mm Hg), N-terminal pro-B-type natriuretic peptide (NT-proBNP) levels (median: -480 ng/L; Q1-Q3: -1,104 to -179 ng/L), and NYHA functional class (67.5% improved by ≥1 class). LVOT gradients and NT-proBNP reductions were sustained through 84 weeks in patients who reached this timepoint. Over 315 patient-years of exposure, 8 patients experienced an adverse event of cardiac failure, and 21 patients had an adverse event of atrial fibrillation, including 11 with no prior history of atrial fibrillation. Twelve patients (5.2%) developed transient reductions in site-read echocardiogram left ventricular ejection fraction of <50%, resulting in temporary treatment interruption; all recovered. Ten patients discontinued treatment due to treatment-emergent adverse events. CONCLUSIONS: Mavacamten treatment showed clinically important and durable improvements in LVOT gradients, NT-proBNP levels, and NYHA functional class, consistent with EXPLORER-HCM. Mavacamten treatment was well tolerated over a median 62-week follow-up.

Recurring homozygous ACTN2 variant (p.Arg506Gly) causes a recessive myopathy.

OBJECTIVE: ACTN2, encoding alpha-actinin-2, is essential for cardiac and skeletal muscle sarcomeric function. ACTN2 variants are a known cause of cardiomyopathy without skeletal muscle involvement. Recently, specific dominant monoallelic variants were reported as a rare cause of core myopathy of variable clinical onset, although the pathomechanism remains to be elucidated. The possibility of a recessively inherited ACTN2-myopathy has also been proposed in a single series. METHODS: We provide clinical, imaging, and histological characterization of a series of patients with a novel biallelic ACTN2 variant. RESULTS: We report seven patients from five families with a recurring biallelic variant in ACTN2: c.1516A>G (p.Arg506Gly), all manifesting with a consistent phenotype of asymmetric, progressive, proximal, and distal lower extremity predominant muscle weakness. None of the patients have cardiomyopathy or respiratory insufficiency. Notably, all patients report Palestinian ethnicity, suggesting a possible founder ACTN2 variant, which was confirmed through haplotype analysis in two families. Muscle biopsies reveal an underlying myopathic process with disruption of the intermyofibrillar architecture, Type I fiber predominance and atrophy. MRI of the lower extremities demonstrate a distinct pattern of asymmetric muscle involvement with selective involvement of the hamstrings and adductors in the thigh, and anterior tibial group and soleus in the lower leg. Using an in vitro splicing assay, we show that c.1516A>G ACTN2 does not impair normal splicing. INTERPRETATION: This series further establishes ACTN2 as a muscle disease gene, now also including variants with a recessive inheritance mode, and expands the clinical spectrum of actinopathies to adult-onset progressive muscle disease.

Large-scale mutational analysis identifies UNC93B1 variants that drive TLR-mediated autoimmunity in mice and humans.

Nucleic acid-sensing Toll-like receptors (TLR) 3, 7/8, and 9 are key innate immune sensors whose activities must be tightly regulated to prevent systemic autoimmune or autoinflammatory disease or virus-associated immunopathology. Here, we report a systematic scanning-alanine mutagenesis screen of all cytosolic and luminal residues of the TLR chaperone protein UNC93B1, which identified both negative and positive regulatory regions affecting TLR3, TLR7, and TLR9 responses. We subsequently identified two families harboring heterozygous coding mutations in UNC93B1, UNC93B1+/T93I and UNC93B1+/R336C, both in key negative regulatory regions identified in our screen. These patients presented with cutaneous tumid lupus and juvenile idiopathic arthritis plus neuroinflammatory disease, respectively. Disruption of UNC93B1-mediated regulation by these mutations led to enhanced TLR7/8 responses, and both variants resulted in systemic autoimmune or inflammatory disease when introduced into mice via genome editing. Altogether, our results implicate the UNC93B1-TLR7/8 axis in human monogenic autoimmune diseases and provide a functional resource to assess the impact of yet-to-be-reported UNC93B1 mutations.

Quality of Life and Exercise Capacity in Early Stage and Subclinical Hypertrophic Cardiomyopathy: A Secondary Analysis of the VANISH Trial.

BACKGROUND: The health-related quality of life (HRQOL) and cardiopulmonary exercise testing (CPET) performance of individuals with subclinical and early stage hypertrophic cardiomyopathy (HCM) have not been systematically studied. Improved understanding will inform the natural history of HCM and factors influencing well-being. METHODS: VANISH trial (Valsartan for Attenuating Disease Evolution in Early Sarcomeric HCM) participants with early stage sarcomeric HCM (primary analysis cohort) and subclinical HCM (sarcomere variant without left ventricular hypertrophy comprising the exploratory cohort) who completed baseline and year 2 HRQOL assessment via the pediatric quality of life inventory and CPET were studied. Metrics correlating with baseline HRQOL and CPET performance were identified. The impact of valsartan treatment on these measures was analyzed in the early stage cohort. RESULTS: Two hundred participants were included: 166 with early stage HCM (mean age, 23±10 years; 40% female; 97% White; and 92% New York Heart Association class I) and 34 subclinical sarcomere variant carriers (mean age, 16±5 years; 50% female; and 100% White). Baseline HRQOL was good in both cohorts, although slightly better in subclinical HCM (composite pediatric quality of life score 84.6±10.6 versus 90.2±9.8; P=0.005). Both cohorts demonstrated mildly reduced functional status (mean percent predicted peak oxygen uptake 73±16 versus 78±12 mL/kg per minute; P=0.18). Percent predicted peak oxygen uptake and peak oxygen pulse correlated with HRQOL. Valsartan improved physical HRQOL in early stage HCM (adjusted mean change in pediatric quality of life score +4.1 versus placebo; P=0.01) but did not significantly impact CPET performance. CONCLUSIONS: Functional capacity can be impaired in young, healthy people with early stage HCM, despite New York Heart Association class I status and good HRQOL. Peak oxygen uptake was similarly decreased in subclinical HCM despite normal left ventricular wall thickness and excellent HRQOL. Valsartan improved physical pediatric quality of life scores but did not significantly impact CPET performance. Further studies are needed for validation and to understand how to improve patient experience. REGISTRATION: URL: https://www.clinicaltrials.gov; Unique identifier: NCT01912534.

The mitochondrial multi-omic response to exercise training across rat tissues.

Mitochondria have diverse functions critical to whole-body metabolic homeostasis. Endurance training alters mitochondrial activity, but systematic characterization of these adaptations is lacking. Here, the Molecular Transducers of Physical Activity Consortium mapped the temporal, multi-omic changes in mitochondrial analytes across 19 tissues in male and female rats trained for 1, 2, 4, or 8 weeks. Training elicited substantial changes in the adrenal gland, brown adipose, colon, heart, and skeletal muscle. The colon showed non-linear response dynamics, whereas mitochondrial pathways were downregulated in brown adipose and adrenal tissues. Protein acetylation increased in the liver, with a shift in lipid metabolism, whereas oxidative proteins increased in striated muscles. Exercise-upregulated networks were downregulated in human diabetes and cirrhosis. Knockdown of the central network protein 17-beta-hydroxysteroid dehydrogenase 10 (HSD17B10) elevated oxygen consumption, indicative of metabolic stress. We provide a multi-omic, multi-tissue, temporal atlas of the mitochondrial response to exercise training and identify candidates linked to mitochondrial dysfunction.

Assessing the Assisted Six-Minute Cycling Test as a Measure of Endurance in Non-Ambulatory Patients with Spinal Muscular Atrophy (SMA).

Assessing endurance in non-ambulatory individuals with Spinal Muscular Atrophy (SMA) has been challenging due to limited evaluation tools. The Assisted 6-Minute Cycling Test (A6MCT) is an upper limb ergometer assessment used in other neurologic disorders to measure endurance. To study the performance of the A6MCT in the non-ambulatory SMA population, prospective data was collected on 38 individuals with SMA (13 sitters; 25 non-sitters), aged 5 to 74 years (mean = 30.3; SD = 14.1). The clinical measures used were A6MCT, Revised Upper Limb Module (RULM), Adapted Test of Neuromuscular Disorders (ATEND), and Egen Klassifikation Scale 2 (EK2). Perceived fatigue was assessed using the Fatigue Severity Scale (FSS), and effort was assessed using the Rate of Perceived Exertion (RPE). Data were analyzed for: (1) Feasibility, (2) Clinical discrimination, and (3) Associations between A6MCT with clinical characteristics and outcomes. Results showed the A6MCT was feasible for 95% of the tested subjects, discriminated between functional groups (p = 0.0086), and was significantly associated with results obtained from RULM, ATEND, EK2, and Brooke (p < 0.0001; p = 0.029; p < 0.001; p = 0.005). These findings indicate the A6MCT's potential to evaluate muscular endurance in non-ambulatory SMA individuals, complementing clinician-rated assessments. Nevertheless, further validation with a larger dataset is needed for broader application.