Clonal diversity and genetic profiling of antibiotic resistance among multidrug/carbapenem-resistant Klebsiella pneumoniae isolates from a tertiary care hospital in Saudi Arabia.

BACKGROUND: The nexus between resistance determinants, plasmid type, and clonality appears to play a crucial role in the dissemination and survival of carbapenem-resistant Klebsiella pneumoniae (CRKP). The incidence of infections involving CRKP in Saudi Arabia is increasing and there is a need for detailed molecular profiling of this pathogen for CRKP surveillance and control. METHODS: The resistance determinants of 71 non-redundant CRKP isolates were investigated by polymerase chain reaction (PCR) and sequencing. Plasmid typing was performed using PCR-based replicon typing and the clonality of isolates was determined by multilocus sequence typing. Capsular polysaccharide synthesis genes and other virulence factors were examined using multiplex PCR. Diversity was calculated using DIVEIN, clonal relationship was determined using eBURST, and phylogenetic analysis was performed using SplitsTree4. RESULTS: A polyclonal OXA-48 gene alone was the most common carbapenemase detected in 48/71 (67.6%) isolates followed by NDM-1 alone in 9/71 (12.7%) isolates. Coproduction of OXA-48 and NDM-1 was observed in 6/71 (8.5%) isolates. Both carbapenemase genes could be transferred into an Escherichia coli recipient. CTX-M-15 was the most abundant extended-spectrum ß-lactamase gene detected in 47/71 (66.2%) isolates, whereas clone-specific CTX-M-14 (ST-199 and -709) was found in 15/71 (21%) isolates. Sixty-seven of 71 isolates were positive for one or more plasmid replicons. The replicons detected were: IncFII; IncFIIK; IncFIA; IncFIB; L/M; IncI1; and IncN. FIIK and L/M were predominant, with 69 and 67% positivity, respectively. All isolates were negative for the magA (K1), rmpA, and K2 genes and presented a non-hypermucoviscous phenotype. CONCLUSION: A polyclonal CRKP reservoir of sequence types (STs)-37, - 199, and - 152 was observed and ST-152 appeared to be a "frequent carrier" of the NDM-1 gene. ST-199, a singleton not previously reported, showed a sequence diversity suggestive of positive selection. A significant association was evident between resistance determinants and the clonal types of K. pneumoniae: all ST-152 isolates were positive for NDM-1 but negative for OXA-48; ST-199 isolates were positive for OXA-48 but negative for NDM-1; and ST-709 and -199 isolates were positive for CTX-M-14. The incidence of certain clonal types in large numbers predicts an outbreak-like situation and warrants stringent surveillance and infection control.

What could be the role of genetic tests and machine learning of AXIN2 variant dominance in non-syndromic hypodontia? A case-control study in orthodontically treated patients.

BACKGROUND: Hypodontia is the most prevalent dental anomaly in humans, and is primarily attributed to genetic factors. Although genome-wide association studies (GWAS) have identified single-nucleotide polymorphisms (SNP) associated with hypodontia, genetic risk assessment remains challenging due to population-specific SNP variants. Therefore, we aimed to conducted a genetic analysis and developed a machine-learning-based predictive model to examine the association between previously reported SNPs and hypodontia in the Saudi Arabian population. Our case-control study included 106 participants (aged 8-50 years; 64 females and 42 males), comprising 54 hypodontia cases and 52 controls. We utilized TaqManTM Real-Time Polymerase Chain Reaction and allelic genotyping to analyze three selected SNPs (AXIN2: rs2240308, PAX9: rs61754301, and MSX1: rs12532) in unstimulated whole saliva samples. The chi-square test, multinomial logistic regression, and machine-learning techniques were used to assess genetic risk by using odds ratios (ORs) for multiple target variables. RESULTS: Multivariate logistic regression indicated a significant association between homozygous AXIN2 rs2240308 and the hypodontia phenotype (ORs [95% confidence interval] 2.893 [1.28-6.53]). Machine-learning algorithms revealed that the AXIN2 homozygous (A/A) genotype is a genetic risk factor for hypodontia of teeth #12, #22, and #35, whereas the AXIN2 homozygous (G/G) genotype increases the risk for hypodontia of teeth #22, #35, and #45. The PAX9 homozygous (C/C) genotype is associated with an increased risk for hypodontia of teeth #22 and #35. CONCLUSIONS: Our study confirms a link between AXIN2 and hypodontia in Saudi orthodontic patients and suggests that combining machine-learning models with SNP analysis of saliva samples can effectively identify individuals with non-syndromic hypodontia.

Mutated neuron navigator 3 as a candidate gene for a rare neurodevelopmental disorder.

BACKGROUND: Neuron navigator 3 (NAV3) is characterized as one of the neuron navigator family (NAV1, NAV2, NAV3) proteins predominantly expressed in the nervous system. The NAV3-encoded protein comprises a conserved AAA and coiled-coil domains characteristic of ATPases, which are associated with different cellular activities. METHODS: We describe a Saudi proband presenting a complex recessive neurodevelopmental disorder (NDD). Whole exome sequencing (WES) followed by Sanger sequencing, 3D protein modeling and RT-qPCR was performed. RESULTS: WES revealed a bi-allelic frameshift variant (c.2604_2605delAG; p.Val870SerfsTer12) in exon 12 of the NAV3 gene. Furthermore, RT-qPCR revealed a significant decrease in the NAV3 mRNA expression in the patient sample, and 3D protein modeling revealed disruption of the overall secondary structure. CONCLUSION: For the time, we associate a bi-allelic variant in the NAV3 gene causing NDD in humans.

Biallelic HMGXB4 loss-of-function variant causes intellectual disability, developmental delay, and dysmorphic features.

Background: HMGXB4 (additionally known as HMG2L1) is a non-histone DNA-binding protein that contains a single HMG-box domain. HMGXB4 was originally described in Xenopus where it was seen to negatively regulate the Wnt/ß-catenin signaling pathway. Materials and methods: In this study, we conducted a genetic and clinical evaluation of a single family with three affected individuals suffering from intellectual disability (ID), global developmental delay (GDD) and dysmorphic facial features.Whole genome sequencing (WGS) and Sanger sequencing were performed on the affected individuals' DNA to identify genetic variations. Additionally, a reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to assess gene expression in both the affected and unaffected individuals in the family. Result: WGS identified a homozygous frameshift variant c.1193_1196del p. (Lys398Argfs × 25) in exon 5 of the HMGXB4 gene (OMIM 604702), which completely segregated the disease phenotype in the family. Furthermore, RT-qPCR revealed a substantial decrease in the HMGXB4 gene expression in the affected individuals as compared to the unaffected individuals of the family. Conclusions: The current study is the first evidence linking a genetic variant in the HMGXB4 gene to ID, GDD, and dysmorphic facial features. Therefore, it is possible that HMGXB4 contributes significantly to developmental milestones and may be responsible for neurodevelopmental disorders in humans.

Congenital disorder of glycosylation with defective fucosylation 2 (FCSK gene defect): The third report in the literature with a mild phenotype.

BACKGROUND: Congenital disorders of glycosylation (CDG) are a group of heterogeneous disorders caused by abnormal lipid or protein glycosylation. Variants in the FCSK gene have been reported to cause CDG. Defective FCSK-induced CDG (FCSK-CDG) has only been reported previously in three unrelated children. METHODS: In this study, we genetically and clinically examined a 3-year-old proband with resolved infantile spasms and normal development. Standard whole-exome sequencing (WES) and Sanger sequencing were performed to identify the functional impact of the variant. RESULTS: WES revealed a rare biallelic missense variant (c.3013G>C; p.Val1005Leu) in FCSK. RT-qPCR showed a significant depletion in FCSK gene expression in the affected individual. Western blotting revealed reduced FCSK expression at the protein level compared to that in the control. Furthermore, 3D protein modeling suggested changes in the secondary structure, which might affect the overall FCSK protein function. CONCLUSION: This study broadens the mutation and phenotypic spectrum of FCSK-associated developmental disorders.

Novel homozygous pathogenic mitochondrial DNAJC19 variant in a patient with dilated cardiomyopathy and global developmental delay.

BACKGROUND: Dilated cardiomyopathy with ataxia syndrome (DCMA) or 3-methylglutaconic aciduria type V is a rare global autosomal recessive mitochondrial syndrome that is clinically and genetically heterogeneous. It is characterized by early-onset dilated cardiomyopathy and increased urinary excretion of 3-methylglutaconic acid. As a result, some patients die due to cardiac failure, while others manifest with growth retardation, microcytic anemia, mild ataxia, and mild muscle weakness. DCMA is caused by variants in the DnaJ heat shock protein family (Hsp40) member C19 gene (DNAJC19), which plays an important role in mitochondrial protein import machinery in the inner mitochondrial membrane. METHODS: We describe a single affected family member who presented with cardiomyopathy, global developmental delay, chest infection, seizures, elevated excretion of 3-methylglutaconic acid, and 3-methylglutaric acid in the urine. RESULTS: Whole-exome sequencing followed by Sanger sequencing revealed a homozygous frameshift variant in the reading frame starting at codon 54 in exon 4 in the DNAJC19 gene (c.159del [Phe54Leufs*5]), which results in a stop codon four positions downstream. Quantitative gene expression analysis revealed that DNAJC19 mRNA expression in this patient was substantially reduced compared to the control. CONCLUSIONS: We present a novel variant in the DNAJC19 gene that causes rare autosomal recessive mitochondrial 3-methylglutaconic aciduria type V. By comparing the current case with previously reported ones, we conclude that the disease is extremely heterogeneous for reasons that are still unknown.

Mutated VWA8 Is Associated With Developmental Delay, Microcephaly, and Scoliosis and Plays a Novel Role in Early Development and Skeletal Morphogenesis in Zebrafish.

Von Willebrand A domain-containing protein 8 (VWA8), also named KIAA0564, is a poorly characterized, mitochondrial matrix-targeted protein having a putative ATPase activity. VWA8 is comprising of ATPase-associated domains and a VWFA domain associated with ATPase activity inside the cell. In the present study, we describe a large consanguineous family of Saudi origin segregating a complex developmental syndrome in an autosomal recessive fashion. All the affected individuals exhibited severe developmental disorders. DNA from three patients was subjected to whole-exome sequencing followed by Sanger sequencing. VWA8 knock-down zebrafish morpholinos were used to study the phenotypic effect of this gene on zebrafish development. A homozygous missense variant [c.947A > G; p.(Asp316Gly)] was identified in exon 8 of the VWA8 gene, which perfectly segregated with the disease phenotype. Using zebrafish morpholino, we observed delayed development at an early stage, lack of movement, light sensitivity, severe skeletal deformity such as scoliosis, and facial dysmorphism. This is the first homozygous variant identified in the VWA8 gene underlying global developmental delay, microcephaly, scoliosis, limbs, and cardiovascular malformations in humans. We provide genetic and molecular evidence using zebrafish morpholino for a homozygous variant in the VWA8 gene, associated with such a complex developmental syndrome in humans.

Pancytopenia, Recurrent Infection, Poor Wound Healing, Heterotopia of the Brain Probably Associated with A Candidate Novel de Novo CDC42 Gene Defect: Expanding the Molecular and Phenotypic Spectrum.

CDC42 (cell division cycle protein 42) belongs to the Rho GTPase family that is known to control the signaling axis that regulates several cellular functions, including cell cycle progression, migration, and proliferation. However, the functional characterization of the CDC42 gene in mammalian physiology remains largely unclear. Here, we report the genetic and functional characterization of a non-consanguineous Saudi family with a single affected individual. Clinical examinations revealed poor wound healing, heterotopia of the brain, pancytopenia, and recurrent infections. Whole exome sequencing revealed a de novo missense variant (c.101C > A, p.Pro34Gln) in the CDC42 gene. The functional assays revealed a substantial reduction in the growth and motility of the patient cells as compared to the normal cells control. Homology three-dimensional (3-D) modeling of CDC42 revealed that the Pro34 is important for the proper protein secondary structure. In conclusion, we report a candidate disease-causing variant, which requires further confirmation for the etiology of CDC42 pathogenesis. This represents the first case from the Saudi population. The current study adds to the spectrum of mutations in the CDC42 gene that might help in genetic counseling and contributes to the CDC42-related genetic and functional characterization. However, further studies into the molecular mechanisms that are involved are needed in order to determine the role of the CDC42 gene associated with aberrant cell migration and immune response.

The effect of the VKORC1 promoter variant on warfarin responsiveness in the Saudi WArfarin Pharmacogenetic (SWAP) cohort.

Warfarin is a frequently prescribed oral anticoagulant with a narrow therapeutic index, requiring careful dosing and monitoring. However, patients respond with significant inter-individual variability in terms of the dose and responsiveness of warfarin, attributed to genetic polymorphisms within the genes responsible for the pharmacokinetics and pharmacodynamics of warfarin. Extensive warfarin pharmacogenetic studies have been conducted, including studies resulting in genotype-guided dosing guidelines, but few large scale studies have been conducted with the Saudi population. In this study, we report the study design and baseline characteristics of the Saudi WArfarin Pharmacogenomics (SWAP) cohort, as well as the association of the VKORC1 promoter variants with the warfarin dose and the time to a stable INR. In the 936 Saudi patients recruited in the SWAP study, the minor allele C of rs9923231 was significantly associated with a 8.45 mg higher weekly warfarin dose (p value = 4.0 × 10-46), as well as with a significant delay in achieving a stable INR level. The addition of the rs9923231 status to the model, containing all the significant clinical variables, doubled the warfarin dose explained variance to 31%. The SWAP cohort represents a valuable resource for future research with the objective of identifying rare and prevalent genetic variants, which can be incorporated in personalized anticoagulation therapy for the Saudi population.

Evolving sequence mutations in the Middle East Respiratory Syndrome Coronavirus (MERS-CoV).

BACKGROUND: Middle East respiratory syndrome coronavirus (MERS-CoV) has continued to cause sporadic outbreaks of severe respiratory tract infection over the last 8 years. METHODS: Complete genome sequencing using next-generation sequencing was performed for MERS-CoV isolates from cases that occurred in Riyadh between 2015 and 2019. Phylogenetic analysis and molecular mutational analysis were carried out to investigate disease severity. RESULTS: A total of eight MERS-CoV isolates were subjected to complete genome sequencing. Phylogenetic analysis resulted in the assembly of 7/8 sequences within lineage 3 and one sequence within lineage 4 showing complex genomic recombination. The isolates contained a variety of unique amino acid substitutions in ORF1ab (41), the N protein (10), the S protein (9) and ORF4b (5). CONCLUSION: Our study shows that MERS-CoV is evolving. The emergence of new variants carries the potential for increased virulence and could impose a challenge to the global health system. We recommend the sequencing every new MERS-CoV isolate to observe the changes in the virus and relate them to clinical outcomes.

Elucidating the Reprograming of Colorectal Cancer Metabolism Using Genome-Scale Metabolic Modeling.

Colorectal cancer is the third most incidental cancer worldwide, and the response rate of current treatment for colorectal cancer is very low. Genome-scale metabolic models (GEMs) are systems biology platforms, and they had been used to assist researchers in understanding the metabolic alterations in different types of cancer. Here, we reconstructed a generic colorectal cancer GEM by merging 374 personalized GEMs from the Human Pathology Atlas and used it as a platform for systematic investigation of the difference between tumor and normal samples. The reconstructed model revealed the metabolic reprogramming in glutathione as well as the arginine and proline metabolism in response to tumor occurrence. In addition, six genes including ODC1, SMS, SRM, RRM2, SMOX, and SAT1 associated with arginine and proline metabolism were found to be key players in this metabolic alteration. We also investigated these genes in independent colorectal cancer patients and cell lines and found that many of these genes showed elevated level in colorectal cancer and exhibited adverse effect in patients. Therefore, these genes could be promising therapeutic targets for treatment of a specific colon cancer patient group.

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.

Multi-drug carbapenem-resistant Klebsiella pneumoniae infection carrying the OXA-48 gene and showing variations in outer membrane protein 36 causing an outbreak in a tertiary care hospital in Riyadh, Saudi Arabia.

OBJECTIVES: To investigate the genes of antibiotic resistance among isolates from the first reported carbapenem-resistant Klebsiella pneumoniae (CRKP) outbreak in a tertiary care hospital, Riyadh, Saudi Arabia. METHODS: Antimicrobial susceptibility testing was performed on bacterial isolates using the Microscan Walkaway system (Siemens, Germany) and was confirmed by Etest (AB Biodisk, Sweden). bla-CTX-M, -SHV, -TEM, -OXA-48, OXA-A,B,C,D, -KPC, -NDM, -VIM, -IMP, integron 1, and outer membrane proteins(Omp)-35 and Omp-36 were investigated by PCR amplification and direct sequencing of PCR products. Isolates were sequence-typed by multilocus sequence typing (MLST). RESULTS: All isolates were resistant to cefotaxime, ceftazidime, cefepime, ciprofloxacin, and piperacillin-tazobactam, and 91% (21 out of 23) were resistant to amikacin and gentamicin. All isolates except two from a single patient were resistant to one of the carbapenems. CTX-M and SHV genes were detected in all isolates, CTX-M-15 and SHV-1 types being predominant among these extended-spectrum beta-lactamases (ESBLs). TEM-1 was found in all except one isolate (isolate 3). Significantly, the OXA-48 gene was also found in all isolates. OXA-D-gene was found in three out of 23 isolates. KPC, NDM, OXA-A, -B, -C, VIM, and IMP genes were absent in all isolates. Disruption of the Omp-36 gene due to insertion of transposon IS903 and/or IS4 was detected in four out of 23 isolates, and some unique variations were also observed in this gene, including an insertion of two amino acids in the L3 region of Omp-36 in one isolate (isolate 3) and a mutation resulting in a premature stop codon in another isolate (isolate 25). MLST revealed ST29 to be the predominant sequence type (17 out of 23 isolates, 74%). Three were ST709 and one each was ST37 and ST111; one isolate had an unknown ST. CONCLUSIONS: This is probably the first reported outbreak of multidrug/carbapenem-resistant Klebsiella infection involving the OXA-48 gene from Saudi Arabia. Although the presence of ESBLs such as OXA, CTX-M, TEM, and SHV are predictable reasons for resistance, variations in the Omp-36 gene might also have precipitated this phenomenon. Disruption of the Omp-36 sequence by large insertional elements, the insertion of two amino acids in a very crucial part of this protein, and the presence of a premature stop codon in one isolate might have rendered this protein incomplete and non-functional. The study also demonstrated that more than one type of clone was responsible for this reported apparent outbreak and that ST29, a clone not reported from this region before, was the major clone responsible.