Reclassifying variations of unknown significance in diseases affecting Saudi Arabia's population reveal new associations.

Introduction: Physicians face diagnostic dilemmas upon reports indicating disease variants of unknown significance (VUS). The most puzzling cases are patients with rare diseases, where finding another matched genotype and phenotype to associate their results is challenging. This study aims to prove the value of updating patient files with new classifications, potentially leading to better assessment and prevention. Methodology: We recruited retrospective phenotypic and genotypic data from King Saud Medical City, Riyadh, Kingdom of Saudi Arabia. Between September 2020 and December 2021, 1,080 patients' genetic profiles were tested in a College of American Pathologists accredited laboratory. We excluded all confirmed pathogenic variants, likely pathogenic variants and copy number variations. Finally, we further reclassified 194 VUS using different local and global databases, employing in silico prediction to justify the phenotype-genotype association. Results: Of the 194 VUS, 90 remained VUS, and the other 104 were reclassified as follows: 16 pathogenic, 49 likely pathogenic, nine benign, and 30 likely benign. Moreover, most of these variants had never been observed in other local or international databases. Conclusion: Reclassifying the VUS adds value to understanding the causality of the phenotype if it has been reported in another family or population. The healthcare system should establish guidelines for re-evaluating VUS, and upgrading VUS should reflect on individual/family risks and management strategies.

Druggable proteins influencing cardiac structure and function: Implications for heart failure therapies and cancer cardiotoxicity.

Dysfunction of either the right or left ventricle can lead to heart failure (HF) and subsequent morbidity and mortality. We performed a genome-wide association study (GWAS) of 16 cardiac magnetic resonance (CMR) imaging measurements of biventricular function and structure. Cis-Mendelian randomization (MR) was used to identify plasma proteins associating with CMR traits as well as with any of the following cardiac outcomes: HF, non-ischemic cardiomyopathy, dilated cardiomyopathy (DCM), atrial fibrillation, or coronary heart disease. In total, 33 plasma proteins were prioritized, including repurposing candidates for DCM and/or HF: IL18R (providing indirect evidence for IL18), I17RA, GPC5, LAMC2, PA2GA, CD33, and SLAF7. In addition, 13 of the 25 druggable proteins (52%; 95% confidence interval, 0.31 to 0.72) could be mapped to compounds with known oncological indications or side effects. These findings provide leads to facilitate drug development for cardiac disease and suggest that cardiotoxicities of several cancer treatments might represent mechanism-based adverse effects.

Cardiac MRI to guide heart failure and atrial fibrillation drug discovery: a Mendelian randomization analysis.

Background: drug development and disease prevention of heart failure (HF) and atrial fibrillation (AF) are impeded by a lack of robust early-stage surrogates. We determined to what extent cardiac magnetic resonance (CMR) measurements act as surrogates for the development of HF or AF in healthy individuals. Methods: Genetic data was sourced on the association with 22 atrial and ventricular CMR measurements. Mendelian randomization was used to determine CMR associations with atrial fibrillation (AF), heart failure (HF), non-ischemic cardiomyopathy (CMP), and dilated cardiomyopathy (DCM). Additionally, for the CMR surrogates of AF and HF, we explored their association with non-cardiac traits. Results: In total we found that 9 CMR measures were associated with the development of HF, 7 with development of non-ischemic CMR, 6 with DCM, and 12 with AF. biventricular ejection fraction (EF), biventricular or end-systolic volumes (ESV) and left-ventricular (LV) end diastolic volume (EDV) were associated with all 4 cardiac outcomes. Increased LV-MVR (mass to volume ratio) affected HF (odds ratio (OR) 0.83, 95%CI 0.79; 0.88), and DCM (OR 0.26, 95%CI 0.20; 0.34. We were able to identify 9 CMR surrogates for HF and/or AF (including LV-MVR, biventricular EDV, ESV, and right-ventricular EF) which associated with non-cardiac traits such as blood pressure, lung function traits, BMI, cardioembolic stroke, and late-onset Alzheimer's disease. Conclusion: CMR measurements may act as surrogate endpoints for the development of HF (including non-ischemic CMP and DCM) or AF. Additionally, we show that changes in cardiac function and structure measured through CMR, may affect diseases of other organs leading to lung disease or late-onset Alzheimer's disease.

LoFTK: a framework for fully automated calculation of predicted Loss-of-Function variants and genes.

BACKGROUND: Loss-of-Function (LoF) variants in human genes are important due to their impact on clinical phenotypes and frequent occurrence in the genomes of healthy individuals. The association of LoF variants with complex diseases and traits may lead to the discovery and validation of novel therapeutic targets. Current approaches predict high-confidence LoF variants without identifying the specific genes or the number of copies they affect. Moreover, there is a lack of methods for detecting knockout genes caused by compound heterozygous (CH) LoF variants. RESULTS: We have developed the Loss-of-Function ToolKit (LoFTK), which allows efficient and automated prediction of LoF variants from genotyped, imputed and sequenced genomes. LoFTK enables the identification of genes that are inactive in one or two copies and provides summary statistics for downstream analyses. LoFTK can identify CH LoF variants, which result in LoF genes with two copies lost. Using data from parents and offspring we show that 96% of CH LoF genes predicted by LoFTK in the offspring have the respective alleles donated by each parent. CONCLUSIONS: LoFTK is a command-line based tool that provides a reliable computational workflow for predicting LoF variants from genotyped and sequenced genomes, identifying genes that are inactive in 1 or 2 copies. LoFTK is an open software and is freely available to non-commercial users at https://github.com/CirculatoryHealth/LoFTK .

Health economic burden of COVID-19 in Saudi Arabia.

Background: The coronavirus disease 2019 (COVID-19) pandemic has added a massive economic burden on health care systems worldwide. Saudi Arabia is one of the numerous countries that have been economically affected by this pandemic. The objective of this study was to provide real-world data on the health economic burden of COVID-19 on the Saudi health sector and assess the direct medical costs associated with the management of COVID-19. Methods: A retrospective cohort study was conducted based on data collected from patients hospitalized with COVID-19 across 10 institutions in eight different regions in Saudi Arabia. The study calculated the direct medical costs of all cases during the study period by using SAS statistical analysis software. These costs included costs directly related to medical services, such as the health care treatment, hospital stays, laboratory investigations, treatment, outcome, and other related care. Results: A total of 5,286 adult patients admitted with COVID-19 during the study period were included in the study. The average age of the patients was 54 years, and the majority were male (79%). Among the COVID-19 patients hospitalized in a general ward, the median hospital length of stay was 5.5 days (mean: 9.18 days), while the ICU stay was 4.2 days (mean: 7.94 days). The total medical costs for general ward and ICU patients were US$ 38,895 and US$ 24,207,296.9, respectively. The total laboratory investigations ranked as the highest-cost services US$ 588,975 followed by treatment US$ 3,886,509.8. Overall, the total cost of all medical services for patients hospitalized with COVID-19 was US$ 51,572,393.4. Conclusion: This national study found that COVID-19 was not only a serious concern for patients but also a serious economic burden on the health care system in Saudi Arabia.

Prevalence and Disease Expression of Pathogenic and Likely Pathogenic Variants Associated With Inherited Cardiomyopathies in the General Population.

BACKGROUND: Pathogenic and likely pathogenic variants associated with arrhythmogenic right ventricular cardiomyopathy (ARVC), dilated cardiomyopathy (DCM), and hypertrophic cardiomyopathy (HCM) are recommended to be reported as secondary findings in genome sequencing studies. This provides opportunities for early diagnosis, but also fuels uncertainty in variant carriers (G+), since disease penetrance is incomplete. We assessed the prevalence and disease expression of G+ in the general population. METHODS: We identified pathogenic and likely pathogenic variants associated with ARVC, DCM and/or HCM in 200 643 UK Biobank individuals, who underwent whole exome sequencing. We calculated the prevalence of G+ and analyzed the frequency of cardiomyopathy/heart failure diagnosis. In undiagnosed individuals, we analyzed early signs of disease expression using available electrocardiography and cardiac magnetic resonance imaging data. RESULTS: We found a prevalence of 1:578, 1:251, and 1:149 for pathogenic and likely pathogenic variants associated with ARVC, DCM and HCM respectively. Compared with controls, cardiovascular mortality was higher in DCM G+ (odds ratio 1.67 [95% CI 1.04; 2.59], P=0.030), but similar in ARVC and HCM G+ (P≥0.100). Cardiomyopathy or heart failure diagnosis were more frequent in DCM G+ (odds ratio 3.66 [95% CI 2.24; 5.81], P=4.9×10-7) and HCM G+ (odds ratio 3.03 [95% CI 1.98; 4.56], P=5.8×10-7), but comparable in ARVC G+ (P=0.172). In contrast, ARVC G+ had more ventricular arrhythmias (P=3.3×10-4). In undiagnosed individuals, left ventricular ejection fraction was reduced in DCM G+ (P=0.009). CONCLUSIONS: In the general population, pathogenic and likely pathogenic variants associated with ARVC, DCM, or HCM are not uncommon. Although G+ have increased mortality and morbidity, disease penetrance in these carriers from the general population remains low (1.2-3.1%). Follow-up decisions in case of incidental findings should not be based solely on a variant, but on multiple factors, including family history and disease expression.

Identification of CSF3R Mutations in B-Lineage Acute Lymphoblastic Leukemia Using Comprehensive Cancer Panel and Next-Generation Sequencing.

B-lineage acute lymphocytic leukemia (B-ALL) is characterized by different genetic aberrations at a chromosomal and gene level which are very crucial for diagnosis, prognosis and risk assessment of the disease. However, there is still controversial arguments in regard to disease outcomes in specific genetic abnormalities, e.g., 9p-deletion. Moreover, in absence of cytogenetic abnormalities it is difficult to predict B-ALL progression. Here, we use the advantage of Next-generation sequencing (NGS) technology to study the mutation landscape of 12 patients with B-ALL using Comprehensive Cancer Panel (CCP) which covers the most common mutated cancer genes. Our results describe new mutations in CSF3R gene including S661N, S557G, and Q170X which might be associated with disease progression.

Targeted SLC19A3 gene sequencing of 3000 Saudi newborn: a pilot study toward newborn screening.

BACKGROUND: Biotin-thiamine-responsive basal ganglia disease (BTBGD) is an autosomal recessive neurometabolic disorder mostly presented in children. The disorder is described as having subacute encephalopathy with confusion, dystonia, and dysarthria triggered by febrile illness that leads to neuroregression and death if untreated. Using biotin and thiamine at an early stage of the disease can lead to significant improvement. METHODS: BTBGD is a treatable disease if diagnosed at an early age and has been frequently reported in Saudi population. Keeping this in mind, the current study screened 3000 Saudi newborns for the SLC19A3 gene mutations using target sequencing, aiming to determine the carrier frequency in Saudi Population and whether BTBGD is a good candidate to be included in the newborn-screened disorders. RESULTS: Using targeted gene sequencing, DNA from 3000 newborns Saudi was screened for the SLC19A3 gene mutations using standard methods. Screening of the SLC19A3 gene revealed a previously reported heterozygous missense mutation (c.1264A>G (p.Thr422Ala) in six unrelated newborns. No probands having homozygous pathogenic mutations were found in the studied cohort. The variant has been frequently reported previously in homozygous state in Saudi population, making it a hot spot mutation. The current study showed that the carrier frequency of SLC19A3 gene mutation is 1 of 500 in Saudi newborns. CONCLUSION: For the first time in the literature, we determined the carrier frequency of SLC19A3 gene mutation in Saudi population. The estimated prevalence is too rare in Saudi population (at least one in million); therefore, the data are not in favor of including such very rare disorders in newborn screening program at population level. However, a larger cohort is needed for a more accurate estimate.

Atypical influenza A(H1N1)pdm09 strains caused an influenza virus outbreak in Saudi Arabia during the 2009-2011 pandemic season.

BACKGROUND: The triple assortment influenza A(H1N1) virus emerged in spring 2009 and disseminated worldwide, including Saudi Arabia. This study was carried out to characterize Saudi influenza isolates in relation to the global strains and to evaluate the potential role of mutated residues in transmission, adaptation, and the pathogenicity of the virus. METHODS: Nasopharyngeal samples (n = 6492) collected between September 2009 to March 2011 from patients with influenza-like illness were screened by PCR for influenza A(H1N1). Phylogenetic and Molecular evolutionary analysis were carried out to place the Saudi strains in relation to the global strains followed by Mutation analysis of surface and internal proteins. RESULTS: Concatenated whole-genome phylogenetic analysis along with hemagglutinin (HA) signature changes, that is, Aspartic Acid (D) at position 187, P83S, S203T, and R223Q confirmed that the Saudi strains belong to the antigenic category of A/California/07/2009. However, phylogenetic analysis revealed unusual strains of A(H1N1) circulating in Saudi Arabia, not belonging to any of known clades, appearing in five distinct groups well supported by group-specific mutations and novel mutation complexes. These cases had characteristic inter- and intragroup substitution patterns while few of their closest matches showed up as sporadic cases the world over. Specific mutation patterns were detected within the functional domains of internal proteins PB2, PB1, PA, NP, NS1, and M2 having a putative role in viral fitness and virulence. Bayesian coalescent MCMC analysis revealed that Saudi strains belonged to cluster 2 of A(H1N1)pdm09 and spread a month later as compared to other strains of this cluster. CONCLUSION: Influenza outbreak in Saudi Arabia during 2009-2011 was caused by atypical strains of influenza A(H1N1)pdm09, probably introduced in this community on multiple occasions. To understand the antigenic significance of these novel point mutations and mutation complexes require functional studies, which will be crucial for risk assessment of emergent strains and defining infection control measures.

Analysis of CCR5 gene polymorphisms in 321 healthy Saudis using Next Generation Sequencing.

AIMS: To investigate the extent of CCR5 polymorphism in the healthy Saudi population. METHOD: A total of 321 healthy Saudi individuals were sequenced using the ion Ampliseq™ Exome kit (Life Technologies, USA) on genomic DNA following manufacturer's protocol. Whole Exome Sequencing (WES) reads were aligned to the human reference genome (hg19 build) with Torrent Suite Software (v5.0.2) and the variants were called using the Torrent Variant Caller plugin (v5.0) and imported into Ion Reporter Server (v5.0) for the annotation. CCR5 coding exons variants were filtered and checked against the NHLBI GO Exome Sequencing Project (NHLBI), NCBI Reference dbSNPs database, 1000 genomes and Exome Aggregation Consortium datasets (ExAC). RESULTS: A total of 475 variants were identified. Table 1 shows polymorphisms/mutations detected within exons that introduced an amino acid change, deletion or copy number variants (CNV). Three mutations are predicted to influence CCR5 function, including the 32bp deletion (Rs333). Four polymorphisms were detected, plus two CNV. CONCLUSIONS: This is the first report on sequencing the full CCR5 gene using NGS in the Saudi population. Here we demonstrate seven polymorphisms/mutations that were reported before. All were detected within very low frequency including the delta 32 mutation. However, we report for the first time copy number variants at two CCR5 gene locations; 45072265 and 38591712.