Genome-wide RNA-seq, DNA methylation and small RNA-seq analysis unraveled complex gene regulatory networks in psoriasis pathogenesis.

Psoriasis is a complex inflammatory skin disease characterized by reversible albeit relapsing red scaly plaques in the skin of a patient. In addition to the genetic predisposition, involvement of epigenetic and non-coding RNAs have also been liked with the disease. Nevertheless, any comprehensive study involving transcriptomic, small-RNA and DNA methylation at the genomic level from same patients is lacking. To investigate the complex regulation of molecular pathways in psoriasis, we carried out multi-omics integrative analysis of RNA-sequencing, small RNA-sequencing and DNA methylation profiling from the psoriatic and adjacent normal skin tissues. Our multi-omics analysis identified the genes and biological processes regulated either independently or in combination by DNA methylation and microRNAs. We identified miRNAs that specifically regulated keratinocyte hyper-proliferation, and cell cycle progression and checkpoint signaling in psoriasis. On contrary, DNA methylation was found to be more predominant in regulating immune and inflammatory responses, another causative factor in psoriasis pathogenesis. Many characteristic pathways in psoriasis e.g., Th17 cell differentiation and JAK-STAT signaling, were found to be regulated by both miRNAs and DNA methylation. We carried out functional characterization of a downregulated miRNA hsa-let-7c-5p, predicted to target upregulated genes in psoriasis involved in cell cycle processes, Th17 cell differentiation and JAK-STAT signaling pathways. Overexpression of hsa-let-7c-5p in keratinocytes caused the downregulation of its target genes, resulting in reduced cell proliferation and migration rates, demonstrating potential of miRNAs in regulating psoriasis pathogenesis. In conclusion, our findings identified distinct and shared gene-networks regulated by DNA methylation and miRNAs of a complex disease with reversible phenotype.

Universal penalized regression (Elastic-net) model with differentially methylated promoters for oral cancer prediction.

BACKGROUND: DNA methylation showed notable potential to act as a diagnostic marker in many cancers. Many studies proposed DNA methylation biomarker in OSCC detection, while most of these studies are limited to specific cohorts or geographical location. However, the generalizability of DNA methylation as a diagnostic marker in oral cancer across different geographical locations is yet to be investigated. METHODS: We used genome-wide methylation data from 384 oral cavity cancer and normal tissues from TCGA HNSCC and eastern India. The common differentially methylated CpGs in these two cohorts were used to develop an Elastic-net model that can be used for the diagnosis of OSCC. The model was validated using 812 HNSCC and normal samples from different anatomical sites of oral cavity from seven countries. Droplet Digital PCR of methyl-sensitive restriction enzyme digested DNA (ddMSRE) was used for quantification of methylation and validation of the model with 22 OSCC and 22 contralateral normal samples. Additionally, pyrosequencing was used to validate the model using 46 OSCC and 25 adjacent normal and 21 contralateral normal tissue samples. RESULTS: With ddMSRE, our model showed 91% sensitivity, 100% specificity, and 95% accuracy in classifying OSCC from the contralateral normal tissues. Validation of the model with pyrosequencing also showed 96% sensitivity, 91% specificity, and 93% accuracy for classifying the OSCC from contralateral normal samples, while in case of adjacent normal samples we found similar sensitivity but with 20% specificity, suggesting the presence of early disease methylation signature at the adjacent normal samples. Methylation array data of HNSCC and normal tissues from different geographical locations and different anatomical sites showed comparable sensitivity, specificity, and accuracy in detecting oral cavity cancer with across. Similar results were also observed for different stages of oral cavity cancer. CONCLUSIONS: Our model identified crucial genomic regions affected by DNA methylation in OSCC and showed similar accuracy in detecting oral cancer across different geographical locations. The high specificity of this model in classifying contralateral normal samples from the oral cancer compared to the adjacent normal samples suggested applicability of the model in early detection.

Identifying the genetic associations among the psoriasis patients in eastern India.

Psoriasis is a multifactorial genetic disorder manifested by hyperproliferation and abnormal differentiation of epidermal keratinocytes, along with the infiltration of inflammatory cells into the skin. Although ~80 genetic susceptibility variants were reported in psoriasis, many loci showed population-specific associations, warranting the need for more population-specific association studies in psoriasis. We determined the association of forty single nucleotide polymorphisms (SNPs) among 2136 psoriasis patients and normal individuals from eastern India. We investigated the expression of corresponding genes and evaluated the protein structure stability for the genes with susceptible coding variants. We found fifteen SNPs significantly associated with psoriasis, while additional three SNPs showed significant association when we classified the patients based on the presence of HLA-Cw6 allele. Epistatic interaction between HLA-Cw6 and other associated loci showed significant association with the SNPs at PSORS1 region, along with other five SNPs outside PSORS1. Three genes showed significant differential expression in psoriatic tissues compared to the adjacent normal skin tissues but were not differential when classified the patients based on their genotypes. SNP rs495337 at SPATA2 (Spermatogenesis Associated 2) showed a 1.2-fold increased risk among the HLA-Cw6 patients compared to combined samples. We found significant downregulation of SPATA2 among the patients with risk genotypes and HLA-Cw6 allele compared to the non-risk genotypes. Protein structure stability analysis showed reduced structural stability for all the mutant residues caused by the associated coding variants. Our study evaluated the genetic associations of psoriasis-susceptible variants in India and evaluated the possible functional significance of these associated variants in psoriasis.

Serum miRNA profiling identified miRNAs associated with disease severity in psoriasis.

Psoriasis vulgaris is a chronic, autoimmune skin disease involving a complex interplay of epidermal keratinocytes, dermal fibroblast and infiltrating immune cells. Differential expressions of miRNAs are observed in psoriasis and the deregulated miRNAs are sometimes associated with disease severity. This study aims to identify miRNAs altered in the serum of psoriasis patients that are associated with the Psoriasis Area and Severity Index (PASI). In order to assess miRNA levels in the serum of psoriasis patients, we selected 24 differentially expressed miRNAs in the psoriatic skin are possibly derived from the skin and immune cells, as well as five miRNAs that are enriched in other tissues. We identified 16 miRNAs that exhibited significantly (p < 0.05) altered levels in the serum of psoriasis patients compared to healthy individuals. Among these, 13 miRNAs showed similar expression pattern in the serum of psoriasis patients as also observed in the psoriatic skin tissues. Ten miRNAs showed an accuracy of greater than 75% in classifying the psoriasis patients from healthy individuals. Further analysis of differential miRNA levels between the low PASI group and the high PASI group identified three miRNAs (miR-147b, miR-3614-5p, and miR-125a-5p) with significantly altered levels between the low severity and the high severity psoriasis patients. Our systematic investigation of skin and immune cell-derived miRNAs in the serum of psoriasis patients revealed alteration in miRNA levels to be associated with disease severity, which may help in monitoring the disease progression and therapeutic response.

Emerging roles of non-coding RNAs in psoriasis pathogenesis.

Psoriasis is a complex genetic skin disorder typically manifested by red, scaly, and itchy plaques most commonly over the scalp, trunk, elbows, and knees. Histopathological features include thickening of the epidermal layer due to hyper-proliferation and abnormal differentiation of epidermal keratinocytes along with infiltration of immune cells in the psoriatic skin. It is a chronic inflammatory relapsing disease, and there is currently no permanent cure for psoriasis. Proper medications can reduce the severity of the disease and improve the quality of life of the patients. While the genetic components of psoriasis pathogenesis are well explored, the full understanding of its epigenetic component remains elusive. Non-coding RNAs (ncRNAs) are documented to regulate various epigenetic processes that lead to the pathogenesis of different diseases including psoriasis. In this review, we have discussed the molecular interplay of different ncRNAs in psoriasis pathogenesis. The roles of microRNAs (miRNAs) in psoriasis are pretty well studied, whereas the roles of long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) are emerging. This review provides ideas covering some of the latest findings of different modes of functions played by those different ncRNAs documented in the literature. As an ever-evolving topic, some works are still ongoing as well as there are several fields that need rigorous scientific ventures. We have proposed the areas which claim more explorations to better understand the roles played by the ncRNAs in psoriasis pathogenesis.

Parkinson's disease-associated 18 bp promoter variant of DJ-1 alters REST binding and regulates its expression.

Parkinson's disease (PD) is a progressive neurodegenerative disorder with a complex etiology. Presence of autosomal mutations in PARK7/DJ-1 gene has been associated with early-onset PD. Growing evidence has suggested that DJ-1 acts as a putative sensor of oxidative stress. Reduced levels of DJ-1 protein have been reported in the cerebrospinal fluid of sporadic PD patients. Several case-control association studies have identified DJ-1 g.168_185del (rs200968609) variants conferring susceptibility towards PD pathogenesis. Similarly, among the PD patients in eastern India, the deletion allele (g.168_185) of this DJ-1 promoter polymorphism was found to be associated with PD. Hence, we aimed to find out the functional contribution of this promoter variant of DJ-1 in PD pathogenesis. The expression of DJ-1 was observed to be significantly reduced in the presence of both deletion and duplication sequences as identified from the luciferase promoter activity assay. The transcription factor binding prediction tool identified DJ-1 promoter 18 bp insertion polymorphism as the only binding partner of REST (RE1 Silencing Transcription Factor). Transient Chromatin Immuno-precipitation (ChIP) assay further confirmed this prediction. Previous reports have highlighted the role of REST in regulating the expression of stress-responsive genes. Our study has identified the functional involvement of DJ-1 promoter variants and REST-mediated regulation of DJ-1 expression in PD pathogenesis.

MiRNAs as non-invasive biomarkers in the serum of Oral Squamous Cell Carcinoma (OSCC) and Oral Potentially Malignant Disorder (OPMD) patients.

OBJECTIVE: Cell-free microRNAs have shown differential levels in the serum of individuals under disease conditions suggesting its potential to act as biomarkers. A population specific miRNA signature in oral cancer is reported in different studies. We aim to identify a set of serum specific miRNAs that may differentiate oral cancer, oral pre-malignant conditions from the healthy individuals. DESIGN: We investigated the levels of 24 miRNAs in the serum of 47 Oral squamous cell carcinoma (OSCC) patients, 20 patients with Oral potentially malignant disorders (OPMD) and 42 healthy controls from Eastern India. Small RNAs were isolated from serum samples followed by cDNA synthesis. Levels of miRNAs were determined using qRT-PCR. The sources of serum specific miRNAs were evaluated using GTEx-RNAseq and TCGA-HNSCC database. RESULTS: Five miRNAs, miR-483-5p, miR-31-5p, Let-7b-5p, miR-486-5p and miR-30e-5p showed significant elevation in OSCC patients. An Elastic-Net model with 4 miRNAs classified OSCC from healthy controls with 80 % sensitivity, 64.3 % specificity, and 72.4 % accuracy. Mir-483-5p and miR-31-5p was significantly overexpressed in OSCC tissues as well as significantly higher in the serum of Leukoplakia and Verrucous carcinoma patients suggesting their potential as early disease markers. MiR-483-5p showed a consistent elevated level in the serum/plasma of oral cancer patients across different population and was found to be tumour specific while, the rest of the miRNAs showed variable results across different studies. CONCLUSIONS: Our study suggested that the serum miRNAs in oral cancer and pre-malignant disorder conditions can be used as a non-invasive marker for screening of these oral conditions.

Effect of deletions in the α-globin gene on the phenotype severity of ß-thalassemia.

Thalassemia is the most common inherited hemoglobinopathy worldwide. Variation of clinical symptoms in this hemoglobinopathy entails differences in disease-onset and transfusion requirements. The aim of this study was to investigate the role of α-globin gene deletions in modulating the clinical heterogeneity of ß-thalassemia (ß-thal) syndromes. A total number 270 ß-thal subjects were enrolled. Hematological parameters were recorded. ß-Globin mutations were determined by amplified refractory mutation system-polymerase chain reaction (ARMS-PCR), gap-PCR and Sanger sequencing. α-Globin gene deletions were determined by multiplex PCR. Out of 270 ß-thal subjects, 19 carried ß+/ß+, 74 had ß0/ß0 and 177 had the ß0/ß+ genotype. When we determined the severity of the different ß-thal subjects in coinherited with the α gene deletion, it was revealed that, 84.2% ß+/ß+ subjects carried a non severe phenotype and did not have an α gene deletion. Of the ß0/ß0 individuals, 95.9% presented a severe phenotype, irrespective of α-globin gene deletions. In cases with the ß0/ß+ genotype, 19.2% subjects also carried a deletion on the α gene. Of these, 61.8% presented a non severe phenotype and 38.2% were severely affected. Only in the ß0/ß+ category did α gene deletions make a significant contribution (p < 0.001) toward alleviation of clinical severity. Therefore, it can be stated that α-globin gene deletions play a role in ameliorating the phenotype in patients with a ß+/ß0 genotype.

Performance assessment of genomic island prediction tools with an improved version of Design-Island.

Genomic Islands (GIs) play an important role in the evolution and adaptation of prokaryotes. The origin and extent of ecological diversity of prokaryotes can be analyzed by comparing GIs across closely or distantly related prokaryotes. Understanding the importance of GI and to study the bacterial evolution, several GI prediction tools have been generated. An unsupervised method, Design-Island, was developed to identify GIs using Monte-Carlo statistical test on randomly selected segments of a chromosome. Here, in the present study Design-Island was modified with the incorporation of majority voting, multiple hypothesis testing correction. The performance of the modified version, Design-Island-II was tested and compared with the existing GI prediction tools. The performance assessment and benchmarking of the GI prediction tools require experimentally validated dataset, which is lacking. So, different datasets, generated or taken from literature were utilized to compare the sensitivity (SN), specificity (SP), precision (PPV) and accuracy (AC) of Design-Island-II. It showed substantial enhancement in term of SN, SP, PPV and AC, and significantly reduced the computation time of the algorithm. The performance of Design-Island-II has also been compared with several GI prediction tools using curated dataset of putative horizontally transferred genes. Design-Island-II showed the highest sensitivity and F1 score, comparable specificity, precision and accuracy in comparison to the other available methods. IslandViewer4 and Islander outperformed all the available methods in terms of AC and PPV respectively. Our study suggested Design-Island-II, IslandViewer4 and GIHunter among the top performing GI prediction tools considering both sensitivity and specificity of the methods.

Comparative genomics analysis of statistically significant genomic islands of Helicobacter pylori strains for better understanding the disease prognosis.

Bacterial virulence factors are often located in their genomic islands (GIs). Helicobacter pylori, a highly diverse organism is reported to be associated with several gastrointestinal diseases like, gastritis, gastric cancer (GC), peptic ulcer, duodenal ulcer (DU) etc. A novel similarity score (Sm)-based comparative analysis with GIs of 50 H. pylori strains revealed clear idea of the various factors which promote disease progression. Two putative pathogenic GIs in some of the H. pylori strains were identified. One GI, having a putative labile enterotoxin and other dynamin-like proteins (DLPs), is predicted to increase the release of toxin by membrane vesicular formation. Another island contains a virulence-associated protein D (vapD) which is a component of a type-II toxin-antitoxin system (TAs), leads to enhance the severity of the H. pylori infection. Besides the well-known virulence factors like Cytotoxin-associated gene A (CagA) and vacA, several GIs have been identified which showed to have direct or indirect impact on H. pylori clinical outcomes. One such GI, containing lipopolysaccharide (LPS) biosynthesis genes was revealed to be directly connected with disease development by inhibiting the immune response. Another collagenase-containing GI worsens ulcers by slowing down the healing process. GI consisted of fliD operon was found to be connected to flagellar assembly and biofilm production. By residing in biofilms, bacteria can avoid antibiotic therapy, resulting in chronic infection. Along with well-studied CagA and vacuolating toxin A (vacA) virulent genes, it is equally important to study these identified virulence factors for better understanding H. pylori-induced disease prognosis.

MICaps: Multi-instance capsule network for machine inspection of Munro's microabscess.

Munro's Microabscess (MM) is the diagnostic hallmark of psoriasis. Neutrophil detection in the Stratum Corneum (SC) of the skin epidermis is an integral part of MM detection in skin biopsy. The microscopic inspection of skin biopsy is a tedious task and staining variations in skin histopathology often hinder human performance to differentiate neutrophils from skin keratinocytes. Motivated from this, we propose a computational framework that can assist human experts and reduce potential errors in diagnosis. The framework first segments the SC layer, and multiple patches are sampled from the segmented regions which are classified to detect neutrophils. Both UNet and CapsNet are used for segmentation and classification. Experiments show that of the two choices, CapsNet, owing to its robustness towards better hierarchical object representation and localisation ability, appears as a better candidate for both segmentation and classification tasks and hence, we termed our framework as MICaps. The training algorithm explores both minimisation of Dice Loss and Focal Loss and makes a comparative study between the two. The proposed framework is validated with our in-house dataset consisting of 290 skin biopsy images. Two different experiments are considered. Under the first protocol, only 3-fold cross-validation is done to directly compare the current results with the state-of-the-art ones. Next, the performance of the system on a held-out data set is reported. The experimental results show that MICaps improves the state-of-the-art diagnosis performance by 3.27% (maximum) and reduces the number of model parameters by 50%.

A combination of circulating microRNA-375-3p and chemokines CCL11, CXCL12, and G-CSF differentiate Crohn's disease and intestinal tuberculosis.

Differentiation of Crohn's disease (CD) from intestinal tuberculosis (ITB) is a big challenge to gastroenterologists because of their indistinguishable features and insensitive diagnostic tools. A non-invasive biomarker is urgently required to distinguish ITB/CD patients particularly in India, a TB endemic region, where CD frequency is increasing rapidly due to urbanization. Among the three differentially expressed miRNAs obtained from small RNA transcriptomic profiling of ileocaecal/terminal ileal tissue of ITB/CD patients (n = 3), only two down-regulated miRNAs, miR-31-5p, and miR-215-5p showed comparable data in qRT-PCR. Out of which, only miR-215-5p was detectable in the patient's plasma, but there was no significant difference in expression between ITB/CD. On the other hand, miR-375-3p, the pulmonary TB specific marker was found in higher amount in the plasma of ITB patients than CD while reverse expression was observed in the ileocaecal/terminal ileal tissues of the same patients. Next, using Bioplex pro-human cytokine 48-plex screening panel, only three chemokines, Eotaxin-1/CCL11, SDF-1α/CXCL12, and G-CSF have noted significantly different levels in the serum of ITB/CD patients. ROC analysis has revealed that compared to a single molecule, a combination of miR-375-3p + Eotaxin-1/CCL11 + SDF-1α /CXCL12 + G-CSF showed a better AUC of 0.83, 95% CI (0.69-0.96) with 100% specificity and positive predictive value while sensitivity, negative predictive value, and accuracy were 56%, 69%, and 78% respectively in distinguishing ITB from CD. This study suggests that a combination of plasma markers shows better potential in differentiating ITB from CD than a single marker and this panel of markers may be used for clinical management of ITB/CD patients.

Co-Regulation of Protein Coding Genes by Transcription Factor and Long Non-Coding RNA in SARS-CoV-2 Infected Cells: An In Silico Analysis.

Altered expression of protein coding gene (PCG) and long non-coding RNA (lncRNA) have been identified in SARS-CoV-2 infected cells and tissues from COVID-19 patients. The functional role and mechanism (s) of transcriptional regulation of deregulated genes in COVID-19 remain largely unknown. In the present communication, reanalyzing publicly available gene expression data, we observed that 66 lncRNA and 5491 PCG were deregulated in more than one experimental condition. Combining our earlier published results and using different publicly available resources, it was observed that 72 deregulated lncRNA interacted with 3228 genes/proteins. Many targets of deregulated lncRNA could also interact with SARS-CoV-2 coded proteins, modulated by IFN treatment and identified in CRISPR screening to modulate SARS-CoV-2 infection. The majority of the deregulated lncRNA and PCG were targets of at least one of the transcription factors (TFs), interferon responsive factors (IRFs), signal transducer, and activator of transcription (STATs), NFκB, MYC, and RELA/p65. Deregulated 1069 PCG was joint targets of lncRNA and TF. These joint targets are significantly enriched with pathways relevant for SARS-CoV-2 infection indicating that joint regulation of PCG could be one of the mechanisms for deregulation. Over all this manuscript showed possible involvement of lncRNA and mechanisms of deregulation of PCG in the pathogenesis of COVID-19.

Temporal variations of country-specific mutational profile of SARS-CoV-2: effect on vaccine efficacy.

Aim: In order to curb the transmission of SARS-CoV-2, nation-wide travel restrictions at different levels were implemented in different countries. Country-specific mutational profile may exist and have an impact on vaccine efficacy. Materials & methods: We identified nonsynonymous mutations in approximately 215,000 SARS-CoV-2 sequences during the 1st year of the pandemic in 35 countries. Mutational profiles on a bimonthly basis were traced over time. We also examined the mutations that overlapped with the spike protein vaccine epitopes. Results: Several new mutations emerged over time and were dominating in specific countries. Many nonsynonymous mutations were within multiple spike protein epitopes that might impact the vaccine efficacy. Conclusion: Our study advocates requirement of active monitoring of country-specific mutations and vaccine efficacies in respective countries.

Functional Mapping of Genetic Interactions between HLA-Cw6 and LCE3A in Psoriasis.

Functional studies to delineate the molecular mechanisms of causal genetic variants are the main focus in the post-GWAS era. Previous GWASs have identified >50 susceptibility loci associated with psoriasis. Functional understanding of the biology underlying the disease risk of most of these associated loci is unclear. In this study, we identified a regulatory SNP at the putative enhancer of the LCE3A gene within the epidermal differentiation complex that showed epistatic interaction with HLA-Cw6. The variant allele disrupted signal transducer and activator of transcription 3 binding to the region, thereby regulating the expression of the downstream LCE3A gene. Electrophoretic mobility shift and pulldown assay confirmed the preferential binding of signal transducer and activator of transcription 3 to the DNA with a wild-type allele compared with the DNA with a variant allele. The reporter assay further validated the IL-6‒stimulated phosphorylated signal transducer and activator of transcription 3‒mediated LCE3A activation in the presence of the wild-type allele. Interestingly, the presence of the HLA-Cw6 allele leads to IL-6‒mediated phosphorylation of signal transducer and activator of transcription 3, followed by its nuclear localization in the epidermal keratinocytes of psoriatic skin, suggesting indirect interaction of the HLA-Cw6 allele and a regulatory SNP upstream of the LCE3A gene. This study reflects an interesting approach to dissecting the molecular mechanism underlying the genetic interaction observed between HLA-Cw6 and LCE3A in psoriasis pathogenesis.

Dysbiosis of Oral Microbiota During Oral Squamous Cell Carcinoma Development.

Infection with specific pathogens and alterations in tissue commensal microbial composition are intricately associated with the development of many human cancers. Likewise, dysbiosis of oral microbiome was also shown to play critical role in the initiation as well as progression of oral cancer. However, there are no reports portraying changes in oral microbial community in the patients of Indian subcontinent, which has the highest incidence of oral cancer per year, globally. To establish the association of bacterial dysbiosis and oral squamous cell carcinoma (OSCC) among the Indian population, malignant lesions and anatomically matched adjacent normal tissues were obtained from fifty well-differentiated OSCC patients and analyzed using 16S rRNA V3-V4 amplicon based sequencing on the MiSeq platform. Interestingly, in contrast to the previous studies, a significantly lower bacterial diversity was observed in the malignant samples as compared to the normal counterpart. Overall our study identified Prevotella, Corynebacterium, Pseudomonas, Deinococcus and Noviherbaspirillum as significantly enriched genera, whereas genera including Actinomyces, Sutterella, Stenotrophomonas, Anoxybacillus, and Serratia were notably decreased in the OSCC lesions. Moreover, we demonstrated HPV-16 but not HPV-18 was significantly associated with the OSCC development. In future, with additional validation, this panel could directly be applied into clinical diagnostic and prognostic workflows for OSCC in Indian scenario.

In silico analysis of altered expression of long non-coding RNA in SARS-CoV-2 infected cells and their possible regulation by STAT1, STAT3 and interferon regulatory factors.

Altered expression of long noncoding RNA (lncRNA), longer than 200 nucleotides without potential for coding protein, has been observed in diverse human diseases including viral diseases. It is largely unknown whether lncRNA would deregulate in SARS-CoV-2 infection, causing ongoing pandemic COVID-19. To identify, if lncRNA was deregulated in SARS-CoV-2 infected cells, we analyzed in silico the data in GSE147507. It was revealed that expression of 20 lncRNA like MALAT1, NEAT1 was increased and 4 lncRNA like PART1, TP53TG1 was decreased in at least two independent cell lines infected with SARS-CoV-2. Expression of NEAT1 was also increased in lungs tissue of COVID-19 patients. The deregulated lncRNA could interact with more than 2800 genes/proteins and 422 microRNAs as revealed from the database that catalogs experimentally determined interactions. Analysis with the interacting gene/protein partners of deregulated lncRNAs revealed that these genes/proteins were associated with many pathways related to viral infection, inflammation and immune functions. To find out whether these lncRNAs could be regulated by STATs and interferon regulatory factors (IRFs), we used ChIPBase v2.0 that catalogs experimentally determined binding from ChIP-seq data. It was revealed that any one of the transcription factors IRF1, IRF4, STAT1, STAT3 and STAT5A had experimentally determined binding at regions within -5kb to +1kb of the deregulated lncRNAs in at least 2 independent cell lines/conditions. Our analysis revealed that several lncRNAs could be regulated by IRF1, IRF4 STAT1 and STAT3 in response to SARS-CoV-2 infection and lncRNAs might be involved in antiviral response. However, these in silico observations are necessary to be validated experimentally.

Characterizations of SARS-CoV-2 mutational profile, spike protein stability and viral transmission.

The recent pandemic of SARS-CoV-2 infection has affected more than 3.0 million people worldwide with more than 200 thousand reported deaths. The SARS-CoV-2 genome has the capability of gaining rapid mutations as the virus spreads. Whole-genome sequencing data offers a wide range of opportunities to study mutation dynamics. The advantage of an increasing amount of whole-genome sequence data of SARS-CoV-2 intrigued us to explore the mutation profile across the genome, to check the genome diversity, and to investigate the implications of those mutations in protein stability and viral transmission. We have identified frequently mutated residues by aligning ~660 SARS-CoV-2 genomes and validated in 10,000 datasets available in GISAID Nextstrain. We further evaluated the potential of these frequently mutated residues in protein structure stability of spike glycoprotein and their possible functional consequences in other proteins. Among the 11 genes, surface glycoprotein, nucleocapsid, ORF1ab, and ORF8 showed frequent mutations, while envelop, membrane, ORF6, ORF7a and ORF7b showed conservation in terms of amino acid substitutions. Combined analysis with the frequently mutated residues identified 20 viral variants, among which 12 specific combinations comprised more than 97% of the isolates considered for the analysis. Some of the mutations across different proteins showed co-occurrences, suggesting their structural and/or functional interaction among different SARS-COV-2 proteins, and their involvement in adaptability and viral transmission. Analysis of protein structure stability of surface glycoprotein mutants indicated the viability of specific variants and are more prone to be temporally and spatially distributed across the globe. A similar empirical analysis of other proteins indicated the existence of important functional implications of several variants. Identification of frequently mutated variants among COVID-19 patients might be useful for better clinical management, contact tracing, and containment of the disease.