Integrative analysis of genomic and transcriptomic data of normal, tumour, and co-occurring leukoplakia tissue triads drawn from patients with gingivobuccal oral cancer identifies signatures of tumour initiation and progression.

A thickened, white patch - leukoplakia - in the oral cavity is usually benign, but sometimes (in ~9% of individuals) it progresses to malignant tumour. Because the genomic basis of this progression is poorly understood, we undertook this study and collected samples of four tissues - leukoplakia, tumour, adjacent normal, and blood - from each of 28 patients suffering from gingivobuccal oral cancer. We performed multiomics analysis of the 112 collected tissues (four tissues per patient from 28 patients) and integrated information on progressive changes in the mutational and transcriptional profiles of each patient to create this genomic narrative. Additionally, we generated and analysed whole-exome sequence data from leukoplakia tissues collected from 11 individuals not suffering from oral cancer. Nonsynonymous somatic mutations in the CASP8 gene were identified as the likely events to initiate malignant transformation, since these were frequently shared between tumour and co-occurring leukoplakia. CASP8 alterations were also shown to enhance expressions of genes that favour lateral spread of mutant cells. During malignant transformation, additional pathogenic mutations are acquired in key genes (TP53, NOTCH1, HRAS) (41% of patients); chromosomal-instability (arm-level deletions of 19p and q, focal-deletion of DNA-repair pathway genes and NOTCH1, amplification of EGFR) (77%), and increased APOBEC-activity (23%) are also observed. These additional alterations were present singly (18% of patients) or in combination (68%). Some of these alterations likely impact immune-dynamics of the evolving transformed tissue; progression to malignancy is associated with immune suppression through infiltration of regulatory T-cells (56%), depletion of cytotoxic T-cells (68%), and antigen-presenting dendritic cells (72%), with a concomitant increase in inflammation (92%). Patients can be grouped into three clusters by the estimated time to development of cancer from precancer by acquiring additional mutations (range: 4-10 years). Our findings provide deep molecular insights into the evolutionary processes and trajectories of oral cancer initiation and progression. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Molecular phylogeny and missense mutations at envelope proteins across coronaviruses.

Envelope (E) protein is one of the structural viroporins (76-109 amino acids long) present in the coronavirus. Sixteen sequentially different E-proteins were observed from a total of 4917 available complete SARS-CoV-2 genomes as on 18th June 2020 in the NCBI database. The missense mutations over the envelope protein across various coronaviruses of the ß-genus were analyzed to know the immediate parental origin of the envelope protein of SARS-CoV-2. The evolutionary origin is also endorsed by the phylogenetic analysis of the envelope proteins comparing sequence homology as well as amino acid conservations.

SARS-CoV2 envelope protein: non-synonymous mutations and its consequences.

In the NCBI database, as on June 6, 2020, total number of available complete genome sequences of SARS-CoV2 across the world is 3617. The envelope (E) protein of SARS-CoV2 possesses several non-synonymous mutations over the transmembrane and C-terminus domains in 15 (0.414%) genomes among 3617 SARS-CoV2 genomes, analyzed. More precisely, 10(0.386%) out of 2588 genomes from the USA, 3(0.806%) from Asia, 1 (0.348%) from Europe and 1 (0.274%) from Oceania contained the missense mutations over the E-protein of SARS-CoV2 genomes. The C-terminus motif DLLV has been to DFLV and YLLV in the proteins from QJR88103 (Australia: Victoria) and QKI36831 (China: Guangzhou) respectively, which might affect the binding of this motif with the host protein PALS1.

Missense mutations in SARS-CoV2 genomes from Indian patients.

The genetic diversity of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV2) in several countries sums up to worldwide genetic diversity. In this present study, variations in terms of missense mutations among the SARS-CoV2 genomes from 128 Indian patients, as of May 2020, are accounted and thereby some key findings with some hypotheses were made. These mutations across various genes of these genomes show wide genetic variations in sequence and rapid evolution of SARS-CoV2 virus. The presence of unique mutations in the studied SARS-CoV2 genomes may led to their attenuation. Few Genes such as ORF6, ORF10 are free from any mutations in the Indian context of 339 genomes observed as of 14th July 2020. Further, E protein contains only one mutation. This may suggest that designing a therapeutic approach against ORF6, ORF10 and E genes may have a beneficial effect in controlling COVID-19 pandemic especially in India.

A novel hotspot and rare somatic mutation p.A138V, at TP53 is associated with poor survival of pancreatic ductal and periampullary adenocarcinoma patients.

BACKGROUND: Pancreatic Ductal Adenocarcinoma (PDAC) is a cancer of the exocrine pancreas and 5-year survival rates remain constant at 7%. Along with PDAC, Periampullary Adenocarcinoma (PAC) accounts for 0.5-2% of all gastrointestinal malignancies. Genomic observations were well concluded for PDAC and PACs in western countries but no reports are available from India till now. METHODS: Targeted Next Generation Sequencing were performed in 8 (5 PDAC and 3 PAC) tumour normal pairs, using a panel of 412 cancer related genes. Primary findings were replicated in 85 tumour samples (31 PDAC and 54 PAC) using the Sanger sequencing. Mutations were also validated by ASPCR, RFLP, and Ion Torrent sequencing. IHC along with molecular dynamics and docking studies were performed for the p.A138V mutant of TP53. Key polymorphisms at TP53 and its associated genes were genotyped by PCR-RFLP method and association with somatic mutations were evaluated. All survival analysis was done using the Kaplan-Meier survival method which revealed that the survival rates varied significantly depending on the somatic mutations the patients harboured. RESULTS: Among the total 114 detected somatic mutations, TP53 was the most frequently mutated (41%) gene, followed by KRAS, SMAD4, CTNNB1, and ERBB3. We identified a novel hotspot TP53 mutation (p.A138V, in 17% of all patients). Low frequency of KRAS mutation (33%) was detected in these samples compared to patients from Western counties. Molecular Dynamics (MD) simulation and DNA-protein docking analysis predicted p.A138V to have oncogenic characteristics. Patients with p.A138V mutation showed poorer overall survival (p = 0.01). So, our finding highlights elevated prevalence of the p53p.A138V somatic mutation in PDAC and pancreatobiliary PAC patients. CONCLUSION: Detection of p.A138V somatic variant in TP53 might serve as a prognostic marker to classify patients. It might also have a role in determining treatment regimes. In addition, low frequency of KRAS hotspot mutation mostly in Indian PDAC patient cohort indicates presence of other early drivers in malignant transformation.

Application of Random Forest and data integration identifies three dysregulated genes and enrichment of Central Carbon Metabolism pathway in Oral Cancer.

BACKGROUND: Studies of epigenomic alterations associated with diseases primarily focus on methylation profiles of promoter regions of genes, but not of other genomic regions. In our past work (Das et al. 2019) on patients suffering from gingivo-buccal oral cancer - the most prevalent form of cancer among males in India - we have also focused on promoter methylation changes and resultant impact on transcription profiles. Here, we have investigated alterations in non-promoter (gene-body) methylation profiles and have carried out an integrative analysis of gene-body methylation and transcriptomic data of oral cancer patients. METHODS: Tumor and adjacent normal tissue samples were collected from 40 patients. Data on methylation in the non-promoter (gene-body) regions of genes and transcriptome profiles were generated and analyzed. Because of high dimensionality and highly correlated nature of these data, we have used Random Forest (RF) and other data-analytical methods. RESULTS: Integrative analysis of non-promoter methylation and transcriptome data revealed significant methylation-driven alterations in some genes that also significantly impact on their transcription levels. These changes result in enrichment of the Central Carbon Metabolism (CCM) pathway, primarily by dysregulation of (a) NTRK3, which plays a dual role as an oncogene and a tumor suppressor; (b) SLC7A5 (LAT1) which is a transporter dedicated to essential amino acids, and is overexpressed in cancer cells to meet the increased demand for nutrients that include glucose and essential amino acids; and, (c) EGFR which has been earlier implicated in progression, recurrence, and stemness of oral cancer, but we provide evidence of epigenetic impact on overexpression of this gene for the first time. CONCLUSIONS: In rapidly dividing cancer cells, metabolic reprogramming from normal cells takes place to enable enhanced proliferation. Here, we have identified that among oral cancer patients, genes in the CCM pathway - that plays a fundamental role in metabolic reprogramming - are significantly dysregulated because of perturbation of methylation in non-promoter regions of the genome. This result compliments our previous result that perturbation of promoter methylation results in significant changes in key genes that regulate the feedback process of DNA methylation for the maintenance of normal cell division.

SNP rs1049430 in the 3'-UTR of SH3GL2 regulates its expression: Clinical and prognostic implications in head and neck squamous cell carcinoma.

Single nucleotide polymorphisms (SNPs) in the 3'-UTR region are emerging cis-regulatory factors associated with the occurrences of several human diseases. SH3GL2, which is located at chromosome 9p21-22, is associated with hyperplastic/mildly dysplastic lesions of the head and neck and has a long 3'-UTR with multiple SNPs. The aim of the present study was to determine the susceptible allele(s) in the 3'-UTR SNPs of SH3GL2 in head and neck squamous cell carcinoma (HNSCC). First, we screened the genotypes of all SNPs located in the 3'-UTR of SH3GL2 in 110 controls and 147 cases in Indian populations by sequencing. A SNP (rs1049430:>G/T) that showed only heterozygosity was further confirmed by genotyping with an Illumina GoldenGate platform in 530 controls and 764 cases. Genotype-specific survival analysis of the HNSCC patients was performed. In addition, genotype-specific mRNA stability, isoform expression and protein expression were analyzed. SNP rs1049430 was not associated with disease occurrence, but it was associated with poor patient outcome. The G allele was associated with decreased SH3GL2 mRNA stability, differential splicing and low protein expression. Thus, our data demonstrate that the presence of the susceptible G allele in SNP rs1049430 is associated with the inactivation of SH3GL2 and could be used as a prognostic marker of HNSCC.

Genome-wide mitochondrial DNA sequence variations and lower expression of OXPHOS genes predict mitochondrial dysfunction in oral cancer tissue.

Several studies reported that mtDNA mutations may play important roles in carcinogenesis although the mechanism is not clear yet. Most of the studies compared mtDNA sequences in a tumor with those in normal tissues from different individuals ignoring inter-individual variations. In this study, 271 SNPs, 7 novel SNPs (or SNVs), and 15 somatic mutations were detected in mtDNA of 8 oral cancer tissues with respect to reference (rCRS) and adjacent normal tissues, respectively, using Ion PGM next generation sequencing method. Most of the sequence variations (76 SNPs and 1 somatic) are present in D-loop region followed by CyB (36 SNPs), ATP6 (24 SNPs), ND5 (17 SNPs and 5 somatic), ND4 (18 coding and 2 somatic) and other non-coding and coding DNA sequences. A total of 53 and 8 non-synonymous SNPs and somatic mutations, respectively, were detected in tumor tissues and some of these variations may have deleterious effects on the protein function as predicted by bioinformatic analysis. Moreover, significantly low mtDNA contents and expression of several mitochondrial genes in tumor compared to adjacent normal tissues may have also affected mitochondrial functions. Taken together, this study suggests that mtDNA mutations as well as low expression of mtDNA coded genes may play important roles in tumor growth. Although the sample size is low, an important aspect of the study is the use of adjacent control tissues to find out somatic mutations and a change in the expression of mitochondrial genes, to rule out inter-individual and inter-tissue variations which are important issues in the study of mitochondrial genomics.

D-loop somatic mutations and ∼5 kb "common" deletion in mitochondrial DNA: important molecular markers to distinguish oral precancer and cancer.

Apart from genomic DNA, mutations at mitochondrial DNA (mtDNA) have been hypothesized to play vital roles in cancer development. In this study, ∼5 kb deletion and D-loop mutations in mtDNA and alteration in mtDNA content were investigated in buccal smears from 104 healthy controls and 74 leukoplakia and 117 cancer tissue samples using Taqman-based quantitative assay and re-sequencing. The ∼5 kb deletion in mtDNA was significantly less (9.8 and 10.5 folds, P < 0.0001) in cancer tissues compared to control and leukoplakia tissues, respectively. On the other hand, somatic mutations in D-loop, investigated in 54 controls, 50 leukoplakias and 56 cancer patients, were found to be significantly more in cancer tissues, but not in leukoplakia tissues, compared to control (Z-score = 5.4). MtDNA contents were observed to be significantly more in leukoplakia (2.1 folds, P = 0.004) and cancer (1.6 folds, P = 0.03) tissues compared to control tissues. So, D-loop somatic mutations and ∼5 kb deletion patterns could be used as distinguishing markers between precancer and cancer tissues. This observation further suggests that somatic mutations in D-loop may facilitate carcinogenesis and cancer cells with less ∼5 kb deletion, i.e., intact mtDNA, may become resistant to apoptosis.

Genetic analysis of familial spontaneous pneumothorax in an Indian family.

Familial spontaneous pneumothorax is one of the phenotypes of Birt-Hogg-Dubé syndrome (BHDS), an autosomal dominant condition associated with folliculin (FLCN). We investigated clinical and genetic data of an Indian family having two patients suffering from spontaneous pneumothorax in the absence of skin lesions or renal tumors. HRCT scan of patient's lung revealed paracardiac cysts, and DNA sequencing of all 14 exons of FLCN from patients showed the presence of heterozygous "C allele" deletion in the poly-cytosine (poly-C) tract of exon 11 leading to truncated folliculin. This mutation was also observed in four asymptomatic members of the family. Our results confirmed the presence of deletion mutation in poly-C tract of FLCN in members of BHDS family. This is the first report of genetic insight in a BHDS family from India but in-depth studies with a larger sample set are necessary to understand mechanism of familial pneumothorax.

Genetic variations at microRNA and processing genes and risk of oral cancer.

Genetic variations at microRNA and microRNA processing genes are known to confer risk of cancer in different populations. Here, we studied variations at eight microRNA (miRNA) and four miRNA processing genes in 452 controls and 451 oral cancer patients by TaqMan genotyping assays. Variant allele-containing genotypes at mir-196a2 and variant allele homozygous genotype at Ran increased the risk of cancer significantly [adjusted odds ratio (OR) (95% confidence interval (CI)) = 1.3 (1-1.7) and 2.3 (1.1-4.6), respectively]. Conversely, variant allele-containing genotypes at mir-34b and variant allele homozygous genotype at Gemin3 reduced the risk of cancer significantly [adjusted OR (95% CI) = 0.7 (0.5-0.9) and 0.6 (0.4-1), respectively]. Cumulative risk was also increased by three times with increase in the number of risk alleles at these four loci. In tobacco stratified analysis, variant allele homozygous genotypes at mir-29a and Ran increased [adjusted OR (95% CI) = 1.5 (1-2.3) and 3 (1.1-8.4) respectively], while variant allele-containing genotypes at mir-34b decreased [adjusted OR (95% CI) = 0.6 (0.4-0.9)] the risk of cancer significantly. Thus, genetic variation at miRNA and processing genes altered the risk of oral cancer in this population thereby corroborating studies in other populations. However, it is necessary to validate this result in different Indian sub populations with larger sample sizes and examine the effect of these variations in tumour tissues to explain the mechanism of risk alteration.

Association between risk of oral precancer and genetic variations in microRNA and related processing genes.

BACKGROUND: MicroRNAs have been implicated in cancer but studies on their role in precancer, such as leukoplakia, are limited. Sequence variations at eight miRNA and four miRNA processing genes were studied in 452 healthy controls and 299 leukoplakia patients to estimate risk of disease. RESULTS: Genotyping by TaqMan assay followed by statistical analyses showed that variant genotypes at Gemin3 and mir-34b reduced risk of disease [OR = 0.5(0.3-0.9) and OR = 0.7(0.5-0.9) respectively] in overall patients as well as in smokers [OR = 0.58(0.3-1) and OR = 0.68(0.5-0.9) respectively]. Among chewers, only mir29a significantly increased risk of disease [OR = 1.8(1-3)]. Gene-environment interactions using MDR-pt program revealed that mir29a, mir34b, mir423 and Xpo5 modulated risk of disease (p < 0.002) which may be related to change in expression of these genes as observed by Real-Time PCR assays. But association between polymorphisms and gene expressions was not found in our sample set as well as in larger datasets from open access platforms like Genevar and 1000 Genome database. CONCLUSION: Variations in microRNAs and their processing genes modulated risk of precancer but further in-depth study is needed to understand mechanism of disease process.

Differential haplotype amplification leads to misgenotyping of heterozygote as homozygote when using single nucleotide mismatch primer.

Mismatches at the 3'end of /or within a primer are reported to affect the efficiency of PCR and cause allele drop. Here, we report preferential amplification of one haplotype and misgenotyping, when double heterozygotes at NAT1 (rs1057126 and rs15561) were genotyped by sequencing and PCR-RFLP methods using mismatch reverse primers located next to the target SNP. Detailed study revealed highest (100%) and lowest (0%) misgenotyping when the mismatch was at the 3rd and 15th nucleotide positions from 3' end of the primer, respectively. But, the same primers, without any mismatch genotyped heterozygotes correctly. Homozygotes can always be detected correctly irrespective of mismatch position in the primer. Similar results were observed for two SNPs (rs12947788 and rs 12951053) at TP53. Using mismatch NAT1 reverse primers, located three nucleotides away from the target SNP, both TaqMan and sequencing methods showed preferential synthesis of one haplotype strand and misgenotyping in heterozygotes, respectively. So, mismatch primer, located next to target SNP, should be avoided to genotype heterozygotes, since, PCR and sequencing based genotyping methods may lead the investigators to report faulty allelic and genotypic frequencies. This study mimics a situation when an unknown variation is present in the primer-binding sites of both chromosomes.

Association between polymorphisms at N-acetyltransferase 1 (NAT1) & risk of oral leukoplakia & cancer.

BACKGROUND & OBJECTIVES: N-acetyltransferases 1 and 2 (NAT1 and NAT2) are important enzymes for metabolism of tobacco carcinogens. Due to polymorphisms, improper activities of these enzymes might lead to the formation of DNA adducts that may modulate risk of tobacco related oral precancer and cancer. Previously, it was shown that NAT2 polymorphisms did not modulate the risk of oral precancer and cancer. We undertook this study to check whether polymorphisms at NAT1 can modulate the risk of oral leukoplakia and cancer either alone or in combination with NAT2. METHODS: Genotypes at four SNPs on NAT1 were determined by TaqMan method in 389 controls, 224 leukoplakia and 310 cancer patients. Genotype data were analyzed to know haplotypes and acetylation status of individuals and, then to estimate the risk of diseases. Using our previously published NAT2 data, combination of NAT1 and NAT2 acetylation genotypes of patients and controls were also analyzed to estimate the risk of diseases. RESULTS: Analysis of NAT1 genotype data revealed that 1088T and 1095C alleles exist in strong linkage disequilibrium (r 2 =0.97, P<0.0001) and SNPs are in Hardy-Weinberg Equilibrium (P=0.1). Wild type or normal acetylating and variant or rapid acetylating alleles were two major alleles (frequencies 0.62 and 0.36, respectively) present in the control population. NAT1 rapid acetylation could not modulate the risk of leukoplakia and cancer (OR=0.9, 95% CI: 0.6-1.3; OR=1.0, 95% CI: 0.7-1.4, respectively). Analysis of combined NAT1 and NAT2 acetylating data also showed no significant enhancement of the risk of diseases. INTERPRETATION & CONCLUSIONS: NAT1 rapid acetylation alone as well as combination of NAT1 rapid-NAT2 slow acetylation did not modulate the risk of oral precancer and cancer in our patient population. So, NAT1/NAT2 metabolized carcinogen products may not be involved in tobacco related oral precancer and cancer. It may be interpreted that large sample size as well as combination of polymorphisms at other candidate loci may be important to estimate the risk of a complex disease like oral cancer.

"Loss of polarity of basal cells" - The term revisited in the light of genomics.

The concept of polarity in development, homeostasis and pathological alteration of tissues has emerged as an interesting aspect of the concerned biology. The epithelial cells exhibit apicobasal polarity which is maintained by Crumbs complex located at apical region of tight junction, 'PAR' complex at sub-apical region of tight junction and Scribble complex at adherens junction. Any functional perturbation of these proteins cause alteration of normal epithelial physiology en-route to epithelial pathology. In this maiden scientific exercise, we have tried to explore the association of expression of cell polarity proteins in OPMD and OSCC. Here, we have chosen DLG1 as a representative of polarity protein. RNA was isolated from the tissue samples. Then cDNA was prepared by RTPCR technique. qPCR was performed on cDNA samples. Expression data was analysed on the basis of Ct values. Paired t-test was performed with normalized Ct values of disease and normal tissue to determine whether there was any significant difference in expression between them. The statistical tests were done using SPSS software. Results of this study reflected increase in DLG1 expression in high grade dysplasia. There was no significant alteration in expression of DLG1 in WDSCC where there is formation of cluster of neoplastic cells ultimately producing epithelial islands and keratin pearls. But, in case of MDSCC, when the same neoplastic cells keep on invading with minimal keratin pearl formation, they again gain the mesenchymal character in full potential. And this phenomenon supports the upregulation of DLG1 in MDSCC in our study.

Homozygous GRHPR C.494G>A mutation is deleterious that causes early onset of nephrolithiasis in West Bengal, India.

Pediatric nephrolithiasis (NL) or Kidney stone disease (KSD) is an untethered topic in Asian population. In Western countries, the annual incidence of paediatric NL is around 6-10%. Here, we present data from West Bengal, India, on lower age (LA, 0-20 years) NL and its prevalence for the first time. To discover the mutations associated with KSD, twenty-four (18 + 6) rare LA-NL patients were selected for Whole Exome Sequencing (WES) and Sanger sequencing, respectively. It was found that GRHPR c. 494G>A mutation (MZ826703) is predominant in our study cohort. This specific homozygous mutation is functionally studied for the first time directly from human peripheral mononuclear cell (PBMC) samples. Using expression study with biochemical activity and computational analysis we assumed that the mutation is pathogenic with loss of function. Moreover, three genes, AGXT, HOGA1 and GRHPR with Novel variants known to cause hyperoxaluria were found frequently in the study cohort. Our study analyses the genes and variations that cause LA-NL, as well as the molecular function of the GRHPR mutation, which may serve as a clinical marker in the population of West Bengal, Eastern India.

Genetic insight into Birt-Hogg-Dubé syndrome in Indian patients reveals novel mutations at FLCN.

BACKGROUND: Birt-Hogg-Dubé syndrome (BHDS) is a rare monogenic condition mostly associated with germline mutations at FLCN. It is characterized by either one or more manifestations of primary spontaneous pneumothorax (PSP), skin fibrofolliculomas and renal carcinoma (chromophobe). Here, we comprehensively studied the mutational background of 31 clinically diagnosed BHDS patients and their 74 asymptomatic related members from 15 Indian families. RESULTS: Targeted amplicon next-generation sequencing (NGS) and Sanger sequencing of FLCN in patients and asymptomatic members revealed a total of 76 variants. Among these variants, six different types of pathogenic FLCN mutations were detected in 26 patients and some asymptomatic family members. Two of the variants were novel mutations: an 11-nucleotide deletion (c.1150_1160delGTCCAGTCAGC) and a splice acceptor mutation (c.1301-1G > A). Two variants were Clinvar reported pathogenic mutations: a stop-gain (c.634C > T) and a 4-nucleotide duplication (c.1329_1332dupAGCC). Two known variants were: hotspot deletion (c.1285delC) and a splice donor mutation (c.1300 + 1G > A). FLCN mutations could not be detected in patients and asymptomatic members from 5 families. All these mutations greatly affected the protein stability and FLCN-FNIP2 interaction as observed by molecular docking method. Family-based association study inferred pathogenic FLCN mutations are significantly associated with BHDS. CONCLUSION: Six pathogenic FLCN mutations were detected in patients from 10 families out of 15 families in the cohort. Therefore, genetic screening is necessary to validate the clinical diagnosis. The pathogenic mutations at FLCN affects the protein-protein interaction, which plays key roles in various metabolic pathways. Since, pathogenic mutations could not be detected in exonic regions of FLCN in 5 families, whole genome sequencing is necessary to detect all mutations at FLCN and/or any undescribed gene/s that may also be implicated in BHDS.

Rare mutations in the accessory proteins ORF6, ORF7b, and ORF10 of the SARS-CoV-2 genomes.

A total number of 3080 SARS-CoV-2 genomes from all continents are considered from the NCBI database. Every accessory protein ORF6, ORF7b, and ORF10 of SARS-CoV-2 possess a single missense mutation in less than 1.5% of the 3080 genomes. It has now been observed that different non-synonymous mutations occurred in these three accessory proteins. Most of these rare mutations are changing the amino acids such as hydrophilic to hydrophobic, acidic or basic to hydrophobic, and vice versa etc. So these highly conserved proteins might play an essential role in virus pathogenicity. This study opens a question whether it carries some messages about the virus rapid replications, and virulence.

Pathogenic perspective of missense mutations of ORF3a protein of SARS-CoV-2.

One of the most important proteins for COVID-19 pathogenesis in SARS-CoV-2 is the ORF3a which is the largest accessory protein among others coded by the SARS-CoV-2 genome. The major roles of the protein include virulence, infectivity, ion channel activity, morphogenesis, and virus release. The coronavirus, SARS-CoV-2 is mutating rapidly, therefore, critical study of mutations in ORF3a is certainly important from the pathogenic perspective. Here, a sum of 175 non-synonymous mutations in the ORF3a of SARS-CoV-2 were identified from 7194 complete genomes of SARS-CoV-2 available from NCBI database. Effects of these mutations on structural stability, and functions of ORF3a were also studied. Broadly, three different classes of mutations, such as neutral, disease, and mixed (neutral and disease) types of mutations were observed. Consecutive phenomena of mutations in ORF3a protein were studied based on the timeline of detection of the mutations. Considering the amino acid compositions of the ORF3a protein, twenty clusters were detected using the K-means clustering method. The present findings on 175 novel mutations of ORF3a proteins will extend our knowledge on ORF3a, a vital accessory protein in SARS-CoV-2, to enlighten the pathogenicity of this life-threatening virus.

Genes and pathways monotonically dysregulated during progression from normal through leukoplakia to gingivo-buccal oral cancer.

Oral squamous cell carcinoma of the gingivo-buccal region (OSCC-GB) accounts for the highest cancer morbidity and mortality among men in India. It has been observed that about one-third of individuals with oral leukoplakia, a dysplastic precancerous lesion in the oral cavity, progress to oral cancer. We aimed to identify systematic transcriptomic changes as a normal tissue in the oral cavity progresses to frank OSCC-GB. Seventy-two OSCC-GB patients, from multiple hospitals, were recruited, and transcriptome analysis of tumor and adjacent normal tissue (of all patients) and adjacent leukoplakia tissue (of a subset of 25 unselected patients with concomitant leukoplakia) was performed. We have identified many differences in the transcriptomic profiles between OSCC-GB and squamous cell carcinoma of the head and neck regions. Compared to the normal/precancerous tissue, significant enrichment of ECM-receptor interaction, PI3K-Akt signaling, cytokine-cytokine receptor interaction, focal adhesion, and cell cycle pathways were observed in OSCC-GB. Using gene set enrichment analysis, we identified a profound role of interferon receptor signaling in tumor growth by activating immune evasion mechanisms. The role of tumor-infiltrating immune cells further supported the growth and immunosuppressive mechanism of tumor tissues. Some immune evasion genes-CD274, CD80, and IDO1-were found to be activated even in the precancerous tissue. Taken together, our findings provide a clear insight into the sequential genetic dysregulation associated with progression to oral cancer. This insight provides a window to the development of predictive biomarkers and therapeutic targets for gingivo-buccal oral cancer.