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.

MAL expression downregulation through suppressive H3K27me3 marks at the promoter in HPV16-related cervical cancers is prognostically relevant and manifested by the interplay of novel MAL antisense long noncoding RNA AC103563.8, E7 oncoprotein and EZH2.

BACKGROUND: MAL (T-lymphocyte maturation-associated protein) is highly downregulated in most cancers, including cervical cancer (CaCx), attributable to promoter hypermethylation. Long noncoding RNA genes (lncGs) play pivotal roles in CaCx pathogenesis, by interacting with human papillomavirus (HPV)-encoded oncoproteins, and epigenetically regulating coding gene expression. Hence, we attempted to decipher the impact and underlying mechanisms of MAL downregulation in HPV16-related CaCx pathogenesis, by interrogating the interactive roles of MAL antisense lncRNA AC103563.8, E7 oncoprotein and PRC2 complex protein, EZH2. RESULTS: Employing strand-specific RNA-sequencing, we confirmed the downregulated expression of MAL in association with poor overall survival of CaCx patients bearing HPV16, along with its antisense long noncoding RNA (lncRNA) AC103563.8. The strength of positive correlation between MAL and AC103563.8 was significantly high among patients compared to normal individuals. While downregulated expression of MAL was significantly associated with poor overall survival of CaCx patients bearing HPV16, AC103563.8 did not reveal any such association. We confirmed the enrichment of chromatin suppressive mark, H3K27me3 at MAL promoter, using ChIP-qPCR in HPV16-positive SiHa cells. Subsequent E7 knockdown in such cells significantly increased MAL expression, concomitant with decreased EZH2 expression and H3K27me3 marks at MAL promoter. In silico analysis revealed that both E7 and EZH2 bear the potential of interacting with AC103563.8, at the same binding domain. RNA immunoprecipitation with anti-EZH2 and anti-E7 antibodies, respectively, and subsequent quantitative PCR analysis in E7-silenced and unperturbed SiHa cells confirmed the interaction of AC103563.8 with EZH2 and E7, respectively. Apparently, AC103563.8 seems to preclude EZH2 and bind with E7, failing to block EZH2 function in patients. Thereby, enhanced EZH2 expression in the presence of E7 could potentially inactivate the MAL promoter through H3K27me3 marks, corroborating our previous results of MAL expression downregulation in patients. CONCLUSION: AC103563.8-E7-EZH2 axis, therefore, appears to crucially regulate the expression of MAL, through chromatin inactivation in HPV16-CaCx pathogenesis, warranting therapeutic strategy development.

A FACS-based novel isolation technique identifies heterogeneous CTCs in oral squamous cell carcinoma.

Purpose: Isolating circulating tumour cells (CTCs) from the blood is challenging due to their low abundance and heterogeneity. Limitations of conventional CTC detection methods highlight the need for improved strategies to detect and isolate CTCs. Currently, the Food and Drug Administration (FDA)-approved CellSearch™ and other RUO techniques are not available in India. Therefore, we wanted to develop a flexible CTC detection/isolation technique that addresses the limitation(s) of currently available techniques and is suitable for various downstream applications. Methods: We developed a novel, efficient, user-friendly CTC isolation strategy combining density gradient centrifugation and immuno-magnetic hematogenous cell depletion with fluorescence-activated cell sorting (FACS)-based positive selection using multiple CTC-specific cell-surface markers. For FACS, a stringent gating strategy was optimised to exclude debris and doublets by side scatter/forward scatter (SSC/FSC) discriminator, remove dead cells by 4',6-diamidino-2-phenylindole (DAPI) staining, and eliminate non-specific fluorescence using a "dump" channel. APC-labelled anti-CD45mAB was used to gate remaining hematogenous cells, while multiple epithelial markers (EpCAM, EGFR, and Pan-Cytokeratin) and an epithelial-mesenchymal transition (EMT) marker (Vimentin) labelled with fluorescein isothiocyanate (FITC) were used to sort cancer cells. The technique was initially developed by spiking Cal 27 cancer cells into the blood of healthy donors and then validated in 95 biopsy-proven oral squamous cell carcinoma (OSCC) patients. CTCs isolated from patients were reconfirmed by Giemsa staining, immuno-staining, and whole transcriptome amplification (WTA), followed by qRT-PCR. In vitro culture and RNA sequencing (RNA-Seq) were also performed to confirm their suitability for various downstream applications. Results: The mean detection efficiency for the Cal 27 tongue cancer cells spiked in the whole blood of healthy donors was 32.82% ± 12.71%. While ~75% of our patients (71/95) had detectable CTCs, the CTC positivity was independent of the TNM staging. The isolated potential cancer cells from OSCC patients were heterogeneous in size. They expressed different CTC-specific markers in various combinations as identified by qRT-PCR after WTA in different patients. Isolated CTCs were also found to be suitable for downstream applications like short-term CTC culture and RNA-Seq. Conclusion: We developed a sensitive, specific, flexible, and affordable CTC detection/isolation technique, which is scalable to larger patient cohorts, provides a snapshot of CTC heterogeneity, isolates live CTCs ready for downstream molecular analysis, and, most importantly, is suitable for developing countries.

Unfurling the functional association between long intergenic noncoding RNAs (lincRNAs) and HPV16-related cervical cancer pathogenesis through weighted gene co-expression network analysis of differentially expressed lincRNAs and coding genes.

Long intergenic noncoding RNAs (lincRNAs) do not overlap annotated coding genes and are located in intergenic regions, as opposed to antisense and sense-intronic lncRNAs, located in genic regions. LincRNAs influence gene expression profiles and are thereby key to disease pathogenesis. In this study, we assessed the association between lincRNAs and HPV16-positive cervical cancer (CaCx) pathogenesis using weighted gene co-expression network analysis (WGCNA) with coding genes, comparing differentially expressed lincRNA and coding genes (DElincGs and DEcGs, respectively) in HPV16-positive patients with CaCx (n = 44) with those in HPV-negative healthy individuals (n = 34). Our analysis revealed five DElincG modules, co-expressing and correlating with DEcGs. We validated a substantial number of such module-specific correlations in the HPV16-positive cancer TCGA-CESC dataset. Four such modules, displayed significant correlations with patient traits, such as HPV16 physical status, lymph node involvement, and overall survival (OS), highlighting a collaborative effect of all genes within specific modules on traits. Using the DAVID bioinformatics knowledgebase, we identified the underlying biological processes associated with these modules as cancer development and progression-associated pathways. Next, we identified the top 10 DElincGs with the highest connectivity within each functional module. Focusing on the prognostic module hub genes, downregulated CTD-2619J13.13 expression was associated with poor patient OS. This lincRNA gene interacted with 25 coding genes of its module and was associated with such biological processes as keratinization loss and keratinocyte differentiation, reflecting severe disease phenotypes. This study has translational relevance in fighting various cancers with high mortality rates in underdeveloped countries.

Correlation of severity & clinical outcomes of COVID-19 with virus variants: A prospective, multicentre hospital network study.

BACKGROUND OBJECTIVES: The clinical course of COVID-19 and its prognosis are influenced by both viral and host factors. The objectives of this study were to develop a nationwide platform to investigate the molecular epidemiology of SARS-CoV-2 (Severe acute respiratory syndrome Corona virus 2) and correlate the severity and clinical outcomes of COVID-19 with virus variants. METHODS: A nationwide, longitudinal, prospective cohort study was conducted from September 2021 to December 2022 at 14 hospitals across the country that were linked to a viral sequencing laboratory under the Indian SARS-CoV-2 Genomics Consortium. All participants (18 yr and above) who attended the hospital with a suspicion of SARS-CoV-2 infection and tested positive by the reverse transcription-PCR method were included. The participant population consisted of both hospitalized as well as outpatients. Their clinical course and outcomes were studied prospectively. Nasopharyngeal samples collected were subjected to whole genome sequencing to detect SARS-CoV-2 variants. RESULTS: Of the 4972 participants enrolled, 3397 provided samples for viral sequencing and 2723 samples were successfully sequenced. From this, the evolution of virus variants of concern including Omicron subvariants which emerged over time was observed and the same reported here. The mean age of the study participants was 41 yr and overall 49.3 per cent were female. The common symptoms were fever and cough and 32.5 per cent had comorbidities. Infection with the Delta variant evidently increased the risk of severe COVID-19 (adjusted odds ratio: 2.53, 95% confidence interval: 1.52, 4.2), while Omicron was milder independent of vaccination status. The independent risk factors for mortality were age >65 yr, presence of comorbidities and no vaccination. INTERPRETATION CONCLUSIONS: The authors believe that this is a first-of-its-kind study in the country that provides real-time data of virus evolution from a pan-India network of hospitals closely linked to the genome sequencing laboratories. The severity of COVID-19 could be correlated with virus variants with Omicron being the milder variant.

Deep sequencing reveals recurrent somatic mutations and distinct molecular subgroups in gastric cancer in Mizo population, North East India.

In India, Mizoram has the highest incidence of gastric cancer (GC) which might be associated with environmental factors such as diet, Helicobacter pylori (H.pylori) and Epstein-Barr virus (EBV) infections, and somatic genomic alterations. We performed PCR cum sequencing and fragment analysis for detection of H. pylori/EBV infection and microsatellite Instability (MSI) in GC patients (N = 68). Somatic mutations were identified by targeted and exome sequencing. We found 87% of GC patients infected with H. pylori and or EBV. Pathogenic infections were mostly mutually exclusive with only 16% of coinfection. TP53, MUC6, and ARID1A were significantly mutated. Two molecular subgroups with distinctive mutational profiles were identified: (1) patients harboring mutations in TP53 and (2) patients harboring mutations in RTK/RAS/PI3-K signaling pathway and chromatin-remodeling genes. Therefore, EBV and H. pylori infections and somatic mutations in the genes involved in RTK/RAS/PI3K signaling pathway, chromatin-remodeling, and TP53 might drive GC development and progression in Mizo patients.

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.

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.

Genetic variants associated with spontaneous preterm birth in women from India: a prospective cohort study.

Background: Despite having the highest number of preterm births globally, no genomic study on preterm birth was previously published from India or other South-Asian countries. Methods: We conducted a genome-wide association (GWA) study of spontaneous preterm birth (sPTB) on 6211 women from India. We used a novel resampling procedure to identify the associated single nucleotide polymorphisms (SNPs) followed by haplotype association analysis and imputation. Findings: We found that 512 maternal SNPs were associated with sPTB (p < 2.51e-3), of which minor allele at 19 SNPs (after Bonferroni correction) had increased genotype relative risk. Haplotypes containing six of the 19 SNPs (rs13011430, rs8179838, rs2327290, rs4798499, rs7629800, and rs13180906) were associated with sPTB (p < 9.9e-4; Bonferroni adjusted p-value <0.05). After imputation in regions around the 19 SNPs, 15 imputed SNPs were found to be associated with sPTB (Bonferroni adjusted p-value <0.05). One of these imputed SNPs, rs35760881, and three other SNPs (rs17307697, rs4308815, and rs10983507) were also reported to be associated with sPTB in women belonging to European ancestry. Moreover, we found that GG genotype at rs1152954, one of the associated SNPs, enhanced risk of sPTB and reduced telomere length. Interpretation: This is the first study from South Asia on the genome-wide identification of maternal SNPs associated with sPTB. These SNPs are known to alter the expression of genes associated with major pathways in sPTB viz. inflammation, apoptosis, cervical ripening, telomere maintenance, selenocysteine biosynthesis, myometrial contraction, and innate immunity. From a public health perspective, the trans-ethnic association of four SNPs identified in our study may help to stratify women with risk of sPTB in most populations. Funding: Department of Biotechnology (India), Grand Challenges India - All Children Thriving Program and Biotechnology Industry Research Assistance Council (BIRAC).

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.

Ocular conjunctival microbiome profiling in dry eye disease: A case control pilot study.

Purpose: Keratoconjunctivitis sicca (KCS) or dry eye disease (DED) is a multifactorial disease that results in discomfort, visual disturbance, and tear film instability with potential damage to the ocular surface. A pilot study was undertaken to determine if there were any major substantial differences in the ocular microbiome in DED patients versus healthy controls. Methods: The bacterial communities residing in the conjunctiva of patients with DED (n = 4) and healthy controls (n = 4) were assessed by 16S ribosomal RNA (rRNA) gene sequencing of the V4-V5 region. Results: The phyla Proteobacteria, Actinobacteria, Bacteroidetes, and Firmicutes were most dominant and accounted for 97% and 94.5% of all bacterial sequences in patients and controls, respectively. At the genus level, 27 bacterial genera were found with more than two-fold difference between patients and controls. Four of these - Acinetobacter, Corynebacterium, Lactobacillus, and Pseudomonas spp. - dominated the ocular microbiome of all subjects, but were proportionately lower in DED (16.5%) compared to controls (37.7%). Several bacterial genera were found to be unique in DED (34) and controls (24). Conclusion: This pilot study is an attempt to profile the ocular microbiome in patients with DED that demonstrated a higher concentration of microbial DNA compared to controls, with Firmicutes phyla dominating the bacterial population in patients with DED.

Association of exosomal miR-96-5p and miR-146a-5p with the disease severity in dengue virus infection.

Exosomes are small extracellular vesicles secreted by cells and have a major role in cell-to-cell signaling. As dengue infection progresses from a mild to a severe form of infection, the exosome's microRNA (miRNA) composition might change, which may contribute to pathogenesis. In this study, a comprehensive analysis of serum exosomal miRNAs was performed and their involvement in dengue virus-induced disease progression in an Indian cohort was assessed. Small RNA-seq showed 50 differentially expressed exosomal miRNAs that were significantly dysregulated during dengue infection. After extensive validation, miR-96-5p was found to be significantly upregulated, whereas miR-146a-5p was significantly downregulated with the progression of disease to severe form. Interestingly, a strong positive correlation was found between the expression levels of miR-96-5p and miR-146a-5p and the platelet levels of the patients. Further, study of miR-146a-5p showed that it regulates the expression of the proteins which are involved in the immune responses. These results suggest that miR-96-5p and miR-146a-5p could be used as diagnostic and prognostic markers for dengue disease progression, in addition to the already available biochemical and pathological parameters.

Emergence of a unique SARS-CoV-2 Delta sub-cluster harboring a constellation of co-appearing non-Spike mutations.

Accumulation of diverse mutations across the structural and nonstructural genes is leading to rapid evolution of SARS-CoV-2, altering its pathogenicity. We performed whole genome sequencing of 239 SARS-CoV-2 RNA samples collected from both adult and pediatric patients across eastern India (West Bengal), during the second pandemic wave in India (April-May 2021). In addition to several common spike mutations within the Delta variant, a unique constellation of eight co-appearing non-Spike mutations was identified, which revealed a high degree of positive mutual correlation. Our results also demonstrated the dynamics of SARS-CoV-2 variants among unvaccinated pediatric patients. 41.4% of our studied Delta strains harbored this signature set of eight co-appearing non-Spike mutations and phylogenetically out-clustered other Delta sub-lineages like 21J, 21A, or 21I. This is the first report from eastern India that portrayed a landscape of co-appearing mutations in the non-Spike proteins, which might have led to the evolution of a distinct Delta subcluster. Accumulation of such mutations in SARS-CoV-2 may lead to the emergence of "vaccine-evading variants." Hence, monitoring of such non-Spike mutations will be significant in the formulation of any future vaccines against those SARS-CoV-2 variants that might evade the current vaccine-induced immunity, among both the pediatric and adult populations.

Contrasting Distribution of SARS-CoV-2 Lineages across Multiple Rounds of Pandemic Waves in West Bengal, the Gateway of East and North-East States of India.

The evolution of viral variants and their impact on viral transmission have been an area of considerable importance in this pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We analyzed the viral variants in different phases of the pandemic in West Bengal, a state in India that is important geographically, and compared the variants with other states like Delhi, Maharashtra, and Karnataka, located in other regions of the country. We have identified 57 pango-lineages in 3,198 SARS-CoV-2 genomes, alteration in their distribution, as well as contrasting profiles of amino acid mutational dynamics across different waves in different states. The evolving characteristics of Delta (B.1.617.2) sublineages and alterations in hydrophobicity profiles of the viral proteins caused by these mutations were also studied. Additionally, implications of predictive host miRNA binding/unbinding to emerging spike or nucleocapsid mutations were highlighted. Our results throw considerable light on interesting aspects of the viral genomic variation and provide valuable information for improved understanding of wave-defining mutations in unfolding the pandemic. IMPORTANCE Multiple waves of infection were observed in many states in India during the coronavirus disease 2019 (COVID19) pandemic. Fine-scale evolution of major SARS-CoV-2 lineages and sublineages during four wave-window categories: Pre-Wave 1, Wave 1, Pre-Wave 2, and Wave 2 in four major states of India: Delhi (North), Maharashtra (West), Karnataka (South), and West Bengal (East) was studied using large-scale virus genome sequencing data. Our comprehensive analysis reveals contrasting molecular profiles of the wave-defining mutations and their implications in host miRNA binding/unbinding of the lineages in the major states of India.

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

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

Transcriptome analysis reveals SALL4 as a prognostic key gene in gastric adenocarcinoma.

BACKGROUND: Stomach adenocarcinoma (STAD) dominates 80-90% of gastric cancer (GC). Over the years, it has been realized that the identification of the genes responsible for gastric carcinogenesis is essential to understand the biomarker discovery. METHODS: This study aims to identify candidate genes for biomarker discovery in STAD. RNA-Seq was performed on three paired tumor-normal and one unpaired tumor samples from four GC patients and investigated for differentially expressed genes (DEGs) using DESeq2. Gene set enrichment analysis were performed. The DEGs were compared with two STAD microarray datasets available on Gene Expression Omnibus (GEO) database. Survival study (OS) were performed using KM-Plotter on the common genes between all the datasets. RESULTS: Totally, 148 DEGs were identified, wherein 55 genes were upregulated and 93 genes were downregulated with |log2foldchange| > 1 and Benjamini-Hochberg (BH) Adjusted P value < 0.01. Cell adhesion molecule (CAM) Pathway was found to be the most significant among the upregulated genes. Gastric acid secretion and mineral absorption pathways were the most significant pathways among the downregulated genes. Comparison with two GEO datasets followed by OS analysis revealed two upregulating genes, APOC1 and SALL4 with prognostic significance. CONCLUSION: Upregulation of APOC1 is associated with marginal overall survival (OS) and SALL4 over-expression was associated with the poor OS using KM-Plotter during 5 years data period. Our study suggests that SALL4 could be a promising biomarker candidate in STAD.

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.

Hypoxic ischemic brain injury with status myoclonus in a case with post-operative massive pulmonary thromboembolism: A diagnostic and therapeutic challenge.

Thromboembolic event is an important and widely regarded risk factor, influencing the outcome and overall survival in post-operative cases with second highest association with gynaecological surgeries. Here we report a 52 years old female who underwent total laparoscopic hysterectomy with bilateral salpingo-oophorectomy (TLH+BSO) presented to our emergency in unconscious and intubated state on her post-operative day-1 with cardiorespiratory failure. Based on the immediate investigations after successful resuscitation, an initial diagnosis of massive pulmonary thromboembolism was made following supportive and conservative management. Interestingly on subsequent days she started to display episodic non-rhythmic symmetric myoclonic jerks persisting for seconds to minutes; resembling 'post-hypoxic' or 'post-anoxic' status myoclonus due to hypoxic ischemic brain injury. Further neuro-electrophysiological examination revealed rhythmic well-formed delta activity on both sides symmetrically with a transformation into theta activity with forward displacement on subsequent days, bilaterally over the same areas. Later on, magnetic resonance imaging with magnetic resonance-spectroscopic findings also corroborated to hypoxic ischemic brain injury related changes over bilateral basal ganglia mostly in caudate nucleus. Post operative venous thromboembolism is a major clinical challenge and a subject of constant intensive care yet the neurological sequalae of such event is given less priority despite of its major contribution on overall prognostic and survival index.

A genomewide association study on individuals with occludable angles identifies potential risk loci for intraocular pressure.

Glaucoma is a heterogeneous group of optic neuropathies and is one of the leading causes of irreversible blindness worldwide. Primary angle closure glaucoma (PACG) is a major subtype, prevalent mostly in east and south Asia, where occludable anterior chamber angle is considered as a primary risk factor, which in turn could be responsible for high intraocular pressure (IOP) and subsequent neurodegeneration of retinal ganglion cells. Clinically, IOP is considered as a major risk factor for glaucoma and viewed as an important endophenotype to promote the disease severity. To investigate the comprehensive genomic insights, we conducted a genomewide association study (GWAS) on IOP in individuals with occludable angle (<15 degrees), thus anatomically predisposed to PACG. After performing GWAS on IOP, we identified 25 genomewide suggestive significant loci (P<1e-05, n = 240) of which, six were in complete linkage disequilibrium with the ABCA4 genic region. We successfully replicated the most significant discovery, SNPs of ABCA4 (rs2065712) in a separate cohort of 89 individuals (P =1.16e-09). We identified multiple SNPs in ABCA4 to be associated with IOP. Also, we obtained genes harbouring significantly associated SNPs, included in relevant biological pathways that could potentially be involved in IOP variation and glaucomatous neurodegeneration.