Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Epigenomics of conventional type-I dendritic cells depicted preferential control of TLR9 versus TLR3 response by NCoR1 through differential IRF3 activation.

Tight control of gene regulation in dendritic cells (DCs) is important to mount pathogen specific immune responses. Apart from transcription factor binding, dynamic regulation of enhancer activity through global transcriptional repressors like Nuclear Receptor Co-repressor 1 (NCoR1) plays a major role in fine-tuning of DC responses. However, how NCoR1 regulates enhancer activity and gene expression in individual or multiple Toll-like receptor (TLR) activation in DCs is largely unknown. In this study, we did a comprehensive epigenomic analysis of murine conventional type-I DCs (cDC1) across different TLR ligation conditions. We profiled gene expression changes along with H3K27ac active enhancers and NCoR1 binding in the TLR9, TLR3 and combined TLR9 + TLR3 activated cDC1. We observed spatio-temporal activity of TLR9 and TLR3 specific enhancers regulating signal specific target genes. Interestingly, we found that NCoR1 differentially controls the TLR9 and TLR3-specific responses. NCoR1 depletion specifically enhanced TLR9 responses as evident from increased enhancer activity as well as TLR9-specific gene expression, whereas TLR3-mediated antiviral response genes were negatively regulated. We validated that NCoR1 KD cDC1 showed significantly decreased TLR3 specific antiviral responses through decreased IRF3 activation. In addition, decreased IRF3 binding was observed at selected ISGs leading to their decreased expression upon NCoR1 depletion. Consequently, the NCoR1 depleted cDC1 showed reduced Sendai Virus (SeV) clearance and cytotoxic potential of CD8+ T cells upon TLR3 activation. NCoR1 directly controls the majority of these TLR specific enhancer activity and the gene expression. Overall, for the first time, we revealed NCoR1 mediates transcriptional control towards TLR9 as compared to TLR3 in cDC1.

Zbtb10 transcription factor is crucial for murine cDC1 activation and cytokine secretion.

Dendritic cell (DC) activation and cytokine production is tightly regulated. In this study, we found that Zbtb10 expression is activation dependent and it is essential for the immunogenic function of cDC1. Zbtb10 knockdown (KD) significantly reduced the expression of co-stimulatory genes CD80 and CD86 along with cytokines including IL-12, IL-6, and IL-10, in activated cDC1 Mutu-DC line. Consequently, the clonal expansion of CD44+ effector T cells in co-cultured CD4+ T cells was drastically reduced owing to significantly reduced IL-2. At the same time, these CD44+ effector T cells were unable to differentiate toward Tbet+ IFNγ+ Th1 subtype. Instead, an increased frequency of Th2 cells expressing GATA3+ and IL-13+ was observed. Interestingly, in Zbtb10 KD condition the co-cultured T cells depicted increased expression of PD1 and LAG3, the T-cell anergic markers. Moreover, the global transcriptome analysis identified that Zbtb10 is pertinent for DC activation and its depletion in cDC1 completely shuts down their immune responses. Mechanistic analysis revealed that Zbtb10 KD enhanced the expression of NKRF (NF-κB repressing factor) leading to drastic suppression of NF-κB related genes. Zbtb10 KD abrogated p65 and RelB nuclear translocation, thereby controlling the activation and maturation of cDC1 and the ensuing adaptive T cell responses.

The show and tell of cross-presentation.

Cross-presentation is the culmination of complex subcellular processes that allow the processing of exogenous proteins and the presentation of resultant peptides on major histocompatibility class I (MHC-I) molecules to CD8 T cells. Dendritic cells (DCs) are a cell type that uniquely specializes in cross-presentation, mainly in the context of viral or non-viral infection and cancer. DCs have an extensive network of endovesicular pathways that orchestrate the biogenesis of an ideal cross-presentation compartment where processed antigen, MHC-I molecules, and the MHC-I peptide loading machinery all meet. As a central conveyor of information to CD8 T cells, cross-presentation allows cross-priming of T cells which carry out robust adaptive immune responses for tumor and viral clearance. Cross-presentation can be canonical or noncanonical depending on the functional status of the transporter associated with antigen processing (TAP), which in turn influences the vesicular route of MHC-I delivery to internalized antigen and the cross-presented repertoire of peptides. Because TAP is a central node in MHC-I presentation, it is targeted by immune evasive viruses and cancers. Thus, understanding the differences between canonical and noncanonical cross-presentation may inform new therapeutic avenues against cancer and infectious disease. Defects in cross-presentation on a cellular and genetic level lead to immune-related disease progression, recurrent infection, and cancer progression. In this chapter, we review the process of cross-presentation beginning with the DC subsets that conduct cross-presentation, the signals that regulate cross-presentation, the vesicular trafficking pathways that orchestrate cross-presentation, the modes of cross-presentation, and ending with disease contexts where cross-presentation plays a role.

NCoR1 controls immune tolerance in conventional dendritic cells by fine-tuning glycolysis and fatty acid oxidation.

Dendritic cells (DCs) undergo rapid metabolic reprogramming to generate signal-specific immune responses. The fine control of cellular metabolism underlying DC immune tolerance remains elusive. We have recently reported that NCoR1 ablation generates immune-tolerant DCs through enhanced IL-10, IL-27 and SOCS3 expression. In this study, we did comprehensive metabolic profiling of these tolerogenic DCs and identified that they meet their energy requirements through enhanced glycolysis and oxidative phosphorylation (OXPHOS), supported by fatty acid oxidation-driven oxygen consumption. In addition, the reduced pyruvate and glutamine oxidation with a broken TCA cycle maintains the tolerogenic state of the cells. Mechanistically, the AKT-mTOR-HIF-1α-axis mediated glycolysis and CPT1a-driven ß-oxidation were enhanced in these tolerogenic DCs. To confirm these observations, we used synthetic metabolic inhibitors and found that the combined inhibition of HIF-1α and CPT1a using KC7F2 and etomoxir, respectively, compromised the overall transcriptional signature of immunological tolerance including the regulatory cytokines IL-10 and IL-27. Functionally, treatment of tolerogenic DCs with dual KC7F2 and etomoxir treatment perturbed the polarization of co-cultured naïve CD4+ T helper (Th) cells towards Th1 than Tregs, ex vivo and in vivo. Physiologically, the Mycobacterium tuberculosis (Mtb) infection model depicted significantly reduced bacterial burden in BMcDC1 ex vivo and in CD103+ lung DCs in Mtb infected NCoR1DC-/-mice. The spleen of these infected animals also showed increased Th1-mediated responses in the inhibitor-treated group. These findings suggested strong involvement of NCoR1 in immune tolerance. Our validation in primary human monocyte-derived DCs (moDCs) showed diminished NCOR1 expression in dexamethasone-derived tolerogenic moDCs along with suppression of CD4+T cell proliferation and Th1 polarization. Furthermore, the combined KC7F2 and etomoxir treatment rescued the decreased T cell proliferative capacity and the Th1 phenotype. Overall, for the first time, we demonstrated here that NCoR1 mediated control of glycolysis and fatty acid oxidation fine-tunes immune tolerance versus inflammation balance in murine and human DCs.

SMRT and NCoR1 fine-tune inflammatory versus tolerogenic balance in dendritic cells by differentially regulating STAT3 signaling.

Dendritic cell (DC) fine-tunes inflammatory versus tolerogenic responses to protect from immune-pathology. However, the role of co-regulators in maintaining this balance is unexplored. NCoR1-mediated repression of DC immune-tolerance has been recently reported. Here we found that depletion of NCoR1 paralog SMRT (NCoR2) enhanced cDC1 activation and expression of IL-6, IL-12 and IL-23 while concomitantly decreasing IL-10 expression/secretion. Consequently, co-cultured CD4+ and CD8+ T-cells depicted enhanced Th1/Th17 frequency and cytotoxicity, respectively. Comparative genomic and transcriptomic analysis demonstrated differential regulation of IL-10 by SMRT and NCoR1. SMRT depletion represses mTOR-STAT3-IL10 signaling in cDC1 by down-regulating NR4A1. Besides, Nfkbia and Socs3 were down-regulated in Ncor2 (Smrt) depleted cDC1, supporting increased production of inflammatory cytokines. Moreover, studies in mice showed, adoptive transfer of SMRT depleted cDC1 in OVA-DTH induced footpad inflammation led to increased Th1/Th17 and reduced tumor burden after B16 melanoma injection by enhancing oncolytic CD8+ T-cell frequency, respectively. We also depicted decreased Ncor2 expression in Rheumatoid Arthritis, a Th1/Th17 disease.

Cancer-Associated circRNA-miRNA-mRNA Regulatory Networks: A Meta-Analysis.

Recent advances in sequencing technologies and the discovery of non-coding RNAs (ncRNAs) have provided new insights in the molecular pathogenesis of cancers. Several studies have implicated the role of ncRNAs, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and recently discovered circular RNAs (circRNAs) in tumorigenesis and metastasis. Unlike linear RNAs, circRNAs are highly stable and closed-loop RNA molecules. It has been established that circRNAs regulate gene expression by controlling the functions of miRNAs and RNA-binding protein (RBP) or by translating into proteins. The circRNA-miRNA-mRNA regulatory axis is associated with human diseases, such as cancers, Alzheimer's disease, and diabetes. In this study, we explored the interaction among circRNAs, miRNAs, and their target genes in various cancers using state-of-the-art bioinformatics tools. We identified differentially expressed circRNAs, miRNAs, and mRNAs on multiple cancers from publicly available data. Furthermore, we identified many crucial drivers and tumor suppressor genes in the circRNA-miRNA-mRNA regulatory axis in various cancers. Together, this study data provide a deeper understanding of the circRNA-miRNA-mRNA regulatory mechanisms in cancers.

Single-Cell Immunogenomic Approach Identified SARS-CoV-2 Protective Immune Signatures in Asymptomatic Direct Contacts of COVID-19 Cases.

The response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is largely impacted by the level of virus exposure and status of the host immunity. The nature of protection shown by direct asymptomatic contacts of coronavirus disease 2019 (COVID-19)-positive patients is quite intriguing. In this study, we have characterized the antibody titer, SARS-CoV-2 surrogate virus neutralization, cytokine levels, single-cell T-cell receptor (TCR), and B-cell receptor (BCR) profiling in asymptomatic direct contacts, infected cases, and controls. We observed significant increase in antibodies with neutralizing amplitude in asymptomatic contacts along with cytokines such as Eotaxin, granulocyte-colony stimulating factor (G-CSF), interleukin 7 (IL-7), migration inhibitory factor (MIF), and macrophage inflammatory protein-1α (MIP-1α). Upon single-cell RNA (scRNA) sequencing, we explored the dynamics of the adaptive immune response in few representative asymptomatic close contacts and COVID-19-infected patients. We reported direct asymptomatic contacts to have decreased CD4+ naive T cells with concomitant increase in CD4+ memory and CD8+ Temra cells along with expanded clonotypes compared to infected patients. Noticeable proportions of class switched memory B cells were also observed in them. Overall, these findings gave an insight into the nature of protection in asymptomatic contacts.

Aurora Kinase B Expression, Its Regulation and Therapeutic Targeting in Human Retinoblastoma.

Purpose: Aurora kinase B (AURKB) plays a pivotal role in the regulation of mitosis and is gaining prominence as a therapeutic target in cancers; however, the role of AURKB in retinoblastoma (RB) has not been studied. The purpose of this study was to determine if AURKB plays a role in RB, how its expression is regulated, and whether it could be specifically targeted. Methods: The protein expression of AURKB was determined using immunohistochemistry in human RB patient specimens and immunoblotting in cell lines. Pharmacological inhibition and shRNA-mediated knockdown were used to understand the role of AURKB in cell viability, apoptosis, and cell cycle distribution. Cell viability in response to AURKB inhibition was also assessed in enucleated RB specimens. Immunoblotting was employed to determine the protein levels of phospho-histone H3, p53, p21, and MYCN. Chromatin immunoprecipitation-qPCR was performed to verify the binding of MYCN on the promoter region of AURKB. Results: The expression of AURKB was found to be markedly elevated in human RB tissues, and the overexpression significantly correlated with optic nerve and anterior chamber invasion. Targeting AURKB with small-molecule inhibitors and shRNAs resulted in reduced cell survival and increased apoptosis and cell cycle arrest at the G2/M phase. More importantly, primary RB specimens showed decreased cell viability in response to pharmacological AURKB inhibition. Additional studies have demonstrated that the MYCN oncogene regulates the expression of AURKB in RB. Conclusions: AURKB is overexpressed in RB, and targeting it could serve as a novel therapeutic strategy to restrict tumor cell growth.

Clinical, Virological, Immunological, and Genomic Characterization of Asymptomatic and Symptomatic Cases With SARS-CoV-2 Infection in India.

Purpose: The current global pandemic of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), led to the investigation with clinical, biochemical, immunological, and genomic characterization from patients to understand the pathophysiology of viral infection. Methods: Samples were collected from six asymptomatic and six symptomatic SARS-CoV-2-confirmed hospitalized patients in Bhubaneswar, Odisha, India. Clinical details, biochemical parameters, and treatment regimen were collected from a hospital; viral load was determined by RT-PCR; and the levels of cytokines and circulating antibodies in plasma were assessed by Bio-Plex and isotyping, respectively. In addition, whole-genome sequencing of viral strains and mutational analysis were carried out. Results: Analysis of the biochemical parameters highlighted the increased levels of C-reactive protein (CRP), lactate dehydrogenase (LDH), serum SGPT, serum SGOT, and ferritin in symptomatic patients. Symptomatic patients were mostly with one or more comorbidities, especially type 2 diabetes (66.6%). The virological estimation revealed that there was no significant difference in viral load of oropharyngeal (OP) samples between the two groups. On the other hand, viral load was higher in plasma and serum samples of symptomatic patients, and they develop sufficient amounts of antibodies (IgG, IgM, and IgA). The levels of seven cytokines (IL-6, IL-1α, IP-10, IL-8, IL-10, IFN-α2, IL-15) were found to be highly elevated in symptomatic patients, while three cytokines (soluble CD40L, GRO, and MDC) were remarkably higher in asymptomatic patients. The whole-genome sequence analysis revealed that the current isolates were clustered with 19B, 20A, and 20B clades; however, 11 additional changes in Orf1ab, spike, Orf3a, Orf8, and nucleocapsid proteins were acquired. The D614G mutation in spike protein is linked with higher virus replication efficiency and severe SARS-CoV-2 infection as three patients had higher viral load, and among them, two patients with this mutation passed away. Conclusions: This is the first comprehensive study of SARS-CoV-2 patients from India. This will contribute to a better understanding of the pathophysiology of SARS-CoV-2 infection and thereby advance the implementation of effective disease control strategies.

BedSect: An Integrated Web Server Application to Perform Intersection, Visualization, and Functional Annotation of Genomic Regions From Multiple Datasets.

A large number of genomic regions, such as transcription factor binding sites (TFBSs) captured from next generation sequencing (NGS) data analyses or those available from the public resource database ENCODE, are generally overlapped to answer a variety of biological questions. Though several command-line tools are available to perform such an analysis, there is a notable lack of an integrated webserver application with which to identify genomic region intersections, generate publication-ready plots depicting subsets of the overlapped regions, and perform functional annotation. Thus, there is an ardent need for a comprehensive and user-friendly webserver application that allows the users to either upload multiple datasets or select from the integrated Gene Transcription Regulation Database (GTRD). We thus introduce BedSect (http://imgsb.org/bedsect/.), which not only fulfils the above criteria but also performs intersection analysis along with visualization of the intersection regions as an UpSet and correlation plot using the integrated Shiny application. Moreover, analyses, including functional annotation, gene ontology, and biological pathways enrichment for the identified unique and intersected genomic regions, can also be performed using the integrated GREAT tool. To view the genomic regions in the genome browser, the inbuilt hyperlink for UCSC can redirect the user to visualize the results as custom tracks.

Analysis of Indian SARS-CoV-2 Genomes Reveals Prevalence of D614G Mutation in Spike Protein Predicting an Increase in Interaction With TMPRSS2 and Virus Infectivity.

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, has emerged as a global pandemic worldwide. In this study, we used ARTIC primers-based amplicon sequencing to profile 225 SARS-CoV-2 genomes from India. Phylogenetic analysis of 202 high-quality assemblies identified the presence of all the five reported clades 19A, 19B, 20A, 20B, and 20C in the population. The analyses revealed Europe and Southeast Asia as two major routes for introduction of the disease in India followed by local transmission. Interestingly, the19B clade was found to be more prevalent in our sequenced genomes (17%) compared to other genomes reported so far from India. Haplotype network analysis showed evolution of 19A and 19B clades in parallel from predominantly Gujarat state in India, suggesting it to be one of the major routes of disease transmission in India during the months of March and April, whereas 20B and 20C appeared to evolve from 20A. At the same time, 20A and 20B clades depicted prevalence of four common mutations 241 C > T in 5' UTR, P4715L, F942F along with D614G in the Spike protein. D614G mutation has been reported to increase virus shedding and infectivity. Our molecular modeling and docking analysis identified that D614G mutation resulted in enhanced affinity of Spike S1-S2 hinge region with TMPRSS2 protease, possibly the reason for increased shedding of S1 domain in G614 as compared to D614. Moreover, we also observed an increased concordance of G614 mutation with the viral load, as evident from decreased Ct value of Spike and the ORF1ab gene.