Impact of Genetic Variations on Thromboembolic Risk in Saudis with Sickle Cell Disease.

BACKGROUND: Sickle cell disease (SCD) is a Mendelian disease characterized by multigenic phenotypes. Previous reports indicated a higher rate of thromboembolic events (TEEs) in SCD patients. A number of candidate polymorphisms in certain genes (e.g., FVL, PRT, and MTHFR) were previously reported as risk factors for TEEs in different clinical conditions. This study aimed to genotype these genes and other loci predicted to underlie TEEs in SCD patients. METHODOLOGY: A multi-center genome-wide association study (GWAS) involving Saudi SCD adult patients with a history of TEEs (n = 65) and control patients without TEE history (n = 285) was performed. Genotyping used the 10× Affymetrix Axiom array, which includes 683,030 markers. Fisher's exact test was used to generate p-values of TEE associations with each single-nucleotide polymorphism (SNP). The haplotype analysis software tool version 1.05, designed by the University of Göttingen, Germany, was used to identify the common inherited haplotypes. RESULTS: No association was identified between the targeted single-nucleotide polymorphism rs1801133 in MTHFR and TEEs in SCD (p = 0.79). The allele frequency of rs6025 in FVL and rs1799963 in PRT in our cohort was extremely low (<0.01); thus, both variants were excluded from the analysis as no meaningful comparison was possible. In contrast, the GWAS analysis showed novel genome-wide associations (p < 5 × 10-8) with seven signals; five of them were located on Chr 11 (rs35390334, rs331532, rs317777, rs147062602, and rs372091), one SNP on Chr 20 (rs139341092), and another on Chr 9 (rs76076035). The other 34 SNPs located on known genes were also detected at a signal threshold of p < 5 × 10-6. Seven of the identified variants are located in olfactory receptor family 51 genes (OR51B5, OR51V1, OR51A1P, and OR51E2), and five variants were related to family 52 genes (OR52A5, OR52K1, OR52K2, and OR52T1P). The previously reported association between rs5006884-A in OR51B5 and fetal hemoglobin (HbF) levels was confirmed in our study, which showed significantly lower levels of HbF (p = 0.002) and less allele frequency (p = 0.003) in the TEE cases than in the controls. The assessment of the haplotype inheritance pattern involved the top ten significant markers with no LD (rs353988334, rs317777, rs14788626882, rs49188823, rs139349992, rs76076035, rs73395847, rs1368823, rs8888834548, and rs1455957). A haplotype analysis revealed significant associations between two haplotypes (a risk, TT-AA-del-AA-ins-CT-TT-CC-CC-AA, and a reverse protective, CC-GG-ins-GG-del-TT-CC-TT-GG-GG) and TEEs in SCD (p = 0.024, OR = 6.16, CI = 1.34-28.24, and p = 0.019, OR = 0.33, CI = 0.13-0.85, respectively). CONCLUSIONS: Seven markers showed novel genome-wide associations; two of them were exonic variants (rs317777 in OLFM5P and rs147062602 in OR51B5), and less significant associations (p < 5 × 10-6) were identified for 34 other variants in known genes with TEEs in SCD. Moreover, two 10-SNP common haplotypes were determined with contradictory effects. Further replication of these findings is needed.

Impact of telemedicine on glycemic control in type 2 diabetes mellitus during the COVID-19 lockdown period.

Background: The lockdown at the start of coronavirus disease 2019 (COVID-19) pandemic in Saudi Arabia (March 2020 to June 2020) shifted routine in-person care for patients with type 2 diabetes mellitus (T2DM) to telemedicine. The aim of this study was to investigate the impact telemedicine had during this period on glycemic control (HbA1c) in patients with T2DM. Methods: 4,266 patients with T2DM were screened from five Ministry of National Guard Health Affairs hospitals in the Kingdom of Saudi Arabia. Age, gender, body mass index (BMI), HbA1c (before and after the COVID-19 lockdown), duration of T2DM, comorbidities and antidiabetic medications data were obtained. Mean and standard deviation of differences in HbA1c were calculated to assess the impact of telemedicine intervention. Correlations between clinically significant variances (when change in the level is ≥0.5%) in HbA1c with demographics and clinical characteristic data were determined using chi square test. Results: Most of the participants were Saudis (97.7%) with 59.7% female and 56.4% ≥60 years of age. Obesity was 63.8%, dyslipidemia 91%, and hypertension 70%. Mean HbA1c of all patients slightly rose from 8.52% ± 1.5% before lockdown to 8.68% ± 1.6% after lockdown. There were n=1,064 patients (24.9%) whose HbA1c decreased by ≥0.5%, n =1,574 patients whose HbA1c increased by ≥0.5% (36.9%), and n =1,628 patients whose HbA1c changed by <0.5% in either direction (38.2%). More males had significant improvements in glycemia compared to females (28.1% vs 22.8%, p<0.0001), as were individuals below the age of 60 years (28.1% vs 22.5%, p<0.0001). Hypertensive individuals were less likely than non-hypertensive to have glycemic improvement (23.7% vs 27.9%, p=0.015). More patients on sulfonylureas had improvements in HbA1c (42.3% vs 37.9%, p=0.032), whereas patients on insulin had higher HbA1c (62.7% vs 56.2%, p=0.001). HbA1c changes were independent of BMI, duration of disease, hyperlipidemia, heart and kidney diseases. Conclusion: Telemedicine was helpful in delivering care to T2DM patients during COVID-19 lockdown, with 63.1% of patients maintaining HbA1c and improving glycemia. More males than females showed improvements. However, the HbA1c levels in this cohort of patients pre- and post-lockdown were unsatisfactorily high, and may be due to in part lifestyle, age, education, and hypertension.

Pharmacogenes that demonstrate high association evidence according to CPIC, DPWG, and PharmGKB.

Background: Different levels of evidence related to the variable responses of individuals to drug treatment have been reported in various pharmacogenomic (PGx) databases. Identification of gene-drug pairs with strong association evidence can be helpful in prioritizing the implementation of PGx guidelines and focusing on a gene panel. This study aimed to determine the pharmacogenes with the highest evidence-based association and to indicate their involvement in drug-gene interactions. Methodology: The publicly available datasets CPIC, DPWG, and PharmGKB were selected to determine the pharmacogenes with the highest drug outcome associations. The upper two levels of evidence rated by the three scoring methods were specified (levels A-B in CPIC, 3-4 in DPWG, or 1-2 levels in PharmGKB). The identified pharmacogenes were further ranked in this study based on the number of medications they interacted with. Results: Fifty pharmacogenes, with high to moderately high evidence of associations with drug response alterations, with potential influence on the therapeutic and/or toxicity outcomes of 152 drugs were identified. CYP2D6, CYP2C9, CYP2C19, G6PD, HLA-B, SLCO1B1, CACNA1S, RYR1, MT-RNR1, and IFNL4 are the top 10 pharmacogenes, where each is predicted to impact patients' responses to ≥5 drugs. Conclusion: This study identified the most important pharmacogenes based on the highest-ranked association evidence and their frequency of involvement in affecting multiple drugs. The obtained data is useful for customizing a gene panel for PGx testing. Identifying the strength of scientific evidence supporting drug-gene interactions aids drug prescribers in making the best clinical decision.

Prevalence of exposure to pharmacogenetic drugs by the Saudis treated at the health care centers of the Ministry of National Guard.

Background: The drugs impacted by genetic variants are known as pharmacogenetic (PGx) drugs. Patients' responses to these drugs may vary according to the variability in patients' genetic makeup. Hence, exploring the pharmacogenes that affect drug treatment is vital to ensure optimal therapy and patients' safety. This study aimed to describe the usage rate of PGx drugs and the frequency of relevant variants in the Saudi population. Methodology: Prescription patterns over seven years (2015-2021) for Saudi patients on PGx drugs treated at the Ministry of National Guard-Health Affairs (MNG-HA) were investigated. Only registered drugs in the MNG-HA formulary (n = 78) were included. The patients were subgrouped into four age groups: ≤24, 25-44, 45-64, and ≥65 years. Further subgrouping was made according to gender and drugs' therapeutic categories following anatomical therapeutic chemical (ATC) classification.Furthermore, an online searching was carried out to identify the pharmacogenes reported in the literature among healthy Saudis. The search included 45 genes that may affect drug outcomes based on evidence rated by either CPIC (A-B levels) or PharmGKB (1-2 levels). Results: The screened patients were 1,483,905. Patients on PGx drugs accounted for 46.7% (n = 693,077 patients). The analgesic group was the most prescribed drug category (47%), which included ibuprofen (20.5%), celecoxib (6.3%), tramadol (5.8%), and others. Cardiovascular agents were the second-most utilized drug class (24.4%). Omeprazole was the second most commonly used medication (11.1%) but ranked third as a class (gastroenterology). Females used PGx drugs more frequently than males (53.5% versus 46.5%) and a higher usage rate by patients aged 45-64 years (31.3%) was noted. The cytochrome P450 genes (CYP2C9, CYP2C19, and CYP2D6) were estimated to impact responses of 54.3% (n = 1,156,113) of the used drugs (27.2% are possibly affected by CYP2C9, 12.8% by CYP2C19, and 14.3% by CYP2D6). Thirty-five pharmacogenes that characterize Saudi population and their variants' allele frequencies were identified from previous reports. This study presents the largest reported number of genes that may affect drug therapies among Saudis. Conclusion: This study confirmed that a high percentage of Saudi patients use PGx drugs and various genotypes of certain pharmacogenes are inherited by the Saudi population.

Impact of SLCO1B1*5 on Flucloxacillin and Co-Amoxiclav-Related Liver Injury.

Background: Idiosyncratic drug-induced liver injury (DILI) is a serious uncommon disease that may develop as a result of the intake of certain drugs such as the antimicrobials flucloxacillin and co-amoxiclav. The reported cases showed significant associations between DILI and various human leukocyte (HLA) markers. The solute carrier organic anion transporter 1B1 (SLCO1B1), a non-HLA candidate gene, was previously reported as a risk factor for liver injury induced by rifampin and methimazole. This study presumed that SLCO1B1 may play a general role in the DILI susceptibility and therefore investigated the association of rs4149056 (SLCO1B1*5, T521C) polymorphism with flucloxacillin- and co-amoxiclav-induced liver injury. Methodology: We recruited 155 and 165 DILI cases of white ancestral origin from various European countries but mainly from the United Kingdom owing to flucloxacillin and co-amoxiclav, respectively. Only adult patients (≥18 years) who were diagnosed with liver injury and who showed i) clinical jaundice or bilirubin >2x the upper limit of normal (ULN), ii) alanine aminotransferase (ALT) >5x ULN or iii) alkaline phosphatase (ALP) >2x ULN and bilirubin > ULN were selected. The population reference sample (POPRES), a European control group (n = 282), was used in comparison with the investigated cases. TaqMan SNP genotyping custom assay designed by Applied Biosystems was used to genotype both DILI cohorts for SLCO1B1 polymorphism (rs4149056). Allelic discrimination analysis was performed using a step one real-time PCR machine. Genotype differences between cases and controls were examined using Fisher's exact test. GraphPad Prism version 5.0 was used to determine the p-value, odds ratio, and 95% confidence interval. Compliance of the control group with Hardy-Weinberg equilibrium was proven using a web-based calculator available at https://wpcalc.com/en/equilibrium-hardy-weinberg/. Results: A small number of cases failed genotyping in each cohort. Thus, only 149 flucloxacillin and 162 co-amoxiclav DILI cases were analyzed. Genotyping of both DILI cohorts did not show evidence of association with the variant rs4149056 (T521C) (OR = 0.71, 95% CI = 0.46-1.12; p = 0.17 for flucloxacillin cases and OR = 0.87, 95% CI = 0.56-1.33; p = 0.58 for co-amoxiclav), although slightly lower frequency (22.8%) of positive flucloxacillin cases was noticed than that of POPRES controls (29.4%). Conclusion: Carriage of the examined allele SLCO1B1*5 is not considered a risk factor for flucloxacillin DILI or co-amoxiclav DILI as presumed. Testing a different allele (SLCO1B1*1B) and another family member gene (SLCO1B3) may still be needed to provide a clearer role of SLCO1B drug transporters in DILI development-related to the chosen antimicrobials.

HIV-Tat Exacerbates the Actions of Atazanavir, Efavirenz, and Ritonavir on Cardiac Ryanodine Receptor (RyR2).

The incidence of sudden cardiac death (SCD) in people living with HIV infection (PLWH), especially those with inadequate viral suppression, is high and the reasons for this remain incompletely characterized. The timely opening and closing of type 2 ryanodine receptor (RyR2) is critical for ensuring rhythmic cardiac contraction-relaxation cycles, and the disruption of these processes can elicit Ca2+ waves, ventricular arrhythmias, and SCD. Herein, we show that the HIV protein Tat (HIV-Tat: 0-52 ng/mL) and therapeutic levels of the antiretroviral drugs atazanavir (ATV: 0-25,344 ng/mL), efavirenz (EFV: 0-11,376 ng/mL), and ritonavir (RTV: 0-25,956 ng/mL) bind to and modulate the opening and closing of RyR2. Abacavir (0-14,315 ng/mL), bictegravir (0-22,469 ng/mL), Rilpivirine (0-14,360 ng/mL), and tenofovir disoproxil fumarate (0-18,321 ng/mL) did not alter [3H]ryanodine binding to RyR2. Pretreating RyR2 with low HIV-Tat (14 ng/mL) potentiated the abilities of ATV and RTV to bind to open RyR2 and enhanced their ability to bind to EFV to close RyR2. In silico molecular docking using a Schrodinger Prime protein-protein docking algorithm identified three thermodynamically favored interacting sites for HIV-Tat on RyR2. The most favored site resides between amino acids (AA) 1702-1963; the second favored site resides between AA 467-1465, and the third site resides between AA 201-1816. Collectively, these new data show that HIV-Tat, ATV, EFV, and RTV can bind to and modulate the activity of RyR2 and that HIV-Tat can exacerbate the actions of ATV, EFV, and RTV on RyR2. Whether the modulation of RyR2 by these agents increases the risk of arrhythmias and SCD remains to be explored.

Two Cases of Recessive Intellectual Disability Caused by NDST1 and METTL23 Variants.

Intellectual disability (ID) is a highly heterogeneous genetic condition with more than a thousand genes described so far. By exome sequencing of two consanguineous families presenting hallmark features of ID, we identified two homozygous variants in two genes previously associated with autosomal recessive ID: NDST1 (c.1966G>A; p.Asp656Asn) and METTL23 (c.310T>C; p.Phe104Leu). The segregation of the variants was validated by Sanger sequencing in all family members. In silico homology modeling of wild-type and mutated proteins revealed substantial changes in the secondary structure of both proteins, indicating a possible effect on function. The identification and validation of new pathogenic NDST1 and METTL23 variants in two cases of autosomal recessive ID further highlight the importance of these genes in proper brain function and development.

Investigation of Oxidative Stress-Related Candidate Genes as Risk Factors for Drug-Induced Liver Injury due to Co-Amoxiclav.

The liver is susceptible to drug toxicity due to its vital role in xenobiotic metabolism and elimination. In addition to human leukocyte antigen (HLA) variants, which were previously determined as risk factors for drug-induced liver injury (DILI) due to co-amoxiclav, other non-HLA genes may contribute to hepatotoxicity risk. In this study, the association between DILI due to co-amoxiclav and several non-HLA genes was investigated. Association of variants in candidate genes (SOD2, GPX1, GSTM1, and GSTT1) with DILI due to various drugs was reported previously in other DILI cohorts. This study examined relevance in a co-amoxiclav-DILI cohort. One hundred sixty-five co-amoxiclav DILI cases were recruited from several European countries by two different studies (DILIGEN and iDILIC). A North-East England population group (n = 334) was used as the control group. PCR assays were used to genotype for the GSTM1 and GSTT1 null alleles with TaqMan SNP genotyping assays used for SOD2 (rs4880) and GPX1 (rs1050450). Fisher's exact test was used to assess differences in significance between cases and controls. None of the studied variants (SOD2 rs4880, GPX1 rs1050450, GSTM1 null allele, and GSTT1 null allele) was significantly associated with co-amoxiclav DILI compared with the control group. No significant differences between cases and controls were seen when combined SOD2/GPX1 genotypes and GST genotypes were considered. Despite the possible functional relevance and the previously reported contribution of the selected genes to DILI, our study failed to confirm associations between the selected genes and liver injury induced by co-amoxiclav.

Correction to: Homozygous missense variant in the TTN gene causing autosomal recessive limb-girdle muscular dystrophy type 10.

Please be advised that following publication of the original article [1], the authors have identified the following errors with the scientific content.

Homozygous missense variant in the TTN gene causing autosomal recessive limb-girdle muscular dystrophy type 10.

BACKGROUND: Limb-girdle muscular dystrophies (LGMDs) are large group of heterogeneous genetic diseases, having a hallmark feature of muscle weakness. Pathogenic mutations in the gene encoding the giant skeletal muscle protein titin (TTN) are associated with several muscle disorders, including cardiomyopathy, recessive congenital myopathies and limb-girdle muscular dystrophy (LGMD) type10. The phenotypic spectrum of titinopathies is expanding, as next generation sequencing (NGS) technology makes screening of this large gene possible. AIM: This study aimed to identify the pathogenic variant in a consanguineous Pakistani family with autosomal recessive LGMD type 10. METHODS: DNA from peripheral blood samples were obtained, whole exome sequencing (WES) was performed and several molecular and bioinformatics analysis were conducted to identify the pathogenic variant. TTN coding and near coding regions were further amplified using PCR and sequenced via Sanger sequencing. RESULTS: Whole exome sequencing analysis revealed a novel homozygous missense variant (c.98807G > A; p.Arg32936His) in the TTN gene in the index patients. No heterozygous individuals in the family presented LGMD features. The variant p.Arg32936His leads to a substitution of the arginine amino acid at position 32,936 into histidine possibly causing LGMD type 10. CONCLUSION: We identified a homozygous missense variant in TTN, which likely explains LGMD type 10 in this family in line with similar previously reported data. Our study concludes that WES is a successful molecular diagnostic tool to identify pathogenic variants in large genes such as TTN in highly inbred population.

Use of Pharmacogenetic Drugs by the Dutch Population.

INTRODUCTION: The Dutch Pharmacogenetics Working Group (DPWG) indicated a list of actionable genotypes that affect patients' response to more 50 drugs; these drugs which show variable effects based on patients' genetic traits were named as pharmacogenetics (PGX) drugs. Preemptive genetic testing before using these drugs may protect certain patients from serious adverse reactions and could help in avoiding treatment failures. The objectives of this study include identifying the rate of PGX drug usage among Dutch population, estimating the level of users who carry the actionable genotypes and determining the main genes involved in drug's effect variability. METHODS: Usage of PGX drugs over 2011-2017 by the insured population (an average of 11.4 million) in outpatient clinics in Netherlands was obtained from the publically available GIP databank. The data of 45 drugs were analyzed and their interactions with selected pharmacogenes were estimated. Frequency of actionable genotypes of 249 Dutch parents was obtained from the public database: Genome of Netherlands (GoNL), to identify the pattern of genetic characteristics of Dutch population. RESULTS: Over a 7 year period, 51.3 million exposures of patients to PGX drugs were reported with an average of 5.3 exposures per each drug user. One quarterof the exposures (12.4 million) are predicted to be experienced by individuals with actionable genotypes (risky exposures). Up to 60% of the risky exposures (around 7.5 million) were related to drugs metabolized by CYP2D6. SLCO1B1, and CYP2C19 were also identified among the top genes affecting response of drugs users (involved in about 22 and 12.4% of the risky exposures, respectively). Cardiovascular medications were the top prescribed PGX drug class (43%), followed by gastroenterology (29%) and psychiatry/neurology medications (15%). Women use more PGX drugs than men (55.8 vs. 44.2%, respectively) with the majority (84%) of users in both sexes are above 45 years. CONCLUSION: PGX drugs are commonly used in Netherlands. Preemptive panel testing for CYP2D6, SLCO1B1, and CYP2C19 only could be useful to predict 95% of vulnerable patients' exposures to PGX drugs. Future studies to assess the economic impact of preemptive panel testing on patients of older age are suggested.

Drug-Induced Liver Injury due to Flucloxacillin: Relevance of Multiple Human Leukocyte Antigen Alleles.

Some patients prescribed flucloxacillin (~ 0.01%) develop drug-induced liver injury (DILI). HLA-B*57:01 is an established genetic risk factor for flucloxacillin DILI. To consolidate this finding, identify additional genetic factors, and assess relevance of risk factors for flucloxacillin DILI in relation to DILI due to other penicillins, we performed a genomewide association study involving 197 flucloxacillin DILI cases and 6,835 controls. We imputed single-nucleotide polymorphism and human leukocyte antigen (HLA) genotypes. HLA-B*57:01 was the major risk factor (allelic odds ratio (OR) = 36.62; P = 2.67 × 10-97 ). HLA-B*57:03 also showed an association (OR = 79.21; P = 1.2 × 10-6 ). Within the HLA-B protein sequence, imputation showed valine97 , common to HLA-B*57:01 and HLA-B*57:03, had the largest effect (OR = 38.1; P = 9.7 × 10-97 ). We found no HLA-B*57 association with DILI due to other isoxazolyl penicillins (n = 6) or amoxicillin (n = 15) and no significant non-HLA signals for any penicillin-related DILI.

A Novel Homozygous Nonsense Mutation p.Cys366* in the WNT10B Gene Underlying Split-Hand/Split Foot Malformation in a Consanguineous Pakistani Family.

Split hand/split foot malformation (SHFM) or ectrodactyly is characterized by a deep median cleft of the hand or foot, hypoplasia or aplasia of the metacarpals, metatarsals, and phalanges. It is a clinically and genetically heterogeneous group of limb malformations. This study aimed to identify the pathogenic variant in a consanguineous Pakistani family with autosomal recessive SHFM. Peripheral blood samples were obtained, DNA was extracted, WNT10B coding and noncoding regions were PCR amplified and Sanger sequencing was performed using workflow suggested by Thermo Fisher Scientific. A novel homozygous nonsense variant (c.1098C>A; p.Cys366*) was identified in the WNT10B gene in the index patients, which probably explains SHFM type 6 in this family in comparison with similar data from the literature.

Druggability of Coronary Artery Disease Risk Loci.

BACKGROUND: Genome-wide association studies have identified multiple loci associated with coronary artery disease and myocardial infarction, but only a few of these loci are current targets for on-market medications. To identify drugs suitable for repurposing and their targets, we created 2 unique pipelines integrating public data on 49 coronary artery disease/myocardial infarction-genome-wide association studies loci, drug-gene interactions, side effects, and chemical interactions. METHODS: We first used publicly available genome-wide association studies results on all phenotypes to predict relevant side effects, identified drug-gene interactions, and prioritized candidates for repurposing among existing drugs. Second, we prioritized gene product targets by calculating a druggability score to estimate how accessible pockets of coronary artery disease/myocardial infarction-associated gene products are, then used again the genome-wide association studies results to predict side effects, excluded loci with widespread cross-tissue expression to avoid housekeeping and genes involved in vital processes and accordingly ranked the remaining gene products. RESULTS: These pipelines ultimately led to 3 suggestions for drug repurposing: pentolinium, adenosine triphosphate, and riociguat (to target CHRNB4, ACSS2, and GUCY1A3, respectively); and 3 proteins for drug development: LMOD1 (leiomodin 1), HIP1 (huntingtin-interacting protein 1), and PPP2R3A (protein phosphatase 2, regulatory subunit b-double prime, α). Most current therapies for coronary artery disease/myocardial infarction treatment were also rediscovered. CONCLUSIONS: Integration of genomic and pharmacological data may prove beneficial for drug repurposing and development, as evidence from our pipelines suggests.