In-silico genomic landscape characterization and evolution of SARS-CoV-2 variants isolated in India shows significant drift with high frequency of mutations

In-silico studies on SARS-CoV-2 genome are considered important to identify the significant pattern of variations and its possible effects on the structural and functional characteristics of the virus. The current study determined such genetic variations and their possible impact among SARS-CoV-2 variants isolated in India. A total of 546 SARS-CoV-2 genomic sequences (India) were retrieved from the gene bank (NCBI) and subjected to alignment against the Wuhan variant (NC_045512.2), the corresponding amino acid changes were analyzed using NCBI Protein-BLAST. These 546 variants revealed 841 mutations; most of these were non-synonymous 464/841 (55.1%), there was no identical variant compared to the original strain. All genes; coding and non-coding showed nucleotide changes, most of the structural genes showed frequent nonsynonymous mutations. The most affected genes were ORF1a/b followed by the S gene which showed 515/841 (61.2%) and 120/841 (14.3%) mutations, respectively. The most frequent non-synonymous mutation 486/546 (89.01%) occurred in the S gene (structural gene) at position 23,403 where A changed to G leading to the replacement of aspartic acid by glycine in position (D614G). Interestingly, four variants also showed deletion. The variants MT800923 and MT800925 showed 12 consecutive nucleotide deletion in position 21982–21993 resulting in 4 consecutive amino acid deletions that were leucine, glycine, valine, and tyrosine in positions 141, 142, 143, and 144 respectively. The present study exhibited a higher mutations rate per variant compared to other studies carried out in India.     

Coronavirus and co-infections: A Saudi Arabian perspective

Mortality due to infectious diseases continues to rise globally, despite advances in antimicrobial therapy and supportive care. This is evident with the occurrence of coronavirus disease 2019 (COVID-19) pandemic, instigated by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Saudi Arabia, an eminent country within the Arab region, has had significant impact during global pandemics, concomitant with the fact that millions of Muslims travel to Saudi Arabia for pilgrimages every year. Herein, we discuss the significance of SARS-CoV-1, SARS-CoV-2, as well as the Middle East respiratory syndrome coronavirus (MERS-CoV) in Saudi Arabia with particular reference to global transmission and/or emergence of new variants due to genetic mixing of different strains. Furthermore, we also discuss the role of Saudi Arabia with reference to novel emerging infectious diseases and re-emerging infections, such as Ebola, zika, and monkeypox, as well as in the context on coinfections. Future strategies to limit the spread of viral infections and the pivotal role of Saudi Arabia, are deliberated upon.     

Computational drug screening against the SARS-CoV-2 Saudi Arabia isolates through a multiple-sequence alignment approach

COVID-19 is a rapidly emerging infectious disease caused by the SARS-CoV-2 virus currently spreading throughout the world. To date, there are no specific drugs formulated for it, and researchers around the globe are racing against the clock to investigate potential drug candidates. The repurposing of existing drugs in the market represents an effective and economical strategy commonly utilized in such investigations. In this study, we used a multiple-sequence alignment approach for preliminary screening of commercially-available drugs on SARS-CoV sequences from the Kingdom of Saudi Arabia (KSA) isolates. The viral genomic sequences from KSA isolates were obtained from GISAID, an open access repository housing a wide variety of epidemic and pandemic virus data. A phylogenetic analysis of the present 164 sequences from the KSA provinces was carried out using the MEGA X software, which displayed high similarity (around 98%). The sequence was then analyzed using the VIGOR4 genome annotator to construct its genomic structure. Screening of existing drugs was carried out by mining data based on viral gene expressions from the ZINC database. A total of 73 hits were generated. The viral target orthologs were mapped to the SARS-CoV-2 KSA isolate sequence by multiple sequence alignment using CLUSTAL OMEGA, and a list of 29 orthologs with purchasable drug information was generated. The results showed that the SARS CoV replicase polyprotein 1a had the highest sequence similarity at 79.91%. Through ZINC data mining, tanshinones were found to have high binding affinities to this target. These compounds could be ideal candidates for SARS-CoV-2. Other matches ranged between 27 and 52%. The results of this study would serve as a significant endeavor towards drug discovery that would increase our chances of finding an effective treatment or prevention against COVID19.     

Prevalence of H1299R polymorphism in the Factor V gene among the Taif-Saudi Arabia population using polymerase chain reaction-reverse hybridization technique

Cardiovascular disease (CVD) remains a major health hazard worldwide. Single nucleotide polymorphisms (SNPs) represent a part of risk factors that contributes to cardiovascular disease. SNP in the coagulation factor V genes have been shown to play a role in the development of cardiovascular disease. Coagulation Factor V is an enzyme cofactor of the coagulation system and contributes to a normal haemostatic balance. The His1299Arg polymorphism in the Factor V gene has been identified and linked to hereditary thrombophilia. The aim of the present study is to determine the prevalence of HR2 haplotype and allele frequency of His1299Arg polymorphism in the Factor V gene among randomly selected healthy individuals from Taif population which belonging to western region of Saudi Arabia. Genotyping of this SNP was carried out via CVD StripAssay, which based on a polymerase chain reaction-reverse hybridization technique. Two hundred healthy unrelated individuals were included in the study. Seventeen out of the studied population (8.5%) had the HR2 haplotype; 14 (7%) were heterozygous (R1/R2), and three (1.5%) were homozygous (R2/R2), with an allelic frequency of 0.05. This is the first report for a Saudi Arabian population that estimates the prevalence of HR2 haplotype and its allele frequencies. In conclusion, the His1299Arg mutant was noticeable within population of western Saudi Arabia. Further larger studies are needed to (1) estimate the prevalence of this mutant among individuals belonging to different KSA locations (2) assess the relative contribution of this mutational event separately and in combination with other thrombophilic polymorphisms in the etiology of cardiovascular disease in KSA.     

Sequence analysis of sub-genotype D hepatitis B surface antigens isolated from Jeddah,Saudi Arabia

Little is known about the prevalence of HBV genotypes/sub-genotypes in Jeddah province, although the hepatitis B virus (HBV) was identified as the most predominant type of hepatitis in Saudi Arabia. To characterize HBV genotypes/sub-genotypes, serum samples from 15 patients with chronic HBV were collected and subjected to HBsAg gene amplification and sequence analysis. Phylogenetic analysis of the HBsAg gene sequences revealed that 11 (48%) isolates belonged to HBV/D while 4 (18%) were associated with HBV/C. Notably, a HBV/D sub-genotype phylogenetic tree identified that eight current isolates (72%) belonged to HBV/D1, whereas three isolates (28%) appeared to be more closely related to HBV/D5, although they formed a novel cluster supported by a branch with 99% bootstrap value. Isolates belonging to D1 were grouped in one branch and seemed to be more closely related to various strains isolated from different countries. For further determination of whether the three current isolates belonged to HBV/D5 or represented a novel sub-genotype, HBV/DA, whole HBV genome sequences would be required. In the present study, we verified that HBV/D1 is the most prevalent HBV sub-genotype in Jeddah, and identified novel variant mutations suggesting that an additional sub-genotype designated HBV/DA should be proposed. Overall, the results of the present HBsAg sequence analyses provide us with insights regarding the nucleotide differences between the present HBsAg/D isolates identified in the populace of Jeddah, Saudi Arabia and those previously isolated worldwide. Additional studies with large numbers of subjects in other areas might lead to the discovery of the specific HBV strain genotypes or even additional new sub-genotypes that are circulating in Saudi Arabia.     

Sequence complementarity between human noncoding RNAs and SARS-CoV-2 genes: What are the implications for human health?

ObjectivesTo investigate in silico the presence of nucleotide sequence complementarity between the RNA genome of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) and human non-coding (nc)RNA genes.MethodsThe FASTA sequence (NC_045512.2) of each of the 11 SARS-CoV-2 isolate Wuhan-Hu-1 genes was retrieved from NCBI.nlm.nih.gov/gene and the Ensembl.org library interrogated for any base-pair match with human ncRNA genes. SARS-CoV-2 gene-matched human ncRNAs were screened for functional activity using bioinformatic analysis. Finally, associations between identified ncRNAs and human diseases were searched in GWAS databases.ResultsA total of 252 matches were found between the nucleotide sequence of SARS-CoV-2 genes and human ncRNAs. With the exception of two small nuclear RNAs, all of them were long non-coding (lnc)RNAs expressed mainly in testis and central nervous system under physiological conditions. The percentage of alignment ranged from 91.30% to 100% with a mean nucleotide alignment length of 17.5 ± 2.4. Thirty-three (13.09%) of them contained predicted R-loop forming sequences, but none of these intersected the complementary sequences of SARS-CoV-2. However, in 31 cases matches fell on ncRNA regulatory sites, whose adjacent coding genes are mostly involved in cancer, immunological and neurological pathways. Similarly, several polymorphic variants of detected non-coding genes have been associated with neuropsychiatric and proliferative disorders.ConclusionThis pivotal in silico study shows that SARS-CoV-2 genes have Watson-Crick nucleotide complementarity to human ncRNA sequences, potentially disrupting ncRNA epigenetic control of target genes. It remains to be elucidated whether this could result in the development of human disease in the long term.     

Molecular characterization of endangered endemic plant Aloe pseudorubroviolacea using chloroplast matK and plastid rbcL gene

An endangered and rare species Aloe pseudorubroviolacea from the plant family Asphodelaceae which is presently recorded as endangered in Saudi Arabia collected from Al-Baha region of Saudi Arabia its GPS Latitude and Longitude coordinates 19.8345, 41.5481. The chloroplast matK and rbcL gene was considered in this study based on molecular identification the size is about 571 and 664 bp respectively. From the sequence analysis the gene matK and rbcL confirm that this species is very much closely related with A. rubroviolacea and also inter related with the species Astroloba rubriflora, Chrysopogon gryllus, Chortolirion angolense shows about 98.7% sequence homology. The partial matK and rbcL gene sequence discriminate Aloe pseudorubroviolacea from the closely related plant species, A. rubroviolacea. The gene sequence of rbcL discriminates the species from Chrysopogon gryllus and Chortolirion angolense, demonstrates the nucleotide variations in 3 different sites (623C/T; 653C/T; 700C/A). This study showed that matK and rbcL sequence region of chloroplast gene used to authenticate the samples of A. pseudorubroviolacea and which provide to help in correct identification and conservation process of this medicinally valuable endangered plant species.     

Polycystic ovary syndrome is linked with the fat mass obesity (FTO) gene variants rs17817449 and rs1421085 in western Saudi Arabia

Polycystic ovary syndrome (PCOS) is characterised by infertility, obesity, insulin resistance and clinical and/or biochemical signs of hyperandrogenism. Obesity is known to be correlated with PCOS causing ovulatory dysfunction and hormone imbalances. Moreover, fat mass and the obesity gene (FTO) were linked with obesity and PCOS. Therefore, it is of interest to determine the genotype and allele frequency for three FTO variants - rs17817449 (G/T), rs1421085 (C/T) and rs8050136 (A/C) -in western Saudi population. 95 PCOS patients and 94 controls were recruited for this study. The genetic variants were assayed using real-time polymerase chain reaction using TaqMan genotyping assays. The chi-squared test was applied to investigate the difference between single nucleotide polymorphisms on PCOS and control subjects, and binary logistic regression was used to determine the association of FTO variants with PCOS symptoms. Variants rs17817449 and rs1421085 were significantly linked with PCOS susceptibility in the study population. Rs17817449 and rs8050136 were significantly associated with hair loss in the PCOS group. Furthermore, rs1421085 and rs8050136 were associated with a high body mass index (BMI>30 kg/m2). Risk alleles in our population associated with hair loss and elevated BMI in women with PCOS were homozygous C for rs8050136. This data will help in defining the genetic predisposition of PCOS among women in western Saudi Arabia.     

The seroprevalence of SARS-CoV-2 IgG antibodies among asymptomatic blood donors in Saudi Arabia

BackgroundIn late 2019, cases of severe pneumonia with unidentified etiology began to emerge in Wuhan, China, before progressively spreading first nationally and then globally.The current study sought to investigate the seroprevalence of immunoglobulin G (IgG) antibodies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) among blood donors in Al-Madinah, Saudi Arabia. To our knowledge, this is the first study in Saudi Arabia to screen blood donors who were not known to be previously infected with SARS-CoV-2.MethodsThis study was a cross-sectional study to assess individuals who donated blood to the central blood bank in Al-Madinah between mid-May and mid-July 2020. An enzyme-linked immunosorbent assay (ELISA) was designed and established to detect antibodies directed against the SARS-CoV-2 spike protein in serum samples. A total of 1,212 healthy blood donors participated in this study. The donors were males and met the requirements for blood donation during the COVID-19 pandemic period in Saudi Arabia.ResultsThe SARS-CoV-2 seroprevalence among blood donors in Al-Madinah was 19.31% (n = 234/1212; 95% confidence interval: 17.12%–21.64%). No statistically significant difference was identified in seropositivity according to age. However, significant differences (p < 0.001) were identified according to ABO blood groups, with those with type A blood presenting the highest rate of seropositivity (29.18%) compared with the other blood groups (12.65% for type B, 16.36% for type AB, and 15.11% for type O).ConclusionA high prevalence of SARS-CoV-2 antibodies was detected among blood donors in Al-Madinah, which indicated a high level of exposure to the virus within the population. This further suggested that as high as one-fifth of the population may have acquired innate immunity against the virus.     

Molecular adaptive evolution of SARS-COV-2 spike protein in Saudi Arabia

The sequences of SARS-CoV-2 spike (S) from Saudi Arabia along with SARS-CoV and bat SARS-like CoVs were obtained. Positive selection analysis and secondary structure investigation of spike sequences were performed. Adaptive molecular evolution was observed in SARS-CoV-2 displayed by positive selection pressure at N-terminal domain (NTD; codons 41, 163, 174 and 218), Receptor binding domain (RBD; codons 378 and 404) and S1/S2 Cleavage site (codon 690). Furthermore, the spike protein secondary structure depicted by the homo-trimer structure showed a high similarity between Saudi SARS-CoV-2 isolate and the parental strain (bat SL-COVZC45). Despite the high similarity depicted in the spike sequence model alignment, it displayed a significant difference when each chain was treated solely owing to 7 motif differences in the three composing chains. In addition, SARS-CoV-2 S trimer model uncovered the presence of N-acetyl glucosamine ligands. Eventually, 3C-like proteinase cleavage site was observed in S2 domain could be used as a site for drug discovery. Genetics and molecular evolutionary facts are useful for assessment of evolution, host adaptation and epidemic patterns ultimately helpful for adaptation of control strategies.     

Systematic screening for mutations in the human serotonin-2A (5-HT2A) receptor gene: Identification of two naturally occurring receptor variants and association analysis in schizophrenia

A statistically significant association between a silent mutation (102T/C) in the serotonin-2A (5-HT_2A) receptor gene and schizophrenia has recently been reported in a sample of Japanese patients and healthy controls. This finding suggests that genetic predisposition to schizophrenia may be affected by a functional 5-HT_2A receptor variant that is in linkage disequilibrium with 102T/C. In the present study, we have sought to identify genetic variation in the 5-HT_2A receptor gene by screening genomic DNA samples from 91 unrelated subjects comprising 45 patients with schizophrenia and 46 healthy controls by using single-strand conformation analysis. We have identified four nucleotide sequence variants. Two sequence changes would result in protein alterations: a substitution of threonine by asparagine at position 25 (Thr25Asn), and a substitution of histidine by tyrosine at position 452 (His452Tyr). In order to test for a possible contribution to the development of schizophrenia, we have determined allele frequencies in extended samples of unrelated patients and healthy controls. The two amino acid substitutions are found with similar frequencies in patients and controls, indicating that the presence of these variants is not causally related to the development of schizophrenia. However, the reported association of the non-coding polymorphism 102T/C with the disease has also been detected in our sample (P = 0.041, odds ratio = 1.28, 95% confidence interval 1.012–1.623).     

Virological surveillance,molecular phylogeny,and evolutionary dynamics of hepatitis C virus subtypes 1a and 4a isolates in patients from Saudi Arabia

Hepatitis C virus (HCV) subtypes are pre-requisite to predict endemicity, epidemiology, clinical pathogenesis, diagnosis, and treatment of chronic hepatitis C infection. HCV genotypes 4 and 1 are the most prevalent in Saudi Arabia, however; less consensus data exist on circulating HCV subtypes in infected individuals. This study was aimed to demonstrate the virological surveillance, phylogenetic analysis, and evolutionary relationship of HCV genotypes 4 and 1 subtypes in the Saudi population with the rest of the world. Fifty-five clinical specimens from different parts of the country were analyzed based on 5′ untranslated region (5′ UTR) amplification, direct sequencing, and for molecular evolutionary genetic analysis. Pair-wise comparison and multiple sequence alignment were performed to determine the nucleotide conservation, nucleotide variation, and positional mutations within the sequenced isolates. The evolutionary relationship of sequenced HCV isolates with referenced HCV strains from the rest of the world was established by computing pairwise genetic distances and generating phylogenetic trees. Twelve new sequences were submitted to GenBank, NCBI database. The results revealed that HCV subtype 4a is more prevalent preceded by 1a in the Saudi population. Molecular phylogeny predicts the descendants’ relationship of subtype 4a isolates very close to Egyptian prototype HCV strains, while 1a isolates were homogeneous and clustering to the European and North American genetic lineages. The implications of this study highlight the importance of HCV subtyping as an indispensable tool to monitor the distribution of viral strains, to determine the risk factors of infection prevalence, and to investigate clinical differences of treatment outcomes among intergenotypic and intragenotypic isolates in the treated population.     

Molecular characterization of the Fe‑hydrogenase gene marker in Trichomonas gallinae isolated from birds in Riyadh,Saudi Arabia

Trichomonas gallinae causes avian oropharyngeal trichomonosis. This pathogen affects a large number of bird species and may cause substantial economic losses to the poultry industry. Al-Azizia poultry market in Riyadh, Saudi Arabia is among the largest poultry markets in the Arabian Gulf. Birds traded in this market may be exposed to a variety of T. gallinae strains. Genetic diversity of T. gallinae among birds in the market was examined using Fe‑hydrogenase gene sequences. These sequences were amplified by PCR for twenty-nine isolates of T. gallinae from four different avian species, including 21 feral pigeons, one common mynah, three chickens, and four turkeys. Sequence analysis showed ten variant gene sequences. Nine sequences comprise a new subtype, including A(KSAF1), C(KSAF1) and C(KSAF3) with 34.48% (n = 10), 6.90% (n = 2), 6.90% (n = 2) of the isolates, respectively. Analyses also showed an additional five new sequences (KSAF1.1., KSAF2, KSAF13, KSAF14, KSAF15), representing 17.24% of the isolates. Subtype II (KSAF) was found in four feral pigeons (13.80%). To our knowledge, this report is the first to describe genotypes of T. gallinae from pigeons in Saudi Arabia using Fe‑hydrogenase gene sequences for subtyping. Subtype analysis infers the presence of multiple genotypes of T. gallinae in Saudi avian populations.     

Evaluation of RDS/Peripherin and ROM1 as candidate genes in generalised progressive retinal atrophy and exclusion of digenic inheritance

Generalised progressive retinal atrophy (gPRA) is a heterogeneous group of hereditary diseases causing degeneration of the retina in dogs and cats. As a combination of mutations in the RDS/Peripherin and the ROM1 genes leads to the phenotype of retinitis pigmentosa in man we first performed mutation analysis to screen these genes for disease causing mutations followed by the investigation of a digenic inheritance in dogs. We cloned the RDS/Peripherin gene and investigated the RDS/Peripherin and ROM1 genes for disease causing mutations in 13 gPRA-affected dog breeds including healthy animals, obligate gPRA carriers and gPRA-affected dogs. We screened for mutations using single strand conformation polymorphism (SSCP) analysis. Sequence analysis revealed several sequence variations. In the coding region of the RDS/Peripherin gene three nucleotide exchanges were identified (A277C; C316T; G1255A), one of which leads to an amino acid substitution (Ala339Thr). Various silent sequence variations were found in the coding region of the ROM1 gene (A536G, G1006A, T1018C, T1111C, C1150T, C1195T), as well as an amino acid substitution (G252T; Ala54Ser). By excluding the respective gene as a cause for gPRA several sequence variations in the intronic regions were investigated. None of these sequence variations cosegregated with autosomal recessively (ar) transmitted gPRA in 11 breeds. The candidate gene RDS/Peripherin obviously does not harbour the critical mutation causing the autosomal recessive form of gPRA because diseased individuals show heterozygous genotypes for sequence variations in the Miniature Poodle, Dachshund, Australian Cattle Dog, Cocker Spaniel, Chesapeake Bay Retriever, Entlebucher Sennenhund, Sloughi, Yorkshire Terrier, Tibet Mastiff, Tibet Terrier and Labrador Retriever breeds. In the following breeds the ROM1 gene was also excluded indirectly for gPRA: Miniature Poodle, Dachshund, Australian Cattle Dog, Sloughi, Collie, Tibet Terrier, Labrador Retriever and Saarloos/Wolfhound. Digenic inheritance for gPRA is practically excluded for both these genes in four breeds: Miniature Poodle, Dachshund, Labrador Retriever and Saarloos/Wolfhound.     

Characterization of Hepatitis C Virus Genotypes by Direct Sequencing of HCV 5′UTR Region of Isolates from Saudi Arabia

The current study was designed to determine the Hepatitis C Virus (HCV) genotypes in a representative sample of HCV chronically infected patients in Saudi Arabia. All HCV isolates were genotyped by sequencing of the 5′UTR region and newly identified HCV isolates were identified. Specific universal primers targeting 5′UTR region were used for both amplification and sequencing of all isolates that resulted in 244 bp fragment which represent about 80% of 5′UTR region. Most of HCV isolates in this study were genotype 4 (76.4%) where only few isolates were recognized as genotype 1 (19.6%). All results were compared to HCV reference sequences from LOS ALAMOS HCV database, considering only the complete full genomes for the main phylogenetic analysis. Sequences that showed maximum identity (98% –100%) were selected. Most isolates were identical with HCV genotype 4 references. Some isolates were similar to different subtypes of HCV genotypes 4, 1 and 6. Phylogenetic analysis showed resemblance of most isolates to similar ones from the Far East, North America and Egypt. Using sequence Weblogo, Alignment analysis of isolated HCV genotypes 4 and 1 showed 92% and 95.5% nucleotide conservation, respectively. There was no predominant nucleotide in the varied sites, in both genotypes. All isolated sequences were submitted to GenBank database.     

High prevalence of SLC6A8 deficiency in X-linked mental retardation

A novel X-linked mental retardation (XLMR) syndrome was recently identified, resulting from creatine deficiency in the brain caused by mutations in the creatine transporter gene, SLC6A8. We have studied the prevalence of SLC6A8 mutations in a panel of 290 patients with nonsyndromic XLMR archived by the European XLMR Consortium. The full-length open reading frame and splice sites of the SLC6A8 gene were investigated by DNA sequence analysis. Six pathogenic mutations, of which five were novel, were identified in a total of 288 patients with XLMR, showing a prevalence of at least 2.1% (6/288). The novel pathogenic mutations are a nonsense mutation (p.Y317X) and four missense mutations. Three missense mutations (p.G87R, p.P390L, and p.P554L) were concluded to be pathogenic on the basis of conservation, segregation, chemical properties of the residues involved, as well as the absence of these and any other missense mutation in 276 controls. For the p.C337W mutation, additional material was available to biochemically prove (i.e., by increased urinary creatine : creatinine ratio) pathogenicity. In addition, we found nine novel polymorphisms (IVS1+26G-->A, IVS7+37G-->A, IVS7+87A-->G, IVS7-35G-->A, IVS12-3C-->T, IVS2+88G-->C, IVS9-36G-->A, IVS12-82G-->C, and p.Y498) that were present in the XLMR panel and/or in the control panel. Two missense variants (p.V629I and p.M560V) that were not highly conserved and were not associated with increased creatine : creatinine ratio, one translational silent variant (p.L472), and 10 intervening sequence variants or untranslated region variants (IVS6+9C-->T, IVS7-151_152delGA, IVS7-99C-->A, IVS8-35G-->A, IVS8+28C-->T, IVS10-18C-->T, IVS11+21G-->A, IVS12+15C-->T, *207G-->C, IVS12+32C-->A) were found only in the XLMR panel but should be considered as unclassified variants or as a polymorphism (p.M560V). Our data indicate that the frequency of SLC6A8 mutations in the XLMR population is close to that of CGG expansions in FMR1, the gene responsible for fragile-X syndrome.     

An alpha-1 antitrypsin genetic variant from India

The highly polymorphic human alpha-1 antitrypsin (AAT) gene codes for the most abundant circulating plasma serine protease inhibitor. Previously, genetic variants of the AAT gene were reported from different regions of the world. In the present study, the AAT gene was characterized in an Indian sample. The AAT gene was isolated and cloned from a liver biopsy sample through RT-PCR and the full-length gene was sequenced. Nucleotide sequence comparison with the human genome and the AAT sequences available in the GenBank (NCBI) demonstrated four unique variations--(i) an A to G variation at position 286 (Thr96Ala), (ii) an A to G variation at position 839 (Asp280Gly), (iii) a T to C variation at position 1182 that did not result in any change in the protein sequence (TTT to TTC both code for Phe) and (iv) an A to C variation at position 1200 (Glu400Asp) that resulted in replacement by an amino acid of similar nature. Other variations found were T to C at position 710 (Val273Ala) and T to C position 863 (Val288Glu), which were also reported earlier. In conclusion, this study reports the entire 1257 bp nucleotide sequence of protein coding region of the human AAT gene from an Indian sample. This preliminary finding is significant, as it reports for the first time the AAT gene sequence in the Indian sample.     

Species-specific molecular signature of Commiphora species of Saudi Arabia inferred from internal transcribed spacer sequences of nuclear ribosomal DNA

The deciduous habit and tendency to produce flowers prior to developing leaves, and a predominantly dioecious system of breeding in the genus Commiphora leads to difficulties in its taxonomic identification at species level. The characteristics of easy amplification by universal primer, shorter length and higher discrimination power at the species level makes the internal transcribed spacer (ITS) sequence of nuclear ribosomal DNA (nrDNA) to a smart gene for generating species-specific phylogenetic inferences in most of the plants groups. The present study deals the ITS sequence of nrDNA based molecular genotyping of seven species of the genus Commiphora of Saudi Arabia. The molecular phylogenetic analysis of ITS sequences of nrDNA of Commiphora species distributed in Saudi Arabia reveals the the occurrence of C. madagascariens in Saudi Arabia.     

SAIVGeM: spreadsheet analysis of immunoglobulin VH gene mutations

The analysis of mutations in immunoglobulin heavy chain variable (IGHV) region genes is a tedious process when performed by hand on multiple sequences. This report describes a set of linked Microsoft Excel files that perform several common analyses on large numbers of IGHV sequences. The spreadsheet analysis of immunoglobulin VH gene mutations (SAIVGeM) package determines the distribution of mutations among each nucleotide, the nature of the mutation at both the nucleotide and amino acid level, the frequency of mutation in the A/G G C/T A/T (RGYW) hotspot motifs of both strand polarity, and the distribution of replacement and silent mutations among the complementarity determining regions (CDRs) and the framework regions (FRs) of the immunoglobulin gene as defined by either the Kabat or IMGT conventions. These parameters are summarized and graphically presented where appropriate. In addition, the SAIVGeM package analyzes those mutations that occur in third positions of redundant codons. Because any nucleotide change in these positions is inherently silent, these positions can be used to study the mutational spectra without biases from the selection of protein structure.