Genetic sex validation for sample tracking in next-generation sequencing clinical testing.

OBJECTIVE: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. RESULTS: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors (49.09%), samples from transgender participants (3.64%) and stem cell or bone marrow transplant patients (7.27%) along with undetermined sample mix-ups (40%) for which sample swaps occurred prior to arrival at genome centers, however the exact cause of the events at the sampling sites resulting in the mix-ups were not able to be determined.

Evaluation of Malignant Hyperthermia Features in Patients with Pathogenic or Likely Pathogenic RYR1 Variants Disclosed through a Population Genomic Screening Program.

BACKGROUND: Malignant hyperthermia (MH) susceptibility is a heritable musculoskeletal disorder that can present as a potentially fatal hypermetabolic response to triggering anesthesia agents. Genomic screening for variants in MH-associated genes RYR1 and CACNA1S provides an opportunity to prevent morbidity and mortality. There are limited outcomes data from disclosing variants in RYR1, the most common MH susceptibility gene, in unselected populations. The authors sought to identify the rate of MH features or fulminant episodes after triggering agent exposure in an unselected population undergoing genomic screening including actionable RYR1 variants. METHODS: The MyCode Community Health Initiative by Geisinger (USA) is an electronic health record-linked biobank that discloses pathogenic and likely pathogenic variants in clinically actionable genes to patient-participants. Available electronic anesthesia and ambulatory records for participants with actionable RYR1 results returned through December 2020 were evaluated for pertinent findings via double-coded chart reviews and reconciliation. Descriptive statistics for observed phenotypes were calculated. RESULTS: One hundred fifty-two participants had an actionable RYR1 variant disclosed during the study period. None had previous documented genetic testing for MH susceptibility; one had previous contracture testing diagnosing MH susceptibility. Sixty-eight participants (44.7%) had anesthesia records documenting triggering agent exposure during at least one procedure. None received dantrolene treatment or had documented muscle rigidity, myoglobinuria, hyperkalemia, elevated creatine kinase, severe myalgia, or tea-colored urine. Of 120 possibly MH-related findings (postoperative intensive care unit admissions, hyperthermia, arterial blood gas evaluation, hypercapnia, or tachycardia), 112 (93.3%) were deemed unlikely to be MH events; 8 (6.7%) had insufficient records to determine etiology. CONCLUSIONS: Results demonstrate a low frequency of classic intraanesthetic hypermetabolic phenotypes in an unselected population with actionable RYR1 variants. Further research on the actionability of screening for MH susceptibility in unselected populations, including economic impact, predictors of MH episodes, and expanded clinical phenotypes, is necessary.

Low adenoma burden in unselected patients with a pathogenic APC variant.

PURPOSE: Genomic screening can improve clinical outcomes, but presentation of individuals with risk for polyposis identified via genomic screening is unknown. To inform assessment of clinical utility of genomic screening for polyposis risk, clinical presentation of individuals in an unselected health care system cohort with an APC pathogenic or likely pathogenic (P/LP) variant causative of familial adenomatous polyposis are described. METHODS: Electronic health records of individuals with an APC P/LP variant identified via the MyCode program (MyCode APC+) were reviewed to assess adenoma burden and compare it among individuals with a clinical diagnosis of familial adenomatous polyposis and matched variant-negative controls. RESULTS: The prevalence of APC P/LP variants in this health care cohort is estimated to be 1 in 2800. Twenty-four MyCode APC+ individuals were identified during the study period. Median age at result disclosure was 53 years. Rate of clinical polyposis was 8%. Two of six participants with a classic region variant and none of those with an attenuated region variant had polyposis. MyCode APC+ participants did not differ from controls in cumulative adenoma count. CONCLUSION: APC P/LP variant prevalence estimate in the MyCode cohort is higher than prior published prevalence rates. Individuals with APC P/LP variants identified via genomic screening had a low adenoma burden.

Yield of Familial Hypercholesterolemia Genetic and Phenotypic Diagnoses After Electronic Health Record and Genomic Data Screening.

Background Data mining of electronic health records to identify patients suspected of familial hypercholesterolemia (FH) has been limited by absence of both phenotypic and genomic data in the same cohort. Methods and Results Using the Geisinger MyCode Community Health Initiative cohort (n=130 257), we ran 2 screening algorithms (Mayo Clinic [Mayo] and flag, identify, network, deliver [FIND] FH) to determine FH genetic and phenotypic diagnostic yields. With 29 243 excluded by Mayo (for secondary causes of hypercholesterolemia, no lipid value in electronic health records), 52 034 excluded by FIND FH (insufficient data to run the model), and 187 excluded for prior FH diagnosis, a final cohort of 59 729 participants was created. Genetic diagnosis was based on presence of a pathogenic or likely pathogenic variant in FH genes. Charts from 180 variant-negative participants (60 controls, 120 identified by FIND FH and Mayo) were reviewed to calculate Dutch Lipid Clinic Network scores; a score ≥5 defined probable phenotypic FH. Mayo flagged 10 415 subjects; 194 (1.9%) had a pathogenic or likely pathogenic FH variant. FIND FH flagged 573; 34 (5.9%) had a pathogenic or likely pathogenic variant, giving a net yield from both of 197 out of 280 (70%). Confirmation of a phenotypic diagnosis was constrained by lack of electronic health record data on physical findings or family history. Phenotypic FH by chart review was present by Mayo and/or FIND FH in 13 out of 120 versus 2 out of 60 not flagged by either (P<0.09). Conclusions Applying 2 recognized FH screening algorithms to the Geisinger MyCode Community Health Initiative identified 70% of those with a pathogenic or likely pathogenic FH variant. Phenotypic diagnosis was rarely achievable due to missing data.

Observational study of population genomic screening for variants associated with endocrine tumor syndromes in a large, healthcare-based cohort.

BACKGROUND: In current care, patients' personal and self-reported family histories are primarily used to determine whether genetic testing for hereditary endocrine tumor syndromes (ETS) is indicated. Population genomic screening for other conditions has increased ascertainment of individuals with pathogenic/likely pathogenic (P/LP) variants, leading to improved management and earlier diagnoses. It is unknown whether such benefits occur when screening broader populations for P/LP ETS variants. This manuscript assesses clinical utility outcomes of a large, unselected, healthcare-based genomic screening program by describing personal and family history of syndrome-related features, risk management behaviors after result disclosure, and rates of relevant post-disclosure diagnoses in patient-participants with P/LP ETS variants. METHODS: Observational study of individuals informed of a P/LP variant in MEN1, RET, SDHAF2, SDHB, SDHC, SDHD, or VHL through Geisinger's MyCode Community Health Initiative between June 2016 and October 2019. Electronic health records (EHRs) of participants were evaluated for a report of pre-disclosure personal and self-reported family histories and post-disclosure risk management and diagnoses. RESULTS: P/LP variants in genes of interest were identified in 199 of 130,490 (1 in 656) adult Geisinger MyCode patient-participants, 80 of which were disclosed during the study period. Eighty-one percent (n = 65) did not have prior evidence of the result in their EHR and, because they were identified via MyCode, were included in further analyses. Five participants identified via MyCode (8%) had a personal history of syndrome-related features; 16 (25%) had a positive self-reported family history. Time from result disclosure to EHR review was a median of 0.7 years. Post-disclosure, 36 (55.4%) completed a recommended risk management behavior; 11 (17%) were diagnosed with a syndrome-related neoplasm after completing a risk management intervention. CONCLUSIONS: Broader screening for pathogenic/likely pathogenic variants associated with endocrine tumor syndromes enables detection of at-risk individuals, leads to the uptake of risk management, and facilitates relevant diagnoses. Further research will be necessary to continue to determine the clinical utility of screening diverse, unselected populations for such variants.

Genomic Screening for Pathogenic Transthyretin Variants Finds Evidence of Underdiagnosed Amyloid Cardiomyopathy From Health Records.

BACKGROUND: New treatments for transthyretin amyloidosis improve survival, but diagnosis remains challenging. Pathogenic or likely pathogenic (P/LP) variants in the transthyretin (TTR) gene are one cause of transthyretin amyloidosis, and genomic screening has been proposed to identify at-risk individuals. However, data on disease features and penetrance are lacking to inform the utility of such population-based genomic screening for TTR. OBJECTIVES: This study characterized the prevalence of P/LP variants in TTR identified through exome sequencing and the burden of associated disease from electronic health records for individuals with these variants from a large (N = 134,753), primarily European-ancestry cohort. METHODS: We compared frequencies of common disease features and cardiac imaging findings between individuals with and without P/LP TTR variants. RESULTS: We identified 157 of 134,753 (0.12%) individuals with P/LP TTR variants (43% male, median age 52 [Q1-Q3: 37-61] years). Seven P/LP variants accounted for all observations, the majority being V122I (p.V142I; 113, 0.08%). Approximately 60% (n = 91) of individuals with P/LP TTR variants (all V122I) had African ancestry. Diagnoses of amyloidosis were limited (2 of 157 patients), although related heart disease diagnoses, including cardiomyopathy and heart failure, were significantly increased in individuals with P/LP TTR variants who were aged >60 years. Fourteen percent (7 of 49) of individuals aged ≥60 or older with a P/LP TTR variant had heart disease and ventricular septal thickness >1.2 cm, only one of whom was diagnosed with amyloidosis. CONCLUSIONS: Individuals with P/LP TTR variants identified by genomic screening have increased odds of heart disease after age 60 years, although amyloidosis is likely underdiagnosed without knowledge of the genetic variant.

The genetic architecture of Plakophilin 2 cardiomyopathy.

PURPOSE: The genetic architecture of Plakophilin 2 (PKP2) cardiomyopathy can inform our understanding of its variant pathogenicity and protein function. METHODS: We assess the gene-wide and regional association of truncating and missense variants in PKP2 with arrhythmogenic cardiomyopathy (ACM), and arrhythmogenic right ventricular cardiomyopathy (ARVC) specifically. A discovery data set compares genetic testing requisitions to gnomAD. Validation is performed in a rigorously phenotyped definite ARVC cohort and non-ACM individuals in the Geisinger MyCode cohort. RESULTS: The etiologic fraction (EF) of ACM-related diagnoses from truncating variants in PKP2 is significant (0.85 [0.80,0.88], p < 2 × 10-16), increases for ARVC specifically (EF = 0.96 [0.94,0.97], p < 2 × 10-16), and is highest in definite ARVC versus non-ACM individuals (EF = 1.00 [1.00,1.00], p < 2 × 10-16). Regions of missense variation enriched for ACM probands include known functional domains and the C-terminus, which was not previously known to contain a functional domain. No regional enrichment was identified for truncating variants. CONCLUSION: This multicohort evaluation of the genetic architecture of PKP2 demonstrates the specificity of PKP2 truncating variants for ARVC within the ACM disease spectrum. We identify the PKP2 C-terminus as a potential functional domain and find that truncating variants likely cause disease irrespective of transcript position.

Clinical Findings and Diagnostic Yield of Arrhythmogenic Cardiomyopathy Through Genomic Screening of Pathogenic or Likely Pathogenic Desmosome Gene Variants.

BACKGROUND: Genomic screening holds great promise for presymptomatic identification of hidden disease, and prevention of dramatic events, including sudden cardiac death associated with arrhythmogenic cardiomyopathy (ACM). Herein, we present findings from clinical follow-up of carriers of ACM-associated pathogenic/likely pathogenic desmosome variants ascertained through genomic screening. METHODS: Of 64 548 eligible participants in Geisinger MyCode Genomic Screening and Counseling program (2015-present), 92 individuals (0.14%) identified with pathogenic/likely pathogenic desmosome variants by clinical laboratory testing were referred for evaluation. We reviewed preresult medical history, patient-reported family history, and diagnostic testing results to assess both arrhythmogenic right ventricular cardiomyopathy and left-dominant ACM. RESULTS: One carrier had a prior diagnosis of dilated cardiomyopathy with arrhythmia; no other related diagnoses or diagnostic family history criteria were reported. Fifty-nine carriers (64%) had diagnostic testing in follow-up. Excluding the variant, 21/59 carriers satisfied at least one arrhythmogenic right ventricular cardiomyopathy task force criterion, 11 (52%) of whom harbored DSP variants, but only 5 exhibited multiple criteria. Six (10%) carriers demonstrated evidence of left-dominant ACM, including high rates of atypical late gadolinium enhancement by magnetic resonance imaging and nonsustained ventricular tachycardia. Two individuals received new cardiomyopathy diagnoses and received defibrillators for primary prevention. CONCLUSIONS: Genomic screening for pathogenic/likely pathogenic variants in desmosome genes can uncover both left- and right-dominant ACM. Findings of overt cardiomyopathy were limited but were most common in DSP-variant carriers and notably absent in PKP2-variant carriers. Consideration of the pathogenic/likely pathogenic variant as a major criterion for diagnosis is inappropriate in the setting of genomic screening.

Developing and Optimizing Innovative Tools to Address Familial Hypercholesterolemia Underdiagnosis: Identification Methods, Patient Activation, and Cascade Testing for Familial Hypercholesterolemia.

BACKGROUND: Familial hypercholesterolemia (FH) is the most common cardiovascular genetic disorder and, if left untreated, is associated with increased risk of premature atherosclerotic cardiovascular disease, the leading cause of preventable death in the United States. Although FH is common, fatal, and treatable, it is underdiagnosed and undertreated due to a lack of systematic methods to identify individuals with FH and limited uptake of cascade testing. METHODS AND RESULTS: This mixed-method, multi-stage study will optimize, test, and implement innovative approaches for both FH identification and cascade testing in 3 aims. To improve identification of individuals with FH, in Aim 1, we will compare and refine automated phenotype-based and genomic approaches to identify individuals likely to have FH. To improve cascade testing uptake for at-risk individuals, in Aim 2, we will use a patient-centered design thinking process to optimize and develop novel, active family communication methods. Using a prospective, observational pragmatic trial, we will assess uptake and effectiveness of each family communication method on cascade testing. Guided by an implementation science framework, in Aim 3, we will develop a comprehensive guide to identify individuals with FH. Using the Conceptual Model for Implementation Research, we will evaluate implementation outcomes including feasibility, acceptability, and perceived sustainability as well as health outcomes related to the optimized methods and tools developed in Aims 1 and 2. CONCLUSIONS: Data generated from this study will address barriers and gaps in care related to underdiagnosis of FH by developing and optimizing tools to improve FH identification and cascade testing.

Clinical outcomes of a genomic screening program for actionable genetic conditions.

PURPOSE: Three genetic conditions-hereditary breast and ovarian cancer syndrome, Lynch syndrome, and familial hypercholesterolemia-have tier 1 evidence for interventions that reduce morbidity and mortality, prompting proposals to screen unselected populations for these conditions. We examined the impact of genomic screening on risk management and early detection in an unselected population. METHODS: Observational study of electronic health records (EHR) among individuals in whom a pathogenic/likely pathogenic variant in a tier 1 gene was discovered through Geisinger's MyCode project. EHR of all eligible participants was evaluated for a prior genetic diagnosis and, among participants without such a diagnosis, relevant personal/family history, postdisclosure clinical diagnoses, and postdisclosure risk management. RESULTS: Eighty-seven percent of participants (305/351) did not have a prior genetic diagnosis of their tier 1 result. Of these, 65% had EHR evidence of relevant personal and/or family history of disease. Of 255 individuals eligible to have risk management, 70% (n = 179) had a recommended risk management procedure after results disclosure. Thirteen percent of participants (41/305) received a relevant clinical diagnosis after results disclosure. CONCLUSION: Genomic screening programs can identify previously unrecognized individuals at increased risk of cancer and heart disease and facilitate risk management and early cancer detection.

Pediatric reporting of genomic results study (PROGRESS): a mixed-methods, longitudinal, observational cohort study protocol to explore disclosure of actionable adult- and pediatric-onset genomic variants to minors and their parents.

BACKGROUND: Exome and genome sequencing are routinely used in clinical care and research. These technologies allow for the detection of pathogenic/likely pathogenic variants in clinically actionable genes. However, fueled in part by a lack of empirical evidence, controversy surrounds the provision of genetic results for adult-onset conditions to minors and their parents. We have designed a mixed-methods, longitudinal cohort study to collect empirical evidence to advance this debate. METHODS: Pediatric participants in the Geisinger MyCode® Community Health Initiative with available exome sequence data will have their variant files assessed for pathogenic/likely pathogenic variants in 60 genes designated as actionable by MyCode. Eight of these genes are associated with adult-onset conditions (Hereditary Breast and Ovarian Cancer Syndrome (HBOC), Lynch syndrome, MUTYH-associated polyposis, HFE-Associated Hereditary Hemochromatosis), while the remaining genes have pediatric onset. Prior to clinical confirmation of results, pediatric MyCode participants and their parents/legal guardians will be categorized into three study groups: 1) those with an apparent pathogenic/likely pathogenic variant in a gene associated with adult-onset disease, 2) those with an apparent pathogenic/likely pathogenic variant in a gene associated with pediatric-onset disease or with risk reduction interventions that begin in childhood, and 3) those with no apparent genomic result who are sex- and age-matched to Groups 1 and 2. Validated and published quantitative measures, semi-structured interviews, and a review of electronic health record data conducted over a 12-month period following disclosure of results will allow for comparison of psychosocial and behavioral outcomes among parents of minors (ages 0-17) and adolescents (ages 11-17) in each group. DISCUSSION: These data will provide guidance about the risks and benefits of informing minors and their family members about clinically actionable, adult-onset genetic conditions and, in turn, help to ensure these patients receive care that promotes physical and psychosocial health. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT03832985. Registered 6 February 2019.

Functional biology of the Steel syndrome founder allele and evidence for clan genomics derivation of COL27A1 pathogenic alleles worldwide.

Previously we reported the identification of a homozygous COL27A1 (c.2089G>C; p.Gly697Arg) missense variant and proposed it as a founder allele in Puerto Rico segregating with Steel syndrome (STLS, MIM #615155); a rare osteochondrodysplasia characterized by short stature, congenital bilateral hip dysplasia, carpal coalitions, and scoliosis. We now report segregation of this variant in five probands from the initial clinical report defining the syndrome and an additional family of Puerto Rican descent with multiple affected adult individuals. We modeled the orthologous variant in murine Col27a1 and found it recapitulates some of the major Steel syndrome associated skeletal features including reduced body length, scoliosis, and a more rounded skull shape. Characterization of the in vivo murine model shows abnormal collagen deposition in the extracellular matrix and disorganization of the proliferative zone of the growth plate. We report additional COL27A1 pathogenic variant alleles identified in unrelated consanguineous Turkish kindreds suggesting Clan Genomics and identity-by-descent homozygosity contributing to disease in this population. The hypothesis that carrier states for this autosomal recessive osteochondrodysplasia may contribute to common complex traits is further explored in a large clinical population cohort. Our findings augment our understanding of COL27A1 biology and its role in skeletal development; and expand the functional allelic architecture in this gene underlying both rare and common disease phenotypes.

Prevalence and Electronic Health Record-Based Phenotype of Loss-of-Function Genetic Variants in Arrhythmogenic Right Ventricular Cardiomyopathy-Associated Genes.

BACKGROUND: Arrhythmogenic right ventricular cardiomyopathy (ARVC) is associated with variants in desmosome genes. Secondary findings of pathogenic/likely pathogenic variants, primarily loss-of-function (LOF) variants, are recommended for clinical reporting; however, their prevalence and associated phenotype in a general clinical population are not fully characterized. METHODS: From whole-exome sequencing of 61 019 individuals in the DiscovEHR cohort, we screened for putative loss-of-function variants in PKP2, DSC2, DSG2, and DSP. We evaluated measures from prior clinical ECG and echocardiograms, manually over-read to evaluate ARVC diagnostic criteria, and performed a PheWAS (phenome-wide association study). Finally, we estimated expected penetrance using Bayesian inference. RESULTS: One hundred forty individuals (0.23%; 59±18 years old at last encounter; 33% male) had an ARVC variant (G+). None had an existing diagnosis of ARVC in the electronic health record, nor significant differences in prior ECG or echocardiogram findings compared with matched controls without variants. Several G+ individuals satisfied major repolarization (n=4) and ventricular function (n=5) criteria, but this prevalence matched controls. PheWAS showed no significant associations of other heart disease diagnoses. Combining our best genetic and disease prevalence estimates yields an estimated penetrance of 6.0%. CONCLUSIONS: The prevalence of ARVC loss-of-function variants is ≈1:435 in a general clinical population of predominantly European descent, but with limited electronic health record-based evidence of phenotypic association in our population, consistent with a low penetrance estimate. Prospective deep phenotyping and longitudinal follow-up of a large sequenced cohort is needed to determine the true clinical relevance of an incidentally identified ARVC loss-of-function variant.

Mcl-1 is a key regulator of apoptosis during CNS development and after DNA damage.

Despite the importance of Mcl-1, an anti-apoptotic Bcl-2 family member, in the regulation of apoptosis, little is known regarding its role in nervous system development and injury-induced neuronal cell death. Because germline deletion of Mcl-1 results in peri-implantation lethality, we address the function of Mcl-1 in the nervous system using two different conditional Mcl-1 mouse mutants in the developing nervous system. Here, we show for the first time that Mcl-1 is required for neuronal development. Neural precursors within the ventricular zone and newly committed neurons in the cortical plate express high levels of Mcl-1 throughout cortical neurogenesis. Loss of Mcl-1 in neuronal progenitors results in widespread apoptosis. Double labeling with active caspase 3 and Tuj1 reveals that newly committed Mcl1 deficient neurons undergo apoptosis as they commence migration away from the ventricular zone. Examination of neural progenitor differentiation in vitro demonstrated that cell death in the absence of Mcl1 is cell autonomous. Although conditional deletion of Mcl-1 in cultured neurons does not trigger apoptosis, loss of Mcl-1 sensitizes neurons to an acute DNA damaging insult. Indeed, the rapid reduction of Mcl-1 mRNA and protein levels are early events after DNA damage in neurons, and maintaining high Mcl-1 levels can protect neurons against death. Together, our results are the first to demonstrate the requirement of Mcl-1, an anti-apoptotic Bcl-2 family protein, for cortical neurogenesis and the survival of neurons after DNA damage.

Matrix metalloproteinase activation and blood-brain barrier breakdown following thrombolysis.

Thrombolysis with tissue plasminogen activator (tPA) is the only pharmacotherapy available for cerebral ischemia. However, the use of tPA can increase the risk of hemorrhage due to blood-brain barrier (BBB) breakdown. Recent evidence suggests that increased activation of matrix metalloproteinases (MMPs) may be involved in this breakdown. This study examines the temporal profile of MMP-2 and -9 following tPA administration to ischemic rats. Male Sprague-Dawley rats were randomly assigned to one of four groups (Sham-tPA; Sham-Saline; Ischemia-tPA; Ischemia-Saline; group n = 6, total N = 120). Focal embolic ischemia was induced by middle cerebral artery occlusion through injection of an autologous clot. One hour post-surgery, tPA (10 mg/kg) or saline was delivered intravenously and animals were euthanized at 3, 6, 12, or 24 h after onset of ischemia. Infarct volume was measured by TTC staining; BBB components examined immunohistochemically; and MMP activation measured by gelatin zymography. Our results show that tPA significantly reduced infarct volumes (overall infarct volume-Sham-tPA: 5.80 +/- 4.55 [mean +/- SE]; Sham-Saline: 5.00 +/- 4.23; Ischemia-tPA: 186.1 +/- 73.45; Ischemia-Saline: 284.8 +/- 88.74; all P < 0.05). Treatment with tPA was also associated with the activation of MMP-9 at 6, 12, and 24 h following ischemia. No temporal changes were observed in MMP-2 activation, although tPA administration increased its activity compared to saline treatment. Analyses of immunohistochemistry showed that destruction of components of the BBB followed MMP-9 activation. Thus, increased MMP-9 activation may, in part, be responsible for the increases in hemorrhagic transformation reported with use of tPA. Our study is the first to demonstrate the temporal profile of MMP activation following thrombolysis with tPA in a model of thrombotic focal cerebral ischemia.