Search

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.

Background Testing strategies is crucial for genetics clinics and testing laboratories. In this study, we tried to compare the hit rate between solo and trio and trio plus testing and between trio and sibship testing. Finally, we studied the impact of extended family analysis, mainly in complex and unsolved cases. Methods Three cohorts were used for this analysis: one cohort to assess the hit rate between solo, trio and trio plus testing, another cohort to examine the impact of the testing strategy of sibship genome vs trio-based analysis, and a third cohort to test the impact of an extended family analysis of up to eight family members to lower the number of candidate variants. Results The hit rates in solo, trio and trio plus testing were 39, 40, and 41%, respectively. The total number of candidate variants in the sibship testing strategy was 117 variants compared to 59 variants in the trio-based analysis. We noticed that the average number of coding candidate variants in trio-based analysis was 1192 variants and 26,454 noncoding variants, and this number was lowered by 50–75% after adding additional family members, with up to two coding and 66 noncoding homozygous variants only, in families with eight family members. Conclusion There was no difference in the hit rate between solo and extended family members. Trio-based analysis was a better approach than sibship testing, even in a consanguineous population. Finally, each additional family member helped to narrow down the number of variants by 50–75%. Our findings could help clinicians, researchers and testing laboratories select the most cost-effective and appropriate sequencing approach for their patients. Furthermore, using extended family analysis is a very useful tool for complex cases with novel genes.

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.

Familial thrombocytosis (FT) is a rare hereditary haematological disorder characterised by increased platelet count, usually caused by germ-line mutations in thrombopoietin (THPO), myeloproliferative leukaemia virus oncogene (MPL) or Janus kinase 2 (JAK2) genes, and can be associated with increased risk of thrombosis. We aimed to determine the yield of diagnostic tests, assess treatment received and describe the clinical course of MPL-associated FT. We retrospectively reviewed all paediatric and adult haematology patients diagnosed with MPL-related FT, who were seen in our clinics from March 2013 to February 2021. Of 64 eligible patients, 26 (41%) were aged

Whole-exome sequencing (WES) and whole-genome sequencing (WGS) are used to diagnose genetic and inherited disorders. However, few studies comparing the detection rates of WES and WGS in clinical settings have been performed. Variant call format files were generated and raw data analysis was performed in cases in which the final molecular results showed discrepancies. We classified the possible explanations for the discrepancies into three categories: the time interval between the two tests, the technical limitations of WES, and the impact of the sequencing system type. This cohort comprised 108 patients with negative array comparative genomic hybridization and negative or inconclusive WES results before WGS was performed. Ten (9%) patients had positive WGS results. However, after reanalysis the WGS hit rate decreased to 7% (7 cases). In four cases the variants were identified by WES but missed for different reasons. Only 3 cases (3%) were positive by WGS but completely unidentified by WES. In this study, we showed that 30% of the positive cases identified by WGS could be identified by reanalyzing the WES raw data, and WGS achieved an only 7% higher detection rate. Therefore, until the cost of WGS approximates that of WES, reanalyzing WES raw data is recommended before performing WGS.

Purpose Whole-exome sequencing (WES) can help identify known and novel pathogenic molecular aberrations. Here, we examined the diagnostic yield of WES in population from consanguineous unions. Methods We preformed retrospective review of multicenter WES results of an unselected cohort of patients with a wide range of phenotypes. Clinical data and WES reports of 454 patients from 5 centers across Saudi Arabia were analyzed. Testing was performed in accredited commercial laboratories, and all the WES laboratory reports were reviewed again using additional clinical information available to the treating physicians. Results Among the 454 probands, we identified highly likely disease-causing variants in 222, thereby achieving a 49% molecular diagnostic yield. The diagnostic yield was 53% in consanguineous unions and 39% in non-consanguineous unions. About 66% of the identified variants in consanguineous families were homozygous, with an autosomal recessive mode of inheritance. Additional clinical data reclassified 11 positive reported results into 4 inconclusive and 7 negative results, and 22 inconclusive results into 17 positive and 5 negative results. Conclusions The diagnostic yield from WES in our unselected cohort is similar to other studies from the same region, which is a higher yield compared to other international regions largely because of the high rate of consanguinity and partly due to simplified variant interpretation and classification in consanguineous unions.

SHFM6 (OMIM 225300) is caused by WNT10B pathogenic variants (12q13.12). It is one of the rarest forms of SHFM; with only seven pathogenic variants described in the world literature. Furthermore, it has not been determined if SHFM6 has specific phenotypic characteristics. In this paper, we present a case series of three unrelated families with SHFM6 caused by three novel WNT10B pathogenic variants. The index patient of the first family was homozygous for the nonsense variant c.676C > T (p.Arg226*) in the WNT10B gene. The index case of the second family had a homozygous splice variant c.338-1G > C in the WNT10B gene. Finally, the index case of the third family carried two different variants in the WNT10B gene: A nonsense variant (p.Arg226*), and a missense variant (p.Gln86Pro). The latter represents the first compound heterozygous pathogenic variant related to SHFM6. We also offer a classification system for the hand/foot defects to illustrate the specific phenotypic characteristics of SHFM6. Based on this classification and a review of all previously reported cases, we demonstrate that SHFM6 caused by WNT10B pathogenic variants have the following characteristics: more severe feet defects (compared to the hand defects), polydactyly, severe flexion digital contractures, and phalangeal dysplasia.

BackgroundDuring critical illness in humans, the effects of caloric restriction on the inflammatory response are not well understood. The aim of this study is to examine the associations of caloric restriction, inflammatory response profiles and outcomes in critically ill patients.MethodsThis is a sub-study of the PermiT trial (Permissive Underfeeding or Standard Enteral Feeding in Critically Ill Adults Trial- ISRCTN68144998). Serum samples were collected on study days 1, 3, 5, 7 and 14 and analyzed for a panel of 29 cytokines. We used principal component analysis to convert possibly correlated variables (cytokine levels) into a limited number of linearly uncorrelated variables (principal components). We constructed repeated measures mixed linear models to assess whether permissive underfeeding compared to standard feeding was associated with difference cytokine levels over time.ResultsA total of 72 critically ill patients were enrolled in this study (permissive underfeeding n = 36 and standard feeding n = 36). Principal component analysis identified 6 components that were responsible for 78% of the total variance. When adjusted to principal components, permissive underfeeding was not associated with 90-day mortality (adjusted odds ratio 1.75, 95% confidence interval 0.44, 6.95, p = 0.43) or with incident renal replacement therapy. The cytokines did not differ with time between permissive underfeeding and standard feeding groups.ConclusionsThe association of permissive underfeeding compared to standard feeding with mortality was not influenced by the inflammatory profile. Permissive underfeeding compared to standard feeding was not associated with differences in the serum levels of cytokines in critically ill patients.

We report our results of 1000 diagnostic WES cases based on 2819 sequenced samples from 54 countries with a wide phenotypic spectrum. Clinical information given by the requesting physicians was translated to HPO terms. WES processes were performed according to standardized settings. We identified the underlying pathogenic or likely pathogenic variants in 307 families (30.7%). In further 253 families (25.3%) a variant of unknown significance, possibly explaining the clinical symptoms of the index patient was identified. WES enabled timely diagnosing of genetic diseases, validation of causality of specific genetic disorders of PTPN23, KCTD3, SCN3A, PPOX, FRMPD4, and SCN1B, and setting dual diagnoses by detecting two causative variants in distinct genes in the same patient. We observed a better diagnostic yield in consanguineous families, in severe and in syndromic phenotypes. Our results suggest that WES has a better yield in patients that present with several symptoms, rather than an isolated abnormality. We also validate the clinical benefit of WES as an effective diagnostic tool, particularly in nonspecific or heterogeneous phenotypes. We recommend WES as a first-line diagnostic in all cases without a clear differential diagnosis, to facilitate personal medical care.

We reviewed our institutional experience from 2011 to 2015 on new cases of Fanconi anemia (FA). Ten unrelated cases were diagnosed during this period. Four patients with severe aplastic anemia (SAA) had c.2392C > T (p.Arg798*) BRIP1/FANCJ mutation. Another child with SAA had novel c.1475T > C (p.Leu492Pro) FANCC mutation. One individual with SAA and acute myeloid leukemia had c.637_643del (p.Tyr213Lysfs*6) FANCG mutation. Three patients presented with early onset of cancer, two had BRCA2 mutation c.7007G > A (p.Arg2336His) and one had a novel c.3425del (p.Leu1142Tyrfs*21) PALB2 mutation. Another infant with c.3425del PALB2 mutation had clonal aberration with partial trisomy of the long arm of chromosome 17. Mutations in FA downstream pathway genes are more frequent in our series than expected. Our preliminary observation will be confirmed in a large multi-institutional study.

The effect of short-term caloric restriction on gene expression in critically ill patients has not been studied. In this sub-study of the PermiT trial (Permissive Underfeeding or Standard Enteral Feeding in Critically Ill Adults Trial- ISRCTN68144998), we examined gene expression patterns in peripheral white blood cells (buffy coat) associated with moderate caloric restriction (permissive underfeeding) in critically ill patients compared to standard feeding. Blood samples collected on study day 1 and 14 were subjected to total RNA extraction and gene expression using microarray analysis. We enrolled 50 patients, 25 in each group. Among 1751 tested genes, 332 genes in 12 pathways were found to be significantly upregulated or downregulated between study day 1 and 14 (global p value for the pathway ≤ 0.05). Using the heatmap, the differential expression of genes from day 1 to 14 in the permissive underfeeding group was compared to the standard feeding group. We further compared gene expression signal intensity in permissive underfeeding compared standard feeding by constructing univariate and multivariate linear regression models on individual patient data. We found differential expression of several genes with permissive underfeeding, most notably those related to metabolism, autophagy and other cellular functions, indicating that moderate differences in caloric intake trigger different cellular pathways.

Background Scientists have previously described an overgrowth syndrome in Saudi patients and named it ‘Upper limb muscle overgrowth with hypoplasia of the index finger’ syndrome. Case presentation We describe a new case and document that the syndrome is caused by the somatic PIK3CA mutation c.3140A>G, p.His1047Arg. We also recruited one of the previously reported cases and found the same somatic mutation in the affected muscles. A wider classification of ‘PIK3CA-related pathology spectrum’ is presented which includes cancer, benign skin lesions/tumors, Cowden syndrome, isolated vascular malformations and various overgrowth syndromes. The latter entity is sub-divided into 3 sub-groups: overgrowth with brain involvement, overgrowth with multiple lipomatosis, and overgrowth without brain involvement/multiple lipomatosis. Conclusion Our literature review indicated that “upper limb muscle overgrowth with hypoplasia of the index finger” is not as rare as previously thought to be. This overgrowth syndrome is unique and is caused by the somatic PIK3CA mutation c.3140A>G.

Background Hepatitis C virus (HCV) genotypes and subtypes are considered an important tool for epidemiological and clinical studies and valuable markers for disease progression and response to antiviral therapy. The aim of this study was to identify the prevalence of HCV genotypes and their relation to socio-demographic factors particularly age and sex, various biochemical profiles and viral load. Methods The records (630) of Saudi patients positive for HCV (2007–2011) reported in the system of the Molecular Pathology Laboratory at a tertiary reference hospital in Riyadh, Saudi Arabia were analyzed. Socio-demographic characteristics, liver biochemical profile, viral load and co-infection with HBV and HIV were retrieved from the hospital database. The associations of continuous and categorical variables with genotypes were analyzed. Result The overall mean age of the surveyed patients was 59 years ±0.5 years (21% were

Background The use of cord blood in the neonatal screening for glucose-6-phosphate dehydrogenase (G6PD) deficiency is being done with increasing frequency but has yet to be adequately evaluated against the use of peripheral blood sample which is usually employed for confirmation. We sought to determine the incidence and gender distribution of G6PD deficiency, and compare the results of cord against peripheral blood in identifying G6PD DEFICIENCY neonates using quantitative enzyme activity assay. Methods We carried out a retrospective and cross-sectional study employing review of primary hospital data of neonates born in a tertiary care center from January to December 2008. Results Among the 8139 neonates with cord blood G6PD assays, an overall incidence of 2% for G6PD deficiency was computed. 79% of these were males and 21% were females with significantly more deficient males (p

We report the results of clinical exome sequencing (CES) on >2,200 previously unpublished Saudi families as a first-tier test. The predominance of autosomal-recessive causes allowed us to make several key observations. We highlight 155 genes that we propose to be recessive, disease-related candidates. We report additional mutational events in 64 previously reported candidates (40 recessive), and these events support their candidacy. We report recessive forms of genes that were previously associated only with dominant disorders and that have phenotypes ranging from consistent with to conspicuously distinct from the known dominant phenotypes. We also report homozygous loss-of-function events that can inform the genetics of complex diseases. We were also able to deduce the likely causal variant in most couples who presented after the loss of one or more children, but we lack samples from those children. Although a similar pattern of mostly recessive causes was observed in the prenatal setting, the higher proportion of loss-of-function events in these cases was notable. The allelic series presented by the wealth of recessive variants greatly expanded the phenotypic expression of the respective genes. We also make important observations about dominant disorders; these observations include the pattern of de novo variants, the identification of 74 candidate dominant, disease-related genes, and the potential confirmation of 21 previously reported candidates. Finally, we describe the influence of a predominantly autosomal-recessive landscape on the clinical utility of rapid sequencing (Flash Exome). Our cohort’s genotypic and phenotypic data represent a unique resource that can contribute to improved variant interpretation through data sharing.

Axial spondylometaphyseal dysplasia (axial SMD) is an autosomal recessive disease characterized by dysplasia of axial skeleton and retinal dystrophy. We conducted whole exome sequencing and identified C21orf2 (chromosome 21 open reading frame 2) as a disease gene for axial SMD. C21orf2 mutations have been recently found to cause isolated retinal degeneration and Jeune syndrome. We found a total of five biallelic C21orf2 mutations in six families out of nine: three missense and two splicing mutations in patients with various ethnic backgrounds. The pathogenic effects of the splicing (splice-site and branch-point) mutations were confirmed on RNA level, which showed complex patterns of abnormal splicing. C21orf2 mutations presented with a wide range of skeletal phenotypes, including cupped and flared anterior ends of ribs, lacy ilia and metaphyseal dysplasia of proximal femora. Analysis of patients without C21orf2 mutation indicated genetic heterogeneity of axial SMD. Functional data in chondrocyte suggest C21orf2 is implicated in cartilage differentiation. C21orf2 protein was localized to the connecting cilium of the cone and rod photoreceptors, confirming its significance in retinal function. Our study indicates that axial SMD is a member of a unique group of ciliopathy affecting skeleton and retina.

Raine syndrome is an autosomal recessive disorder caused by mutations in the FAM20C gene that is characterized by generalized osteosclerosis with periosteal new bone formation and distinctive craniofacial dysmorphism. We report on a child who is homozygous for a 487-kb deletion in 7p22.3 that contains FAM20C. Both parents were heterozygous for the deletion. Our patient had the common craniofacial features as well as, uncommon features such as protruding tongue, short stature, and hypoplastic distal phalanges. In addition, he had wormian bones and pyriform aperture stenosis, features that are usually under diagnosed. It is clear that Raine syndrome has a wide range of expression and may not be lethal in the neonatal period. Furthermore, Raine cases due to whole gene deletion do not seem to have a major difference in the phenotype over those caused by various mutations.

Prior work has suggested that loss of expression of one or more of the many C/D box small nucleolar RNAs (snoRNAs) encoded within the complex, paternally expressed SNRPN (small nuclear ribonuclear protein N) locus may result in the phenotype of Prader-Willi syndrome (PWS). We suggest that the minimal critical region for PWS is approximately 121 kb within the460-kb SNRPN locus, bordered by a breakpoint cluster region identified in three individuals with PWS who have balanced reciprocal translocations and by the proximal deletion breakpoint of a familial deletion found in an unaffected mother, her three children with Angelman syndrome, and her father. The subset of SNRPN-encoded snoRNAs within this region comprises the PWCR1/HBII-85 cluster of snoRNAs and the single HBII-438A snoRNA. These are the only known genes within this region, which suggests that loss of their expression may be responsible for much or all of the phenotype of PWS. This hypothesis is challenged by findings in two individuals with PWS who have balanced translocations with breakpoints upstream of the proposed minimal critical region but whose cells were reported to express transcripts within it, adjacent to these snoRNAs. By use of real-time quantitative reverse-transcriptase polymerase chain reaction, we reassessed expression of these transcripts and of the snoRNAs themselves in fibroblasts of one of these patients. We find that the transcripts reported to be expressed in lymphoblast-somatic cell hybrids are not expressed in fibroblasts, and we suggest that the original results were misinterpreted. Most important, we show that the PWCR1/HBII-85 snoRNAs are not expressed in fibroblasts of this individual. These results are consistent with the hypothesis that loss of expression of the snoRNAs in the proposed minimal critical region confers much or all of the phenotype of PWS.

BackgroundPrader-Willi syndrome (MIM #176270; PWS) is caused by lack of the paternally-derived copies, or their expression, of multiple genes in a 4 Mb region on chromosome 15q11.2. Known mechanisms include large deletions, maternal uniparental disomy or mutations involving the imprinting center.De novobalanced reciprocal translocations in 5 reported individuals had breakpoints clustering inSNRPNintron 2 or exon 20/intron 20. To further dissect the PWS phenotype and define the minimal critical region for PWS features, we have studied a 22 year old male with a milder PWS phenotype and ade novotranslocation t(4;15)(q27;q11.2).MethodsWe used metaphase FISH to narrow the breakpoint region and molecular analyses to map the breakpoints on both chromosomes at the nucleotide level. The expression of genes on chromosome 15 on both sides of the breakpoint was determined by RT-PCR analyses.ResultsPertinent clinical features include neonatal hypotonia with feeding difficulties, hypogonadism, short stature, late-onset obesity, learning difficulties, abnormal social behavior and marked tolerance to pain, as well as sticky saliva and narcolepsy. Relative macrocephaly and facial features are not typical for PWS. The translocation breakpoints were identified withinSNRPNintron 17 and intron 10 of a spliced non-coding transcript in band 4q27. LINE and SINE sequences at the exchange points may have contributed to the translocation event. By RT-PCR of lymphoblasts and fibroblasts, we find that upstream SNURF/SNRPN exons and snoRNAs HBII-437 and HBII-13 are expressed, but the downstream snoRNAs PWCR1/HBII-85 and HBII-438A/B snoRNAs are not.ConclusionAs part of the PWCR1/HBII-85 snoRNA cluster is highly conserved between human and mice, while no copy of HBII-438 has been found in mouse, we conclude that PWCR1/HBII-85 snoRNAs is likely to play a major role in the PWS- phenotype.