SARS-CoV-2 ORF8 Mediates Signals in Macrophages and Monocytes through MyD88 Independently of the IL-17 Receptor.

SARS-CoV-2 has caused an estimated 7 million deaths worldwide to date. A secreted SARS-CoV-2 accessory protein, known as open reading frame 8 (ORF8), elicits inflammatory pulmonary cytokine responses and is associated with disease severity in COVID-19 patients. Recent reports proposed that ORF8 mediates downstream signals in macrophages and monocytes through the IL-17 receptor complex (IL-17RA, IL-17RC). However, generally IL-17 signals are found to be restricted to the nonhematopoietic compartment, thought to be due to rate-limiting expression of IL-17RC. Accordingly, we revisited the capacity of IL-17 and ORF8 to induce cytokine gene expression in mouse and human macrophages and monocytes. In SARS-CoV-2-infected human and mouse lungs, IL17RC mRNA was undetectable in monocyte/macrophage populations. In cultured mouse and human monocytes and macrophages, ORF8 but not IL-17 led to elevated expression of target cytokines. ORF8-induced signaling was fully preserved in the presence of anti-IL-17RA/RC neutralizing Abs and in Il17ra-/- cells. ORF8 signaling was also operative in Il1r1-/- bone marrow-derived macrophages. However, the TLR/IL-1R family adaptor MyD88, which is dispensable for IL-17R signaling, was required for ORF8 activity yet MyD88 is not required for IL-17 signaling. Thus, we conclude that ORF8 transduces inflammatory signaling in monocytes and macrophages via MyD88 independently of the IL-17R.

Enhanced complement activation and MAC formation accelerates severe COVID-19.

Emerging evidence indicates that activation of complement system leading to the formation of the membrane attack complex (MAC) plays a detrimental role in COVID-19. However, their pathogenic roles have never been experimentally investigated before. We used three knock out mice strains (1. C3-/-; 2. C7-/-; and 3. Cd59ab-/-) to evaluate the role of complement in severe COVID-19 pathogenesis. C3 deficient mice lack a key common component of all three complement activation pathways and are unable to generate C3 and C5 convertases. C7 deficient mice lack a complement protein needed for MAC formation. Cd59ab deficient mice lack an important inhibitor of MAC formation. We also used anti-C5 antibody to block and evaluate the therapeutic potential of inhibiting MAC formation. We demonstrate that inhibition of complement activation (in C3-/-) and MAC formation (in C3-/-. C7-/-, and anti-C5 antibody) attenuates severe COVID-19; whereas enhancement of MAC formation (Cd59ab-/-) accelerates severe COVID-19. The degree of MAC but not C3 deposits in the lungs of C3-/-, C7-/- mice, and Cd59ab-/- mice as compared to their control mice is associated with the attenuation or acceleration of SARS-CoV-2-induced disease. Further, the lack of terminal complement activation for the formation of MAC in C7 deficient mice protects endothelial dysfunction, which is associated with the attenuation of diseases and pathologic changes. Our results demonstrated the causative effect of MAC in severe COVID-19 and indicate a potential avenue for modulating the complement system and MAC formation in the treatment of severe COVID-19.

Meta-i6mA: an interspecies predictor for identifying DNA N6-methyladenine sites of plant genomes by exploiting informative features in an integrative machine-learning framework.

DNA N6-methyladenine (6mA) represents important epigenetic modifications, which are responsible for various cellular processes. The accurate identification of 6mA sites is one of the challenging tasks in genome analysis, which leads to an understanding of their biological functions. To date, several species-specific machine learning (ML)-based models have been proposed, but majority of them did not test their model to other species. Hence, their practical application to other plant species is quite limited. In this study, we explored 10 different feature encoding schemes, with the goal of capturing key characteristics around 6mA sites. We selected five feature encoding schemes based on physicochemical and position-specific information that possesses high discriminative capability. The resultant feature sets were inputted to six commonly used ML methods (random forest, support vector machine, extremely randomized tree, logistic regression, naïve Bayes and AdaBoost). The Rosaceae genome was employed to train the above classifiers, which generated 30 baseline models. To integrate their individual strength, Meta-i6mA was proposed that combined the baseline models using the meta-predictor approach. In extensive independent test, Meta-i6mA showed high Matthews correlation coefficient values of 0.918, 0.827 and 0.635 on Rosaceae, rice and Arabidopsis thaliana, respectively and outperformed the existing predictors. We anticipate that the Meta-i6mA can be applied across different plant species. Furthermore, we developed an online user-friendly web server, which is available at http://kurata14.bio.kyutech.ac.jp/Meta-i6mA/.

Deficiency of Caspase-1 Attenuates HIV-1-Associated Atherogenesis in Mice.

Within arterial plaque, HIV infection creates a state of inflammation and immune activation, triggering NLRP3/caspase-1 inflammasome, tissue damage, and monocyte/macrophage infiltration. Previously, we documented that caspase-1 activation in myeloid cells was linked with HIV-associated atherosclerosis in mice and people with HIV. Here, we mechanistically examined the direct effect of caspase-1 on HIV-associated atherosclerosis. Caspase-1-deficient (Casp-1-/-) mice were crossed with HIV-1 transgenic (Tg26+/-) mice with an atherogenic ApoE-deficient (ApoE-/-) background to create global caspase-1-deficient mice (Tg26+/-/ApoE-/-/Casp-1-/-). Caspase-1-sufficient (Tg26+/-/ApoE-/-/Casp-1+/+) mice served as the controls. Next, we created chimeric hematopoietic cell-deficient mice by reconstituting irradiated ApoE-/- mice with bone marrow cells transplanted from Tg26+/-/ApoE-/-/Casp-1-/- (BMT Casp-1-/-) or Tg26+/-/ApoE-/-/Casp-1+/+ (BMT Casp-1+/+) mice. Global caspase-1 knockout in mice suppressed plaque deposition in the thoracic aorta, serum IL-18 levels, and ex vivo foam cell formation. The deficiency of caspase-1 in hematopoietic cells resulted in reduced atherosclerotic plaque burden in the whole aorta and aortic root, which was associated with reduced macrophage infiltration. Transcriptomic analyses of peripheral mononuclear cells and splenocytes indicated that caspase-1 deficiency inhibited caspase-1 pathway-related genes. These results document the critical atherogenic role of caspase-1 in chronic HIV infection and highlight the implication of this pathway and peripheral immune activation in HIV-associated atherosclerosis.

i6mA-Fuse: improved and robust prediction of DNA 6 mA sites in the Rosaceae genome by fusing multiple feature representation.

DNA N6-methyladenine (6 mA) is one of the most vital epigenetic modifications and involved in controlling the various gene expression levels. With the avalanche of DNA sequences generated in numerous databases, the accurate identification of 6 mA plays an essential role for understanding molecular mechanisms. Because the experimental approaches are time-consuming and costly, it is desirable to develop a computation model for rapidly and accurately identifying 6 mA. To the best of our knowledge, we first proposed a computational model named i6mA-Fuse to predict 6 mA sites from the Rosaceae genomes, especially in Rosa chinensis and Fragaria vesca. We implemented the five encoding schemes, i.e., mononucleotide binary, dinucleotide binary, k-space spectral nucleotide, k-mer, and electron-ion interaction pseudo potential compositions, to build the five, single-encoding random forest (RF) models. The i6mA-Fuse uses a linear regression model to combine the predicted probability scores of the five, single encoding-based RF models. The resultant species-specific i6mA-Fuse achieved remarkably high performances with AUCs of 0.982 and 0.978 and with MCCs of 0.869 and 0.858 on the independent datasets of Rosa chinensis and Fragaria vesca, respectively. In the F. vesca-specific i6mA-Fuse, the MBE and EIIP contributed to 75% and 25% of the total prediction; in the R. chinensis-specific i6mA-Fuse, Kmer, MBE, and EIIP contribute to 15%, 65%, and 20% of the total prediction. To assist high-throughput prediction for DNA 6 mA identification, the i6mA-Fuse is publicly accessible at https://kurata14.bio.kyutech.ac.jp/i6mA-Fuse/.

ProIn-Fuse: improved and robust prediction of proinflammatory peptides by fusing of multiple feature representations.

A proinflammatory peptide (PIP) is a type of signaling molecules that are secreted from immune cells, which contributes to the first line of defense against invading pathogens. Numerous experiments have shown that PIPs play an important role in human physiology such as vaccines and immunotherapeutic drugs. Considering high-throughput laboratory methods that are time consuming and costly, effective computational methods are great demand to timely and accurately identify PIPs. Thus, in this study, we proposed a computational model in conjunction with a multiple feature representation, called ProIn-Fuse, to improve the performance of PIPs identification. Specifically, a feature representation learning model was utilized to generate the probabilistic scores by using the random forest models employing eight sequence encoding schemes. Finally, the ProIn-Fuse was constructed by linearly combining the resultant eight probabilistic scores. Evaluated through independent test, the ProIn-Fuse yielded an accuracy of 0.746, which was 10% higher than those obtained by the state-of-the-art PIP predictors. The proposed ProIn-Fuse can facilitate faster and broader applications of PIPs in drug design and development. The web server, datasets and online instruction are freely accessible at http://kurata14.bio.kyutech.ac.jp/ProIn-Fuse .

Evolution of Sequence-based Bioinformatics Tools for Protein-protein Interaction Prediction.

Protein-protein interactions (PPIs) are the physical connections between two or more proteins via electrostatic forces or hydrophobic effects. Identification of the PPIs is pivotal, which contributes to many biological processes including protein function, disease incidence, and therapy design. The experimental identification of PPIs via high-throughput technology is time-consuming and expensive. Bioinformatics approaches are expected to solve such restrictions. In this review, our main goal is to provide an inclusive view of the existing sequence-based computational prediction of PPIs. Initially, we briefly introduce the currently available PPI databases and then review the state-of-the-art bioinformatics approaches, working principles, and their performances. Finally, we discuss the caveats and future perspective of the next generation algorithms for the prediction of PPIs.

NTyroSite: Computational Identification of Protein Nitrotyrosine Sites Using Sequence Evolutionary Features.

Nitrotyrosine is a product of tyrosine nitration mediated by reactive nitrogen species. As an indicator of cell damage and inflammation, protein nitrotyrosine serves to reveal biological change associated with various diseases or oxidative stress. Accurate identification of nitrotyrosine site provides the important foundation for further elucidating the mechanism of protein nitrotyrosination. However, experimental identification of nitrotyrosine sites through traditional methods are laborious and expensive. In silico prediction of nitrotyrosine sites based on protein sequence information are thus highly desired. Here, we report a novel predictor, NTyroSite, for accurate prediction of nitrotyrosine sites using sequence evolutionary information. The generated features were optimized using a Wilcoxon-rank sum test. A random forest classifier was then trained using these features to build the predictor. The final NTyroSite predictor achieved an area under a receiver operating characteristics curve (AUC) score of 0.904 in a 10-fold cross-validation test. It also significantly outperformed other existing implementations in an independent test. Meanwhile, for a better understanding of our prediction model, the predominant rules and informative features were extracted from the NTyroSite model to explain the prediction results. We expect that the NTyroSite predictor may serve as a useful computational resource for high-throughput nitrotyrosine site prediction. The online interface of the software is publicly available at https://biocomputer.bio.cuhk.edu.hk/NTyroSite/.

Precision-cut lung slices as an ex vivo model to study Pneumocystis murina survival and antimicrobial susceptibility.

IMPORTANCE: Our study reveals the potential of precision-cut lung slices as an ex vivo platform to study the growth/survival of Pneumocystis spp. that can facilitate the development of new anti-fungal drugs.

Deficiency of Tlr7 and Irf7 in mice increases the severity of COVID-19 through the reduced interferon production.

Toll-like receptor 7 (Tlr7) deficiency-accelerated severe COVID-19 is associated with reduced production of interferons (IFNs). However, the underlying mechanisms remain elusive. To address these questions, we utilize Tlr7 and Irf7 deficiency mice, single-cell RNA analysis together with bone marrow transplantation approaches. We demonstrate that at the early phase of infection, SARS-CoV-2 causes the upregulation of Tlr7, Irf7, and IFN pathways in the lungs of the infected mice. The deficiency of Tlr7 and Irf7 globally and/or in immune cells in mice increases the severity of COVID-19 via impaired IFN activation in both immune and/or non-immune cells, leading to increased lung viral loads. These effects are associated with reduced IFN alpha and gamma levels in the circulation. The deficiency of Tlr7 tends to cause the reduced production and nuclear translocation of interferon regulatory factor 7 (IRF7) in the lungs of the infected mice, indicative of reduced IRF7 activation. Despite higher amounts of lung viral antigen, Tlr7 or Irf7 deficiency resulted in substantially reduced production of antibodies against SARS-CoV-2, thereby delaying the viral clearance. These results highlight the importance of the activation of TLR7 and IRF7 leading to IFN production on the development of innate and adaptive immunity against COVID-19.

Upregulation of inflammatory genes and pathways links obesity to severe COVID-19.

Obesity is a risk factor for developing severe COVID-19. However, the mechanism underlying obesity-accelerated COVID-19 remains unclear. Here, we report results from a study in which 2-3-month-old K18-hACE2 (K18) mice were fed a western high-fat diet (WD) or normal chow (NC) over 3 months before intranasal infection with a sublethal dose of SARS-CoV2 WA1 (a strain ancestral to the Wuhan variant). After infection, the WD-fed K18 mice lost significantly more body weight and had more severe lung inflammation than normal chow (NC)-fed mice. Bulk RNA-seq analysis of lungs and adipose tissue revealed a diverse landscape of various immune cells, inflammatory markers, and pathways upregulated in the infected WD-fed K18 mice when compared with the infected NC-fed control mice. The transcript levels of IL-6, an important marker of COVID-19 disease severity, were upregulated in the lung at 6-9 days post-infection in the WD-fed mice when compared to NC-fed mice. Transcriptome analysis of the lung and adipose tissue obtained from deceased COVID-19 patients found that the obese patients had an increase in the expression of genes and the activation of pathways associated with inflammation as compared to normal-weight patients (n = 2). The K18 mouse model and human COVID-19 patient data support a link between inflammation and an obesity-accelerated COVID-19 disease phenotype. These results also indicate that obesity-accelerated severe COVID-19 caused by SARS-CoV-2 WA1 infection in the K18 mouse model would be a suitable model for dissecting the cellular and molecular mechanisms underlying pathogenesis.

Cell-intrinsic and -extrinsic effects of SARS-CoV-2 RNA on pathogenesis: single-cell meta-analysis.

Single-cell RNA-seq has been used to characterize human COVID-19. To determine if preclinical models successfully mimic the cell-intrinsic and -extrinsic effects of severe disease, we conducted a meta-analysis of single-cell data across five model species. To assess whether dissemination of viral RNA in lung cells tracks pathology and results in cell-intrinsic and -extrinsic transcriptomic changes in COVID-19. We conducted a meta-analysis by analyzing six publicly available, scRNA-seq data sets. We used dual mapping (host and virus) and differential gene expression analyses to compare viral+ and viral- cell populations. We conducted a principal component analysis to identify successful models of human COVID-19. We found expression of viral RNA in many non-epithelial cell types. Fibroblasts, macrophages, and endothelial cells exhibit clear evidence of viral-intrinsic and -extrinsic effects on host gene expression. Using viral RNA expression, we found that K18-hACE2 mice most closely modeled severe human COVID-19, followed by hamsters. Ferrets and macaques are poor models of human disease due to the low presence of viral RNA. Moreover, we found that increased transcripts of certain key inflammatory genes such as IL1B, IL18, and CXCL10 are not restricted to virally infected cells, suggesting these genes are regulated in a paracrine or autocrine fashion. These data affirm widespread dissemination of viral RNA in the lung, which may be key in the pathogenesis of severe COVID-19 and demonstrate ferrets and Rhesus macaques are poor models of human COVID-19. IMPORTANCE We conducted a high-resolution meta-analysis of scRNA-seq data from humans and five animal models of COVID-19. This study reports viral RNA dissemination in several cell types in human data as well as in some of the pre-clinical models. Using this metric, the K18-hACE2 mouse model, followed by the hamster model, most closely resembled human COVID-19. We observed clear evidence of viral-intrinsic effects within cells (e.g., IRF5 expression) as well as viral-extrinsic cytokine modulation (e.g., IL1B, IL18, CXCL10). We observed proinflammatory chemokine expression in cells devoid of viral RNA expression, suggesting autocrine/paracrine interferon regulation. This report serves as a resource-synthesizing data from COVID-19 humans and animal models and suggesting improvements for relevant pre-clinical models that may aid future diagnostic and therapeutic development projects.

Novel Pneumocystis Antigens for Seroprevalence Studies.

Pneumocystis jirovecii is the most common cause of fungal pneumonia in children under the age of 2 years. However, the inability to culture and propagate this organism has hampered the acquisition of a fungal genome as well as the development of recombinant antigens to conduct seroprevalence studies. In this study, we performed proteomics on Pneumocystis-infected mice and used the recent P. murina and P. jirovecii genomes to prioritize antigens for recombinant protein expression. We focused on a fungal glucanase due to its conservation among fungal species. We found evidence of maternal IgG to this antigen, followed by a nadir in pediatric samples between 1 and 3 months of age, followed by an increase in prevalence over time consistent with the known epidemiology of Pneumocystis exposure. Moreover, there was a strong concordance of anti-glucanase responses and IgG against another Pneumocystis antigen, PNEG_01454. Taken together, these antigens may be useful tools for Pneumocystis seroprevalence and seroconversion studies.

COVID-19 and influenza infections mediate distinct pulmonary cellular and transcriptomic changes.

SARS-CoV-2 infection can cause persistent respiratory sequelae. However, the underlying mechanisms remain unclear. Here we report that sub-lethally infected K18-human ACE2 mice show patchy pneumonia associated with histiocytic inflammation and collagen deposition at 21 and 45 days post infection (DPI). Transcriptomic analyses revealed that compared to influenza-infected mice, SARS-CoV-2-infected mice had reduced interferon-gamma/alpha responses at 4 DPI and failed to induce keratin 5 (Krt5) at 6 DPI in lung, a marker of nascent pulmonary progenitor cells. Histologically, influenza- but not SARS-CoV-2-infected mice showed extensive Krt5+ "pods" structure co-stained with stem cell markers Trp63/NGFR proliferated in the pulmonary consolidation area at both 7 and 14 DPI, with regression at 21 DPI. These Krt5+ "pods" structures were not observed in the lungs of SARS-CoV-2-infected humans or nonhuman primates. These results suggest that SARS-CoV-2 infection fails to induce nascent Krt5+ cell proliferation in consolidated regions, leading to incomplete repair of the injured lung.

A Comprehensive Comparative Review of Protein Sequence-Based Computational Prediction Models of Lysine Succinylation Sites.

Lysine succinylation is a post-translational modification (PTM) of protein in which a succinyl group (-CO-CH2-CH2-CO2H) is added to a lysine residue of protein that reverses lysine's positive charge to a negative charge and leads to the significant changes in protein structure and function. It occurs on a wide range of proteins and plays an important role in various cellular and biological processes in both eukaryotes and prokaryotes. Beyond experimentally identified succinylation sites, there have been a lot of studies for developing sequence-based prediction using machine learning approaches, because it has the promise of being extremely time-saving, accurate, robust, and cost-effective. Despite these benefits for computational prediction of lysine succinylation sites for different species, there are a number of issues that need to be addressed in the design and development of succinylation site predictors. In spite of the fact that many studies used different statistical and machine learning computational tools, only a few studies have focused on these bioinformatics issues in depth. Therefore, in this comprehensive comparative review, an attempt is made to present the latest advances in the prediction models, datasets, and online resources, as well as the obstacles and limits, to provide an advantageous guideline for developing more suitable and effective succinylation site prediction tools.

Recent Development of Bioinformatics Tools for microRNA Target Prediction.

MicroRNAs (miRNAs) are central players that regulate the post-transcriptional processes of gene expression. Binding of miRNAs to target mRNAs can repress their translation by inducing the degradation or by inhibiting the translation of the target mRNAs. Highthroughput experimental approaches for miRNA target identification are costly and timeconsuming, depending on various factors. It is vitally important to develop bioinformatics methods for accurately predicting miRNA targets. With the increase of RNA sequences in the post-genomic era, bioinformatics methods are being developed for miRNA studies especially for miRNA target prediction. This review summarizes the current development of state-of-the-art bioinformatics tools for miRNA target prediction, points out the progress and limitations of the available miRNA databases, and their working principles. Finally, we discuss the caveat and perspectives of the next-generation algorithms for the prediction of miRNA targets.

Prediction of serine phosphorylation sites mapping on Schizosaccharomyces Pombe by fusing three encoding schemes with the random forest classifier.

Serine phosphorylation is one type of protein post-translational modifications (PTMs), which plays an essential role in various cellular processes and disease pathogenesis. Numerous methods are used for the prediction of phosphorylation sites. However, the traditional wet-lab based experimental approaches are time-consuming, laborious, and expensive. In this work, a computational predictor was proposed to predict serine phosphorylation sites mapping on Schizosaccharomyces pombe (SP) by the fusion of three encoding schemes namely k-spaced amino acid pair composition (CKSAAP), binary and amino acid composition (AAC) with the random forest (RF) classifier. So far, the proposed method is firstly developed to predict serine phosphorylation sites for SP. Both the training and independent test performance scores were used to investigate the success of the proposed RF based fusion prediction model compared to others. We also investigated their performances by 5-fold cross-validation (CV). In all cases, it was observed that the recommended predictor achieves the largest scores of true positive rate (TPR), true negative rate (TNR), accuracy (ACC), Mathew coefficient of correlation (MCC), Area under the ROC curve (AUC) and pAUC (partial AUC) at false positive rate (FPR) = 0.20. Thus, the prediction performance as discussed in this paper indicates that the proposed approach may be a beneficial and motivating computational resource for predicting serine phosphorylation sites in the case of Fungi. The online interface of the software for the proposed prediction model is publicly available at http://mollah-bioinformaticslab-stat.ru.ac.bd/PredSPS/ .

iLBE for Computational Identification of Linear B-cell Epitopes by Integrating Sequence and Evolutionary Features.

Linear B-cell epitopes are critically important for immunological applications, such as vaccine design, immunodiagnostic test, and antibody production, as well as disease diagnosis and therapy. The accurate identification of linear B-cell epitopes remains challenging despite several decades of research. In this work, we have developed a novel predictor, Identification of Linear B-cell Epitope (iLBE), by integrating evolutionary and sequence-based features. The successive feature vectors were optimized by a Wilcoxon-rank sum test. Then the random forest (RF) algorithm using the optimal consecutive feature vectors was applied to predict linear B-cell epitopes. We combined the RF scores by the logistic regression to enhance the prediction accuracy. iLBE yielded an area under curve score of 0.809 on the training dataset and outperformed other prediction models on a comprehensive independent dataset. iLBE is a powerful computational tool to identify the linear B-cell epitopes and would help to develop penetrating diagnostic tests. A web application with curated datasets for iLBE is freely accessible at http://kurata14.bio.kyutech.ac.jp/iLBE/.

i4mC-ROSE, a bioinformatics tool for the identification of DNA N4-methylcytosine sites in the Rosaceae genome.

One of the most important epigenetic modifications is N4-methylcytosine, which regulates many biological processes including DNA replication and chromosome stability. Identification of N4-methylcytosine sites is pivotal to understand specific biological functions. Herein, we developed the first bioinformatics tool called i4mC-ROSE for identifying N4-methylcytosine sites in the genomes of Fragaria vesca and Rosa chinensis in the Rosaceae, which utilizes a random forest classifier with six encoding methods that cover various aspects of DNA sequence information. The i4mC-ROSE predictor achieves area under the curve scores of 0.883 and 0.889 for the two genomes during cross-validation. Moreover, the i4mC-ROSE outperforms other classifiers tested in this study when objectively evaluated on the independent datasets. The proposed i4mC-ROSE tool can serve users' demand for the prediction of 4mC sites in the Rosaceae genome. The i4mC-ROSE predictor and utilized datasets are publicly accessible at http://kurata14.bio.kyutech.ac.jp/i4mC-ROSE/.