Deleterious missense variants in the aflatoxin biosynthesis genes explain the low toxicity of Aspergillus flavus from infected rice.

Aspergillus flavus is one of the most natural contaminants of the improperly stored rice grains. It produces several secondary metabolites, like aflatoxins, which are well known hepatotoxic, hepatocarcinogenic and mutagenic agents. This study describes the in silico consequences of the missense mutations identified in several genes of aflatoxins biosynthesis in rice-contaminating A. flavus isolates. In the in vitro portion of the study, aflatoxins production profile was measured, and PCR-single strand-conformation polymorphism (SSCP)-sequencing method was used to genotype the studied genetic loci: aflP, aflM, aflR, PEP, and cob. Results showed aflatoxigenic potential in 79 out of 109 A. flavus isolates. Twenty-two missense and fifty-five synonymous mutations were found to be distributed variably on the studied loci. In the in silico portion of this study, several computations were utilized to predict the effect of each observed missense mutation on proteins structure, function, and stability. Seven mutations (O-methyl transferase: p.G256C; ver-1 dehydrogenase: p.K179 N and p.V183L; aspergillopepsin-1: p.P137L, p.S138F, p.G154C, and p.S158C) were found to be highly deleterious among the missense variants with damaging effects on their proteins' structure and function. In contrast to these detected variations in the aflatoxigenic loci, all missense mutations in the control non-aflatoxigenic cob gene were found to be neutral. These findings indicated that the observed mutations may reduce the concomitant biohazard of their biosynthesized aflatoxins. The current findings suggest that the naturally available variants may reduce or eliminates the dangerous consequences of aflatoxins upon ingesting the rice infected with A. flavus. To the best of our knowledge, this study is the first comprehensive report to analyze the missense mutations on the aflatoxin biosynthesis genes using in vitro and the state-of-art bio-computational tools.

Risk Y-haplotypes and pathogenic variants of Arab-ancestry boys with autism by an exome-wide association study.

Autism is heterogeneous multifactorial neurodevelopmental and neuropsychiatric disorder with repetitive and limited behaviors as well as communication deficits. Prevalence of autism in males is predominant than females, but their genetic association is unclear. The study was performed to investigate Y-chromosome haplotypes, significant risk variants and susceptibility genes associated with autistic Saudi young males with autism. Exome genotyping microarray analysis was performed in Saudi young boys with autism (cases, n = 47) and without autism and other genetic or neurodevelopmental disorders (control, n = 43) to identify the functional exonic risk variations among 243,345 exonic variations. The most significant single nucleotide polymorphisms (SNPs) of protein coding associated with autism in Saudi young boys were studied for functional enrichment. Y-chromosome haplotyping analysis of 6 SNPs such as rs1865680, rs2032624, rs2032658, rs2032631, rs9786153 and rs13447352 uncovered the most significant protective (ACGACA p = 2.94 × 10-9) among the controls and the high risk Y-haplotype (GAAGTC p = 6.85 × 10-6) among autistic boys. Exome association study revealed 6 susceptible genes, MCC, AUTS2, VSX1, SETBP1, CNTN3, and PCDH11Y that were known for autistic disorder. The significant predisposed genes with functional variants of Y-chromomere are strongly connected with spermatogenic failure (p = 8.02 × 10-8), azoospermia (p = 6.32 × 10-7), partial chromosome Y deletion (p = 7.66 × 10-6), HDMs demethylate histones pathway (p = 3.55 × 10-4) and immune system diseases (p = 4.11 × 10-3). Y-haplotypes and highly significant pathogenic exonic variants in MCC, AUTS2, VSX1, SETBP1, CNTN3 and PCDH11Y genes are more influential genetic factors for developing autism in boys of Arab origin.

Co-inheritance of alpha globin gene deletion lowering serum iron level in female beta thalassemia patients.

In the Eastern province of Saudi Arabia, thalassemia is highly common. Data on the effect of alpha globin gene variation on the concentration of iron on transfusion dependent Saudis are scanty. A total of 166 transfusions dependent ß-thalassemia were included in this study to understand association between the alpha globin gene variation and concentration of iron. Using multiplex PCR, the alpha globin gene deletions were identified. Also, HBA1 and HBA2 genes were sequenced by Sanger sequencing. Saudi transfusion dependent female ß-thalassemia patients with wild alpha globin genotype (αα/αα) were observed with iron level beyond the normal range. However, normal range of iron was observed in transfusion dependent Saudi female beta thalassemia patients co-inherited with double (-α3.7/-α3.7, or --Fil/αα or --MED/αα or - (α) 20.5/αα) or double heterozygosity (- -/-α3.7) alpha globin gene deletions, which is significantly (p < 0.0001) less compared to the Saudi transfused female with wild alpha globin genotype (αα/αα). The co-inheritance alpha globin gene deletions in female beta thalassemia patients were significantly lowering serum iron. Detailed studies can be taken forward to identify the molecular pathways involved in globin gene deletion as modulator.

Optimal pcr primers for rapid and accurate detection of Aspergillus flavus isolates.

Aspergillus flavus is among the most devastating opportunistic pathogens of several food crops including rice, due to its high production of carcinogenic aflatoxins. The presence of these organisms in economically important rice strip farming is a serious food safety concern. Several polymerase chain reaction (PCR) primers have been designed to detect this species; however, a comparative assessment of their accuracy has not been conducted. This study aims to identify the optimal diagnostic PCR primers for the identification of A. flavus, among widely available primers. We isolated 122 A. flavus native isolates from randomly collected rice strips (N = 300). We identified 109 isolates to the genus level using universal fungal PCR primer pairs. Nine pairs of primers were examined for their PCR diagnostic specificity on the 109 isolates. FLA PCR was found to be the optimal PCR primer pair for specific identification of the native isolates, over aflP(1), aflM, aflA, aflD, aflP(3), aflP(2), and aflR. The PEP primer pair was found to be the most unsuitable for A. flavus identification. In conclusion, the present study indicates the powerful specificity of the FLA PCR primer over other commonly available diagnostic primers for accurate, rapid, and large-scale identification of A. flavus native isolates. This study provides the first simple, practical comparative guide to PCR-based screening of A. flavus infection in rice strips.

Association of MBL2 Gene Polymorphism with Dental Caries in Saudi Children.

The high prevalence of dental caries in children worldwide is a major oral health problem which requires early intervention. Dental caries is mainly caused by the action of acids produced by bacteria in addition to many other factors. Recent genetic studies have reported that a number of genes are associated with the susceptibility to dental caries. The majority of these genes are associated with inflammation, increased susceptibility to infection, and dentine matrix formation. Using the TaqMan assay and direct DNA sequencing, the prevalence of 6 single-nucleotide polymorphisms (SNPs) in MMP9, MBL2, MMP2, and TIMP2 genes was determined in 102 children with caries and in 100 age-matched caries-free controls. Out of the 6 SNPs tested in the 4 selected genes, only rs11003125 in the MBL2 gene was shown to be associated with a high prevalence of caries in our cohort. In addition, haplotype analysis of the 6 SNPs tested revealed that certain haplotypes, namely GT of rs11003125G and rs7501477T and GT of rs7096206G and rs7501477T, were found to be associated with a high prevalence of dental caries in our cohort, while haplotype AG of rs17576A and rs7501477G was found to have a protective effect against dental caries. In conclusion, the data indicate that rs11003125 in the MBL2 gene was shown to be associated with a high prevalence of caries in our cohort, and 2 haplotypes are also involved in the increased susceptibility to dental caries.

Type 2 diabetes associated variants of KCNQ1 strongly confer the risk of cardiovascular disease among the Saudi Arabian population.

Genome-wide association studies have identified several loci associated with an increased risk for cardiovascular disease (CVD) and type 2 diabetes (T2D). Polymorphisms within the KCNQ1 (potassium voltage-gated channel, KQT-like subfamily, member 1) gene are consistently associated with T2D in a number of populations. The current study was undertaken to evaluate the association of 3 polymorphisms of KCNQ1 (rs2237892, rs151290 and rs2237895) with T2D and/or CVD. Patients diagnosed with either T2D (320 patients), CVD (250 patients) or both (60 patients) and 516 healthy controls were genotyped by TaqMan assay run on a real time PCR thermocycler. A statistically significant association was found for SNPs rs151290 (OR = 1.76; 95%CI = 1.02-3.05; p = 0.0435) and rs2237895 (OR = 2.49; 95%CI = 1.72-3.61; p < 0.0001) with CVD. SNP rs151290 (OR = 7.43; 95%CI = 1.00-55.22; p = 0.0499) showed a strong association in patients with both T2D and CVD. None of the SNPs showed any significant association with T2D. Haploview analysis showed that the ACC (rs151290, rs2237892 and rs2237895) haplotype is the most significant risk allele combination for CVD, while CCA is the most significant risk haplotype for co-morbidity with T2D. KCNQ1 polymorphism at SNPs rs151290 and rs2237895 is strongly associated with CVD in this population, but presented no association with T2D.

Intronic Polymorphisms in the CDKN2B-AS1 Gene Are Strongly Associated with the Risk of Myocardial Infarction and Coronary Artery Disease in the Saudi Population.

Recent genome-wide association studies identified single nucleotide polymorphisms (SNPs) on the chromosome 9p21.3 conferring the risk for CAD (coronary artery disease) in individuals of Caucasian ancestry. We performed a genetic association study to investigate the effect of 12 candidate SNPs within 9p21.3 locus on the risk of CAD in the Saudi population of the Eastern Province of Saudi Arabia. A total of 250 Saudi CAD patients who had experienced an myocardial infarction (MI) and 252 Saudi age-matched healthy controls were genotyped using TaqMan assay. Controls with evidenced lack of CAD provided 90% of statistical power at the type I error rate of 0.05. Five percent of the results were rechecked for quality control using Sanger sequencing, the results of which concurred with the TaqMan genotyping results. Association analysis of 12 SNPs indicated a significant difference in the genotype distribution for four SNPs between cases and controls (rs564398 p = 0.0315, χ² = 4.6, odds ratio (OD) = 1.5; rs4977574 p = 0.0336, χ² = 4.5, OD = 1.4; rs2891168 p = 1.85 × 10 - 10, χ² = 40.6, OD = 2.1 and rs1333042 p = 5.14 × 10 - 9, χ² = 34.1, OD = 2.2). The study identified three protective haplotypes (TAAG p = 1.00 × 10 - 4; AGTA p = 0.022 and GGGCC p = 0.0175) and a risk haplotype (TGGA p = 2.86 × 10 - 10) for the development of CAD. This study is in line with others that indicated that the SNPs located in the intronic region of the CDKN2B-AS1 gene are associated with CAD.

Co-inheritance of novel ATRX gene mutation and globin (α & ß) gene mutations in transfusion dependent beta-thalassemia patients.

α-Thalassemia X-linked mental retardation syndrome is a rare inherited intellectual disability disorder due to mutations in the ATRX gene. In our previous study of the prevalence of ß-thalassemia mutations in the Eastern Province of Saudi Arabia, we confirmed the widespread coinheritance of α-thalassemia mutation. Some of these subjects have a family history of mental retardation, the cause of which is unknown. Therefore, we investigated the presence or absence of mutations in the ATRX gene in these patients. Three exons of the ATRX gene and their flanking regions were directly sequenced. Only four female transfusion dependent ß-thalassemia patients were found to be carriers of a novel mutation in the ATRX gene. Two of the ATRX gene mutations, c.623delA and c.848T>C were present in patients homozygous for IVS I-5(G→C) and homozygous for Cd39(C → T) ß-thalassemia mutation, respectively. While the other two that were located in the intronic region (flanking regions), were present in patients homozygous for Cd39(C → T) ß-thalassemia mutation. The two subjects with the mutations in the coding region had family members with mental retardation, which suggests that the novel frame shift mutation and the missense mutation at coding region of ATRX gene are involved in ATRX syndrome.

Phenomics and genomic features of Enterococcus avium IRMC1622a isolated from a clinical sample of hospitalized patient.

BACKGROUND: Enterococcus avium (E. avium) is a Gram-positive nosocomial pathogen that is commonly isolated from the alimentary tract. The objective of this functional genomics study was to identify the resistant genes by analyzing the genome of E. avium IRMC1622a, a type of bacteria found in feces collected from a patient at a Saudi Arabian tertiary hospital. METHODS: The bacterial strain IRMC1622a was identified by 16 S rRNA sequencing as Enterococcus sp. The resistance phenomics were performed using VITEK® 2, and morphological analysis was achieved using a scanning electron microscope (SEM). Finally, the whole bacterial genome of the bacterial strain IRMC1622a was subjected to sequencing during October 2023 using Oxford Nanopore long-read sequencing technology, and mining for resistant genes. RESULTS: The results of antimicrobial resistant phenomics indicated that the IRMC1622a strain was sensitive to all tested antimicrobial agents except for erythromycin, and the same result was confirmed by genomic analysis in addition to other classes of antibiotics. SEM showed E. avium IRMC1622a is ovoid shape, in single cells (L 1.2797 ± 0.1490 µm), in pairs (L 1.7333 ± 0.1054 µm), and in chains (L 2.44033 ± 0.1978 µm). The E. avium IRMC1622a genome has 14 (in CARD) antimicrobial resistance genes that were identified with several mechanisms of antimicrobial resistance, such as the efflux pump and conferring antibiotic resistance. The present study revealed that the E. avium IRMC1622a genome contains a high number of genes associated with virulence factors, and 14 matched pathogenic protein families and predicted as human pathogen (probability score 0.855). We report two (ISEnfa4 and ISEfa5) mobile genetic elements for the first time in the E. avium genome. CONCLUSIONS: The study concludes that E. avium IRMC1622a is susceptible to all tested antibacterials except erythromycin. The IRMC1622a has 14 genes encoding antimicrobial resistance mechanisms, including the efflux pump and conferring antibiotic resistance. This could indicate a potential rise in E. avium resistance in healthcare facilities. These observations may raise concerns regarding E. avium resistance in healthcare. We need more research to understand the pathophysiology of E. avium, which leads to hospital-acquired infections.

The Potential of Dutasteride for Treating Multidrug-Resistant Candida auris Infection.

Novel antifungal drugs are urgently needed to treat candidiasis caused by the emerging fungal multidrug-resistant pathogen Candida auris. In this study, the most cost-effective drug repurposing technology was adopted to identify an appropriate option among the 1615 clinically approved drugs with anti-C. auris activity. High-throughput virtual screening of 1,3-beta-glucanosyltransferase inhibitors was conducted, followed by an analysis of the stability of 1,3-beta-glucanosyltransferase drug complexes and 1,3-beta-glucanosyltransferase-dutasteride metabolite interactions and the confirmation of their activity in biofilm formation and planktonic growth. The analysis identified dutasteride, a drug with no prior antifungal indications, as a potential medication for anti-auris activity in seven clinical C. auris isolates from Saudi Arabian patients. Dutasteride was effective at inhibiting biofilm formation by C. auris while also causing a significant reduction in planktonic growth. Dutasteride treatment resulted in disruption of the cell membrane, the lysis of cells, and crushed surfaces on C. auris, and significant (p-value = 0.0057) shrinkage in the length of C. auris was noted at 100,000×. In conclusion, the use of repurposed dutasteride with anti-C. auris potential can enable rapid recovery in patients with difficult-to-treat candidiasis caused by C. auris and reduce the transmission of nosocomial infection.

Genomic Insights into Virulence Factors and Multi-Drug Resistance in Clostridium perfringens IRMC2505A.

Clostridium perfringens is a spore-forming, Gram-positive anaerobic pathogen that causes several disorders in humans and animals. A multidrug-resistant Clostridium strain was isolated from the fecal sample of a patient who was clinically suspected of gastrointestinal infection and had a recent history of antibiotic exposure and diarrhea. The strain was identified by 16s rRNA sequencing as Clostridium perfringens. The strain's pathogenesis was analyzed through its complete genome, specifically antimicrobial resistance-related genes. The Clostridium perfringens IRMC2505A genome contains 19 (Alr, Ddl, dxr, EF-G, EF-Tu, folA, Dfr, folP, gyrA, gyrB, Iso-tRNA, kasA, MurA, rho, rpoB, rpoC, S10p, and S12p) antibiotic-susceptible genetic species according to the k-mer-based detection of antimicrobial resistance genes. Genome mapping using CARD and VFDB databases revealed significant (p-value = 1 × 10-26) genes with aligned reads against antibiotic-resistant genes or virulence factors, including phospholipase C, perfringolysin O, collagenase, hyaluronidase, alpha-clostripain, exo-alpha-sialidase, and sialidase activity. In conclusion, this is the first report on C. perfringens from Saudi Arabia that conducted whole genome sequencing of IRMC2505A and confirmed the strain as an MDR bacterium with several virulence factors. Developing control strategies requires a detailed understanding of the epidemiology of C. perfringens, its virulence factors, and regional antimicrobial resistance patterns.

Exome-wide analysis identify multiple variations in olfactory receptor genes (OR12D2 and OR5V1) associated with autism spectrum disorder in Saudi females.

Background: Autism Spectrum Disorder (ASD) is a multifactorial, neurodevelopmental disorder, characterized by deficits in communication, restricted and repetitive behaviors. ASD is highly heritable in Saudi Arabia; indecencies of affected individuals are increasing. Objectives: To identify the most significant genes and SNPs associated with the increased risk of ASD in Saudi females to give an insight for early diagnosis. Methods: Pilot case-control study mostly emphasized on the significant SNPs and haplotypes contributing to Saudi females with ASD patients (n = 22) compared to controls (n = 51) without ASD. With the use of allelic association analysis tools, 243,345 SNPs were studied systematically and classified according to their significant association. The significant SNPs and their genes were selected for further investigation for mapping of ASD candidate causal variants and functional impact. Results: In females, five risk SNPs at p ≤ 2.32 × 10-05 was identified in association with autism. The most significant exonic variants at chromosome 6p22.1 with olfactory receptor genes (OR12D2 and OR5V1) clustered with high linkage disequilibrium through haplotyping analysis. Comparison between highly associated genes (56 genes) of male and female autistic patients with female autistic samples revealed that 39 genes are unique biomarkers for Saudi females with ASD. Conclusion: Multiple variations in olfactory receptor genes (OR5V1 and OR12D2) and single variations on SPHK1, PLCL2, AKAP9 and LOC107984893 genes are contributing to ASD in females of Arab origin. Accumulation of these multiple predisposed coding SNPs can increase the possibility of developing ASD in Saudi females.

Mining the nanotube-forming Bacillus amyloliquefaciens MR14M3 genome for determining anti-Candida auris and anti-Candida albicans potential by pathogenicity and comparative genomics analysis.

There is a global health concern associated with the emergence of the multidrug-resistant (MDR) fungus Candida auris, which has significant mortality rates. Finding innovative and distinctive anti-Candida compounds is essential for treating infections caused by MDR C. auris. A bacterial strain with anti-Candida activity was isolated and identified using 16 S rRNA gene sequencing. The whole genome was sequenced to identify biosynthesis-related gene clusters. The pathogenicity and cytotoxicity of the isolate were analyzed in Candida and HFF-1 cell lines, respectively. This study set out to show that whole-genome sequencing, cytotoxicity testing, and pathogenicity analysis combined with genome mining and comparative genomics can successfully identify biosynthesis-related gene clusters in native bacterial isolates that encode antifungal natural compounds active against Candida albicans and C. auris. The native isolate MR14M3 has the ability to inhibit C. auris (zone of inhibition 25 mm) and C. albicans (zone of inhibition 25 mm). The 16 S rRNA gene sequence of MR14M3 aligned with Bacillus amyloliquefaciens with similarity (100%). Bacillus amyloliquefaciens MR14M3 establishes bridges of intercellular nanotubes (L 258.56 ± 35.83 nm; W 25.32 ± 6.09 nm) connecting neighboring cells. Candida cell size was reduced significantly, and crushed phenotypes were observed upon treatment with the defused metabolites of B. amyloliquefaciens MR14M3. Furthermore, the pathogenicity of B. amyloliquefaciens MR14M3 on Candida cells was observed through cell membrane disruption and lysed yeast cells. The whole-genome alignment of the MR14M3 genome (3981,643 bp) using 100 genes confirmed its affiliation with Bacillus amyloliquefaciens. Genome mining analysis revealed that MR14M3-coded secondary metabolites are involved in the biosynthesis of polyketides (PKs) and nonribosomal peptide synthases (NRPSs), including 11 biosynthesis-related gene clusters with one hundred percent similarity. Highly conserved biosynthesis-related gene clusters with anti-C. albicans and anti-C. auris potentials and cytotoxic-free activity of B. amyloliquefaciens MR14M3 proposes the utilization of Bacillus amyloliquefaciens MR14M3 as a biofactory for an anti-Candida auris and anti-C. albicans compound synthesizer.

Genomic Landscape of Multidrug Resistance and Virulence in Enterococcus faecalis IRMC827A from a Long-Term Patient.

We report on a highly virulent, multidrug-resistant strain of Enterococcus faecalis IRMC827A that was found colonizing a long-term male patient at a tertiary hospital in Khobar, Saudi Arabia. The E. faecalis IRMC827A strain carries several antimicrobial drug resistance genes and harbours mobile genetic elements such as Tn6009, which is an integrative conjugative element that can transfer resistance genes between bacteria and ISS1N via an insertion sequence. Whole-genome-sequencing-based antimicrobial susceptibility testing on strains from faecal samples revealed that the isolate E. faecalis IRMC827A is highly resistant to a variety of antibiotics, including tetracycline, doxycycline, minocycline, dalfopristin, virginiamycin, pristinamycin, chloramphenicol, streptomycin, clindamycin, lincomycin, trimethoprim, nalidixic acid and ciprofloxacin. The isolate IRMC827A carries several virulence factors that are significantly associated with adherence, biofilm formation, sortase-assembled pili, manganese uptake, antiphagocytosis, and spreading factor of multidrug resistance. The isolate also encompasses two mutations (G2576T and G2505A) in the 23S rRNA gene associated with linezolid resistance and three more mutations (gyrA p.S83Y, gyrA p.D759N and parC p.S80I) of the antimicrobial resistance phenotype. The findings through next-generation sequencing on the resistome, mobilome and virulome of the isolate in the study highlight the significance of monitoring multidrug-resistant E. faecalis colonization and infection in hospitalized patients. As multidrug-resistant E. faecalis is a serious pathogen, it is particularly difficult to treat and can cause fatal infections. It is important to have quick and accurate diagnostic tests for multidrug-resistant E. faecalis, to track the spread of multidrug-resistant E. faecalis in healthcare settings, and to improve targeted interventions to stop its spread. Further research is necessary to develop novel antibiotics and treatment strategies for multidrug-resistant E. faecalis infections.

The Effectiveness of Nigella sativa and Ginger as Appetite Suppressants: An Experimental Study on Healthy Wistar Rats.

Background: Obesity is a global pandemic that is associated with high morbidity and mortality. Natural herbs are commonly used for weight reduction and appetite suppression. Therefore, we aim to investigate the role and mechanism of Nigella sativa (NS) and ginger on weight reduction and appetite regulation. Methods: This experimental study was performed at Imam Abdulrahman Bin Faisal University. Twenty-five female rats were distributed into 5 groups: NS (oral 1000mg/kg), Ginger (500 mg/kg), NS-ginger (both interventions), a positive control (intraperitoneal 50 µg/kg Liraglutide), and a negative control. Each intervention was given for 9 weeks. Food intake and body weight were assessed weekly. Serum lipid profile and peptides involved in appetite control (cholecystokinin (CCK), glucagon-like peptide 1(GLP-1), gastric inhibitory polypeptide (GIP), ghrelin, peptide YY, and orexin) were assayed at the end of the experiment. Results: None of the interventions showed a statistically significant difference regarding food consumption or weight gain (p > 0.05). However, the three interventions significantly reduced total cholesterol (TC), NS and NS-ginger significantly increased HDL, NS increased ghrelin and ginger increased orexin. Conclusion: The present dose and duration of NS, ginger, or in combination did not demonstrate a significant change in body weight or food consumption in comparison to the negative or positive controls. However, NS or ginger has improved the lipid profile by reducing TC and increasing HDL. In addition, NS or ginger can influence some of the peptides involved in appetite regulation such as the increase in ghrelin induced by NS and the reduction of orexin induced by ginger. We believe that these latter effects are novel and might indicate a promising effect of these natural products on appetite regulation.

Emerging of composition variations of SARS-CoV-2 spike protein and human ACE2 contribute to the level of infection: in silico approaches.

SARS-CoV-2 is causative of pandemic COVID-19. There is a sequence similarity between SARS-CoV-2 and SARS-CoV; however, SARS-CoV-2 RBDs (receptor-binding domain) binds 20-fold strongly with human angiotensin-converting enzyme 2 (hACE2) than SARS-CoV. The study aims to investigate protein-protein interactions (PPI) of hACE2 with SARS-CoV-2 RBD between wild and variants to detect the most influential interaction. Variants of hACE2 were retrieved from NCBI and subjected to determine the most pathogenic nsSNPs. Probability of PPIs determines the binding affinity of hACE2 genetic variants with RBD was investigated. Composition variations at the hACE2 and RBD were processed for PatchDock and refined by FireDock for the PPIs. Twelve nsSNPs were identified as the top pathogenic from SNPs (n = 7489) in hACE2 using eight bioinformatics tools. Eight RBD variants were complexed with 12 nSNPS of hACE2, and the global energy scores (Kcal/mol) were calculated and classified as very weak (-3.93 to -18.43), weak (-18.42 to -32.94), moderate (-32.94 to -47.44), strong (-47.44 to -61.95) and very strong (-61.95 to -76.46) zones. Seven composition variants in the very strong zone [G726R-G476S; R768W-V367F; Y252N-V483A; Y252N-V367F; G726R-V367F; N720D-V367F and N720D-F486L], and three in very weak [P263S-S383C; RBD-H378R; G726R-A348T] are significantly (p < 0.00001) varied for global energy score. Zonation of the five zones was established based on the scores to differentiate the effect of hACE2 and RBD variants on the binding affinity. Moreover, our findings support that the combination of hACE2 and RBD is key players for the risk of infection that should be done by further laboratory studies.Communicated by Ramaswamy H. Sarma.

Safety of Tocilizumab in COVID-19 Patients and Benefit of Single-Dose: The Largest Retrospective Observational Study.

Severe acute respiratory coronavirus-2 (SARS-CoV-2) still presents a public threat and puts extra strain on healthcare facilities. Without an effective antiviral drug, all available treatment options are considered supportive. Tocilizumab as a treatment option has to date shown variable results. In this retrospective study, we aimed to assess predictors of mortality of COVID-19 patients (n = 300) on tocilizumab and the clinical effectiveness of this drug. The results showed that ICU admission OR = 64.6 (95% CI: 8.2, 507.4); age of the patient OR = 1.1 (95% CI: 1.0, 1.1); and number of tocilizumab doses administered by the patient OR(two doses) = 4.0 (95% CI: 1.5, 10.9), OR(three doses) = 1.5 (95% CI: 0.5, 5.1), and OR(four doses or more) = 7.2 (95% CI: 2.0, 25.5) presented strong correlation factors that may be linked to COVID-19 mortality. Furthermore, our study showed the beneficial effects of early administration of tocilizumab OR = 1.2 (95% CI: 1.1, 1.4) and longer hospital length of stay OR = 0.974 (95% CI: 0.9, 1.0) in reducing COVID-19 mortalities. High blood D-dimer concentration OR = 1.1 (95% CI: 1.0, 1.2) and reciprocal blood phosphate concentration OR = 0.008 (95% CI: 0.0, 1.2) were correlated to high mortality under SARS-CoV-2 infection. The short-term effect of a single dose of tocilizumab was a significant increase in blood BUN and liver enzymes (ALT, AST, and LDH) above their normal ranges. Furthermore, it significantly reduced CRP blood concentration, but not to normal levels (13.90 to 1.40 mg/dL, p < 0.001). Assessing the effect of different doses of tocilizumab (in terms of the number of doses, total mg, and total mg/kg administered by the patients) indicated that administering more than one dose may lead to increases in ICU length of stay and hospital length of stay of up to 14 and 22 days after the last dose of tocilizumab (6 to 14, p = 0.06, and 10 to 22, p < 0.001), with no improvement in 28- and 90-day mortality, as confirmed by Kaplan−Meier analysis. There were also clear correlations and trends between the number of doses of tocilizumab and increased blood CO2, MCV, RDW, and D-dimer concentrations and between number of doses of tocilizumab and decreased CRP, AST, and hemoglobin concentrations. Microbiology analysis showed a significant increase in the incidence of infection after tocilizumab administration (28 to 119, p < 0.001) with a median time of incidence within 6 days of the first dose of tocilizumab. A significant correlation was also found between the number of tocilizumab doses and the number of incidences of infections after tocilizumab administration r (298) = 0.396, p = 1.028 × 10−12. Based on these results and depending on the pharmacokinetic parameters of the drug, we recommend single-dose administration of tocilizumab as the optimal dosage for COVID-19 patients who do not have active bacterial infection or liver diseases, to be administered as soon as the patient is admitted to the hospital.

Novel Feather Degrading Keratinases from Bacillus cereus Group: Biochemical, Genetic and Bioinformatics Analysis.

In this study, five keratinolytic bacteria were isolated from poultry farm waste of Eastern Province, Saudi Arabia. The highest keratinase activity was obtained at 40-45 °C, pH 8-9, feather concentration 0.5-1%, and using white chicken feather as keratin substrate for 72 h. Enhancement of keratinase activity through physical mutagen UV radiation and/or chemical mutagen ethyl methanesulfonate (EMS) resulted in five mutants with 1.51-3.73-fold increased activity over the wild type. When compared with the wild type, scanning electron microscopy validated the mutants' effectiveness in feather degradation. Bacterial isolates are classified as members of the S8 family peptidase Bacillus cereus group based on sequence analysis of the 16S rRNA and keratinase genes. Interestingly, keratinase KerS gene shared 95.5-100% identity to keratinase, thermitase alkaline serine protease, and thermophilic serine protease of the B. cereus group. D137N substitution was observed in the keratinase KerS gene of the mutant strain S13 (KerS13uv+ems), and also seven substitution variations in KerS26 and KerS26uv of strain S26 and its mutant S26uv. Functional analysis revealed that the subtilisin-like serine protease domain containing the Asp/His/Ser catalytic triad of KerS gene was not affected by the predicted substitutions. Prediction of physicochemical properties of KerS gene showed instability index between 17.5-19.3 and aliphatic index between 74.7-75.7, which imply keratinase stability and significant thermostability. The docking studies revealed the impact of substitutions on the superimposed structure and an increase in binding of mutant D137N of KerS13uv+ems (affinity: -7.17; S score: -6.54 kcal/mol) and seven mutants of KerS26uv (affinity: -7.43; S score: -7.17 kcal/mol) compared to the wild predicted structure (affinity: -6.57; S score: -6.68 kcal/mol). Together, the keratinolytic activity, similarity to thermostable keratinases, and binding affinity suggest that keratinases KerS13uv+ems and KerS26uv could be used for feather processing in the industry.

The association between two genetic polymorphisms in ITGB3 and increase risk of venous thromboembolism in cancer patients in Eastern Province of Saudi Arabia.

Venous thromboembolism (VTE) is one of the major complications in most cancer patients leading to poor prognosis and short survival. Several common clinical risk factors coexist in cancer patients are used as risk predictive biomarkers to help in the management and prevention of VTE. These include cancer site and stage, chemotherapy regimen and elevated biological markers. However, Genetic polymorphisms in genes controlling coagulation and fibrinolysis are significantly associated with VTE if detected, then they might be more sensitive individual predictive biomarkers for VTE risk assessment. This study was conducted to evaluate the association between ITGB3 rs3809865 and rs5918 with VTE risk as well as monitor the effect of VTE on overall survival of these cancer patients. In this retrospective case-control study, 195 cancer patients' formalin-fixed paraffin embedded tissue (FFPE) samples were collected (controls n = 157, case n = 38) using the stored data through Jan 2010 to Sep 2018 from King Fahad Specialist Hospital in Dammam. Samples were genotyped using TaqMan genotyping assay, then logistic regression analysis and Chi-square were used to predict the association between risk factors and VTE. Survival Comparison was tested by the log-rank test. Genetic polymorphisms in ITGB3 (rs3809865 and rs5918) found not to be associated with VTE increasing risk in cancer patients (p>0.05). While the advanced stage was potentially increasing the risk of VTE events (OR 5.1 CI 2.01-12.9p = 0.001). Patients with VTE showed a poor overall survival reflected by the median survival rate of only three years compared to seven years for cancer patients without VTE. This study highlighted the potential influence of VTE on prognosis and survival of cancer patients and raised the importance of exploring risk predictive biomarkers in our population. This will improve the risk prediction biomarkers leading to implementing safe and effective thrombosis prophylaxis strategies.

A Pilot Mitochondrial Genome-Wide Association on Migraine Among Saudi Arabians.

Background: Mitochondrial DNA (mtDNA) mutations have been reported in multiple neurological diseases and helped to explain the pathophysiology of these diseases. Similarly, variations in mtDNA might exist in migraine and can explain the effect of low ATP production in the neurons on the initiation of migraine attack. Therefore, in the current study we aim to explore the association of mtDNA mutations on migraine in the Saudi population. Subjects and Methods: Over 1950 young Saudi female students were screened for migraine, among that a total of 103 satisfied the ICHD-3 criteria. However,  20 migraine cases confirmed in the neurology clinic and gave consent to participate in the study. Another 20 age-matched healthy controls were also recruited. Mitochondrial sequence variations were filtered from exome sequencing using NCBI GenBank Reference Sequence: NC_012920.1 and analysed using MITOMAP. Genes with significant single nucleotide polymorphisms (SNPs) were investigated by the gene functional classification tool DAVID and functional enrichment analysis of protein-protein interaction networks through STRING 11.5 for the most significant associated genes. Results: Genome wide analysis of the mitochondrial sequence variations between the patients with migraine and control revealed the association of 30 SNPs (p < 0.05) in the mitochondrial genome. The highest significance (p = 0.001033) was observed in a coding SNP (rs1603225278) in the CYTB gene and rs386829281 in the region of origin of replication. Twenty-four significant SNPs were in the coding region of nine (ND5, ND4, COX2, COX1, ND3, CYTB, COX3, ND2 and ND1) genes. Conclusion: This is the first study to demonstrate the association of mtDNA variations with migraine in the Saudi population. The current findings will help to highlight the significance of mtDNA mutations to migraine pathophysiology and will serve as a reference data for larger national and international studies.