Single-cell transcriptomic analysis of gingivo-buccal oral cancer reveals two dominant cellular programs.

Oral squamous cell carcinoma of the gingivo-buccal region (OSCC-GB) is the most common cancer among men in India, and is associated with poor prognosis and frequent recurrence. Cellular heterogeneity in OSCC-GB was investigated by single-cell RNA sequencing of tumors derived from the oral cavity of 12 OSCC-GB patients, 3 of whom had concomitant presence of a precancerous lesion (oral submucous fibrosis [OSMF]). Unique malignant cell types, features, and phenotypic shifts in the stromal cell population were identified in oral tumors with associated submucous fibrosis. Expression levels of FOS, ATP1A, and DUSP1 provided robust discrimination between tumors with or without the concomitant presence of OSMF. Malignant cell populations shared between tumors with and without OSMF were enriched with the expression of partial epithelial-mesenchymal transition (pEMT) or fetal cell type signatures indicative of two dominant cellular programs in OSCC-GB-pEMT and fetal cellular reprogramming. Malignant cells exhibiting fetal cellular and pEMT programs were enriched with the expression of immune-related pathway genes known to be involved in antitumor immune response. In the tumor microenvironment, higher infiltration of immune cells than the stromal cells was observed. The T cell population was large in tumors and diverse subtypes of T cells with varying levels of infiltration were found. We also detected double-negative PLCG2+ T cells and cells with intermediate M1-M2 macrophage polarization. Our findings shed light on unique aspects of cellular heterogeneity and cell states in OSCC-GB.

Genomic associations with antibody response to an oral cholera vaccine.

Oral cholera vaccine is one of the key interventions used in our fight to end the longest pandemic of our time, cholera. The immune response conferred by the currently available cholera vaccines, as measured by serum antibody levels, is variable amongst its recipients. We undertook a genome wide association study (GWAS) on antibody response to the cholera vaccine; globally, the first GWAS on cholera vaccine response. We identified three clusters of bi-allelic SNPs, in high within-cluster linkage disequilibrium that were moderately (p < 5 × 10-6) associated with antibody response to the cholera vaccine and mapped to chromosomal regions 4p14, 4p16.1 and 6q23.3. Intronic SNPs of TBC1D1 comprised the cluster on 4p14, intronic SNPs of TBC1D14 comprised that on 4p16.1 and SNPs upstream of TNFAIP3 formed the cluster on 6q23.3. SNPs within and around these clusters have been implicated in immune cell function and immunological aspects of autoimmune or infectious diseases (e.g., diseases caused by Helicobacter pylori and malarial parasite). 6q23.3 is a prominent region harbouring many loci associated with immune related diseases, including multiple sclerosis, rheumatoid arthritis and systemic lupus erythematosus, as well as IL2 and INFα response to a smallpox vaccine. The gene clusters identified in this study play roles in vesicle-mediated pathway, autophagy and NF-κB signaling. No significant effect of O blood group on antibody response to the cholera vaccine was observed in this study.

Identification of HPV16 positive cervical cancer subsets characterized by divergent immune and oncogenic phenotypes with potential implications for immunotherapy.

BACKGROUND: Cervical cancers (CaCx), like many other cancer types, portray high molecular heterogeneity that affects response to therapy, including immunotherapy. In India and other developing countries, CaCx mortality rates are very high because women report to the clinics with advanced cancers in absence of organized screening programs. This calls for implementation of newer therapeutic regimens for CaCx, like immunotherapy, which is again not used commonly in such countries. OBJECTIVE: Therefore, we focused on dissecting tumour immune heterogeneity, if any, identify immune gene-based biomarkers of heterogeneity and subsets of such cancers with the potential for immunotherapy. We also attempted to characterize the cancer-associated phenotypes of such subsets, including viral load, to decipher the relationship of tumour immunogenicity with oncogenicity. METHODS: Employing RNA-seq analysis of 44 HPV16 positive CaCx patients, immune subtypes were identified by unsupervised hierarchical clustering of global immune-gene expression profiles. Proportions of tumor infiltrating immune cells in the tumor milieu were estimated, employing Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT), using gene expression data from RNA-seq. The oncogenic phenotypes of the immune subtypes of CaCx were deciphered through differential gene expression (DEGs) and pathway enrichment analysis. Viral load was estimated through TaqMan-based qRT-PCR analysis. RESULTS: Analysis revealed the presence of two immune subtypes of CaCx, A (26/44; 59.09%) and B (18/44; 40.90%). Compared to Subtype-A, Subtype-B portrayed overexpression of immune genes and high infiltration of immune cells, specifically CD8+ T cells (p < 0.0001). Besides, a significant correlation between PD-1 and PD-L1 co-expression among Subtype-B, as opposed to Subtype-A, confirmed the interactive roles of these immune checkpoint molecules in Subtype B. Stepwise discriminant analysis pin-pointed ten immune-genes that could classify 100% of the patients significantly (p < 0.0001) into the two immune subtypes and serve as potential biomarkers of CaCx immunity. Differential gene expression analysis between the subtypes unveiled that Subtype-B was more biologically aggressive than Subtype-A, reflecting loss of structural integrity and promotion of cancer progression. The viral load was significantly lower in Subtype-B (average viral load = 10.74/100 ng of genomic DNA) compared to Subtype-A (average viral load = 14.29/100 ng of genomic DNA). Thus viral load and the ten-gene panel underscore their association with immunogenicity and oncogenicity. CONCLUSION: Our study provides strong evidence that only a subset, about 41% of HPV16 positive CaCx patients in India, portray immune enrichment of the tumor milieu coupled with aggressive phenotypes. Such subtypes are therefore likely to benefit through checkpoint molecule-based or tumor infiltrating lymphocyte-based immunotherapy, which could be a leap forward in tackling aggressive forms of such CaCx in India and other developing countries.

Overexpression of CD73 is associated with recurrence and poor prognosis of gingivobuccal oral cancer as revealed by transcriptome and deep immune profiling of paired tumor and margin tissues.

BACKGROUND: For various cancers, differences in response to treatment and subsequent survival period have been reported to be associated with variation in immune contextures. AIM: We sought to identify whether such association exists in respect of gingivobuccal oral cancer. MATERIALS AND METHODS: We performed deep immune profiling of tumor and margin tissues collected from 46 treatment naïve, Human Papillomavirus (HPV) negative, patients. Each patient was followed for 24 months and prognosis (recurrence/death) noted. Key findings were validated by comparing with TCGA-HNSC cohort data. RESULTS: About 28% of patients showed poor post-treatment prognosis. These patients exhibited a high probability of recurrence even within 1 year and death within 2 years. There was restricted immune cell infiltration in tumor, but not in margin, among these patients. Reduced expression of eight immune-related genes (IRGs) (NT5E, THRA, RBP1, TLR4, ITGA6, BMPR1B, ITGAV, SSTR1) in tumor strongly predicted better quality of prognosis, both in our patient cohort and in TCGA-HNSC cohort. Tumors of patients with better prognosis were associated with (a) lower CD73+ cells with concomitant lower expression level of NT5E/CD73, (b) higher proportions of CD4+ and CD8+ T cells, B cells, NK cells, M1 macrophages, (c) higher %Granzyme+ cells, (d) higher TCR and BCR repertoire diversities. CD73 expression in tumor was associated with low CD8+ and CD4+ T cells, low immune repertoire diversity, and advanced cancer stage. DISCUSSION AND CONCLUSION: High infiltration of anti-tumor immune cells in both tumors and margins results in good prognosis, while in patients with minimal infiltration in tumors in spite of high infiltration in margins results in poor prognosis. Targeted CD73 immune-checkpoint inhibition may improve clinical outcome.

Genome-wide temporal landscaping of DNA methylation in pregnant women delivering at term: a GARBH-InI study.

Background: We performed an epigenome-wide longitudinal DNA methylation study on an Indian cohort of pregnant women, GARBH-Ini, at three time points during pregnancy and at delivery. Aim & objective: Our aim was to identify temporal DNA methylation changes in maternal peripheral blood during the period of gestation and assess their impact on biological pathways critical for term delivery. Results: Significantly differentially methylated CpGs were identified by linear mixed model analysis (Bonferroni p < 0.01) and classified into two distinct temporal methylation trends: increasing and decreasing during gestation. Genes with upward methylation trend were enriched for T-cell activity, while those with a downward trend were enriched for solute transport and cell structure organization functions. Conclusion: Consistent trends of DNA methylation in maternal peripheral blood point to the sentinel function of T cells in the maintenance of pregnancy, and the importance of coordinated cellular remodeling to facilitate term delivery.

DNA methylation is the addition of a methyl group to the molecular structure of DNA, which then alters the gene expression. The goal of the study was to find out how DNA methylation patterns change over time during pregnancy and how these changes are related to the biological processes that are important for the delivery of a healthy baby at full term. Using statistical modeling, we identified specific patterns of DNA methylation changes during pregnancy and classified them into two groups based on the direction of the changes. The genes associated with increasing methylation levels were related to the activities of T cells, which are important for the immune system. The genes associated with decreasing methylation levels were related to processes like transporting substances and organizing cell structures. In conclusion, our findings suggest that T cells play an important role in maintaining a healthy pregnancy, and the study highlights the importance of coordinated changes in cells to support a successful delivery of a baby at term.

Multiple migrations from East Asia led to linguistic transformation in NorthEast India and mainland Southeast Asia.

NorthEast India, with its unique geographic location in the midst of the Himalayas and Bay of Bengal, has served as a passage for the movement of modern humans across the Indian subcontinent and East/Southeast Asia. In this study we look into the population genetics of a unique population called the Khasi, speaking a language (also known as the Khasi language) belonging to the Austroasiatic language family and residing amidst the Tibeto-Burman speakers as an isolated population. The Khasi language belongs to one of the three major broad classifications or phyla of the Austroasiatic language and the speakers of the three sub-groups are separated from each other by large geographical distances. The Khasi speakers are separated from their nearest Austroasiatic language-speaking sub-groups: the "Mundari" sub-family from East and peninsular India and the "Mon-Khmers" in Mainland Southeast Asia. We found the Khasi population to be genetically distinct from other Austroasiatic speakers, i.e. Mundaris and Mon-Khmers, but relatively similar to the geographically proximal Tibeto Burmans. The possible reasons for this genetic-linguistic discordance lie in the admixture history of different migration events that originated from East Asia and proceeded possibly towards Southeast Asia. We found at least two distinct migration events from East Asia. While the ancestors of today's Tibeto-Burman speakers were affected by both, the ancestors of Khasis were insulated from the second migration event. Correlating the linguistic similarity of Tibeto-Burman and Sino-Tibetan languages of today's East Asians, we infer that the second wave of migration resulted in a linguistic transition while the Khasis could preserve their linguistic identity.

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.

dbGENVOC: database of GENomic Variants of Oral Cancer, with special reference to India.

Oral cancer is highly prevalent in India and is the most frequent cancer type among Indian males. It is also very common in southeast Asia. India has participated in the International Cancer Genome Consortium (ICGC) and some national initiatives to generate large-scale genomic data on oral cancer patients and analyze to identify associations and systematically catalog the associated variants. We have now created an open, web-accessible database of these variants found significantly associated with Indian oral cancer patients, with a user-friendly interface to enable easy mining. We have value added to this database by including relevant data collated from various sources on other global populations, thereby providing opportunities of comparative geographical and/or ethnic analyses. Currently, no other database of similar nature is available on oral cancer. We have developed Database of GENomic Variants of Oral Cancer, a browsable online database framework for storage, retrieval and analysis of large-scale data on genomic variants and make it freely accessible to the scientific community. Presently, the web-accessible database allows potential users to mine data on ∼24 million clinically relevant somatic and germline variants derived from exomes (n = 100) and whole genomes (n = 5) of Indian oral cancer patients; all generated by us. Variant data from The Cancer Genome Atlas and data manually curated from peer-reviewed publications were also incorporated into the database for comparative analyses. It allows users to query the database by a single gene, multiple genes, multiple variant sites, genomic region, patient ID and pathway identities. Database URL: http://research.nibmg.ac.in/dbcares/dbgenvoc/.

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.

SARS-CoV-2 mutation 614G creates an elastase cleavage site enhancing its spread in high AAT-deficient regions.

SARS-CoV-2 was first reported from China. Within three months, it evolved to 10 additional subtypes. Two evolved subtypes (A2 and A2a) carry a non-synonymous Spike protein mutation (D614G). We conducted phylodynamic analysis of over 70,000 SARS-CoV-2 coronaviruses worldwide, sequenced until July2020, and found that the mutant subtype (614G) outcompeted the pre-existing type (614D), significantly faster in Europe and North-America than in East Asia. Bioinformatically and computationally, we identified a novel neutrophil elastase (ELANE) cleavage site introduced in the G-mutant, near the S1-S2 junction of the Spike protein. We hypothesised that elevation of neutrophil elastase level at the site of infection will enhance the activation of Spike protein thus facilitating host cell entry for 614G, but not the 614D, subtype. The level of neutrophil elastase in the lung is modulated by its inhibitor α1-antitrypsin (AAT). AAT prevents lung tissue damage by elastase. However, many individuals exhibit genotype-dependent deficiency of AAT. AAT deficiency eases host-cell entry of the 614G virus, by retarding inhibition of neutrophil elastase and consequently enhancing activation of the Spike protein. AAT deficiency is highly prevalent in European and North-American populations, but much less so in East Asia. Therefore, the 614G subtype is able to infect and spread more easily in populations of the former regions than in the latter region. Our analyses provide a molecular biological and evolutionary model for the higher observed virulence of the 614G subtype, in terms of causing higher morbidity in the host (higher infectivity and higher viral load), than the non-mutant 614D subtype.

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.

Dysbiotic Lesional Microbiome With Filaggrin Missense Variants Associate With Atopic Dermatitis in India.

Background: Atopic Dermatitis (AD) has been associated with the loss of function (LoF) mutations in Filaggrin (FLG) gene and increase in relative abundance of specific microbes in the lesional skin, predominantly in Caucasians. Our study aims to determine, in Indian AD patients, (a) the prevalence of FLG LoF and missense mutations, and (b) the nature and extent of dysbiosis and altered microbial pathways with and without mutations in FLG. AD patients (n = 34) and healthy controls (n = 54) were recruited from India in this study and shotgun sequencing was carried out in a subset of samples with adequate microbiome DNA concentration. Host DNA from the same subset of samples was subjected to FLG coding region sequencing and host-microbiome association was estimated. Results: The prevalence of FLG LoFs that are associated with AD globally were significantly lesser in our cases and controls (8.6%, 0%) than those reported in Europeans (27%, 2.6%). Staphylococcus aureus was present only on AD skin [abundance in Pediatric AD: 32.86%; Adult AD: 22.17%], but not on healthy skin on which Staphylococcus hominis (Adult controls: 16.43%, Adult AD: 0.20%; p = 0.002), Cutibacterium acnes (Adult controls:10.84%, Adult AD: 0.90%; p = 0.02), and Malassezia globosa (Adult controls: 8.89%, Adult AD: 0.005%; p = 0.001) were significantly more abundant. Microbial pathways mostly associated with skin barrier permeability, ammonia production and inflammation (Arginine and Proline metabolism, Histidine Metabolism and Staphylococcus aureus infection) were significantly enriched on AD skin metagenome. These pathways are also reported to impair antimicrobial peptide activity. Among AD patients with missense single nucleotide polymorphisms harboring "potentially damaging" alleles in FLG gene, damaging allele dosage was significantly (p < 0.02) positively correlated with relative abundance of phylum_Proteobacteria up to order_Pseudomonadales and negatively correlated with phylum_Firmicutes up to species_Staphylococcus aureus. Conclusion: Our study has provided evidence that host DNA profile is significantly associated with microbiome composition in the development of AD. Species and strain level analysis showed that the microbial pathways enriched in AD cases were mostly found in MRSA strains. These evidences can be harnessed to control AD by modulating the microbiome using a personalized strategy. Our findings on the association of FLG genotypes with the microbiome dysbiosis may pave the way for a personalized strategy to provide a more effective control of AD.

A framework for pathway knowledge driven prioritization in genome-wide association studies.

Many variants with low frequencies or with low to modest effects likely remain unidentified in genome-wide association studies (GWAS) because of stringent genome-wide thresholds for detection. To improve the power of detection, variant prioritization based on their functional annotations and epigenetic landmarks has been used successfully. Here, we propose a novel method of prioritization of a GWAS by exploiting gene-level knowledge (e.g., annotations to pathways and ontologies) and show that it further improves power. Often, disease associated variants are found near genes that are coinvolved in specific biological pathways relevant to disease process. Utilization of this knowledge to conduct a prioritized scan increases the power to detect loci that map to genes clustered in a few specific pathways. We have developed a computationally scalable framework based on penalized logistic regression (termed GKnowMTest-Genomic Knowledge-guided Multiplte Testing) to enable a prioritized pathway-guided GWAS scan with a very large number of gene-level annotations. We demonstrate that the proposed strategy improves overall power and maintains the Type 1 error globally. Our method works on genome-wide summary level data and a user-specified list of pathways (e.g., those extracted from large pathway databases without reference to biology of a specific disease). It automatically reweights the input p values by incorporating the pathway enrichments as "adaptively learned" from the data using a cross-validation technique to avoid overfitting. We used whole-genome simulations and some publicly available GWAS data sets to illustrate the application of our method. The GKnowMTest framework has been implemented as a user-friendly open-source R package.

Hepatic transcriptome signature correlated with HOMA-IR explains early nonalcoholic fatty liver disease pathogenesis.

INTRODUCTION AND OBJECTIVES: Non-alcoholic fatty liver disease (NAFLD) is multistage with heterogeneous outcomes. We studied the influence of insulin resistance (IR) on the hepatic transcriptome of early NAFLD stages, to understand disease development. MATERIALS AND METHODS: In this cross-sectional study, possible clinicopathological risk factors were compared between mild-NAFL (N = 72) and non-alcoholic steatohepatitis (NASH; N = 51) patients. Liver tissue-transcriptome difference was studied between a subset of 25 mild-NAFL and 20 NASH biopsies and validated on another subset of 12 mild-NAFL and 13 NASH biopsies, using RT-PCR. The relationship between IR driven gene expression changes with fibrosis in NASH was investigated. RESULTS: Significantly higher weight (p = 0.005) and elevated levels of HbA1c (p = 0.009), FBG (p = 0.03) and HOMA-IR (p = 0.009) were found in NASH patients. Five differentially expressed genes (DEGs, fold change > 1.5) were identified in NASH-FABP4, FABP5L2, CD24, PRAP1, and SPP1. The DEGs were positively associated with disease severity and HOMA-IR, and were found to be efficient classifiers of mild-NAFL and NASH. Additional 1218 genes identified related to IR (IrCGs), which can classify NASH-with-fibrosis patients separately from mild-NAFL and NASH patients. IrCGs can promote intra-hepatic fat accumulation, dysregulation of the lipid metabolism, lipotoxicity, and activation of cell survival pathways including activation of cell proliferation and differentiation pathways. CONCLUSIONS: Hepatic expression of genes associated with insulin resistance may drive NAFLD development and progression.

Study of Caspase 8 mutation in oral cancer and adjacent precancer tissues and implication in progression.

It is hypothesized that same driver gene mutations should be present in both oral leukoplakia and cancer tissues. So, we attempted to find out mutations at one of the driver genes, CASP8, in cancer and adjacent leukoplakia tissues. Patients (n = 27), affected by both of cancer and adjacent leukoplakia, were recruited for the study. Blood and tissue DNA samples were used to identify somatic mutations at CASP8 by next generation sequencing method. In total, 56% (15 out of 27) cancer and 30% (8 out of 27) leukoplakia tissues had CASP8 somatic mutations. In 8 patients, both cancer and adjacent leukoplakia tissues, located within 2-5 cm of tumor sites, had identical somatic mutations. But, in 7 patients, cancer samples had somatic mutations but none of the leukoplakia tissues, located beyond 5cm of tumor sites, had somatic mutations. Mutated allele frequencies at CASP8 were found to be more in cancer compared to adjacent leukoplakia tissues. This study provides mutational evidence that oral cancer might have progressed from previously grown leukoplakia lesion. Leukoplakia tissues, located beyond 5cm of cancer sites, were free from mutation. The study implies that CASP8 mutation could be one of the signatures for some of the leukoplakia to progress to oral cancer.

Analysis of RNA sequences of 3636 SARS-CoV-2 collected from 55 countries reveals selective sweep of one virus type.

BACKGROUND & OBJECTIVES: SARS-CoV-2 (Severe acute respiratory syndrome coronavirus-2) is evolving with the progression of the pandemic. This study was aimed to investigate the diversity and evolution of the coronavirus SARS-CoV-2 with progression of the pandemic over time and to identify similarities and differences of viral diversity and evolution across geographical regions (countries). METHODS: Publicly available data on type definitions based on whole-genome sequences of the SARS-CoV-2 sampled during December and March 2020 from 3636 infected patients spread over 55 countries were collected. Phylodynamic analyses were performed and the temporal and spatial evolution of the virus was examined. RESULTS: It was found that (i) temporal variation in frequencies of types of the coronavirus was significant; ancestral viruses of type O were replaced by evolved viruses belonging to type A2a; (ii) spatial variation was not significant; with the spread of SARS-CoV-2, the dominant virus was the A2a type virus in every geographical region; (iii) within a geographical region, there was significant micro-level variation in the frequencies of the different viral types, and (iv) the evolved coronavirus of type A2a swept rapidly across all continents. INTERPRETATION & CONCLUSIONS: SARS-CoV-2 belonging to the A2a type possesses a non-synomymous variant (D614G) that possibly eases the entry of the virus into the lung cells of the host. This may be the reason why the A2a type has an advantage to infect and survive and as a result has rapidly swept all geographical regions. Therefore, large-scale sequencing of coronavirus genomes and, as required, of host genomes should be undertaken in India to identify regional and ethnic variation in viral composition and its interaction with host genomes. Further, careful collection of clinical and immunological data of the host can provide deep learning in relation to infection and transmission of the types of coronavirus genomes.

Epigenomic dysregulation-mediated alterations of key biological pathways and tumor immune evasion are hallmarks of gingivo-buccal oral cancer.

BACKGROUND: Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. Although genomic drivers of OSCC-GB have been identified by whole exome and whole genome sequencing, no epigenome-wide study has been conducted in OSCC-GB; our study has filled this gap, and has discovered and validated epigenomic hallmarks of gingivobuccal oral cancer. METHODS: We have carried out integrative analysis of epigenomic (n = 87) and transcriptomic (n = 72) profiles of paired tumor-normal tissues collected from OSCC-GB patients from India. Genome-wide DNA methylation assays and RNA-sequencing were performed on high-throughput platforms (Illumina) using a half-sample of randomly selected patients to discover significantly differentially methylated probes (DMPs), which were validated on the remaining half-sample of patients. RESULTS: About 200 genes showed significant inverse correlation between promoter methylation and expression, of which the most significant genes included genes that act as transcription factors and genes associated with other cancer types. Novel findings of this study include identification of (a) potential immunosuppressive effect in OSCC-GB due to significant promoter hypomethylation driven upregulation of CD274 and CD80, (b) significant dysregulation by epigenetic modification of DNMT3B (upregulation) and TET1 (downregulation); and (c) known drugs that can reverse the direction of dysregulation of gene expression caused by promoter methylation. CONCLUSIONS: In OSCC-GB patients, there are significant alterations in expression of key genes that (a) regulate normal cell division by maintenance of balanced DNA methylation and transcription process, (b) maintain normal physiological signaling (PPAR, B cell receptor) and metabolism (arachidonic acid) pathways, and (c) provide immune protection against antigens, including tumor cells. These findings indicate novel therapeutic targets, including immunotherapeutic, for treatment of OSCC-GB.

Lymph node metastasis in oral cancer is strongly associated with chromosomal instability and DNA repair defects.

Oral squamous cell carcinoma (OSCC) is highly prevalent in south and southeast Asia. Many (30-50%) OSCC patients develop lymph node metastasis (LNM), which is the most important prognostic factor in OSCC. To identify genomic correlates of LNM, we compared exome sequences and copy number variation data of blood and tumor DNA from highly contrasting subgroups of patients to reduce false inferences-(i) patients with LNM and (ii) patients with late stage disease but without LNM. We found that LNM is associated with (i) specific hotspot somatic mutations in TP53 and CASP8; (ii) rare nonsilent germline mutations in BRCA2 and FAT1; (iii) mutations in mito-G2/M and nonhomologous end joining (NHEJ) pathways; (iv) recurrent deletion of genes for DNA repair by homologous recombination; and (v) chromosomal instability. LN+ patients with NHEJ pathway mutations have longer disease-free survival. Five genomic features have a high predictive value of LNM.