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.

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.

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.

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.

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.

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.

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.

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/

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.

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.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

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.

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.

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.

Mutation Spectrum of GNE Myopathy in the Indian Sub-Continent

BACKGROUND

GNE myopathy is an adult onset recessive genetic disorder that affects distal muscles sparing the quadriceps. GNE gene mutations have been identified in GNE myopathy patients all over the world. Homozygosity is a common feature in GNE myopathy patients worldwide.

OBJECTIVES

The major objective of this study was to investigate the mutation spectrum of GNE myopathy in India in relation to the population diversity in the country.

MATERIALS AND METHODS

We have collated GNE mutation data of Indian GNE myopathy patients from published literature and from recently identified patients. We also used data of people of Indian subcontinent from 1000 genomes database, South Asian Genome database and Strand Life Science database to determine frequency of GNE mutations in the general population.

RESULTS

A total of 67 GNE myopathy patients were studied, of whom 21% were homozygous for GNE variants, while the rest were compound heterozygous. Thirty-five different mutations in the GNE gene were recorded, of which 5 have not been reported earlier. The most frequent mutation was p.Val727Met (65%) found mainly in the heterozygous form. Another mutation, p.Ile618Thr was also common (16%) but was found mainly in patients from Rajasthan, while p.Val727Met was more widely distributed. The latter was also seen at a high frequency in general population of Indian subcontinent in all the databases. It was also present in Thailand but was absent in general population elsewhere in the world.

CONCLUSION

p.Val727Met is likely to be a founder mutation of Indian subcontinent.

Discerning the Origins of the Negritos, First Sundaland People: Deep Divergence and Archaic Admixture

Human presence in Southeast Asia dates back to at least 40,000 years ago, when the current islands formed a continental shelf called Sundaland. In the Philippine Islands, Peninsular Malaysia, and Andaman Islands, there exist indigenous groups collectively called Negritos whose ancestry can be traced to the "First Sundaland People." To understand the relationship between these Negrito groups and their demographic histories, we generated genome-wide single nucleotide polymorphism data in the Philippine Negritos and compared them with existing data from other populations. Phylogenetic tree analyses show that Negritos are basal to other East and Southeast Asians, and that they diverged from West Eurasians at least 38,000 years ago. We also found relatively high traces of Denisovan admixture in the Philippine Negritos, but not in the Malaysian and Andamanese groups, suggesting independent introgression and/or parallel losses involving Denisovan introgressed regions. Shared genetic loci between all three Negrito groups could be related to skin pigmentation, height, facial morphology and malarial resistance. These results show the unique status of Negrito groups as descended from the First Sundaland People.

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.

Detecting cognizable trends of gene expression in a time series RNA-sequencing experiment: a bootstrap approach

Study of temporal trajectory of gene expression is important. RNA sequencing is popular in genome-scale studies of transcription. Because of high expenses involved, many time-course RNA sequencing studies are challenged by inadequacy of sample sizes. This poses difficulties in conducting formal statistical tests of significance of null hypotheses. We propose a bootstrap algorithm to identify 'cognizable' 'time-trends' of gene expression. Properties of the proposed algorithm are derived using a simulation study. The proposed algorithm captured known 'time-trends' in the simulated data with a high probability of success, even when sample sizes were small (n < 10). The proposed statistical method is efficient and robust to capture 'cognizable' 'time-trends' in RNA sequencing data.

Dissecting host factors that regulate the early stages of tuberculosis infection

Incomplete understanding of mechanisms involved in the host-pathogen interactions constrains our efforts to eliminate tuberculosis. In many individuals, resulting from immune response to mycobacterial infection organised structures called granulomas are formed. To identify host responses that may control at least the early stages of infection, we employed an in vitro granuloma model. Here, human PBMCs were infected with live Mycobacterium tuberculosis in culture, and the appearance of granuloma-like structures was monitored over the next several days. Production of cytokines and chemokines in culture supernatants was monitored at various times, and the resulting temporal profiles were examined for possible correlations with either granuloma formation, or bacterial growth. While a positive association of TNF-α and IFN-γ secretion levels with extent of granuloma formation could clearly be identified, we were, however, unable to detect any statistically significant relationship between any cytokine/chemokine and bacterial growth. Examination of specific host cellular biochemical pathways revealed that either modulation of neutral lipid homeostasis through inhibition of the Gi-protein coupled receptor GPR109A, or regulation of host metabolic pathways through addition of vitamin D, provided a more effective means of controlling infection. A subsequent genotypic analysis for a select subset of genes belonging to pathways known to be significant for TB pathology revealed associations of polymorphisms with cytokine secretions and bacterial growth independently. Collectively therefore, the present study supports that key metabolic pathways of the host cell, rather than levels of relevant cytokines/chemokines might be more critical for regulating the intracellular mycobacterial load, in the context of granuloma formation.

Comparative analyses of genetic risk prediction methods reveal extreme diversity of genetic predisposition to nonalcoholic fatty liver disease (NAFLD) among ethnic populations of India

Nonalcoholic fatty liver disease (NAFLD) is a distinct pathologic condition characterized by a disease spectrum ranging from simple steatosis to steato-hepatitis, cirrhosis and hepatocellular carcinoma. Prevalence of NAFLD varies in different ethnic groups, ranging from 12% in Chinese to 45% in Hispanics. Among Indian populations, the diversity in prevalence is high, ranging from 9% in rural populations to 32% in urban populations, with geographic differences as well. Here, we wished to find out if this difference is reflected in their genetic makeup. To date, several candidate genes and a few genomewide association studies (GWAS) have been carried out, and many associations between single nucleotide polymorphisms (SNPs) and NAFLD have been observed. In this study, the risk allele frequencies (RAFs) of NAFLD-associated SNPs in 20 Indian ethnic populations (376 individuals) were analysed. We used two different measures for calculating genetic risk scores and compared their performance. The correlation of additive risk scores of NAFLD for three Hapmap populations with their weighted mean prevalence was found to be high (R(2) = 0.93). Later we used this method to compare NAFLD risk among ethnic Indian populations. Based on our observation, the Indian caste populations have high risk scores compared to Caucasians, who are often used as surrogate and similar to Indian caste population in disease gene association studies, and is significantly higher than the Indian tribal populations.

Genomics of immune response to typhoid and cholera vaccines

Considerable variation in antibody response (AR) was observed among recipients of an injectable typhoid vaccine and an oral cholera vaccine. We sought to find whether polymorphisms in genes of the immune system, both innate and adaptive, were associated with the observed variation in response. For both vaccines, we were able to discover and validate several polymorphisms that were significantly associated with immune response. For the typhoid vaccines, these polymorphisms were on genes that belonged to pathways of polysaccharide recognition, signal transduction, inhibition of T-cell proliferation, pro-inflammatory signalling and eventual production of antimicrobial peptides. For the cholera vaccine, the pathways included epithelial barrier integrity, intestinal homeostasis and leucocyte recruitment. Even though traditional wisdom indicates that both vaccines should act as T-cell-independent antigens, our findings reveal that the vaccines induce AR using different pathways.

First International Workshopson Provocative Questions (PQ) in Cancer Research, October-November 2014, New Delhi, Bengaluru, and Thiruvananthapuram, India

In 2011, the National Cancer Institute (NCI, USA) introduced the Provocative Questions (PQ) Initiative, a new approach allowing active researchers to define major unsolved or neglected problems in oncology unaddressed by existing funding. Last year, the U.S. NCI teamed up with the Indian Department of Biotechnology (DBT) to pilot the PQ approach in three cities in India. Workshop outcomes includedthe generation of fundable "PQs" (perplexing questions understudied by the international scientific community), as well as the identification of several non-PQ projects and research-related issues of importance to DBT and other Indian funding groups. The workshops clearly indicated the need to expand beyond crafting "PQs" when considering the best areas for research funding in international settings. Nonetheless, the first set of PQ workshops provided a forum to discuss key issues regarding cancer research in India, including the paucity of cancer research funding, and the lack of relevant human resource training and technology sharing platforms. Continued open debate between researchers, funders and policymakers will be essential to effectively strengthen the cancer research portfolio in India.

Health & Demographic Surveillance System Profile: The Birbhum population project (Birbhum HDSS)

The Birbhum HDSS was established in 2008 and covers 351 villages in four administrative blocks in rural areas of Birbhum district of West Bengal, India. The project currently follows 54 585 individuals living in 12557 households. The population being followed up is economically underprivileged and socially marginalized. The HDSS, a prospective longitudinal cohort study, has been designed to study changes in population demographic, health and healthcare utilization. In addition to collecting data on vital statistics and antenatal and postnatal tracking, verbal autopsies are being performed. Moreover, periodic surveys capturing socio-demographic and economic conditions have been conducted twice. Data on nutritional status (children as well as adults), non-communicable diseases, smoking etc. have also been collected in special surveys. Currently, intervention studies on anaemia, undernutrition and common preschool childhood morbidities through behavioural changes are under way. For access to the data, a researcher needs to send a request to the Data Manager [suri.shds@gmail.com]. Data are shared in common formats like comma-separated files (csv) or Microsoft Excel (xlsx) or Microsoft Access Database (mdb).The HDSS will soon upgrade its data management system to a more integrated platform, coordinated and guided by INDEPTH data sharing policy.

Somatic mutations in arachidonic acid metabolism pathway genes enhance oral cancer post-treatment disease-free survival

The arachidonic acid metabolism (AAM) pathway promotes tumour progression. Chemical inhibitors of AAM pathway prolong post-treatment survival of cancer patients. Here we test whether non-synonymous somatic mutations in genes of this pathway, acting as natural inhibitors, increase post-treatment survival. We identify loss-of-function somatic mutations in 15 (18%) of 84 treatment-naïve oral cancer patients by whole-exome sequencing, which we map to genes of AAM pathway. Patients (n = 53) who survived ≥ 12 months after surgery without recurrence have significantly (P = 0.007) higher proportion (26% versus 3%) of mutations than those who did not (n = 31). Patients with mutations have a significantly (P = 0.003) longer median disease-free survival (24 months) than those without (13 months). Compared with the presence of a mutation, absence of any mutation increases the hazard ratio for death (11.3) significantly (P = 0.018). The inferences are strengthened when we pool our data with The Cancer Genome Atlas (TCGA) data. In patients with AAM pathway mutations, some downstream pathways, such as the PI3K-Akt pathway, are downregulated.

Global footprints of purifying selection on Toll-like receptor genes primarily associated with response to bacterial infections in humans

Toll-like receptors (TLRs) are directly involved in host–pathogen interactions. Polymorphisms in these genes are associated with susceptibility to infectious diseases. To understand the influence of environment and pathogen diversity on the evolution of TLR genes, we have undertaken a large-scale population-genetic study. Our study included two hunter–gatherer tribal populations and one urbanized nontribal population from India with distinct ethnicities (n = 266) and 14 populations inhabiting four different continents (n = 1,092). From the data on DNA sequences of cell-surface TLR genes, we observed an excess of rare variants and a large number of low frequency haplotypes in each gene. Nonsynonymous changes were few in every population and the commonly used statistical tests for detecting natural selection provided evidence of purifying selection. The evidence of purifying selection acting on the cell-surface TLRs of the innate immune system is not consistent with Haldane’s theory of coevolution of immunity genes, at least of innate immunity genes, with pathogens. Our study provides evidence that genes of the cell-surface TLRs, that is, TLR2 and TLR4, have been so optimized to defend the host against microbial infections that new mutations in these genes are quickly eliminated.

A genomic view of the peopling and population structure of India

Recent advances in molecular and statistical genetics have enabled the reconstruction of human history by studying living humans. The ability to sequence and study DNA by calibrating the rate of accumulation of changes with evolutionary time has enabled robust inferences about how humans have evolved. These data indicate that modern humans evolved in Africa about 150,000 years ago and, consistent with paleontological evidence, migrated out of Africa. And through a series of settlements, demographic expansions, and further migrations, they populated the entire world. One of the first waves of migration from Africa was into India. Subsequent, more recent, waves of migration from other parts of the world have resulted in India being a genetic melting pot. Contemporary India has a rich tapestry of cultures and ecologies. There are about 400 tribal groups and more than 4000 groups of castes and subcastes, speaking dialects of 22 recognized languages belonging to four major language families. The contemporary social structure of Indian populations is characterized by endogamy with different degrees of porosity. The social structure, possibly coupled with large ecological heterogeneity, has resulted in considerable genetic diversity and local genetic differences within India. In this essay, we provide genetic evidence of how India may have been peopled, the nature and extent of its genetic diversity, and genetic structure among the extant populations of India.

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 correlates of variability in immune response to an oral cholera vaccine

Cholera is endemic to many countries. Recent major outbreaks of cholera have prompted World Health Organization to recommend oral cholera vaccination as a public-health strategy. Variation in percentage of seroconversion upon cholera vaccination has been recorded across populations. Vaccine-induced responses are influenced by host genetic differences. We have investigated association between single-nucleotide polymorphic (SNP) loci in and around 296 immunologically relevant genes and total anti-lipopolysaccharide (LPS) antibody response to a killed whole-cell vaccine, comprising LPS from multiple strains of Vibrio cholerae. Titers derived from standard vibriocidal assays were also analyzed to gain further insights on validated SNP associations. Vaccination was administered to 1000 individuals drawn from India. Data on two independent random subsets, each comprising ∼500 vaccinees, were used for discovery of genomic associations and validation, respectively. Significant associations of four SNPs and haplotypes in three genes (MARCO, TNFAIP3 and CXCL12) with AR were discovered and validated, of which two in TNFAIP3 and CXCL12 were also significantly associated with immunity (fourfold increase in vibriocidal titers). CXCL12 is a neutrophil and lymphocyte chemoattractant that is upregulated in response to V. cholerae infection. LPS in the vaccine possibly provides signals that mimic those of the live bacterium. TNFAIP3 promotes intestinal epithelial barrier integrity and provides tight junction protein regulation; possible requirements for adequate response to the vaccine. LPS is a potent activator of innate immune responses and a ligand of MARCO. Variants in this gene have been found to be associated with LPS response, but not with high vibriocidal titer level.

Genetic variants of TNFα, IL10, IL1β, CTLA4 and TGFβ1 modulate the indices of alcohol-induced liver injury in East Indian population

Alcohol induced liver disease or alcoholic liver disease (ALD), a complex trait, encompasses a gamut of pathophysiological alterations in the liver due to continuous exposure to a toxic amount of alcohol (more than 80 g per day). Of all chronic heavy drinkers, only 15-20% develops hepatitis or cirrhosis concomitantly or in succession. Several studies revealed that inter-individual as well as inter-ethnic genetic variation is one of the major factors that predispose to ALD. The role of genetic factors in ALD has long been sought for in ethnically distinct population groups. ALD is fast emerging as an important cause of chronic liver disease in India; even in populations such as "Bengalis" who were "culturally immune" earlier. While the genetic involvement in the pathogenesis of ALD is being sought for in different races, the complex pathophysiology of ALD as well as the knowledge of population level diversity of the relevant alcohol metabolizing and inflammatory pathways mandates the need for well designed studies of genetic factors in ethnically distinct population groups. An array of cytokines plays a critical role as mediators of injury, inflammation, fibrosis and cirrhosis in ALD. We, therefore, studied the association of polymorphisms in five relevant cytokine genes with "clinically significant" ALD in an ethnic "Bengali" population in Eastern India. Compared with "alcoholic" controls without liver disease (n=110), TNFα -238AA genotype, IL1β -511CC genotype, TGFβ1 -509CC genotype and IL10 -592AA genotype were significantly overrepresented in ALD patients (n=181; OR=2.4 and 95% CI 1.2-5.5, P(genotype)=0.042, P(allelic)=0.008; OR=2.7 and 95% CI 1.2-5.9, P(genotype)=0.018, P(allelic)=0.023; OR=4.7 and 95% CI 1.7-13.1, P(genotype)=0.003, P(allelic)=0.014; and OR=2.2 and 95% CI 1.1-4.8, P(genotype)=0.04, P(allelic)=0.039 respectively). Moreover a cumulative genetic risk analysis revealed a significant trend for developing ALD with an increase in the number of risk alleles on IL10 and TGFβ1 loci among alcoholics. The risk genotype of IL1β and TGFβ1 also influences the total bilirubin, albumin and alanine aminotransferase levels among alcoholic "Bengalis". The present study is the first case-control study from Eastern India that comprehensively identified polymorphic markers in TNFα, IL10, IL1β and TGFβ1 genes to be associated with ALD in the Bengali population, accentuating the significance of genetic factors in clinical expressions of ALD.

Association of His1085His INSR gene polymorphism with type 2 diabetes in South Indians

BACKGROUND AND OBJECTIVE

The INSR gene, which encodes the insulin receptor, is a candidate gene for type 2 diabetes (T2D). The objective of the present study was to sequence some of the crucial exons in the INSR gene such as exon 2, which encodes the insulin-binding domain of the INSR protein, and exons 17-21, which encode the protein tyrosine kinase domain for mutations/polymorphisms, and to study their association with T2D in the South Indian population.

SUBJECTS AND METHODS

The INSR gene was sequenced in 25 normal glucose-tolerant (NGT) and 25 T2D subjects, and the variant found was genotyped by polymerase chain reaction-restriction fragment length polymorphism in 1,016 NGT and 1,010 T2D subjects, randomly selected from the Chennai Urban Rural Epidemiology Study.

RESULTS

Only one previously reported polymorphism, His1085His [rs1799817, (C→T)], in exon 17 was detected by sequencing. The frequency of the "T" allele of the His1085His polymorphism was significantly lower in the T2D subjects (31%) compared with the NGT subjects (35%) and showed significant protection against diabetes (odds ratio 0.85, 95% confidence interval 0.75-0.97, P=0.019). Regression analysis according to a recessive model taking the CC+CT genotype as the reference showed that the TT genotype was protective against diabetes (odds ratio 0.71, 95% confidence interval 0.50-0.99, P=0.048). Adjusting this P value by the number of competing models (three) using Bonferroni's correction, we found that the association finding did not remain significant.

CONCLUSIONS

The "T" allele of the His1085His polymorphism in the INSR gene shows significant protection against diabetes. This study gains importance because there are no data available to date on the role of INSR variants in T2D in the Indian population.