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.

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.

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.

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.

A polymorphism in the CYP1B1 promoter is functionally associated with primary congenital glaucoma.

Primary congenital glaucoma (PCG) is a childhood autosomal-recessive disorder caused by developmental defects in the trabecular meshwork and anterior chamber angle. These defects cause raised intraocular pressure (IOP) that damages the optic nerve and if left untreated, results in irreversible blindness. Mutations in CYP1B1 gene at the GLC3A locus (2p21) are associated with PCG. However, there has been very limited exploration of its promoter region. We resequenced the CYP1B1 promoter in a large cohort (n = 835) that included patients with PCG (n = 301), other primary glaucomas (primary open-angle glaucoma: n = 115 and primary angle closure glaucoma: n = 100) and unaffected controls (n = 319). We functionally characterized one associated variant by luciferase reporter assay using the trabecular meshwork (TM3) cell line. We found evidence of strong (P = 6.01 × 10(-4)) association of rs2567206 (T2805C) SNP in PCG and not in other primary glaucomas. Luciferase assay indicated a ∼90% reduction in CYP1B1 promoter activity in the risk-allele (C) compared to the other allele (T). The association of the risk allele was stronger in cases harboring homozygous CYP1B1 mutations (P = 3.42 × 10(-12)). The risk haplotype 'C-C-G' in the promoter had a strong non-random association to the previously characterized risk haplotype 'C-G-G-T-A' in the coding region. The independent effect of genotype at the promoter T2805C locus (P = 0.001), and the interaction effect of genotypes at the promoter and coding region mutations loci (P = 0.001) were significant for the presenting IOP of the worst affected eye. This is the first study that unequivocally shows the functional involvement of a CYP1B1 promoter variant in PCG.