Using high-resolution variant frequencies to empower clinical genome interpretation

Association of Pathogenic/Likely Pathogenic Genetic Variants for Cardiomyopathies With Clinical Outcomes: A Multiancestry Analysis in the All of Us Research Program

BACKGROUND

This study aimed to evaluate the prevalence of pathogenic/likely pathogenic cardiomyopathy variant carriers in a multiancestry US population and examine the risk of adverse clinical outcomes.

METHODS

This retrospective cohort study included multiancestry US adults aged ≥18 years with sequencing data from the All of Us Research Program. Pathogenic/likely pathogenic variants in cardiomyopathy genes were identified using the ClinVar database. The primary outcome was heart failure. Secondary outcomes included cardiomyopathy and arrhythmia. Outcomes were identified from electronic health records. Interval-censored Cox models, taking age on the timescale, were used to assess the risk of outcomes in pathogenic/likely pathogenic variant carriers with noncarriers as the reference group.

RESULTS

Among 167 435 individuals (median age, 55.2 [39.5-66.3] years; 61.7% female; 40.7% non-European ancestry) included, the proportion of pathogenic/likely pathogenic cardiomyopathy variant carriers was 0.7% of the overall population and 0.8%, 0.8%, 0.5%, and 1.2% of European, African, East Asian, and South Asian ancestry individuals, respectively. Over the lifetime, there were 12 867 heart failure events (205 in carriers and 12 662 in noncarriers), with an incidence rate of 3.05 (95% CI, 2.66-3.49) per 1000 person-years in carriers and 1.37 (95% CI, 1.35-1.40) in noncarriers (HRadj, 2.30 [95% CI, 2.04-2.60]). Cardiomyopathy occurred in 5164 (161 in carriers and 5003 in noncarriers), with an incidence rate of 2.38 (95% CI, 2.04-2.78) per 1000 person-years among carriers and 0.54 (95% CI, 0.53-0.56) in noncarriers (HRadj, 4.31 [95% CI, 3.73-4.97]). There were 19 405 arrhythmia events (263 in carriers and 19 142 in noncarriers), with an incidence rate of 3.93 (95% CI, 3.48-4.44) per 1000 person-years among carriers and 2.09 (95% CI, 2.06-2.12) in noncarriers (HRadj, 2.12 [95% CI, 1.78-2.53]).

CONCLUSIONS

Pathogenic/likely pathogenic cardiomyopathy variant carriers have an increased risk of heart failure, cardiomyopathy, and arrhythmias. Despite the modest overall prevalence, the associated risks suggest potential benefits of targeted genetic screening for early detection and management.

Comprehensive functional splicing analysis of non-canonical CNGB3 variants using in vitro minigene splice assays

Variants in the CNGB3 gene, encoding the B3-subunit of the cone photoreceptor cyclic nucleotide gated channel, are a major cause of autosomal recessive achromatopsia, a rare inherited retinal disease. The mutation spectrum of achromatopsia-associated CNGB3 variants comprises all types of mutations, including those that are straightforward to evaluate in molecular genetic diagnostics, such as frame-shifting, nonsense, and canonical splice site variants. Additionally, variants have been identified within splice regions outside the conserved ±1,2 splice site dinucleotides, making their potential impact on disease association challenging to interpret. This poses a major hurdle for clinical interpretation of causality between the patient's genotype and the proposed clinical diagnosis, but also for the inclusion of such patients into clinical trials for gene augmentation therapy, for which only patients with confirmed (likely) pathogenic CNGB3 variants are eligible. We here performed comprehensive genetic functional analysis of 21 candidate spliceogenic CNGB3 variants-15 reported and 6 novel variants-by means of in vitro minigene splice assays and cDNA analysis, and characterization of spliceogenic events by subcloning, Sanger-sequencing, and capillary fragment analysis. For 16 variants, an impact on splicing was confirmed, supporting the reclassification of 86% of variants of uncertain significance as likely pathogenic or pathogenic according to the ACMG/AMP guidelines. This reclassification enables the confirmation of patients' genotypes, both retrospectively and prospectively. © 2025 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Reclassifying NOBOX variants in primary ovarian insufficiency cases with a corrected gene model and a novel quantitative framework

STUDY QUESTION

How updated expression and genomic data combined with a disease/disorder-specific classification system can be used to correct a gene model for a better evaluation of the pathogenicity of variants found in patients?

SUMMARY ANSWER

By combining available genomic and transcriptomic data from several species and a quantitative classification framework with primary ovarian insufficiency (POI)-adjusted parameters, we correct the human NOBOX (newborn ovary homeobox) gene model and provide a reclassification of variants previously reported in POI cases.

WHAT IS KNOWN ALREADY

The NOBOX gene, encoding a gonad-specific transcription factor with a crucial role in early folliculogenesis and considered a major gene involved in POI, is currently described as being expressed as four transcripts, the longest one considered canonical. All the variants identified in POI cases have been evaluated according to this canonical transcript, and the various functional tests have been performed using the corresponding predicted protein.

STUDY DESIGN, SIZE, DURATION

We refined and corrected the NOBOX gene model using available genomic and RNAseq data in human and 16 other mammalian species. Expression data were selected for tissue specificity, strand specificity, and coverage. The analysis of RNAseq data from different ovarian fetal stages allows for a time-course description of NOBOX isoforms. Literature was scanned to retrieve NOBOX variants reported in POI cases, and NOBOX variants present in ClinVar and GnomAD 4 databases were also retrieved.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Strand-specific RNAseq data from human fetal ovaries and human adult testes were analysed to infer the correct human NOBOX gene isoforms. The conservation of the gene structure was verified by combining the aligned genomic sequences from 17 mammalian species covering a wide phylogenetic range and the relevant RNAseq data. As changing a gene model implies a reclassification of variants, we set up a quantitative framework with updated variant frequencies from GnomAD4 and POI-adjusted parameters following the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines. Using this framework, we reclassified 44 NOBOX variants reported in POI patients and families, 117 NOBOX variants reported in ClinVar, and 2613 NOBOX variants present in GnomAD4.

MAIN RESULTS AND THE ROLE OF CHANCE

The corrected NOBOX gene model proposes the invalidation of two transcripts, including the canonical one. The two correct isoforms were present in fetal ovarian samples, and only one was detected in adult testes. Only 14 variants remained as possibly causative for POI. Furthermore, this re-evaluation strongly suggests that NOBOX biallelic variants are the most likely cause of POI.

LARGE SCALE DATA

Large tables are provided as supplementary data sets on the Zenodo repository.

LIMITATIONS, REASONS FOR CAUTION

The proposed gene model is robust but relies on available transcriptomic data covering a range of time points and tissues. Our scoring system was manually adjusted and other laboratories can implement it with different parameters.

WIDER IMPLICATIONS OF THE FINDINGS

For the NOBOX variants that cannot be considered pathogenic or causative anymore, the genome/exome sequencing data of the corresponding patients should be reanalysed. Furthermore, the functional studies performed using the obsolete coding sequence should be reconsidered. The corrected gene model should be taken into account when evaluating novel NOBOX variants identified in POI patients. Our results highlight the importance of the careful assessment of the most updated expression data for validating a gene model, enabling a correct evaluation of the pathogenicity of variants found in patients. The proposed quantitative framework developed here can be used for the classification of variants in other genes underlying POI. Furthermore, the global approach based on quantitatively adjusting the ACMG/AMP guidelines could be extended to other inherited pathologies.

STUDY FUNDING/COMPETING INTEREST(S)

This project was not funded. All the authors have no conflict of interest to disclose.

Low Penetrance Sarcomere Variants Contribute to Additive Risk in Hypertrophic Cardiomyopathy

BACKGROUND

Classically, hypertrophic cardiomyopathy (HCM) has been viewed as a single-gene (monogenic) disease caused by pathogenic variants in sarcomere genes. Pathogenic sarcomere variants are individually rare and convey high risk for developing HCM (highly penetrant). Recently, important polygenic contributions have also been characterized. Low penetrance sarcomere variants (LowSVs) at intermediate frequencies and effect sizes have not been systematically investigated. We hypothesize that LowSVs may be common in HCM with substantial influence on disease risk and severity.

METHODS

Among all sarcomere variants observed in the Sarcomeric Human Cardiomyopathy Registry (SHaRe), we identified putative LowSVs defined by (1) population frequency greater than expected for highly penetrant (monogenic) HCM (allele frequency >5×10-5 in the Genome Aggregation Database, gnomAD) and (2) moderate enrichment (>2×) in patients with HCM compared with gnomAD. LowSVs were examined for their association with disease severity and clinical outcomes. Functional effects of selected LowSVs were assessed using induced pluripotent stem cell-derived cardiomyocytes. Association of LowSVs with HCM-adjacent traits in the general population was tested using UK Biobank cardiac magnetic resonance imaging data.

RESULTS

Among 6045 patients and 1159 unique variants in sarcomere genes, 12 LowSVs were identified. LowSVs were collectively common in the general population (1:350) and moderately enriched in HCM (aggregate odds ratio, 14.9 [95% CI, 12.5-17.9]). Isolated LowSVs were associated with an older age of HCM diagnosis and fewer adverse events. However, LowSVs in combination with a pathogenic sarcomere variant conferred higher morbidity (eg, composite adverse event hazard ratio, 5.4 [95% CI, 3.0-9.8] versus single pathogenic sarcomere variant, 2.0 [95% CI, 1.8-2.2]; P<0.001). An intermediate functional impact was validated for 2 specific LowSVs-MYBPC3 c.442G>A (partial splice gain) and TNNT2 c.832C>T (intermediate effect on contractile mechanics). Cardiac magnetic resonance imaging analysis of the general population revealed 5 of 12 LowSVs were significantly associated with HCM-adjacent traits without overt HCM.

CONCLUSIONS

This study establishes a new class of low penetrance sarcomere variants that are relatively common in the population. When penetrant, isolated LowSVs cause mild HCM. In combination with pathogenic sarcomere variants, LowSVs markedly increase disease severity, supporting a clinically significant additive effect. Last, LowSVs also contribute to age-related remodeling even in the absence of overt HCM.

Optimizing gene panels for equitable reproductive carrier screening: The Goldilocks approach

PURPOSE

Professional organizations recommend pan-ancestry carrier screening for autosomal recessive and X-linked conditions. Advances in DNA sequencing have allowed the analysis of hundreds of genes; however, the optimal number of genes for carrier screening remains unclear. The American College of Medical Genetics and Genomics (ACMG) has proposed a tiered approach recommending screening for 113 genes.

METHODS

We analyzed ClinVar and gnomAD v4.1.0, for genes associated with serious autosomal recessive and X-linked conditions and modeled screening performance across panels of varying compositions and sizes in diverse genetic ancestries. We also reevaluated the ACMG gene list using the updated gnomAD data.

RESULTS

We identified potential inconsistencies in the ACMG gene lists, particularly in the carrier test performance (defined as a positive yield) for underrepresented genetic ancestry groups. Modeling of the population data for 1310 genes revealed that the screening of 152, 248, 531, and 725 genes achieved 90%, 95%, 99%, and 99.7% positive yields, respectively, in couples. Real-world data from the screening of more than 60,000 couples were used to validate the model.

CONCLUSION

Our methodology optimizes the gene content of carrier screening panels for diverse ancestry groups, provides a mechanism for continually updating guidelines, ensures consistency with genomic population data, and improves equity across populations.

Study of the impact of ClinGen Revisions on ACMG/AMP variant semi-automatic classification for Rare Diseases diagnosis

With the rapid development of massive sequencing technologies, the analysis of genetic variants for clinical diagnosis has exponentially escalated, particularly in the context of Rare Diseases (RDs). Diagnosing them involves identifying the genetic variants responsible for the underlying pathology development. In 2015, the American College of Medical Genetics (ACMG) established a set of recommendations to assess the evidence associated with each variant, aiming to achieve a standardized five tier classification. Over the past 5 years, ClinGen, the NIH-funded Clinical Genome Resource, has reviewed these criteria in order to make variant classification a more reproducible and rigorous process. This paper examines the impact of ClinGen-Rev modifications on variant classification, comparing them with the ACMG-2015 original recommendations. After analyzing sets of genetic variants, extracted from VCFs samples, using both criteria, we observed a change in 8.0 % of the clinical verdicts for these variants. ClinGen-Rev modifications correctly categorized 89.2 % of the curated variants, representing a significant improvement compared to the 65.6 % achieved by ACMG-2015. We also analyzed the modifications impact in a real like clinical setting, showing a significant overall reduction of VUS variants and thus potential reduction in analysis time. Finally, we discuss the underlying reasons for the most relevant changes in terms of specific labels and present their implications on the prioritization and selection process of variants, identifying some recommendations of key significant importance.

Differential inclusion of NEB exons 143 and 144 provides insight into NEB-related myopathy variant interpretation and disease manifestation

Biallelic pathogenic variants in the gene encoding nebulin (NEB) are a known cause of congenital myopathy. We present two brothers with congenital myopathy and compound heterozygous variants (NC_000002.12:g.151692086G>T; NM_001271208.2: c.2079C>A; p.(Cys693Ter) and NC_000002.12:g.151533439T>C; NM_001271208.2:c.21522+3A>G) in NEB. Transcriptomic sequencing on affected individual muscles revealed that the extended splice variant c.21522+3A>G causes exon 144 skipping. Nebulin isoforms containing exon 144 are known to be mutually exclusive with isoforms containing exon 143, and these isoforms are differentially expressed during development and in adult skeletal muscles. Affected individuals' MRI patterns of muscle involvement were compared with the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these affected individuals better fits the distribution of exon 144-containing isoforms in muscle than with previously published MRI findings in NEB-related disease due to other variants. Our report introduces disease pathogenesis and manifestation as a result of alteration of isoform distributions in muscle.

Advances in genetic diagnosis and therapy of hereditary heart disease: a bibliometric review from 2004 to 2024

Hereditary heart disease (HHD) is a series of cardiac disorders associated with monogenic or polygenic abnormalities and is one of the leading causes of sudden death, particularly in young adults. The updated European Cardiology guideline for cardiomyopathies provides the first comprehensive summary of genotyping, imaging, and therapy recommendations for inherited cardiomyopathies, but still lacks a comprehensive discussion of research advances and future trends in genetic diagnosis and therapy of HHD. Our research aims to fill this gap. Bibliometric analysis software (CiteSpace 6.3.R1, VOSviewer 1.6.18, and Scimago Graphica) was used to analyze the general information, trends, and emerging foci of HHD in the past 20 years, including author, country, institution, keyword, and so on. There were 5,757 publications were screened and aggregated in the database, including 1876 reviews and 3,881 articles. Hypertrophic cardiomyopathy (HCM), arrhythmogenic cardiomyopathy (ACM), Brugada syndrome (BrS), myocardial amyloidosis, and Fabry disease (FD) were the main types of HHD that were explored in greater depth. Moreover, new diagnostic methods, clinical cohorts, and genetically targeted therapies for HHD patients are key research hotspots. The relationship between the pathogenicity of genes and prognosis will become increasingly important for therapy.

Development and validation of animal variant classification guidelines to objectively evaluate genetic variant pathogenicity in domestic animals

Assessing the pathogenicity of a disease-associated genetic variant in animals accurately is vital, both on a population and individual scale. At the population level, breeding decisions based on invalid DNA tests can lead to the incorrect inclusion or exclusion of animals and compromise the long-term health of a population, and at the level of the individual animal, lead to incorrect treatment and even life-ending decisions. Criteria to determine pathogenicity are not standardized, i.e., no guidelines for animal variants are available. Here, we aimed to develop and validate guidelines to be used by the community for Mendelian disorders in domestic animals to classify variants in categories based on standardized criteria. These so-called animal variant classification guidelines (AVCG) were based on those developed for humans by The American College of Medical Genetics and Genomics (ACMG). In a direct comparison, 83% of the pathogenic variants were correctly classified with ACMG, while this increased to 92% with AVCG. We described methods to develop datasets for benchmarking the criteria and identified the most optimal in silico variant effect predictor tools. As the reproducibility was high, we classified 72 known disease-associated variants in cats and 40 other disease-associated variants in eight additional species.

Prequalification of genome-based newborn screening for severe childhood genetic diseases through federated training based on purifying hyperselection

Genome-sequence-based newborn screening (gNBS) has substantial potential to improve outcomes in hundreds of severe childhood genetic disorders (SCGDs). However, a major impediment to gNBS is imprecision due to variants classified as pathogenic (P) or likely pathogenic (LP) that are not SCGD causal. gNBS with 53,855 P/LP variants, 342 genes, 412 SCGDs, and 1,603 therapies was positive in 74% of UK Biobank (UKB470K) adults, suggesting 97% false positives. We used the phenomenon of purifying hyperselection, which acts to decrease the frequency of SCGD causal diplotypes, to reduce false positives. Training of gene-disease-inheritance mode-diplotype tetrads in 618,290 control and affected subjects identified 293 variants or haplotypes and seven genes with variable inheritance contributing higher positive diplotype counts than consistent with purifying hyperselection and with little or no evidence of SCGD causality. With these changes, 2.0% of UKB470K adults were positive. In contrast, gNBS was positive in 7.2% of 3,118 critically ill children with suspected SCGDs and 7.9% of 705 infant deaths. When compared with rapid diagnostic genome sequencing (RDGS), gNBS had 99.1% recall. In eight true-positive children, gNBS was projected to decrease time to diagnosis by a median of 121 days and avoid life-threatening disease presentations in four children, organ damage in six children, ∼$1.25 million in healthcare cost, and ten (1.4%) infant deaths. Federated training predicated on purifying hyperselection provides a general framework to attain high precision in population screening. Federated training across many biobanks and clinical trials can provide a privacy-preserving mechanism for qualification of gNBS in diverse genetic ancestries.

Genetic and phenotypic architecture of human myocardial trabeculation

Cardiac trabeculae form a network of muscular strands that line the inner surfaces of the heart. Their development depends on multiscale morphogenetic processes and, while highly conserved across vertebrate evolution, their role in the pathophysiology of the mature heart is not fully understood. Here we report variant associations across the allele frequency spectrum for trabecular morphology in 47,803 participants of the UK Biobank using fractal dimension analysis of cardiac imaging. We identified an association between trabeculation and rare variants in 56 genes that regulate myocardial contractility and ventricular development. Genome-wide association studies identified 68 loci in pathways that regulate sarcomeric function, differentiation of the conduction system and cell fate determination. We found that trabeculation-associated variants were modifiers of cardiomyopathy phenotypes with opposing effects in hypertrophic and dilated cardiomyopathy. Together, these data provide insights into mechanisms that regulate trabecular development and plasticity, and identify a potential role in modifying monogenic disease expression.

Pharmacogenetic DPYD allele variant frequencies: A comprehensive analysis across an ancestrally diverse Iranian population

BACKGROUND

Cancer treatment has improved over the past decades, but many cancer patients still experience adverse drug reactions (ADRs). Pharmacogenomics (PGx), known as personalized treatment, is a pillar of precision medicine that aims to optimize the efficacy and safety of medications by studying the germline variations. Germline variations in the DPYD lead to significant ADRs. The present cross-sectional study aims to evaluate the allele frequency of the DPYD gene variations in the Iranian population to provide insights into personalized treatment decisions in the Iranian population.

METHODS

The allele frequency of 51 pharmacogenetic variations in the clinically relevant DPYD was assessed in a representative sample set of 1142 unrelated Iranian individuals and subpopulations of different ethnic groups who were genotyped using the Infinium Global Screening Array-24 BeadChip.

RESULTS

The genotyping assay revealed eight pharmacogenetic variants including DPYD rs1801265 (c.85T > C; DPYD*9A), rs2297595 (c.496A > G), rs1801158 (c.1601G > A; DPYD*4), rs1801159 (c.1627A > G; DPYD*5), rs1801160 (c.2194G > A; DPYD*6), rs17376848 (c.1896T > C), rs56038477 (c.1236G > A; HapB3), and rs75017182 (c.1129-5923C > G; HapB3) with minor allele frequency (MAF) ≥ 1%.

CONCLUSION

The results of the study reveal significant genetic variations among Iranian population that could significantly influence clinical decision-making. These variants, with their potential to explain the substantial variability in drug response phenotypes among different populations, shed light on a crucial aspect of pharmacogenomics. These findings not only provide valuable insights but also inspire the design and implementation of future pharmacogenomic clinical trials, motivating further research in this crucial area.

Analysis of 1386 epileptogenic brain lesions reveals association with DYRK1A and EGFR

Lesional focal epilepsy (LFE) is a common and severe seizure disorder caused by epileptogenic lesions, including malformations of cortical development (MCD) and low-grade epilepsy-associated tumors (LEAT). Understanding the genetic etiology of these lesions can inform medical and surgical treatment. We conducted a somatic variant enrichment mega-analysis in brain tissue from 1386 individuals who underwent epilepsy surgery, including 599 previously unpublished individuals with ultra-deep ( > 1600x) targeted panel sequencing. Here we confirm four known associations (BRAF, SLC35A2, MTOR, PTPN11), support eight associations without prior statistical support (FGFR1, PIK3CA, AKT3, NF1, PTEN, RHEB, KRAS, NRAS), and identify novel associations for two genes, DYRK1A and EGFR. Both novel genes show specific histopathological phenotypes, interact with LFE genes and pathways, and may represent promising candidates as biomarkers and potentially druggable targets.

Classification of PTEN germline non-truncating variants: a new approach to interpretation

BACKGROUND

PTEN hamartoma tumour syndrome (PHTS) encompasses distinct syndromes, including Cowden syndrome resulting from PTEN pathogenic variants. Missense variants account for 30% of PHTS cases, but their classification remains challenging. To address these difficulties, guidelines were published by the Clinical Genome Resource PTEN Variant Curation Expert Panel.

METHODS

Between 2010 and 2020, the Bergonie Institute reference laboratory identified 76 different non-truncating PTEN variants in 166 patients, 17 of which have not previously been reported. Variants were initially classified following the current guidelines. Subsequently, a new classification method was developed based on four main criteria: functional exploration, phenotypic features and familial segregation, in silico modelling, and allelic frequency.

RESULTS

This new method of classification is more discriminative and reclassifies 25 variants, including 8 variants of unknown significance.

CONCLUSION

This report proposes a revision of the current PTEN variant classification criteria which at present rely on functional tests evaluating only the phosphatase activity of PTEN and apply a particularly stringent clinical PHTS score.The classification of non-truncating variants of PTEN is facilitated by taking into consideration protein stability for variants with intact phosphatase activity, clinical and segregation criteria adapted to the phenotypic variability of PHTS and by specifying the allelic frequency of variants in the general population. This novel method of classification remains to be validated in a prospective cohort.

Specifications of the ACMG/AMP variant curation guidelines for the analysis of germline ATM sequence variants

The ClinGen Hereditary Breast, Ovarian, and Pancreatic Cancer (HBOP) Variant Curation Expert Panel (VCEP) is composed of internationally recognized experts in clinical genetics, molecular biology, and variant interpretation. This VCEP made specifications for the American College of Medical Genetics and Association for Molecular Pathology (ACMG/AMP) guidelines for the ataxia telangiectasia mutated (ATM) gene according to the ClinGen protocol. These gene-specific rules for ATM were modified from the ACMG/AMP guidelines and were tested against 33 ATM variants of various types and classifications in a pilot curation phase. The pilot revealed a majority agreement between the HBOP VCEP classifications and the ClinVar-deposited classifications. Six pilot variants had conflicting interpretations in ClinVar, and re-evaluation with the VCEP's ATM-specific rules resulted in four that were classified as benign, one as likely pathogenic, and one as a variant of uncertain significance (VUS) by the VCEP, improving the certainty of interpretations in the public domain. Overall, 28 of the 33 pilot variants were not VUS, leading to an 85% classification rate. The ClinGen-approved, modified rules demonstrated value for improved interpretation of variants in ATM.

Identification of Novel Genetic Risk Variants Associated with Hidradenitis Suppurativa in an Exome Sequencing Cohort of 92,455 Individuals

INTRODUCTION

Hidradenitis suppurativa (HS) is a prevalent and persistent inflammatory skin disorder, lacking a known cure or effective biomarkers for early diagnosis at present. The genetic determinants of HS have not been fully documented, but it is believed to result from a combination of genetic and environmental factors.

METHODS

To identify relevant HS gene variants in sporadic HS patients, this study utilized longitudinal electronic health records (EHRs) and whole-exome sequencing. DNA exome sequencing data from 92,455 participant samples in the MyCode biobank, linked to Geisinger's EHR, were analyzed. This cohort included 1,092 HS cases and 91,363 healthy controls. The MyCode EHR has a median longitudinal follow-up of 15 years per participant, with an average of 87 clinical encounters, 687 laboratory tests, and 7 procedures.

RESULTS

There were 1,092 (901 females and 191 males) participants aged 14-89 years (median 47 years) with HS (L73.2), indicating a 1.18% prevalence and accounting for a 4.7:1 female-to-male ratio among the individuals presenting for clinical care. γ-secretase complex, syndromic, and autoinflammatory gene variants were assessed. Potential pathogenic variants were identified among 66 individuals in the HS genes studied. Molecularly, the estimated HS variant prevalence was 1:1,400 in the cohort, 12.3% of variant carriers had HS diagnosis in EHR.

CONCLUSIONS

Using longitudinal EHR data, genomic screening identified HS-associated gene variants in a defined group of sporadic HS patients to augment the clinical diagnosis, particularly in cases of ambiguity. Based on this study, the field of skin disorders can benefit from a personalized approach to HS diagnosis using large-scale sequencing.

High-resolution functional mapping of RAD51C by saturation genome editing

Pathogenic variants in RAD51C confer an elevated risk of breast and ovarian cancer, while individuals homozygous for specific RAD51C alleles may develop Fanconi anemia. Using saturation genome editing (SGE), we functionally assess 9,188 unique variants, including >99.5% of all possible coding sequence single-nucleotide alterations. By computing changes in variant abundance and Gaussian mixture modeling (GMM), we functionally classify 3,094 variants to be disruptive and use clinical truth sets to reveal an accuracy/concordance of variant classification >99.9%. Cell fitness was the primary assay readout allowing us to observe a phenomenon where specific missense variants exhibit distinct depletion kinetics potentially suggesting that they represent hypomorphic alleles. We further explored our exhaustive functional map, revealing critical residues on the RAD51C structure and resolving variants found in cancer-segregating kindred. Furthermore, through interrogation of UK Biobank and a large multi-center ovarian cancer cohort, we find significant associations between SGE-depleted variants and cancer diagnoses.

The Genetic Factors Influencing Cardiomyopathies and Heart Failure across the Allele Frequency Spectrum

Heart failure (HF) remains a major cause of mortality and morbidity worldwide. Understanding the genetic basis of HF allows for the development of disease-modifying therapies, more appropriate risk stratification, and personalised management of patients. The advent of next-generation sequencing has enabled genome-wide association studies; moving beyond rare variants identified in a Mendelian fashion and detecting common DNA variants associated with disease. We summarise the latest GWAS and rare variant data on mixed and refined HF aetiologies, and cardiomyopathies. We describe the recent understanding of the functional impact of titin variants and highlight FHOD3 as a novel cardiomyopathy-associated gene. We describe future directions of research in this field and how genetic data can be leveraged to improve the care of patients with HF.

Genetic disease risks of under-represented founder populations in New York City

The detection of founder pathogenic variants, those observed in high frequency only in a group of individuals with increased inter-relatedness, can help improve delivery of health care for that community. We identified 16 groups with shared ancestry, based on genomic segments that are shared through identity by descent (IBD), in New York City using the genomic data of 25,366 residents from the All Of Us Research Program and the Mount Sinai BioMe biobank. From these groups we defined 8 as founder populations, mostly communities currently under-represented in medical genomics research, such as Puerto Rican, Garifuna and Filipino/Pacific Islanders. The enrichment analysis of ClinVar pathogenic or likely pathogenic (P/LP) variants in each group identified 202 of these damaging variants across the 8 founder populations. We confirmed disease-causing variants previously reported to occur at increased frequencies in Ashkenazi Jewish and Puerto Rican genetic ancestry groups, but most of the damaging variants identified have not been previously associated with any such founder populations, and most of these founder populations have not been described to have increased prevalence of the associated rare disease. Twenty-five of 51 variants meeting Tier 2 clinical screening criteria (1/100 carrier frequency within these founder groups) have never previously been reported. We show how population structure studies can provide insights into rare diseases disproportionately affecting under-represented founder populations, delivering a health care benefit but also a potential source of stigmatization of these communities, who should be part of the decision-making about implementation into health care delivery.

ParSE-seq: a calibrated multiplexed assay to facilitate the clinical classification of putative splice-altering variants

Interpreting the clinical significance of putative splice-altering variants outside canonical splice sites remains difficult without time-intensive experimental studies. To address this, we introduce Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed assay to quantify variant effects on RNA splicing. We first apply this technique to study hundreds of variants in the arrhythmia-associated gene SCN5A. Variants are studied in 'minigene' plasmids with molecular barcodes to allow pooled variant effect quantification. We perform experiments in two cell types, including disease-relevant induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). The assay strongly separates known control variants from ClinVar, enabling quantitative calibration of the ParSE-seq assay. Using these evidence strengths and experimental data, we reclassify 29 of 34 variants with conflicting interpretations and 11 of 42 variants of uncertain significance. In addition to intronic variants, we show that many synonymous and missense variants disrupted RNA splicing. Two splice-altering variants in the assay also disrupt splicing and sodium current when introduced into iPSC-CMs by CRISPR-Cas9 editing. ParSE-seq provides high-throughput experimental data for RNA-splicing to support precision medicine efforts and can be readily adopted to study other loss-of-function genotype-phenotype relationships.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Clinical exome sequencing unravels the diverse spectrum of genetic heterogeneity and genotype-phenotype correlations in hypertrophic cardiomyopathy

BACKGROUND

Catalogues of pathogenic genetic mutations in hypertrophic cardiomyopathy (HCM) are disproportionately small when compared to that of the size of the population with South Asian ancestry and their collective increased risk of heart disease.

METHODS

We conducted clinical exome sequencing of 200 HCM patients to identified cardiomyopathy-associated genetic mutations. The clinical and echocardiographic characteristics of genotype-positive and genotype-negative patients were compared, and the likelihood of detecting a positive genetic test result was evaluated. Allelic burden analysis was done to compare the minor allele frequencies (MAF) of the pathogenic or likely pathogenic (P/LP) variants and variants of uncertain significance (VUSs) identified in the cohort against various population genomics databases.

RESULTS

The genetic yield was 40% for P/LP variants, with MYBPC3 and MYH7 as the predominant sarcomere genes. Younger age-at-diagnosis, family history of HCM, asymmetric hypertrophic (ASH) pattern, the ratio of the interventricular septum to posterior wall thickness (IVS/PW ratio), left atrial (LA) dimensions, severe mitral regurgitation grade (MR grade), late gadolinium enhancement (LGE) detected fibrosis and absence of hypertension were associated with an increased likelihood of HCM-associated variants. Patients who experienced ventricular tachycardia and premature cardiovascular death were significantly likely to carry MYBPC3 or loss-of-function variants. LA and interventricular septal (IVS) dimensions were associated with MYH7 variants. The rare variant burden for P/LP variants and VUSs was significantly enriched in HCM cases compared to population controls.

CONCLUSION

Our study provides a comprehensive evaluation of HCM-associated genetic mutations from an Indian population. The identified genotype-phenotype associations could improve the yield of targeted genetic testing in HCM.

Evidence-based recommendations for gene-specific ACMG/AMP variant classification from the ClinGen ENIGMA BRCA1 and BRCA2 Variant Curation Expert Panel

The ENIGMA research consortium develops and applies methods to determine clinical significance of variants in hereditary breast and ovarian cancer genes. An ENIGMA BRCA1/2 classification sub-group, formed in 2015 as a ClinGen external expert panel, evolved into a ClinGen internal Variant Curation Expert Panel (VCEP) to align with Food and Drug Administration recognized processes for ClinVar contributions. The VCEP reviewed American College of Medical Genetics and Genomics/Association of Molecular Pathology (ACMG/AMP) classification criteria for relevance to interpreting BRCA1 and BRCA2 variants. Statistical methods were used to calibrate evidence strength for different data types. Pilot specifications were tested on 40 variants and documentation revised for clarity and ease of use. The original criterion descriptions for 13 evidence codes were considered non-applicable or overlapping with other criteria. Scenario of use was extended or re-purposed for eight codes. Extensive analysis and/or data review informed specification descriptions and weights for all codes. Specifications were applied to pilot variants with pre-existing ClinVar classification as follows: 13 uncertain significance or conflicting, 14 pathogenic and/or likely pathogenic, and 13 benign and/or likely benign. Review resolved classification for 11/13 uncertain significance or conflicting variants and retained or improved confidence in classification for the remaining variants. Alignment of pre-existing ENIGMA research classification processes with ACMG/AMP classification guidelines highlighted several gaps in the research processes and the baseline ACMG/AMP criteria. Calibration of evidence strength was key to justify utility and strength of different data types for gene-specific application. The gene-specific criteria demonstrated value for improving ACMG/AMP-aligned classification of BRCA1 and BRCA2 variants.

Prevalence and Reclassification of Genetic Variants in South African Populations with Breast Cancer

Concurrent testing of numerous genes for hereditary breast cancer (BC) is available but can result in management difficulties. We evaluated use of an expanded BC gene panel in women of diverse South African ancestries and assessed use of African genomic data to reclassify variants of uncertain significance (VUS). A total of 331 women of White, Black African, or Mixed Ancestry with BC had a 9-gene panel test, with an additional 75 genes tested in those without a pathogenic/likely pathogenic (P/LP) variant. The proportion of VUS reclassified using ClinGen gene-specific allele frequency (AF) thresholds or an AF > 0.001 in nonguidelines genes in African genomic data was determined. The 9-gene panel identified 58 P/LP variants, but only two of the P/LP variants detected using the 75-gene panel were in confirmed BC genes, resulting in a total of 60 (18.1%) in all participants. P/LP variant prevalence was similar across ancestry groups, but VUS prevalence was higher in Black African and Mixed Ancestry than in White participants. In total, 611 VUS were detected, representing 324 distinct variants. 10.8% (9/83) of VUS met ClinGen AF thresholds in genomic data while 10.8% (26/240) in nonguideline genes had an AF > 0.001. Overall, 27.0% of VUS occurrences could potentially be reclassified using African genomic data. Thus, expanding the gene panel yielded few clinically actionable variants but many VUS, particularly in participants of Black African and Mixed Ancestry. However, use of African genomic data has the potential to reclassify a significant proportion of VUS.

Diagnostic yield from cardiac gene testing for inherited cardiac conditions and re-evaluation of pre-ACMG variants of uncertain significance

BACKGROUND

Inherited cardiomyopathies (HCM, DCM, ACM) and cardiac ion channelopathies (long QT/Brugada syndromes, CPVT) are associated with significant morbidity and mortality; however, diagnosis of a familial pathogenic variant in a proband allows for subsequent cascade screening of their at-risk relatives.

AIMS

We investigated the diagnostic yield from cardiac gene panel testing and reviewed variants of uncertain significance from patients attending three specialist cardiogenetics services in Ireland in the years 2002 to 2020.

RESULTS

Reviewing molecular genetic diagnostic reports of 834 patients from 820 families, the initial diagnostic yield of pathogenic/likely pathogenic variants was 237/834 patients (28.4%), increasing to 276/834 patients (33.1%) following re-evaluation of cases with variant(s) of uncertain significance. Altogether, 42/85 patients with VUS reviewed (49.4%) had a re-classification that could change their clinical management. Females were more likely to carry pathogenic/likely pathogenic variants than males (139/374, 37.2% vs 137/460, 29.8%, respectively, p = 0.03), and the diagnostic yields were highest in the 0 to < 2 years age group (6/12, 50.0%) and amongst those tested for cardiomyopathy gene panels (13/35, 37.1%). Variants in the MYBPC3/MYH7 (87/109, 79.8%) and KCNQ1/KCNH2 (91/100, 91.0%) genes were the predominant genetic causes for hypertrophic cardiomyopathy and long QT syndrome, respectively.

CONCLUSION

Our study highlights the importance of collation and review of pre-ACMG genetic variants to increase diagnostic utility of genetic testing for inherited heart disease. Almost half of patients with pre-ACMG VUS reviewed had their variant re-classified to likely pathogenic/likely benign which resulted in a positive clinical impact for patients and their families.

Multisite Validation of a Functional Assay to Adjudicate SCN5A Brugada Syndrome-Associated Variants

BACKGROUND

Brugada syndrome is an inheritable arrhythmia condition that is associated with rare, loss-of-function variants in SCN5A. Interpreting the pathogenicity of SCN5A missense variants is challenging, and ≈79% of SCN5A missense variants in ClinVar are currently classified as variants of uncertain significance. Automated patch clamp technology enables high-throughput functional studies of ion channel variants and can provide evidence for variant reclassification.

METHODS

An in vitro SCN5A-Brugada syndrome automated patch clamp assay was independently performed at Vanderbilt University Medical Center and Victor Chang Cardiac Research Institute. The assay was calibrated according to ClinGen Sequence Variant Interpretation recommendations using high-confidence variant controls (n=49). Normal and abnormal ranges of function were established based on the distribution of benign variant assay results. Odds of pathogenicity values were derived from the experimental results according to ClinGen Sequence Variant Interpretation recommendations. The calibrated assay was then used to study SCN5A variants of uncertain significance observed in 4 families with Brugada syndrome and other arrhythmia phenotypes associated with SCN5A loss-of-function.

RESULTS

Variant channel parameters generated independently at the 2 research sites showed strong correlations, including peak INa density (R2=0.86). The assay accurately distinguished benign controls (24/25 concordant variants) from pathogenic controls (23/24 concordant variants). Odds of pathogenicity values were 0.042 for normal function and 24.0 for abnormal function, corresponding to strong evidence for both American College of Medical Genetics and Genomics/Association for Molecular Pathology benign and pathogenic functional criteria (BS3 and PS3, respectively). Application of the assay to 4 clinical SCN5A variants of uncertain significance revealed loss-of-function for 3/4 variants, enabling reclassification to likely pathogenic.

CONCLUSIONS

This validated high-throughput assay provides clinical-grade functional evidence to aid the classification of current and future SCN5A-Brugada syndrome variants of uncertain significance.

Unraveling the genetic tapestry of pediatric sarcomeric cardiomyopathies and masquerading phenocopies in Jordan

Pediatric cardiomyopathies are mostly attributed to variants in sarcomere-related genes. Unfortunately, the genetic architecture of pediatric cardiomyopathies has never been previously studied in Jordan. We sought to uncover the genetic landscape of 14 patients from nine families with several subtypes of pediatric cardiomyopathies in Jordan using Exome sequencing (ES). Our investigation identified pathogenic and likely pathogenic variants in seven out of nine families (77.8%), clustering in sarcomere-related genes. Surprisingly, phenocopies of sarcomere-related hypertrophic cardiomyopathies were evident in probands with glycogen storage disorder and mitochondrial-related disease. Our study underscored the significance of streamlining ES or expanding cardiomyopathy-related gene panels to identify plausible phenocopies of sarcomere-related cardiomyopathies. Our findings also pointed out the need for genetic testing in patients with cardiomyopathy and their at-risk family members. This can potentially lead to better management strategies, enabling early interventions, and ultimately enhancing their prognosis. Finally, our findings provide an initial contribution to the currently absent knowledge about the molecular underpinnings of cardiomyopathies in Jordan.

The dominant findings of a recessive man: from Mendel's kid pea to kidney

The research of Mendel, born two centuries ago, still has many direct implications for our everyday clinical work. He introduced the terms "dominant" and "recessive" characters and determined their 3:1 ratio in the offspring of heterozygous "hybrid" plants. This distribution allowed calculation of the number of the phenotype-determining "elements," i.e., the alleles, and has been used ever since to prove the monogenic origin of a disorder. The Mendelian inheritance of monogenic kidney disorders is still of great help in distinguishing them from those with multifactorial origin in clinical practice. Inheritance of most monogenic kidney disorders fits to Mendel's observations: the equal contribution of the two parents and the complete penetrance or the direct correlation between the frequency of the recessive character and the degree of inbreeding. Nevertheless, beyond the truth of these basic concepts, several observations have expanded their genetic characteristics. The extreme genetic heterogeneity, the pleiotropy of the causal genes and the role of modifiers in ciliopathies, the digenic inheritance and parental imprinting in some tubulopathies, and the incomplete penetrance and eventual interallelic interactions in podocytopathies, reflect this expansion. For all these reasons, the transmission pattern in a natural setting may depend not only on the "character" but also on the causal gene and the variant. Mendel's passion for research combined with his modest personality and meticulous approach can still serve as an example in the work required to understand the non-Mendelian universe of genetics.

The role of GPD1L, a sodium channel interacting gene, in the pathogenesis of Brugada Syndrome

Background

Brugada Syndrome (BrS) is an inherited arrhythmia syndrome in which mutations in the cardiac sodium channel SCN5A (NaV1.5) account for approximately 20% of cases. Mutations in sodium channel-modifying genes may account for additional BrS cases, though BrS may be polygenic given common SNPs associated with BrS have been identified. Recent analysis, however, has suggested that SCN5A should be regarded as the sole monogenic cause of BrS.

Objective

We sought to re-assess the genetic underpinnings of BrS in a large mutligenerational family with a putative mutation in GPD1L that affects surface membrane expression of NaV1.5 in vitro.

Methods

Fine linkage mapping was performed in the family using the Illumina Global Screening Array. Whole exome sequencing of the proband was performed to identify rare variants and mutations, and Sanger sequencing was used to assay previously-reported risk single nucleotide polymorphsims (SNPs) for BrS.

Results

Linkage analysis decreased the size of the previously-reported microsatellite linkage region to approximately 3 Mb. GPD1L-A280V was the only coding non-synonymous variation present at less than 1% allele frequency in the proband within the linkage region. No rare non-synonymous variants were present outside the linkage area in affected individuals in genes associated with BrS. Risk SNPs known to predispose to BrS were overrepresented in affected members of the family.

Conclusion

Together, our data suggest GPD1L-A280V remains the most likely cause of BrS in this large multigenerational family. While care should be taken in interpreting variant pathogenicity given the genetic uncertainty of BrS, our data support inclusion of other putative BrS genes in clinical genetic panels.

Mutation Patterns Predict Drug Sensitivity in Acute Myeloid Leukemia

PURPOSE

The inherent genetic heterogeneity of acute myeloid leukemia (AML) has challenged the development of precise and effective therapies. The objective of this study was to elucidate the genomic basis of drug resistance or sensitivity, identify signatures for drug response prediction, and provide resources to the research community.

EXPERIMENTAL DESIGN

We performed targeted sequencing, high-throughput drug screening, and single-cell genomic profiling on leukemia cell samples derived from patients with AML. Statistical approaches and machine learning models were applied to identify signatures for drug response prediction. We also integrated large public datasets to understand the co-occurring mutation patterns and further investigated the mutation profiles in the single cells. The features revealed in the co-occurring or mutual exclusivity pattern were further subjected to machine learning models.

RESULTS

We detected genetic signatures associated with sensitivity or resistance to specific agents, and identified five co-occurring mutation groups. The application of single-cell genomic sequencing unveiled the co-occurrence of variants at the individual cell level, highlighting the presence of distinct subclones within patients with AML. Using the mutation pattern for drug response prediction demonstrates high accuracy in predicting sensitivity to some drug classes, such as MEK inhibitors for RAS-mutated leukemia.

CONCLUSIONS

Our study highlights the importance of considering the gene mutation patterns for the prediction of drug response in AML. It provides a framework for categorizing patients with AML by mutations that enable drug sensitivity prediction.

A Comprehensive Analysis of Non-Desmosomal Rare Genetic Variants in Arrhythmogenic Cardiomyopathy: Integrating in Padua Cohort Literature-Derived Data

Arrhythmogenic cardiomyopathy (ACM) is an inherited myocardial disease at risk of sudden death. Genetic testing impacts greatly in ACM diagnosis, but gene-disease associations have yet to be determined for the increasing number of genes included in clinical panels. Genetic variants evaluation was undertaken for the most relevant non-desmosomal disease genes. We retrospectively studied 320 unrelated Italian ACM patients, including 243 cases with predominant right-ventricular (ARVC) and 77 cases with predominant left-ventricular (ALVC) involvement, who did not carry pathogenic/likely pathogenic (P/LP) variants in desmosome-coding genes. The aim was to assess rare genetic variants in transmembrane protein 43 (TMEM43), desmin (DES), phospholamban (PLN), filamin c (FLNC), cadherin 2 (CDH2), and tight junction protein 1 (TJP1), based on current adjudication guidelines and reappraisal on reported literature data. Thirty-five rare genetic variants, including 23 (64%) P/LP, were identified in 39 patients (16/243 ARVC; 23/77 ALVC): 22 FLNC, 9 DES, 2 TMEM43, and 2 CDH2. No P/LP variants were found in PLN and TJP1 genes. Gene-based burden analysis, including P/LP variants reported in literature, showed significant enrichment for TMEM43 (3.79-fold), DES (10.31-fold), PLN (117.8-fold) and FLNC (107-fold). A non-desmosomal rare genetic variant is found in a minority of ARVC patients but in about one third of ALVC patients; as such, clinical decision-making should be driven by genes with robust evidence. More than two thirds of non-desmosomal P/LP variants occur in FLNC.

Specifications of the ACMG/AMP variant curation guidelines for the analysis of germline ATM sequence variants

The ClinGen Hereditary Breast, Ovarian and Pancreatic Cancer (HBOP) Variant Curation Expert Panel (VCEP) is composed of internationally recognized experts in clinical genetics, molecular biology and variant interpretation. This VCEP made specifications for ACMG/AMP guidelines for the ataxia telangiectasia mutated (ATM) gene according to the Food and Drug Administration (FDA)-approved ClinGen protocol. These gene-specific rules for ATM were modified from the American College of Medical Genetics and Association for Molecular Pathology (ACMG/AMP) guidelines and were tested against 33 ATM variants of various types and classifications in a pilot curation phase. The pilot revealed a majority agreement between the HBOP VCEP classifications and the ClinVar-deposited classifications. Six pilot variants had conflicting interpretations in ClinVar and reevaluation with the VCEP's ATM-specific rules resulted in four that were classified as benign, one as likely pathogenic and one as a variant of uncertain significance (VUS) by the VCEP, improving the certainty of interpretations in the public domain. Overall, 28 the 33 pilot variants were not VUS leading to an 85% classification rate. The ClinGen-approved, modified rules demonstrated value for improved interpretation of variants in ATM.

Improving the Yield of Genetic Diagnosis through Additional Genetic Panel Testing in Hereditary Ophthalmic Diseases

Numerous hereditary ophthalmic diseases display significant genetic diversity. Consequently, the utilization of gene panel sequencing allows a greater number of patients to receive a genetic diagnosis for their clinical manifestations. We investigated how to improve the yield of genetic diagnosis through additional gene panel sequencing in hereditary ophthalmic diseases. A gene panel sequencing consisting of a customized hereditary retinopathy panel or hereditary retinitis pigmentosa (RP) panel was prescribed and referred to a CAP-accredited clinical laboratory. If no significant mutations associated with hereditary retinopathy and RP were detected in either panel, additional gene panel sequencing was requested for research use, utilizing the remaining panel. After additional gene panel sequencing, a total of 16 heterozygous or homozygous variants were identified in 15 different genes associated with hereditary ophthalmic diseases. Of 15 patients carrying any candidate variants, the clinical symptoms could be tentatively accounted for by genetic mutations in seven patients. However, in the remaining eight patients, given the in silico mutation predictive analysis, variant allele frequency in gnomAD, inheritance pattern, and genotype-phenotype correlation, fully elucidating the clinical manifestations with the identified rare variant was challenging. Our study highlights the utility of gene panel sequencing in achieving accurate diagnoses for hereditary ophthalmic diseases and enhancing the diagnostic yield through additional gene panel sequencing. Thus, gene panel sequencing can serve as a primary tool for the genetic diagnosis of hereditary ophthalmic diseases, even in cases where a single genetic cause is suspected. With a deeper comprehension of the genetic mechanisms underlying these diseases, it becomes feasible.

Development of a comprehensive cardiovascular disease genetic risk assessment test

Background

Despite monogenic and polygenic contributions to cardiovascular disease (CVD), genetic testing is not widely adopted, and current tests are limited by the breadth of surveyed conditions and interpretation burden.

Methods

We developed a comprehensive clinical genome CVD test with semi-automated interpretation. Monogenic conditions and risk alleles were selected based on the strength of disease association and evidence for increased disease risk, respectively. Non-CVD secondary findings genes, pharmacogenomic (PGx) variants and CVD polygenic risk scores (PRS) were assessed for inclusion. Test performance was modeled using 2,594 genomes from the 1000 Genomes Project, and further investigated in 20 previously tested individuals.

Results

The CVD genome test is composed of a panel of 215 CVD gene-disease pairs, 35 non-CVD secondary findings genes, 4 risk alleles or genotypes, 10 PGx genes and a PRS for coronary artery disease. Modeling of test performance using samples from the 1000 Genomes Project revealed ~6% of individuals with a monogenic finding in a CVD-associated gene, 6% with a risk allele finding, ~1% with a non-CVD secondary finding, and 93% with CVD-associated PGx variants. Assessment of blinded clinical samples showed complete concordance with prior testing. An average of 4 variants were reviewed per case, with interpretation and reporting time ranging from 9-96 min.

Conclusions

A genome sequencing based CVD genetic risk assessment can provide comprehensive genetic disease and genetic risk information to patients with CVD. The semi-automated and limited interpretation burden suggest that this testing approach could be scaled to support population-level initiatives.

Cardiac arrhythmia and epilepsy genetic variants in sudden unexpected death in epilepsy

Introduction

Sudden Unexpected Death in Epilepsy (SUDEP) is the leading epilepsy-related cause of death, affecting approximately 1 per 1,000 individuals with epilepsy per year. Genetic variants that affect autonomic function, such as genes associated with cardiac arrhythmias, may predispose people with epilepsy to greater risk of both sudden cardiac death and SUDEP. Advances in next generation sequencing allow for the exploration of gene variants as potential biomarkers.

Methods

Genetic testing for the presence of cardiac arrhythmia and epilepsy gene variants was performed via genetic panels in 39 cases of SUDEP identified via autopsy by the Ontario Forensic Pathology Service. Variants were summarized by in-silico evidence for pathogenicity from 4 algorithms (SIFT, PolyPhen-2, PROVEAN, Mutation Taster) and allele frequencies in the general population (GnomAD). A maximum credible population allele frequency of 0.00004 was calculated based on epilepsy prevalence and SUDEP incidence to assess whether a variant was compatible with a pathogenic interpretation.

Results

Median age at the time of death was 33.3 years (range: 2, 60). Fifty-nine percent (n=23) were male. Gene panels detected 62 unique variants in 45 genes: 19 on the arrhythmia panel and 26 on the epilepsy panel. At least one variant was identified in 28 (72%) of decedents. Missense mutations comprised 57 (92%) of the observed variants. At least three in silico models predicted 12 (46%) cardiac arrhythmia panel missense variants and 20 (65%) epilepsy panel missense variants were pathogenic. Population allele frequencies were <0.00004 for 11 (42%) of the cardiac variants and 10 (32%) of the epilepsy variants. Together, these metrics identified 13 SUDEP variants of interest.

Discussion

Nearly three-quarters of decedents in this SUDEP cohort carried variants in comprehensive epilepsy or cardiac arrhythmia gene panels, with more than a third having variants in both panels. The proportion of decedents with cardiac variants aligns with recent studies of the disproportionate cardiac burden the epilepsy community faces compared to the general population and suggests a possible cardiac contribution to epilepsy mortality. These results identified 13 priority targets for future functional studies of these genes potential role in sudden death and demonstrates the necessity for further exploration of potential genetic contributions to SUDEP.

ClinGen variant curation expert panel recommendations for classification of variants in GAMT, GATM and SLC6A8 for cerebral creatine deficiency syndromes

Cerebral creatine deficiency syndromes (CCDS) are inherited metabolic phenotypes of creatine synthesis and transport. There are two enzyme deficiencies, guanidinoacetate methyltransferase (GAMT), encoded by GAMT and arginine-glycine amidinotransferase (AGAT), encoded by GATM, which are involved in the synthesis of creatine. After synthesis, creatine is taken up by a sodium-dependent membrane bound creatine transporter (CRTR), encoded by SLC6A8, into all organs. Creatine uptake is very important especially in high energy demanding organs such as the brain, and muscle. To classify the pathogenicity of variants in GAMT, GATM, and SLC6A8, we developed the CCDS Variant Curation Expert Panel (VCEP) in 2018, supported by The Clinical Genome Resource (ClinGen), a National Institutes of Health (NIH)-funded resource. We developed disease-specific variant classification guidelines for GAMT-, GATM-, and SLC6A8-related CCDS, adapted from the American College of Medical Genetics/Association of Molecular Pathology (ACMG/AMP) variant interpretation guidelines. We applied specific variant classification guidelines to 30 pilot variants in each of the three genes that have variants associated with CCDS. Our CCDS VCEP was approved by the ClinGen Sequence Variant Interpretation Working Group (SVI WG) and Clinical Domain Oversight Committee in July 2022. We curated 181 variants including 72 variants in GAMT, 45 variants in GATM, and 64 variants in SLC6A8 and submitted these classifications to ClinVar, a public variant database supported by the National Center for Biotechnology Information. Missense variants were the most common variant type in all three genes. We submitted 32 new variants and reclassified 34 variants with conflicting interpretations. We report specific phenotype (PP4) using a points system based on the urine and plasma guanidinoacetate and creatine levels, brain magnetic resonance spectroscopy (MRS) creatine level, and enzyme activity or creatine uptake in fibroblasts ranging from PP4, PP4_Moderate and PP4_Strong. Our CCDS VCEP is one of the first panels applying disease specific variant classification algorithms for an X-linked disease. The availability of these guidelines and classifications can guide molecular genetics and genomic laboratories and health care providers to assess the molecular diagnosis of individuals with a CCDS phenotype.

The importance of variant reinterpretation in inherited cardiovascular diseases: Establishing the optimal timeframe

Inherited cardiovascular diseases are rare diseases that are difficult to diagnose by non-expert professionals. Genetic analyses play a key role in the diagnosis of these diseases, in which the identification of a pathogenic genetic variant is often a diagnostic criterion. Therefore, genetic variant classification and routine reinterpretation as data become available represent one of the main challenges associated with genetic analyses. Using the genetic variants identified in an inherited cardiovascular diseases unit during a 10-year period, the objectives of this study were: 1) to evaluate the impact of genetic variant reinterpretation, 2) to compare the reclassification rates between different cohorts of cardiac channelopathies and cardiomyopathies, and 3) to establish the most appropriate periodicity for genetic variant reinterpretation. All the evaluated cohorts (full cohort of inherited cardiovascular diseases, cardiomyopathies, cardiac channelopathies, hypertrophic cardiomyopathy, dilated cardiomyopathy, arrhythmogenic cardiomyopathy, Brugada syndrome, long QT syndrome and catecholaminergic polymorphic ventricular tachycardia) showed reclassification rates above 25%, showing even higher reclassification rates when there is definitive evidence of the association between the gene and the disease in the cardiac channelopathies. Evaluation of genetic variant reclassification rates based on the year of the initial classification showed that the most appropriate frequency for the reinterpretation would be 2 years, with the possibility of a more frequent reinterpretation if deemed convenient. To keep genetic variant classifications up to date, genetic counsellors play a critical role in the reinterpretation process, providing clinical evidence that genetic diagnostic laboratories often do not have at their disposal and communicating changes in classification and the potential implications of these reclassifications to patients and relatives.

Clinical interpretation of KCNH2 variants using a robust PS3/BS3 functional patch-clamp assay

Long QT syndrome (LQTS), caused by the dysfunction of cardiac ion channels, increases the risk of sudden death in otherwise healthy young people. For many variants in LQTS genes, there is insufficient evidence to make a definitive genetic diagnosis. We have established a robust functional patch-clamp assay to facilitate classification of missense variants in KCNH2, one of the key LQTS genes. A curated set of 30 benign and 30 pathogenic missense variants were used to establish the range of normal and abnormal function. The extent to which variants reduced protein function was quantified using Z scores, the number of standard deviations from the mean of the normalized current density of the set of benign variant controls. A Z score of -2 defined the threshold for abnormal loss of function, which corresponds to 55% wild-type function. More extreme Z scores were observed for variants with a greater loss-of-function effect. We propose that the Z score for each variant can be used to inform the application and weighting of abnormal and normal functional evidence criteria (PS3 and BS3) within the American College of Medical Genetics and Genomics variant classification framework. The validity of this approach was demonstrated using a series of 18 KCNH2 missense variants detected in a childhood onset LQTS cohort, where the level of function assessed using our assay correlated to the Schwartz score (a scoring system used to quantify the probability of a clinical diagnosis of LQTS) and the length of the corrected QT (QTc) interval.

Challenges and approaches to calibrating patient phenotype as evidence for cancer gene variant classification under ACMG/AMP guidelines

Since first publication of the American College of Medical Genetics and Genomics/Association for Medical Pathology (ACMG/AMP) variant classification guidelines, additional recommendations for application of certain criteria have been released (https://clinicalgenome.org/docs/), to improve their application in the diagnostic setting. However, none have addressed use of the PS4 and PP4 criteria, capturing patient presentation as evidence towards pathogenicity. Application of PS4 can be done through traditional case-control studies, or "proband counting" within or across clinical testing cohorts. Review of the existing PS4 and PP4 specifications for Hereditary Cancer Gene Variant Curation Expert Panels revealed substantial differences in the approach to defining specifications. Using BRCA1, BRCA2 and TP53 as exemplar genes, we calibrated different methods proposed for applying the "PS4 proband counting" criterion. For each approach, we considered limitations, non-independence with other ACMG/AMP criteria, broader applicability, and variability in results for different datasets. Our findings highlight inherent overlap of proband-counting methods with ACMG/AMP frequency codes, and the importance of calibration to derive dataset-specific code weights that can account for potential between-dataset differences in ascertainment and other factors. Our work emphasizes the advantages and generalizability of logistic regression analysis over simple proband-counting approaches to empirically determine the relative predictive capacity and weight of various personal clinical features in the context of multigene panel testing, for improved variant interpretation. We also provide a general protocol, including instructions for data formatting and a web-server for analysis of personal history parameters, to facilitate dataset-specific calibration analyses required to use such data for germline variant classification.

Role of TBX20 Truncating Variants in Dilated Cardiomyopathy and Left Ventricular Noncompaction

BACKGROUND

Less than 40% of patients with dilated cardiomyopathy (DCM) have a pathogenic/likely pathogenic genetic variant identified. TBX20 has been linked to congenital heart defects; although an association with left ventricular noncompaction (LVNC) and DCM has been proposed, it is still considered a gene with limited evidence for these phenotypes. This study sought to investigate the association between the TBX20 truncating variant (TBX20tv) and DCM/LVNC.

METHODS

TBX20 was sequenced by next-generation sequencing in 7463 unrelated probands with a diagnosis of DCM or LVNC, 22 773 probands of an internal comparison group (hypertrophic cardiomyopathy, channelopathies, or aortic diseases), and 124 098 external controls (individuals from the gnomAD database). Enrichment of TBX20tv in DCM/LVNC was calculated, cosegregation was determined in selected families, and clinical characteristics and outcomes were analyzed in carriers.

RESULTS

TBX20tv was enriched in DCM/LVNC (24/7463; 0.32%) compared with internal (1/22 773; 0.004%) and external comparison groups (4/124 098; 0.003%), with odds ratios of 73.23 (95% CI, 9.90-541.45; P<0.0001) and 99.76 (95% CI, 34.60-287.62; P<0.0001), respectively. TBX20tv was cosegregated with DCM/LVNC phenotype in 21 families for a combined logarythm of the odds score of 4.53 (strong linkage). Among 57 individuals with TBX20tv (49.1% men; mean age, 35.9±20.8 years), 41 (71.9%) exhibited DCM/LVNC, of whom 14 (34.1%) had also congenital heart defects. After a median follow-up of 6.9 (95% CI, 25-75:3.6-14.5) years, 9.7% of patients with DCM/LVNC had end-stage heart failure events and 4.8% experienced malignant ventricular arrhythmias.

CONCLUSIONS

TBX20tv is associated with DCM/LVNC; congenital heart defect is also present in around one-third of cases. TBX20tv-associated DCM/LVNC is characterized by a nonaggressive phenotype, with a low incidence of major cardiovascular events. TBX20 should be considered a definitive gene for DCM and LVNC and routinely included in genetic testing panels for these phenotypes.

Prevalence of common autosomal recessive mutation carriers in women in the Southern Vietnam following the application of expanded carrier screening

The common autosomal recessive (AR) mutation carrier is still unknown in Vietnam. This study aims to identify the most common AR gene mutation carriers in women of reproductive age to build a Vietnamese-specific carrier screening panel for AR and X-linked disorders in the preconception and prenatal healthcare program. A cross-sectional study was conducted at University Medical Center-Branch 2 in Ho Chi Minh City from December 1st, 2020, to June 30th, 2023. 338 women have consented to take a 5 mL blood test to identify 540 recessive genes. The carrier screening panel was designed based on the American College of Medical Genetics and Genomics (ACMG)-recommended genes and suggestions from 104 clinical experts in Vietnam. Obstetricians and genetic experts counseled all positive testing results to discuss the possibility of recessive diseases in their offspring. The most common recessive disorders were defined at a prevalence of 1 in 60 or greater, and those were added to a Vietnamese-specific carrier screening panel. 338 non-pregnant and pregnant women underwent the expanded carrier screening (ECS). The carrier frequency was 63.6%, in which 215 women carried at least one AR gene mutation. GJB2 hearing impairment was identified as the most common chronic condition (1 in 5). The second most common AR disorder was beta-thalassemia (1 in 16), followed by cystic fibrosis (1 in 23), G6PD deficiency (1 in 28), Wilson's disease (1 in 31), Usher's syndrome (1 in 31), and glycogen storage disease (1 in 56). Seven common recessive genes were added in ethnic-based carrier screening. Women in the South of Vietnam have been carried for many recessive conditions at high frequency, such as hearing impairment, genetic anemia, and cystic fibrosis. It is necessary to implement a preconception and prenatal screening program by using seven widely popular AR genes in a Vietnamese-specific carrier screening panel to reduce the burden related to AR and X-linked disorders.

Genome Sequencing of Idiopathic Speech Delay

Genetic investigations of people with speech and language disorders can provide windows into key aspects of human biology. Most genomic research into impaired speech development has so far focused on childhood apraxia of speech (CAS), a rare neurodevelopmental disorder characterized by difficulties with coordinating rapid fine motor sequences that underlie proficient speech. In 2001, pathogenic variants of FOXP2 provided the first molecular genetic accounts of CAS aetiology. Since then, disruptions in several other genes have been implicated in CAS, with a substantial proportion of cases being explained by high-penetrance variants. However, the genetic architecture underlying other speech-related disorders remains less well understood. Thus, in the present study, we used systematic DNA sequencing methods to investigate idiopathic speech delay, as characterized by delayed speech development in the absence of a motor speech diagnosis (such as CAS), a language/reading disorder, or intellectual disability. We performed genome sequencing in a cohort of 23 children with a rigorous diagnosis of idiopathic speech delay. For roughly half of the sample (ten probands), sufficient DNA was also available for genome sequencing in both parents, allowing discovery of de novo variants. In the thirteen singleton probands, we focused on identifying loss-of-function and likely damaging missense variants in genes intolerant to such mutations. We found that one speech delay proband carried a pathogenic frameshift deletion in SETD1A, a gene previously implicated in a broader variable monogenic syndrome characterized by global developmental problems including delayed speech and/or language development, mild intellectual disability, facial dysmorphisms, and behavioural and psychiatric symptoms. Of note, pathogenic SETD1A variants have been independently reported in children with CAS in two separate studies. In other probands in our speech delay cohort, likely pathogenic missense variants were identified affecting highly conserved amino acids in key functional domains of SPTBN1 and ARF3. Overall, this study expands the phenotype spectrum associated with pathogenic SETD1A variants, to also include idiopathic speech delay without CAS or intellectual disability, and suggests additional novel potential candidate genes that may harbour high-penetrance variants that can disrupt speech development.

Differential inclusion of NEB exons 143 and 144 provides insight into NEB-related myopathy variant interpretation and disease manifestation

Biallelic pathogenic variants in the gene encoding nebulin (NEB) are a known cause of congenital myopathy. We present two individuals with congenital myopathy and compound heterozygous variants (NM_001271208.2: c.2079C>A; p.(Cys693Ter) and c.21522+3A>G ) in NEB. Transcriptomic sequencing on patient muscle revealed that the extended splice variant c.21522+3A>G causes exon 144 skipping. Nebulin isoforms containing exon 144 are known to be mutually exclusive with isoforms containing exon 143, and these isoforms are differentially expressed during development and in adult skeletal muscles. Patients MRIs were compared to the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these patients better fits the distribution of exon 144-containing isoforms in muscle than with previously published MRI findings in NEB-related disease due to other variants. To our knowledge this is the first report hypothesizing disease pathogenesis through the alteration of isoform distributions in muscle.

A Comprehensive Analysis of 3 Moroccan Genomes Revealed Contributions From Both African and European Ancestries

Genetic variations in the human genome represent the differences in DNA sequence within individuals. This highlights the important role of whole human genome sequencing which has become the keystone for precision medicine and disease prediction. Morocco is an important hub for studying human population migration and mixing history. This study presents the analysis of 3 Moroccan genomes; the variant analysis revealed 6 379 606 single nucleotide variants (SNVs) and 1 050 577 small InDels. Of those identified SNVs, 219 152 were novel, with 1233 occurring in coding regions, and 5580 non-synonymous single nucleotide variants (nsSNP) variants were predicted to affect protein functions. The InDels produced 1055 coding variants and 454 non-3n length variants, and their size ranged from -49 and 49 bp. We further analysed the gene pathways of 8 novel coding variants found in the 3 genomes and revealed 5 genes involved in various diseases and biological pathways. We found that the Moroccan genomes share 92.78% of African ancestry, and 92.86% of Non-Finnish European ancestry, according to the gnomAD database. Then, population structure inference, by admixture analysis and network-based approach, revealed that the studied genomes form a mixed population structure, highlighting the increased genetic diversity in Morocco.

Gene-specific ACMG/AMP classification criteria for germline APC variants: Recommendations from the ClinGen InSiGHT Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel

PURPOSE

The Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel (VCEP) was established by the International Society for Gastrointestinal Hereditary Tumours and the Clinical Genome Resource, who set out to develop recommendations for the interpretation of germline APC variants underlying Familial Adenomatous Polyposis, the most frequent hereditary polyposis syndrome.

METHODS

Through a rigorous process of database analysis, literature review, and expert elicitation, the APC VCEP derived gene-specific modifications to the ACMG/AMP (American College of Medical Genetics and Genomics and Association for Molecular Pathology) variant classification guidelines and validated such criteria through the pilot classification of 58 variants.

RESULTS

The APC-specific criteria represented gene- and disease-informed specifications, including a quantitative approach to allele frequency thresholds, a stepwise decision tool for truncating variants, and semiquantitative evaluations of experimental and clinical data. Using the APC-specific criteria, 47% (27/58) of pilot variants were reclassified including 14 previous variants of uncertain significance (VUS).

CONCLUSION

The APC-specific ACMG/AMP criteria preserved the classification of well-characterized variants on ClinVar while substantially reducing the number of VUS by 56% (14/25). Moving forward, the APC VCEP will continue to interpret prioritized lists of VUS, the results of which will represent the most authoritative variant classification for widespread clinical use.

SIGMA leverages protein structural information to predict the pathogenicity of missense variants

Leveraging protein structural information to evaluate pathogenicity has been hindered by the scarcity of experimentally determined 3D protein. With the aid of AlphaFold2 predictions, we developed the structure-informed genetic missense mutation assessor (SIGMA) to predict missense variant pathogenicity. In comparison with existing predictors across labeled variant datasets and experimental datasets, SIGMA demonstrates superior performance in predicting missense variant pathogenicity (AUC = 0.933). We found that the relative solvent accessibility of the mutated residue contributed greatly to the predictive ability of SIGMA. We further explored combining SIGMA with other top-tier predictors to create SIGMA+, proving highly effective for variant pathogenicity prediction (AUC = 0.966). To facilitate the application of SIGMA, we pre-computed SIGMA scores for over 48 million possible missense variants across 3,454 disease-associated genes and developed an interactive online platform (https://www.sigma-pred.org/). Overall, by leveraging protein structure information, SIGMA offers an accurate structure-based approach to evaluating the pathogenicity of missense variants.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Rare disease gene association discovery from burden analysis of the 100,000 Genomes Project data

To discover rare disease-gene associations, we developed a gene burden analytical framework and applied it to rare, protein-coding variants from whole genome sequencing of 35,008 cases with rare diseases and their family members recruited to the 100,000 Genomes Project (100KGP). Following in silico triaging of the results, 88 novel associations were identified including 38 with existing experimental evidence. We have published the confirmation of one of these associations, hereditary ataxia with UCHL1 , and independent confirmatory evidence has recently been published for four more. We highlight a further seven compelling associations: hypertrophic cardiomyopathy with DYSF and SLC4A3 where both genes show high/specific heart expression and existing associations to skeletal dystrophies or short QT syndrome respectively; monogenic diabetes with UNC13A with a known role in the regulation of β cells and a mouse model with impaired glucose tolerance; epilepsy with KCNQ1 where a mouse model shows seizures and the existing long QT syndrome association may be linked; early onset Parkinson's disease with RYR1 with existing links to tremor pathophysiology and a mouse model with neurological phenotypes; anterior segment ocular abnormalities associated with POMK showing expression in corneal cells and with a zebrafish model with developmental ocular abnormalities; and cystic kidney disease with COL4A3 showing high renal expression and prior evidence for a digenic or modifying role in renal disease. Confirmation of all 88 associations would lead to potential diagnoses in 456 molecularly undiagnosed cases within the 100KGP, as well as other rare disease patients worldwide, highlighting the clinical impact of a large-scale statistical approach to rare disease gene discovery.

Ethnicity, consanguinity, and genetic architecture of hypertrophic cardiomyopathy

AIMS

Hypertrophic cardiomyopathy (HCM) is characterized by phenotypic heterogeneity that is partly explained by the diversity of genetic variants contributing to disease. Accurate interpretation of these variants constitutes a major challenge for diagnosis and implementing precision medicine, especially in understudied populations. The aim is to define the genetic architecture of HCM in North African cohorts with high consanguinity using ancestry-matched cases and controls.

METHODS AND RESULTS

Prospective Egyptian patients (n = 514) and controls (n = 400) underwent clinical phenotyping and genetic testing. Rare variants in 13 validated HCM genes were classified according to standard clinical guidelines and compared with a prospective HCM cohort of majority European ancestry (n = 684). A higher prevalence of homozygous variants was observed in Egyptian patients (4.1% vs. 0.1%, P = 2 × 10-7), with variants in the minor HCM genes MYL2, MYL3, and CSRP3 more likely to present in homozygosity than the major genes, suggesting these variants are less penetrant in heterozygosity. Biallelic variants in the recessive HCM gene TRIM63 were detected in 2.1% of patients (five-fold greater than European patients), highlighting the importance of recessive inheritance in consanguineous populations. Finally, rare variants in Egyptian HCM patients were less likely to be classified as (likely) pathogenic compared with Europeans (40.8% vs. 61.6%, P = 1.6 × 10-5) due to the underrepresentation of Middle Eastern populations in current reference resources. This proportion increased to 53.3% after incorporating methods that leverage new ancestry-matched controls presented here.

CONCLUSION

Studying consanguineous populations reveals novel insights with relevance to genetic testing and our understanding of the genetic architecture of HCM.