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.

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.

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.

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.

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.

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.

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.

A powerful method to integrate genotype and gene expression data for dissecting the genetic architecture of a disease.

To decipher the genetic architecture of human disease, various types of omics data are generated. Two common omics data are genotypes and gene expression. Often genotype data for a large number of individuals and gene expression data for a few individuals are generated due to biological and technical reasons, leading to unequal sample sizes for different omics data. Unavailability of standard statistical procedure for integrating such datasets motivates us to propose a two-step multi-locus association method using latent variables. Our method is powerful than single/separate omics data analysis and it unravels comprehensively deep-seated signals through a single statistical model. Extensive simulation confirms that it is robust to various genetic models as its power increases with sample size and number of associated loci. It provides p-values very fast. Application to real dataset on psoriasis identifies 17 novel SNPs, functionally related to psoriasis-associated genes, at much smaller sample size than standard GWAS.

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.

A Pregnancy Cohort to Study Multidimensional Correlates of Preterm Birth in India: Study Design, Implementation, and Baseline Characteristics of the Participants.

Globally, preterm birth is a major public health problem. In India, 3.6 million of the 27 million infants born annually are preterm. Risk stratification of women based on multidimensional risk factors assessed during pregnancy is critical for prevention of preterm birth. A cohort study of pregnant women was initiated in May 2015 at the civil hospital in Gurugram, Haryana, India. Women are enrolled within 20 weeks of gestation and are followed until delivery and once postpartum. The objectives are to identify clinical, epidemiologic, genomic, epigenomic, proteomic, and microbial correlates; discover molecular-risk markers by using an integrative -omics approach; and generate a risk-prediction algorithm for preterm birth. We describe here the longitudinal study design, methodology of data collection, and the repositories of data, biospecimens, and ultrasound images being created. A total of 4,326 pregnant women, with documented evidence of recruitment before 20 weeks of gestation, have been enrolled through March 2018. We report baseline characteristics and outcomes of the first 2,000 enrolled participants. A high frequency of preterm births (14.9% among 1,662 live births) is noteworthy. The cohort database and the repositories will become global resources to answer critical questions on preterm birth and other birth outcomes.

Genomic reconstruction of the history of extant populations of India reveals five distinct ancestral components and a complex structure.

India, occupying the center stage of Paleolithic and Neolithic migrations, has been underrepresented in genome-wide studies of variation. Systematic analysis of genome-wide data, using multiple robust statistical methods, on (i) 367 unrelated individuals drawn from 18 mainland and 2 island (Andaman and Nicobar Islands) populations selected to represent geographic, linguistic, and ethnic diversities, and (ii) individuals from populations represented in the Human Genome Diversity Panel (HGDP), reveal four major ancestries in mainland India. This contrasts with an earlier inference of two ancestries based on limited population sampling. A distinct ancestry of the populations of Andaman archipelago was identified and found to be coancestral to Oceanic populations. Analysis of ancestral haplotype blocks revealed that extant mainland populations (i) admixed widely irrespective of ancestry, although admixtures between populations was not always symmetric, and (ii) this practice was rapidly replaced by endogamy about 70 generations ago, among upper castes and Indo-European speakers predominantly. This estimated time coincides with the historical period of formulation and adoption of sociocultural norms restricting intermarriage in large social strata. A similar replacement observed among tribal populations was temporally less uniform.

Pre-operative progesterone benefits operable breast cancer patients by modulating surgical stress.

PURPOSE: We have reported a survival benefit of single injection of hydroxyprogesterone prior to surgery for primary tumour in patients with node-positive operable breast cancer. Hydroxyprogesterone was meant to recapitulate the luteal phase of menstrual cycle in these women. We wanted to understand the molecular basis of action of hydroxyprogesterone on primary breast tumours in a peri-operative setting. METHODS: We performed whole transcriptome sequencing (RNA-Seq) of primary breast tumour samples collected from patients before and after hydroxyprogesterone exposure and controls. Paired breast cancer samples were obtained from patients who were given hydroxyprogesterone before surgery and a group of patients who were subjected to only surgery. RESULTS: A test of significance between the two groups revealed 207 significantly altered genes, after correction for multiple hypothesis testing. We found significantly contrasting gene expression patterns in exposed versus unexposed groups; 142 genes were up-regulated post-surgery among exposed patients, and down-regulated post-surgery among unexposed patients. Significantly enriched pathways included genes that respond to progesterone, cellular stress, nonsense-mediated decay of proteins and negative regulation of inflammatory response. These results suggest that cellular stress is modulated by hydroxyprogesterone. Network analysis revealed that UBC, a mediator of stress response, to be a major node to which many of the significantly altered genes connect. CONCLUSIONS: Our study suggests that pre-operative exposure to progesterone favourably modulates the effect of surgical stress, and this might underlie its beneficial effect when administered prior to surgery.

Y-chromosomal sequences of diverse Indian populations and the ancestry of the Andamanese.

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

Exome sequencing reveals the likely involvement of SOX10 in uveal melanoma.

PURPOSE: To identify the spectrum of somatic mutations in an Asian Indian patient with uveal melanoma (UM) without metastasis using exome sequencing. CASE REPORT: A 49-year-old man from India was diagnosed as having cilio-choroidal (uveal) melanoma (UM), without metastasis, in his right eye with the help of magnetic resonance imaging. This was later confirmed by histopathological evaluation. Two individuals from India with non-neoplastic blind eyes were recruited as controls. The affected eyes from the UM patient and the two control individuals were enucleated, and uveal tissues were collected. DNA was extracted from uveal tissue, and the matched blood sample from each of the three individuals was followed by exome sequencing. Statistical and bioinformatic analyses were done to identify somatic mutations and their putative associations with UM. Thirty-one somatic mutations (25 amino acid altering) in protein-coding (exonic) regions were detected in the UM patient. Of the amino acid-altering somatic mutations, 16 mutations were predicted to be candidate mutations relevant to UM. Somatic mutations, putatively causal for UM, were identified in GNAQ, SF3B1, and SOX10. CONCLUSIONS: Somatic mutations in GNAQ and SF3B1 genes were probable drivers of UM in the Indian patient; these were also reported earlier in some White patients. In addition, a frameshift deletion of 20 base pairs has been identified in SOX10 in the UM patient. Somatic mutations in SOX10, a transcription factor, which acts upstream of microphthalmia-associated transcription factor and synergizes with microphthalmia-associated transcription factor, was identified in some melanoma cell lines. The transcription factor SOX10 was found to have an essential role in melanocyte development and pigmentation. Our finding of the frameshift deletion (p.H387fs) in exon 4 of SOX10 in UM provides an important insight and complements earlier findings of mutations in GNAQ and SF3B1 on the genomic basis of UM.

Association of IL28B single nucleotide polymorphism rs8099917 with response to treatment in genotype 3 HCV-infected patients from India.

BACKGROUND AND AIM: IL28B gene polymorphisms have been associated with treatment-response (sustained virological response, SVR) in genotype 1 hepatitis C virus (HCV)-infected patients, but only with early phase of viral decline (rapid virological response, RVR) with genotype 3 HCV-infected patients. Association between IL28B variants and SVR in genotype 3 HCV- infected patients is unclear. Our study aimed to replicate the association of IL28Bsingle nucleotide polymorphism (SNP) rs8099917 with SVR and to validate its association with RVR in genotype 3 HCV-infected patients. METHODS: 72 patients receiving combination therapy (interferon-alpha and ribavirin) at different Indian centers were retrospectively recruited and their genotype atrs8099917 was determined. The association with RVR and SVR was tested taking in to account the variation in relevant covariates such as age, gender, baseline HCV RNA copy number and liver enzymes. RESULTS: The minor allele frequency (MAF) in the pooled samples was 0.17 at rs8099917 (G allele). 68% had TT, 29% had GT and 3% had the GG genotype. SVR was achieved in 71% of patients. A significant association ofrs8099917 with both RVR (p = 0.026) and SVR (p = 0.016) was observed with none of the covariates showing any significant association. The relapse rate was high (20%) but no association of rs8099917 was observed with relapse (p = 0.420). CONCLUSION: An IL28B SNP associates with both early phase of viral decline and sustained response in a cohort of genotype 3 HCV-infected patients from India.

Genomic determinants of antibody response to a typhoid vaccine in Indian recipients.

Typhoid is endemic in India and has high global incidence. There were large outbreaks of typhoid in India between 1990 and 2018. Available typhoid vaccines induce variable levels of protective antibodies among recipients; thus, there is variability in response to the vaccine. Interindividual genomic differences is hypothesized to be a determinant of the variability in response. We studied the antibody response of ~1000 recipients of the Vi-polysaccharide typhoid vaccine from Kolkata, India, who showed considerable variability of antibody response, i.e., anti-Vi-polysaccharide antibody level 28 days postvaccination relative to prevaccination. For each vaccinee, wholegenome genotyping was performed using the Infinium Global Screening Array (Illumina). We identified 39 SNPs that mapped to 13 chromosomal regions to be associated with antibody response to the vaccine; these included SNPs on genes LRRC28 (15q26.3), RGS7 (1q43), PTPRD (9p23), CERKL (2q31.3), DGKB (7p21.2), and TCF4 (18q21.2). Many of these loci are known to be associated with various blood cell traits, autoimmune traits and responses to other vaccines; these genes are involved in immune related functions, including TLR response, JAK-STAT signalling, phagocytosis and immune homeostasis.

Uracil as a biomarker for spatial pyrimidine metabolism in the development of gingivobuccal oral squamous cell carcinoma.

No biomarker has yet been identified that allows accurate diagnosis and prognosis of oral cancers. In this study, we investigated the presence of key metabolites in oral cancer using proton nuclear magnetic resonance (NMR) spectroscopy to identify metabolic biomarkers of gingivobuccal oral squamous cell carcinoma (GB-OSCC). NMR spectroscopy revealed that uracil was expressed in 83.09% of tumor tissues and pyrimidine metabolism was active in GB-OSCC; these results correlated well with immunohistochemistry (IHC) and RNA sequencing data. Based on further gene and protein analyses, we proposed a pathway for the production of uracil in GB-OSCC tissues. Uridinetriphosphate (UTP) is hydrolyzed to uridine diphosphate (UDP) by CD39 in the tumor microenvironment (TME). We hypothesized that UDP enters the cell with the help of the UDP-specific P2Y6 receptor for further processing by ENTPD4/5 to produce uracil. As the ATP reserves diminish, the weakened immune cells in the TME utilize pyrimidine metabolism as fuel for antitumor activity, and the same mechanism is hijacked by the tumor cells to promote their survival. Correspondingly, the differential expression of ENTPD4 and ENTPD5 in immune and tumor cells, respectively, indicatedtheir involvement in disease progression. Furthermore, higher uracil levels were detected in patients with lymph node metastasis, indicating that metastatic potential is increased in the presence of uracil. The presence of uracil and/or expression patterns of intermediate molecules in purine and pyrimidine pathways, such asCD39, CD73, and P2Y6 receptors together with ENTPD4 and ENTPD5, hold promise as biomarker(s) for oral cancer diagnosis and prognosis.

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.