Compositional diversity and evolutionary pattern of coronavirus accessory proteins.

Accessory proteins play important roles in the interaction between coronaviruses and their hosts. Accordingly, a comprehensive study of the compositional diversity and evolutionary patterns of accessory proteins is critical to understanding the host adaptation and epidemic variation of coronaviruses. Here, we developed a standardized genome annotation tool for coronavirus (CoroAnnoter) by combining open reading frame prediction, transcription regulatory sequence recognition and homologous alignment. Using CoroAnnoter, we annotated 39 representative coronavirus strains to form a compositional profile for all of the accessary proteins. Large variations were observed in the number of accessory proteins of 1-10 for different coronaviruses, with SARS-CoV-2 and SARS-CoV having the most (9 and 10, respectively). The variation between SARS-CoV and SARS-CoV-2 accessory proteins could be traced back to related coronaviruses in other hosts. The genomic distribution of accessory proteins had significant intra-genus conservation and inter-genus diversity and could be grouped into 1, 4, 2 and 1 types for alpha-, beta-, gamma-, and delta-coronaviruses, respectively. Evolutionary analysis suggested that accessory proteins are more conservative locating before the N-terminal of proteins E and M (E-M), while they are more diverse after these proteins. Furthermore, comparison of virus-host interaction networks of SARS-CoV-2 and SARS-CoV accessory proteins showed that they share multiple antiviral signaling pathways, those involved in the apoptotic process, viral life cycle and response to oxidative stress. In summary, our study provides a tool for coronavirus genome annotation and builds a comprehensive profile for coronavirus accessory proteins covering their composition, classification, evolutionary pattern and host interaction.

Genomic annotation and molecular evolution of monkeypox virus outbreak in 2022.

Monkeypox virus (MPXV) has generally circulated in West and Central Africa since its emergence. Recently, sporadic MPXV infections in several nonendemic countries have attracted widespread attention. Here, we conducted a systematic analysis of the recent outbreak of MPXV-2022, including its genomic annotation and molecular evolution. The phylogenetic analysis indicated that the MPXV-2022 strains belong to the same lineage of the MPXV strain isolated in 2018. However, compared with the MPXV strain in 2018, in total 46 new consensus mutations were observed in the MPXV-2022 strains, including 24 nonsynonymous mutations. By assigning mutations to 187 proteins encoded by the MPXV genome, we found that 10 proteins in the MPXV are more prone to mutation, including D2L-like, OPG023, OPG047, OPG071, OPG105, OPG109, A27L-like, OPG153, OPG188, and OPG210 proteins. In the MPXV-2022 strains, four and three nucleotide substitutions are observed in OPG105 and OPG210, respectively. Overall, our studies illustrated the genome evolution of the ongoing MPXV outbreak and pointed out novel mutations as a reference for further studies.

sitePath: a visual tool to identify polymorphism clades and help find fixed and parallel mutations.

BACKGROUND: Identifying polymorphism clades on phylogenetic trees could help detect punctual mutations that are associated with viral functions. With visualization tools coloring the tree, it is easy to visually find clades where most sequences have the same polymorphism state. However, with the fast accumulation of viral sequences, a computational tool to automate this process is urgently needed. RESULTS: Here, by implementing a branch-and-bound-like search method, we developed an R package named sitePath to identify polymorphism clades automatically. Based on the identified polymorphism clades, fixed and parallel mutations could be inferred. Furthermore, sitePath also integrated visualization tools to generate figures of the calculated results. In an example with the influenza A virus H3N2 dataset, the detected fixed mutations coincide with antigenic shift mutations. The highly specificity and sensitivity of sitePath in finding fixed mutations were achieved for a range of parameters and different phylogenetic tree inference software. CONCLUSIONS: The result suggests that sitePath can identify polymorphism clades per site. The clustering of sequences on a phylogenetic tree can be used to infer fixed and parallel mutations. High-quality figures of the calculated results could also be generated by sitePath.

Identification of dynamic signatures associated with smoking-related squamous cell lung cancer and chronic obstructive pulmonary disease.

Chronic obstructive pulmonary disease (COPD) is a risk factor for the development of lung cancer. The aim of this study was to identify early diagnosis biomarkers for lung squamous cell carcinoma (SQCC) in COPD patients and to determine the potential pathogenetic mechanisms. The GSE12472 data set was downloaded from the Gene Expression Omnibus database. Differentially co-expressed links (DLs) and differentially expressed genes (DEGs) in both COPD and normal tissues, or in both SQCC + COPD and COPD samples were used to construct a dynamic network associated with high-risk genes for the SQCC pathogenetic process. Enrichment analysis was performed based on Gene Ontology annotations and Kyoto Encyclopedia of Genes and Genomes pathway analysis. We used the gene expression data and the clinical information to identify the co-expression modules based on weighted gene co-expression network analysis (WGCNA). In total, 205 dynamic DEGs, 5034 DLs and one pathway including CDKN1A, TP53, RB1 and MYC were found to have correlations with the pathogenetic progress. The pathogenetic mechanisms shared by both SQCC and COPD are closely related to oxidative stress, the immune response and infection. WGCNA identified 11 co-expression modules, where magenta and black were correlated with the "time to distant metastasis." And the "surgery due to" was closely related to the brown and blue modules. In conclusion, a pathway that includes TP53, CDKN1A, RB1 and MYC may play a vital role in driving COPD towards SQCC. Inflammatory processes and the immune response participate in COPD-related carcinogenesis.

Azithromycin Protects against Zika virus Infection by Upregulating virus-induced Type I and III Interferon Responses.

Azithromycin (AZM) is a widely used antibiotic, with additional antiviral and anti-inflammatory properties that remain poorly understood. Although Zika virus (ZIKV) poses a significant threat to global health, there are currently no vaccines or effective therapeutics against it. Herein, we report that AZM effectively suppresses ZIKV infection in vitro by targeting a late stage in the viral life cycle. Besides that, AZM upregulates the expression of host type I and III interferons and several of their downstream interferon-stimulated genes (ISGs) in response to ZIKV infection. In particular, we found that AZM upregulates the expression of MDA5 and RIG-I, pathogen recognition receptors (PRRs) induced by ZIKV infection, and increases the levels of phosphorylated TBK1 and IRF3. Interestingly, AZM treatment upregulates phosphorylation of TBK1, without inducing phosphorylation of IRF3 by itself. These findings highlight the potential use of AZM as a broad antiviral agent to combat viral infection and prevent ZIKV associated devastating clinical outcomes, such as congenital microcephaly.

Identification of the mulberry genes involved in ethylene biosynthesis and signaling pathways and the expression of MaERF-B2-1 and MaERF-B2-2 in the response to flooding stress.

The phytohormone ethylene is essential to plant growth and development. It plays crucial roles in responses to biotic and abiotic stress. The mulberry tree is an important crop plant in countries in which people rear silkworms for silk production. The availability of the mulberry genome has made it possible to identify mulberry genes involved in ethylene biosynthesis and signal pathways. A total of 145 mulberry genes were identified by both homology-based and hidden Markov model (HMM) search, including 29 genes associated with ethylene biosynthesis and 116 genes in the AP2/ERF family. Studies on gene structure have provided a genetic basis for understanding the functions of these genes. The differences in gene expression were also observed in different tissues. The expression of two mulberry genes in the AP2/ERF family, MaERF-B2-1 and MaERF-B2-2, was found to be associated with the response to flooding stress.

Suppression of the METTL3-m6A-integrin ß1 axis by extracellular acidification impairs T cell infiltration and antitumor activity.

The acidic metabolic byproducts within the tumor microenvironment (TME) hinder T cell effector functions. However, their effects on T cell infiltration remain largely unexplored. Leveraging the comprehensive The Cancer Genome Atlas dataset, we pinpoint 16 genes that correlate with extracellular acidification and establish a metric known as the "tumor acidity (TuAci) score" for individual patients. We consistently observe a negative association between the TuAci score and T lymphocyte score (T score) across various human cancer types. Mechanistically, extracellular acidification significantly impedes T cell motility by suppressing podosome formation. This phenomenon can be attributed to the reduced expression of methyltransferase-like 3 (METTL3) and the modification of RNA N6-methyladenosine (m6A), resulting in a subsequent decrease in the expression of integrin ß1 (ITGB1). Importantly, enforced ITGB1 expression leads to enhanced T cell infiltration and improved antitumor activity. Our study suggests that modulating METTL3 activity or boosting ITGB1 expression could augment T cell infiltration within the acidic TME, thereby improving the efficacy of cell therapy.

sc-ImmuCC: hierarchical annotation for immune cell types in single-cell RNA-seq.

Accurately identifying immune cell types in single-cell RNA-sequencing (scRNA-Seq) data is critical to uncovering immune responses in health or disease conditions. However, the high heterogeneity and sparsity of scRNA-Seq data, as well as the similarity in gene expression among immune cell types, poses a great challenge for accurate identification of immune cell types in scRNA-Seq data. Here, we developed a tool named sc-ImmuCC for hierarchical annotation of immune cell types from scRNA-Seq data, based on the optimized gene sets and ssGSEA algorithm. sc-ImmuCC simulates the natural differentiation of immune cells, and the hierarchical annotation includes three layers, which can annotate nine major immune cell types and 29 cell subtypes. The test results showed its stable performance and strong consistency among different tissue datasets with average accuracy of 71-90%. In addition, the optimized gene sets and hierarchical annotation strategy could be applied to other methods to improve their annotation accuracy and the spectrum of annotated cell types and subtypes. We also applied sc-ImmuCC to a dataset composed of COVID-19, influenza, and healthy donors, and found that the proportion of monocytes in patients with COVID-19 and influenza was significantly higher than that in healthy people. The easy-to-use sc-ImmuCC tool provides a good way to comprehensively annotate immune cell types from scRNA-Seq data, and will also help study the immune mechanism underlying physiological and pathological conditions.

Psoriasis and Leprosy: An Arcane Relationship.

Purpose: Psoriasis (Ps) and leprosy are chronic inflammatory skin disorders, characterised by enhanced innate and adaptive immunity. Ps and leprosy rarely coexist. The molecular immune mechanism of the Ps and leprosy rarely coexistence is unclear. Patients and Methods: RNA-sequencing (RNA-seq) was performed on 20 patients with Ps, 5 adults with lepromatous leprosy (L-lep), and 5 patients with tuberculoid leprosy (T-lep) to analyse the differentially expressed genes (DEGs) between them. Moreover, the biological mechanism of Ps and leprosy was explored by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, Gene Ontology (GO) analysis, Gene Set Enrichment Analysis analysis, and protein-protein interaction (PPI) analyses. Finally, 13 DEGs of 10 skin biopsies of Ps patients, 6 samples of L-lep patients, 6 samples of T-lep patients and 5 healthy controls were confirmed by quantitative real-time polymerase chain reaction (qRT-PCR). Results: The PPI network was constructed and primarily associated with immune response, IL-17 signalling, and Toll-like receptor pathway between Ps and leprosy. Th17 markers (interleukin (IL)-19, IL-20, IL-36A, IL-36G, IL-22, IL-17A, and lipocalin-2 (LCN2) had higher expression in Ps than in L-lep and T-lep, whereas macrophage biomarkers (CLEC4E and TREM2), SPP1, and dendritic cell (DC)-related hallmarks (ITGAX) and TNF-a had significantly lower expression across Ps and T-lep than in L-lep. Conclusion: To put it simply, Ps patients with IL-17A, IL-19, IL-20, IL-36A, IL-36G, and IL-22 in conjunction with LCN2 with up-graduated expression might be not susceptible to L-lep. However, high levels of CLEC4E, TREM2, and SPP1 in L-lep patients indicated that they unlikely suffered from Ps.

Conserved Pattern and Potential Role of Recurrent Deletions in SARS-CoV-2 Evolution.

SARS-CoV-2 continues adapting to human hosts during the current worldwide pandemic since 2019. This virus evolves through multiple means, such as single nucleotide mutations and structural variations, which has brought great difficulty to disease prevention and control of COVID-19. Structural variation, including multiple nucleotide changes like insertions and deletions, has a greater impact relative to single nucleotide mutation on both genome structures and protein functions. In this study, we found that deletion occurred frequently in not only SARS-CoV-2 but also in other SARS-related coronaviruses. These deletions showed obvious location bias and formed 45 recurrent deletion regions in the viral genome. Some of these deletions showed proliferation advantages, including four high-frequency deletions (nsp6 Δ106-109, S Δ69-70, S Δ144, and Δ28271) that were detected in around 50% of SARS-CoV-2 genomes and other 19 median-frequency deletions. In addition, the association between deletions and the WHO reported variants of concern (VOC) and variants of interest (VOI) of SARS-CoV-2 indicated that these variants had a unique combination of deletion patterns. In the spike (S) protein, the deletions in SARS-CoV-2 were mainly in the N-terminal domain. Some deletions, such as S Δ144/145 and S Δ243-244, have been confirmed to block the binding sites of neutralizing antibodies. Overall, this study revealed a conservative regional pattern and the potential effect of some deletions in SARS-CoV-2 over the whole genome, providing important evidence for potential epidemic control and vaccine development. IMPORTANCE Mutations in SARS-CoV-2 were studied extensively, while only the structure variations on the spike protein were discussed well in previous studies. To study the role of structural variations in virus evolution, we described the distribution of structure variations on the whole genome. Conserved patterns were found of deletions among SARS-CoV-2, SARS-CoV-2-like, and SARS-CoV-like viruses. There were 45 recurrent deletion regions (RDRs) in SARS-CoV-2 generated through the integration of deleted positions. In these regions, four high-frequency deletions parallelly appeared in multiple strains. Furthermore, in the spike protein, the deletions in SARS-CoV-2 were mainly in the N-terminal domain, blocking the binding sites of some neutralizing antibodies, while the structural variations in SARS-related coronavirus were mainly in the N-terminal domain and receptor binding domain. The receptor binding domain is highly related to hosting recognition. The deletions in the receptor binding domain may play a role in host adaption.

Insight into the Molecular Characteristics of Langhans Giant Cell by Combination of Laser Capture Microdissection and RNA Sequencing.

PURPOSE: The presence of Langhans giant cell (LGC) is a hallmark of mycobacterium-induced granuloma. The molecular characteristics and functions of LGC remain unclear to date. The study aimed to systematically characterize the molecular characteristics of LGC and reveal the potential functions. METHODS: Human LGCs were purified through laser capture microdissection (LCM) in vitro. RNA sequencing and in-depth transcriptome analysis were performed for purified LGCs and macrophages in the same system. Skin samples from mycobacterial infection patients were used to confirm some of the transcriptional expression. RESULTS: Human LGCs have different expression pattern from macrophages in the same in vitro system. A total of 967 differentially expressed genes were found. Bioinformatics analysis showed that LGCs are is characterized by active cell shape regulation, increased cytoskeletal components, weakened energy metabolism level, and reduced immune response. CCL7 may be a specific molecular for LGC to communicate with CCR1-expression cells in granuloma. CONCLUSION: LGCs have unique molecular characteristics different from that of macrophages. They may play a role in maintaining the hemostasis in granuloma.

Genetic differentiation and diversity of SARS-CoV-2 Omicron variant in its early outbreak.

The recently emerged Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has quickly spread around the world. Although many consensus mutations of the Omicron variant have been recognized, little is known about its genetic variation during its transmission in the population. Here, we comprehensively analyzed the genetic differentiation and diversity of the Omicron variant during its early outbreak. We found that Omicron achieved more structural variations, especially deletions, on the SARS-CoV-2 genome than the other four variants of concern (Alpha, Beta, Gamma, and Delta) in the same timescale. In addition, the Omicron variant acquired, except for 50 consensus mutations, seven great new non-synonymous nucleotide substitutions during its spread. Three of them are on the S protein, including S_A701V, S_L1081V, and S_R346K, which belong to the receptor-binding domain (RBD). The Omicron BA.1 branch could be divided into five divergent groups spreading across different countries and regions based on these seven novel mutations. Furthermore, we found that the Omicron variant possesses more mutations related to a faster transmission rate than the other SARS-CoV-2 variants by assessing the relationship between the genetic diversity and transmission rate. The findings indicated that more attention should be paid to the significant genetic differentiation and diversity of the Omicron variant for better disease prevention and control.

Immune profiles in mouse brain and testes infected by Zika virus with variable pathogenicity.

The Zika virus is responsible for neurological diseases such as microcephaly, Guillain-Barré syndrome, neuropathy, and myelitis in human adults and children. Previous studies have shown that the Zika virus can infect nerve progenitor cells and interfere with neural development. However, it is unclear how the immune system responds to infection with Zika viruses with variable pathogenicity. Here, we used two Zika strains with relatively different pathogenicity, the Asian ancestral strain CAM/2010 and the America pandemic strain GZ01/2016, to infect the brains of mice. We found that both strains elicited a strong immune response. Notably, the strain with relatively high pathogenicity, GZ01/2016, caused more intense immune regulation, with stronger CD8+ T cell and macrophage activation at 14 days post infection (dpi), as well as a greater immune gene disturbance. Notably, several TNF family genes were upregulated at 14 dpi, including Tnfrsf9, Tnfsf13, Tnfrsf8, Cd40, and Tnfsf10. It was notable that GZ01/2016 could maintain the survival of nerve cells at 7dpi but caused neurological disorders at 14dpi. These results indicate that Zika viruses with high pathogenicity may induce sustained activation of the immune system leading to nerve tissue damage.

Detecting Potentially Adaptive Mutations from the Parallel and Fixed Patterns in SARS-CoV-2 Evolution.

Early identification of adaptive mutations could provide timely help for the control and prevention of the COVID-19 pandemic. The fast accumulation of SARS-CoV-2 sequencing data provides important support, while also raising a great challenge for the recognition of adaptive mutations. Here, we proposed a computational strategy to detect potentially adaptive mutations from their fixed and parallel patterns in the phylogenetic trajectory. We found that the biological meanings of fixed substitution and parallel mutation are highly complementary, and can reasonably be integrated as a fixed and parallel (paraFix) mutation, to identify potentially adaptive mutations. Tracking the dynamic evolution of SARS-CoV-2, 37 sites in spike protein were identified as having experienced paraFix mutations. Interestingly, 70% (26/37) of them have already been experimentally confirmed as adaptive mutations. Moreover, most of the mutations could be inferred as paraFix mutations one month earlier than when they became regionally dominant. Overall, we believe that the concept of paraFix mutations will help researchers to identify potentially adaptive mutations quickly and accurately, which will provide invaluable clues for disease control and prevention.

Oncolytic Zika virus promotes intratumoral T cell infiltration and improves immunotherapy efficacy in glioblastoma.

Glioblastoma (GBM) is the deadliest primary brain tumor and is generally resistant to immunotherapy because of severe dysfunction of T cells. Novel treatment options are critically needed to overcome the immunotherapy resistance of GBM. Here we demonstrate that Zika virus (ZIKV) treatment improves the efficacy of anti-PD ligand 1 (PD-L1) immunotherapy in GBM. We found that ZIKV induces a strong pro-inflammatory response and increases CD4+ and CD8+ T cell intratumoral infiltration and activation in GBM mouse models. ZIKV treatment of mice bearing GBM tumors inhibits tumor growth and prolongs survival. These therapeutic effects of ZIKV on GBM tumors are negated in mice depleted of T cells. Moreover, ZIKV dramatically promotes activation of the type I interferon signaling pathway in GBM cells. ZIKV treatment potently sensitizes GBM to PD-L1 blockade and provides significant and durable survival benefits. Our findings reveal that ZIKV overcomes the resistance of GBM to immune checkpoint blockade, which may lead to therapeutic applications of ZIKV in individuals with GBM receiving immunotherapy.

Transcriptome Analysis Reveals the Molecular Immunological Characteristics of Lesions in Patients with Halo Nevi When Compared to Stable Vitiligo, Normal Nevocytic Nevi and Cutaneous Melanoma.

BACKGROUND: Given their similar appearance and histology, halo nevi (HN) were considered as a type of vitiligo. However, whether HN have stronger immune response than stable vitiligo (VL) remains unclear. In addition, the molecular alterations in HN compared with normal nevocytic nevi (NN) and primary cutaneous melanoma (MM) must be determined. This study aimed to systematically characterize the molecular immunological features of HN. METHODS: Skin samples from patients with HN, VL, NN, and MM were obtained with informed consent. Each of the four groups underwent transcriptome sequencing and data analysis were for pairwise comparison. Quantitative real-time PCR (RT-qPCR) was conducted to confirm the transcriptional expression of some differentially expressed genes (DEGs) that were closely related to immunity. RESULTS: A total of 441 and 1507 DEGs were found in the HN/NN and HN/MM groups, respectively. Compared with those of VL, HN lesions contained 162 up-regulated DEGs and 12 down-regulated DEGs. Bioinformatics analysis showed that the up-regulated genes in HN were substantially enriched in immune response, immune deficiency, and immune rejection; biological stimulation (virus, bacteria); and proliferation and activation of immune cells. Immune cell composition analysis also confirmed high expression levels of multiple immunocytes in HN. CONCLUSION: The molecular immune mechanisms of HN and VL were similar, but the immune activity of HN was stronger than that of VL. Innate and adaptive immunity were involved in the pathogenesis and progression of HN and VL.

seq-ImmuCC: Cell-Centric View of Tissue Transcriptome Measuring Cellular Compositions of Immune Microenvironment From Mouse RNA-Seq Data.

The RNA sequencing approach has been broadly used to provide gene-, pathway-, and network-centric analyses for various cell and tissue samples. However, thus far, rich cellular information carried in tissue samples has not been thoroughly characterized from RNA-Seq data. Therefore, it would expand our horizons to better understand the biological processes of the body by incorporating a cell-centric view of tissue transcriptome. Here, a computational model named seq-ImmuCC was developed to infer the relative proportions of 10 major immune cells in mouse tissues from RNA-Seq data. The performance of seq-ImmuCC was evaluated among multiple computational algorithms, transcriptional platforms, and simulated and experimental datasets. The test results showed its stable performance and superb consistency with experimental observations under different conditions. With seq-ImmuCC, we generated the comprehensive landscape of immune cell compositions in 27 normal mouse tissues and extracted the distinct signatures of immune cell proportion among various tissue types. Furthermore, we quantitatively characterized and compared 18 different types of mouse tumor tissues of distinct cell origins with their immune cell compositions, which provided a comprehensive and informative measurement for the immune microenvironment inside tumor tissues. The online server of seq-ImmuCC are freely available at http://wap-lab.org:3200/immune/.

Zika virus shedding in the stool and infection through the anorectal mucosa in mice.

Zika virus (ZIKV) has elicited global concern due to its unique biological features, unusual transmission routes, and unexpected clinical outcomes. Although ZIKV transmission through anal intercourse has been reported in humans, it remains unclear if ZIKV is detectable in the stool, if it can infect the host through the anal canal mucosa, and what the pathogenesis of such a route of infection might be in the mouse model. Herein, we demonstrate that ZIKV RNA can be recovered from stools in multiple mouse models, as well as from the stool of a ZIKV patient. Remarkably, intra-anal (i.a.) inoculation with ZIKV leads to efficient infection in both Ifnar1-/- and immunocompetent mice, characterized by extensive viral replication in the blood and multiple organs, including the brain, small intestine, testes, and rectum, as well as robust humoral and innate immune responses. Moreover, i.a. inoculation of ZIKV in pregnant mice resulted in transplacental infection and delayed fetal development. Overall, our results identify the anorectal mucosa as a potential site of ZIKV infection in mice, reveal the associated pathogenesis of i.a. infection, and highlight the complexity of ZIKV transmission through anal intercourse.

Draft genome sequence of the mulberry tree Morus notabilis.

Human utilization of the mulberry-silkworm interaction started at least 5,000 years ago and greatly influenced world history through the Silk Road. Complementing the silkworm genome sequence, here we describe the genome of a mulberry species Morus notabilis. In the 330-Mb genome assembly, we identify 128 Mb of repetitive sequences and 29,338 genes, 60.8% of which are supported by transcriptome sequencing. Mulberry gene sequences appear to evolve ~3 times faster than other Rosales, perhaps facilitating the species' spread worldwide. The mulberry tree is among a few eudicots but several Rosales that have not preserved genome duplications in more than 100 million years; however, a neopolyploid series found in the mulberry tree and several others suggest that new duplications may confer benefits. Five predicted mulberry miRNAs are found in the haemolymph and silk glands of the silkworm, suggesting interactions at molecular levels in the plant-herbivore relationship. The identification and analyses of mulberry genes involved in diversifying selection, resistance and protease inhibitor expressed in the laticifers will accelerate the improvement of mulberry plants.