It is commonly accepted that genetic risks for psychiatric disorders, such as attention-deficit hyperactivity disorder (ADHD), schizophrenia, obsessive-compulsive disorder (OCD), bipolar disorder and autism, are conferred by mutations in multiple genes. Attention-deficit hyperactivity disorder (ADHD) is one of the major public health related problems and is a neuro-developmental disorder of childhood . Worldwide about 5% of children have ADHD , and in the Kingdom of Saudi Arabia (KSA) the reported cases of ADHD are 16.4% . According to Saudi ADHD society (https://adhd.org.sa/en/), an estimated 1.6 million children suffer from ADHD in the kingdom. ADHD is defined by inappropriate levels of inattention, hyperactivity, and impulsivity and it affects boys more than girls . ADHD is associated with hypo-activation of executive functions, including neurocognitive processes, such as inhibition, has widespread learning and memory problems, task switching or set-shifting, and the inability to plan, and pays prolonged attention . The aetiology of ADHD is not clear, having multifactorial aspects ; however, recent studies have revealed genetic factors are strongly implicated in the aetiology of ADHD . However, the exact cause of ADHD is not clear; a number of factors can contribute to the development of ADHD, such as genetics, smoking, alcohol abuse, or drug use during pregnancy; exposure to environmental toxins during pregnancy, exposure to environmental toxins during growth and development, low birth weight, and brain injuries are associated with ADHD development. A recent study has indicated that there is a modest association between higher frequency of digital media use and the development of ADHD symptoms; however, further research is required to determine whether this association is causal . Executive function disorder is implicated inmost symptoms of ADHD and is involved in the frontal-subcortical circuits of the brain . On the other hand, the biochemical changes in dopaminergic, adrenergic, and serotonergic systems are associated with the appearance of the symptoms of ADHD . Therefore, the changes in the brain area and neurotransmitter levels contribute to the explanation of the neural dysfunction of ADHD . ADHD is considered one of the familial disorders .
Numerous candidate genes for ADHD were implicated by genome-wide association studies. Several studies suggest the involvement of the dopamine neurotransmitter system in the aetiology of ADHD. Specifically, Ribasés et al. (2012) have shown the association between DRD1 and childhood ADHD . According to previous genetic studies, many genes rather than a single gene were found to play a role in the development of ADHD. Polymorphisms in genes, such as VNTR, DAT1, DRD2, DRD3, DRD4, DRD5, DBH, 5-HTT, HTR1B and SNAP2, were reported to be associated with ADHD . Furthermore, genome-wide association studies recently identified common and rare variants and implicated about 391 SNPs with this disease (https://www.ebi.ac.uk/gwas/). Exome analysis has been used to successfully identify diseases causing protein variants in majority of cancers and other diseases; also recently researchers are identifying potential causative mutations associated with major psychiatric diseases such as schizophrenia (SCZ), bipolar disorder, obsessive-compulsive disorder (OCD), major depression disorder (MDD), anxiety disorders (AD), autism, and ADHD. Hawi et al. (2018) have shown that novel DNA variants in the brain derived neurotrophic factor (BDNF) gene increase risks for ADHD using next generation sequencing (NGS) conducted on an Illumina HiSeq2000 platform . Another recent study also has shown that rare and more frequent variants in multiple genes act together in contributing to ADHD risk by combining linkage analysis and whole exome sequencing (WES) .
ADHD is highly heritable in children with an estimated mean heritability of 76% . Several studies have suggested that the risk for ADHD may be higher (2 - 8 times) among the first-degree relatives . Consanguineous marriages are common in many Middle Eastern countries, including Saudi Arabia, and around 57.7% prevalence of consanguinity with first-cousin marriages is the most frequent . This may be one of the reasons for the high incidence of genetic disorders including ADHD, too. The new finding of this study is the screening of Saudi Arabian school children (age > 14 years) in the Makkah (Mecca) region for molecular identification of polymorphisms associated with ADHD response. This will provide a genetic screening map for the molecular genetic determinants of the disease in Saudi children, with the goal of reducing the impact of chronic diseases on health and the economy this kingdom.
2.1. Study Design and Clinical Presentation
This study is a retrospective case-study approved by the Medical Ethical Committee of the Faculty of Medicine, Umm Al-Qura University (Ethical approval No HAPO-02-K-012-2017-01-204) and conducted in accordance with the principles of the Declaration of Helsinki. Consent to participation in this study was obtained from the legal guardians of all patients. Four Saudi Arabian paediatric patients with ADHD were recruited in this study; the cases were selected from those scoring a significant degree of ADHD according to the Electronic Observation Questionnaire. The total number of these cases was 106, with variations in several parameters including age (6 - 14 years old), and gender (58 males 54.7% and 48 females 45.3%). Inclusion criteria of being at ages 6 - 14, and having met the criteria of primary diagnosis of ADHD based on Diagnostic and Statistical Manual of Mental Disorders (DSM-V) were performed. Any patient of abnormal body mass according to age and gender, and/or having any psychotic or other neurological disorders, and/or having any medical condition conflicting the assessment (e.g. hearing or cardiac abnormalities) were considered as an exclusion criteria.
They were either diagnosed in the clinic by a specialist in the research team, or were previously diagnosed by a specialist. All the cases have met the criteria of Attention-Deficit Hyperactivity Disorder Clinical List (see attached supplementary Table 1) that was derived from the Diagnostic and Statistical Manual of Mental Disorders―fifth edition (DSM-5) .
An electronic questionnaire that contained 22 questions was created, it included questions about the previous mental and medical history of the subjects also  . Then the questionnaire was distributed to health facilities and centres of developmental disorders in general education. One hundred and six contributors have filled the questionnaire. The parents of the 30 cases out of the 106, matched the required criteria from the electronic questionnaire and were seen by the consultant and answered questions about their child to confirm the diagnosis and determine the severity of the case. Then, the patient was observed to better determine the ADHD symptoms. The confirmation of the diagnosis was done according to the criteria of the Clinical Symptoms of ADHD list  where, 10 patients were classified as mild, 15 patients were moderate, and 5 were severe. Each case took about 45 minutes to confirm, and was categorized into one of three groups: ADHD cases, ADHD associated with Autism cases, and ADHD associated with developmental delay cases. Among the 30 cases that were seen in the clinic, 18 of the cases met the criteria of just ADHD excluding any other form of associated diseases.
Table 1. Internal consistency validity for a list of clinical symptoms of attention-deficit hyperactivity disorder.
The protocol was designed for the attention-deficit and hyperactivity disorder. It consists of 18 parameters divided into two parameters: one for attention-deficit which included 9 phrases, and another is for hyperactivity and impulsivity which included 9 phrases. All the parameters were derived from the Diagnostic and Statistical Manual of Mental Disorders-fifth edition (DSM-5). The list was designed as a 3-point scale as follows (parameter does not exist = 0, parameter is partly present = 1, parameter is fully present = 2). The highest score that may be obtained is 36, and the lowest score is 0. A high score of the protocol reflects the presence of the parameter in two places or more, which hinders the performance of the general situation. Our protocol showed high psychometric characteristics, based on the participants’ responses to a question on a 3-point frequency scale, values of the Pearson correlation coefficient as a measure of internal consistency validity were calculated between the participants’ scores on each dimension and the total score of the scale, as will shown in Table 1.
2.2. DNA Isolation and Whole Exome Sequencing
Genomic DNA was extracted from EDTA anticoagulated peripheral venous blood samples using PureLinkTM Genomic DNA Mini Kit by Applied Bio-systems (Catalog number: K182001). Genomic DNA was quantified using Qubit-dsDNA HS (High Sensitivity) Assay Kit (Thermo Scientific). Whole exome sequencing was performed for patient DNA samples using Ion Torrent technology from the Thermo Scientific. 50 - 100 ng of the genomic DNA was subjected to enriched library preparation with an Ion AmpliSeqTM Library Kit Plus kit (Cat No. 4488990) with mean amplicon size of 215 bpand ion AmpliSeqTM Exome RDY Panel 1x8 containing primer pools (Cat. No. 4489838, Life Technologies, USA) kit following the manufacturer’s protocol, and then quantified using qPCR using Ion Library TaqMan quantitation kit (Cat. No. 4468802). Each sample was assigned a distinct barcode using the Ion-Xpress barcode kit (Cat. No. 4471250). Barcoded libraries were diluted to 100 pmol, pooled samples libraries and subjected to template preparation on Ion sphere particles (ISP) using the Ion One Touch 2 System and Ion PI HI-Q Template OT2 200 kit (Cat. No. A26434). Each templated ISP was loaded on the Ion PI v3 chip (Cat. No. A26770) and sequenced on the Ion Proton instrument (Life Technologies, USA) using the Ion PI Hi-Q Sequencing Kit (Cat. No. A26772). The alignments were performed by a Torrent Suite v4.2. Filtering low quality reads, adapter removing followed alignment against the human reference genome (hg19 build); bio-informatic analysis was performed using Torrent Suite Software (v220.127.116.11), and the variants were called using the Torrent Variant Caller plugin (v5.10), and then imported into the Ion Reporter software (v5.10) for the annotation. Gene variants located within exons, introns, predictive splice sites and UTRs were identified with less than 1% minor allele frequency (MAF). Functional exons variants were also identified as frame-shift, synonymous, missense, nonsense and stop loss.
2.3. Biostatistical Analysis and Bio-Informatics
The total number of genes is 271 which were selected form the ADHD the genome-wide association study (GWAS) data base. However, only 222 genes were used in filtration in this study because 49 of the genes which did not match the ion torrent library kit (attached supplementary Table 2). ADHD included exonic variants nonsense, missense, frame-shift insertion, frame-shift deletion, and using the Ion Reporter v5.6 annotation tool for functional consequences for genes and variants.
Table 2. Summary of ADHD patients’ characteristics.
Table 1 shows values of the Pearson correlation coefficient between the Attention-Deficit dimension, hyperactivity and impulsivity dimension and the total scores, which are 0.83 and 0.85, respectively. The values are statistically significant at 0.01 p-value.
To examine the reliability of the clinical symptoms of the ADHD list, Cronbach’s coefficient between each dimension of the list and the list total score were calculated. Cronbach’s value for the instrument is 0.85 while it is 0.80 for the Attention-Deficit dimension and 0.82 for the hyperactivity and impulsivity dimension, which are statistically significant.
Selection of samples for whole exome analysis: 4 patient samples (the most severe cases) from our previously mentioned 18 cases of just ADHD were selected for exome sequencing for this study. The 4 patients’ selection was based on diagnosis, severity, and age. The criteria and patient phenotypic characteristics are summarized in Table 2.
In this study, DNA samples from four children with ADHD were sequenced for whole exome analysis (Sample 1 - 4; Table 2) samples were diagnosed for ADHD according to the Diagnostic and Statistical Manual of Mental Disorders―fifth edition (DSM-5 Diagnostic Criteria). Half of cases were male (50%) in an age range of 6 - 11 years. The average read depth for the targeted base coverage was 82×, with 94% of targeted regions covered at greater than 20×.
Almost, more than 160,000 variants per VCF files were annotated from all the samples. In this study a total of 51,920, 56,215, 55,747 and 54,664 variants were found in the children with ADHD. In Table 3 we have shown the filtered variants applying the filters to retrieve pathogenic variants, the filters including exonicvariants, sorting intolerant from tolerant (SIFT) score < 0.05 and PolyPhen score of >0.5.
Thirty-three variants were identified following filtering for potential functionality, 31 of which had been identified previously in the dbSNP and 2 were novel. The 2 novelmissense variants identified in the study were found to be of ITGA1 (chr5:52243247) c.3451G > T; p. (Ala1151Ser) in one patient, and in SPATA13 (spermatogenesis-associated protein 13) at chr13:24798055 locus c. 988G > A; p. (Ala330Thr) in another patient. Six missense variants (rs16841277, rs2228209, rs2230283, rs3741883, rs1716, rs2272606) were found in 3 patients respectively. Other seven missense variants (rs6558702, rs920829, rs149000560, rs2279574, rs11109968, rs9935113, rs7191351) were found in only two of the four patients. Eighteen other variants were identified in only one patient each (Table 3). The quality statistics of the variants is shown in Table 4. All variants show very low SIFT score indicating these substitutions in the variants have damaging effect on the protein function. Many variants have high PolyPhen score, suggesting variants with score in this range are more confidently predicted to be damaging.
SIFT score predicts protein function deleterious effects based on amino acid substitution. Scores ranges from 0.0 - 1.0, variants with scores less than 0.05 are considered as deleterious, a SIFT value of 1.0 is considered as benign. The heat map generated 56 variants from 27 genes (Figure 1).
Clustering Algorithms to visually identify patterns based upon the scoring of variant impact, the rows and columns of a heat map are often sorted by hierarchical clustering trees. Heat maps of different SNPs analysed using SIFT score and hierarchical clustering of genes in 4 ADHD patients is shown in Figure 1.
In our study we found 33 variants in 222 commonly associated genes with ADHD listed in GWAS database. Out of these 31 SNPs had been reported previously in the dbSNP and 2 SNPs were novel. All the target genes we have selected from the GWAS data base for the analysis in this study (supplementary Table 2) are known to be linked to ADHD (http://adhd.psych.ac.cn); however, we have detected many additional SNPs which are new, because previously these SNPs were not linked to ADHD. This indicates our finding is a new to demonstrate the association of known SNPs with the ADHD in Saudi children. There was no association was demonstrated with the ADHD and our newly detected SNPs in the previous study, information about 14 previously reported variants (rs13166360, rs150179494, rs6558702, rs920829, rs149000560, rs2279574, rs2230283, rs3741883, rs928661, rs2306022, rs45551636, rs9935113, rs1716, and
Table 3. Rare and novel variants identified by exome sequencing in Saudi Arabian ADHD patients.
Table 4. Quality Statistics of Rare and novel variants identified by exome sequencing in ADHD patients.
Figure 1. The details of SIFT scores.
rs2272606). The 2 novel missense variants identified in the study were found to be encoded for ITGA1 (chr5:52243247) c.3451G > T; p. Ala1151Ser in one patient, and in SPATA13 (chr13:24798055) gene c. 988G > A; p. Ala330Thr in one patient, this gene also known as ASEF2 and serves as a GEF for both CDC42 and RAC1 in cell culture and that by binding to APC (adenomatosis polyposis coli) increases its GEF activity, and modulates the cell migration . Also, we have found a missense5’-UTR variant G > A (rs41287016) in SPATA13 gene which has been reported previously. One SNP in ADCY2 (Adenylate cyclase-2) rs13166360 is associated with bipolar disorder . Another SNP we found in the gene DUSP6 (rs2279574) also known to be associated with bipolar disorder .
The identified polymorphisms in this study have high variant frequencies (0.25 - 1). In addition, all amino acid substitutions for the mutant nucleotide are predicted deleterious according to the SIFT scores (≥0.05). One variant in each gene, viz., MTUS2 (rs928661) c.2855A > C; p. Gln952Pro; and another variant in PKD1L2 (rs11150370) c.901C > G; p. Pro301Ala; were found in all the 4 cases respectively. Six missense variants (rs16841277, rs2228209, rs2230283, rs3741883, rs1716, rs2272606) were found in 3 patients. Seven missense variants (rs6558702, rs920829, rs149000560, rs2279574, rs11109968, rs9935113, rs7191351) were found in two of the four patients.
Four different SNPs were identified in our study for the NCK Associated Protein 5 (NCKAP5) gene on ch:2; the diseases associated with NCKAP5 include neurological disorder, hypersomnia which affiliates many tissues including brain, causing behaviour and neurological disorders . The National Sleep Foundation states that ADHD is linked with a variety of sleep problems. For example, children with ADHD had increased rates of daytime sleepiness, had shorter sleep durations, more sleep problems on weekdays, anxiety or depression aggravated these sleep complications in the ADHD children  . Another study found that children with ADHD had signs of sleep-disordered breathing such as habitual snoring, habitual snoring may be three times higher in children with ADHD (Sleep and neurobehavioral characteristics of 5- to 7-year-old children with parentally reported symptoms of ADHD .
The ITGAE gene or integrin subunit alpha E belongs to the Integrin alpha subunits (ITGA) family. It functions as a receptor for E-cadherin. This protein mediates adhesion of intra-epithelial T-lymphocytes to epithelial cell monolayers and may be involved in the activation of intestinal intraepithelial lymphocytes (IEL) . It is expressed in multiple tissues including IEL, bone marrow and brain. ITGA11 (integrin subunit alpha 11) is another member of the Integrin alpha subunits (ITGA) family. It is a collagen receptor that is expressed in muscle tissues, and it plays an important role in attaching muscle tissue to the extracellular matrix . Both ITGAE and ITGA11 are involved in cell adhesion. SNP rs220470 in the ITGAE gene and SNP rs7164335 in ITGA11 were present in the 30 top hit single SNPs in a study conducted using independent DNA pools from 343 ADHD-affected adults and 304 controls for association analyses of the ADHD diagnostic phenotype .
However, in this study we identified two SNPs in ITGAE (integrin, alpha-E) and one SNP in ITGA11 (integrin, alpha-11) that may have a role in causing ADHD, which are: rs1716, and rs2272606 in ITGAE gene and one SNP rs2306022 in ITGA11.However, these three SNPs are not implicated in ADHD development previously. The SNP rs1716 was previously associated with a gastrointestinal stromal tumour subtype  and increased risk of melanoma along with SNP rs2272606 , while SNP rs2306022 of ITGA11 was implicated in increasing risk for lung cancer and fibrosis . In three of the patients we have detected a SNP rs3741883 located on the NAUK1 (AMPK-related protein kinase 5; ARK5 or KIAA0537) gene. The NUAK family, SNF1-like kinase gene is highly expressed in heart and brain and involved in several processes that include tumour progression, cell proliferation, cell adhesion, senescence and regulation of cell ploidy . Two SNPs located on NUAK1 (rs6539247 and rs2569973) have been detected in a recent study done in ADHD adult patients, however, in the present study in children with ADHD these two SNPs were not found , evidence suggests that the NAUK1 gene was associated with cortical amyloid-β protein which plays an essential role in Alzheimer’s disease , it might also be involved in attention performance in persistent ADHD since Alzheimer’s disease is relevant in the context of ADHD.
Additionally, previous studies have identified SNP rs753390486 in KIF21Bgene (kinesin family member 21b; KIF21B or KIAA0449), which was also found in the present study in one ADHD patient. This gene is expressed in high level in the CNS, mostly in the dendrites of neurons . The KIF21B gene is a member of the kinesin superfamily that contains 45 genes known as kinesins, divided into 15 different families. Genetic variation in some members of the kinesin family has been linked to neurodegenerative diseases such as Alzheimer’s disease, amyotrophic lateral sclerosis (ALS) and Huntington’s disease . The KIF21B gene was previously associated with multiple sclerosis . Our results also found five SNPs located on the PKD1L2 gene which are: rs9935113, rs34276551, rs11150370, rs7185774 and rs7191351, this gene belongs to the polycystin family and is expressed in tissues including placenta, mammary gland, liver, lung, skeletal muscle, brain and testis  . PKD1L2 may function as an ion-channel regulator and as a G-protein-coupled receptor  .
Genetic studies have established that several genes have been associated with ADHD such as ACOXL, PADR3B, HIVEP1, TRPA1 and SYNE1. Among these genes, SYNE1 gene contributes to multiple genetic disorders, specifically neuromuscular disorders  . The SYNE1 gene is mainly expressed in skeletal muscle cells, therefore, mutations in SYNE1 could alter muscular functions . Mutations in SYNE1 have been linked with numerous cases of autism, bipolar disorder and ARCA1 [42 , 44]. It has been identified as one of the top risk genes in causing ADHD according to GWAS studies . Our findings from this study revealed a SNP (rs150179494) in the SYNE1 gene. These results suggest the association between these SNPs and ADHD occurrence. Also, there are no current studies on the effects of the other SNPs in the SYNE1 gene on ADHD development. Furthermore, several studies have linked the ACOXL gene to certain types of tumours. ACOXL showed no levels of gene expression in almost 72% of prostate tumours . The exact function of the ACOXL gene is largely undefined; however, loci in this gene are thought to participate in causing alopecia areata . Our study has suggested the contribution of four SNPs (rs779762950, rs77331476, rs1554005, rs17041850) in the ACOXL gene that could cause ADHD. Further studies are required to explore the role of the ACOXL gene and the effects of these SNPs in ADHD and other mental disorders.
Another gene, which is thought to be related to ADHD, is the PARD3 (partitioning-defective protein 3) gene codes for a protein that is involved in exon splicing and neuronal polarity . Among the PARD gene family, the PARD3B gene, located on the 2q33.33 region , has been linked with nicotine dependence defined by the FTND (Fagerström Test for Nicotine Dependence)  . Our findings revealed a novel SNP (rs80119103) to be responsible for increasing the susceptibility of ADHD. A recent study reported an association between PARD3 and schizophrenia and ADHD , and it has been observed in adults with various forms of disinhibitory psychopathology and antisocial behaviour as well . TRPA1 is a member of the transit receptor potential family (TRPA1), and it was suggested to have a role in several sensing reactions. A study that was performed on knockout mice has shown a relationship between TRPA1 and the transduction of chemical and mechanical stimuli through nociceptor sensory neurons  . According to a study that has investigated the role of TRPA1, it is suggested that this gene has an essential role in inflammation and neuropathic pain . Our study identified one of the SNPs in TRPA1 (rs920829) as a risk factor for ADHD. This SNP has been previously presented in another study as a key player in causing asthma in children . The HIVEP1 gene encodes for one of the transcription factors that belong to the ZAS family. It binds to motifs that regulate HIV gene expression . A study has shown that the HIVEP1 gene was one of the highest ranked-genes to be related to mental disorders such as schizophrenia , bipolar disease (BP) and post-traumatic stress disorder (PTSD)  . The findings from our study have revealed an SNP (rs2228209) in the HIVEP1 gene that could potentially be linked to ADHD.A body of evidence has shown that SNPs in the integrin alpha 1 ITGA1 gene may cause ADHD. A study conducted on 1033 ADHD patients, (age ranging from 6 to 16 years). Four SNPs have been reported in the ITGA1 gene located on chromosome 5 (rs1979398, rs16880453, rs1531545 and rs4074793  . ITGA1 gene has a role in cell attachment and neurite outgrowth as well as in peripheral nerve regeneration  however, the exact mechanism ofITGA1 polymorphisms is still not clearly understood in ADHD. Consistent with this study, our data explored the association of the ITGA1 gene at chromosome 5 in position 5:2243247 with ADHD. It has been reported that in ITGAE a SNP rs220470 (intronic variant) was associated with ADHD in Caucasians . To our knowledge, this finding is novel and these SNPs have not been reported in the literature. Two Integrin subunits including alpha and beta have been known to have a role in the development of the neurological system by mediating many neural cell functions .
In addition, our study has identified SNP on the Fermitin family homolog 3FERMT3 gene on chromosome 11 (rs149000560), another study done in 465 trioshas identified different SNP on the same gene (rs12575642) . This missense variant (rs149000560) is responsible for the leukocyte adhesion deﬁciency-1/variant syndrome . Taken together, the data suggested that these SNPs might have an association with ADHD susceptibility. On the other hand, our study has identified three SNPs in the NCK Associated Protein 5; NCKAP5 gene (rs776162189, rs16841277, and rs142329411) which may be related to the ADHD phenotype. Indeed, previous study has demonstrated a set of risk genes or copy number variations (CNVs) involved in the autism spectrum, disorder as in human subjects, one of which was the NCKAP5 gene . Our exome sequencing has shown SNPs in the POC1B-GALNT4 and FAM71C-ANKS1Bgenes, which may have a role in causing ADHD, these SNPs were rs2230283 and rs11109968, respectively. However, to our knowledge the significant association of these two SNPs with ADHD has not been previously reported.
We have identified two novel missense mutations in the ITGA1 and SPATA13 genes respectively in two different ADHD children with exome sequencing analysis. These two novel mutations, c.988G > A, p. Ala330Thr, in SPATA13, and c.3451G > T, p. Ala1151Ser, in ITGA1, have not been reported so far in literature. Our current results revealed mutations which had not previously been linked with ADHD within the database. However, the three documented variants are rs13166360 with bipolar disorder, rs920829 with neuropathic pain, and rs6558702 with schizophrenia. This study recommends screening these in ADHD cases to confirm the prevalence in the Saudi population. It is recommended for future studies to examine the four variants in detail. The strength of our current study is that we have selected satisfactory diagnosis severity ADHD cases in the same region. However, conducting of further studies with larger number of patients is necessary to confirm the ADHD associated SNPs in Saudi Children.
This work was supported financially by the Deanship of Scientific Researchat Umm Al-Qura University (Grant Code: 15-MED-3-1-0062) to Dr. Neda Bogari. We are thankful to Mr. Ehab M. Melibary and Mr. Sami M. Kalantan for their technical in this study, Miss Renad Y. Natto, Miss Dema N. Bogari, and Mr Safar A. Almalki for the continuous support. We all are deeply appreciative for Mr. Abdulmonim Gowda, of Integrated Gulf Biosystems. We would like to acknowledge the Albaseera Center and Successful Communication Center, Khawla Bint Khuwailad Center and Center of Optimism.
Source of Fund
This work was supported financially by the Deanship of Scientific Researchat Umm Al-Qura University (Grant Code: 15-MED-3-1-0062) to Dr. Neda Bogari
We confirm that the manuscript has been read and approved by all named authors and that there are no other persons who satisfied the criteria for authorship but are not listed. We further confirm that the order of authors listed in the manuscript has been approved by all of us.
We confirm that we have given due consideration to the protection of intellectual property associated with this work and that there are no impediments to publication, including the timing of publication, with respect to intellectual property. In so doing we confirm that we have followed the regulations of our institutions concerning intellectual property.