Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection.

Identification of host genes essential for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may reveal novel therapeutic targets and inform our understanding of coronavirus disease 2019 (COVID-19) pathogenesis. Here we performed genome-wide CRISPR screens in Vero-E6 cells with SARS-CoV-2, Middle East respiratory syndrome CoV (MERS-CoV), bat CoV HKU5 expressing the SARS-CoV-1 spike, and vesicular stomatitis virus (VSV) expressing the SARS-CoV-2 spike. We identified known SARS-CoV-2 host factors, including the receptor ACE2 and protease Cathepsin L. We additionally discovered pro-viral genes and pathways, including HMGB1 and the SWI/SNF chromatin remodeling complex, that are SARS lineage and pan-coronavirus specific, respectively. We show that HMGB1 regulates ACE2 expression and is critical for entry of SARS-CoV-2, SARS-CoV-1, and NL63. We also show that small-molecule antagonists of identified gene products inhibited SARS-CoV-2 infection in monkey and human cells, demonstrating the conserved role of these genetic hits across species. This identifies potential therapeutic targets for SARS-CoV-2 and reveals SARS lineage-specific and pan-CoV host factors that regulate susceptibility to highly pathogenic CoVs.

A tryptophan metabolite made by a gut microbiome eukaryote induces pro-inflammatory T cells.

The large intestine harbors microorganisms playing unique roles in host physiology. The beneficial or detrimental outcome of host-microbiome coexistence depends largely on the balance between regulators and responder intestinal CD4+ T cells. We found that ulcerative colitis-like changes in the large intestine after infection with the protist Blastocystis ST7 in a mouse model are associated with reduction of anti-inflammatory Treg cells and simultaneous expansion of pro-inflammatory Th17 responders. These alterations in CD4+ T cells depended on the tryptophan metabolite indole-3-acetaldehyde (I3AA) produced by this single-cell eukaryote. I3AA reduced the Treg subset in vivo and iTreg development in vitro by modifying their sensing of TGFß, concomitantly affecting recognition of self-flora antigens by conventional CD4+ T cells. Parasite-derived I3AA also induces over-exuberant TCR signaling, manifested by increased CD69 expression and downregulation of co-inhibitor PD-1. We have thus identified a new mechanism dictating CD4+ fate decisions. The findings thus shine a new light on the ability of the protist microbiome and tryptophan metabolites, derived from them or other sources, to modulate the adaptive immune compartment, particularly in the context of gut inflammatory disorders.

KDM5B promotes immune evasion by recruiting SETDB1 to silence retroelements.

Tumours use various strategies to evade immune surveillance1,2. Immunotherapies targeting tumour immune evasion such as immune checkpoint blockade have shown considerable efficacy on multiple cancers3,4 but are ineffective for most patients due to primary or acquired resistance5-7. Recent studies showed that some epigenetic regulators suppress anti-tumour immunity2,8-12, suggesting that epigenetic therapies could boost anti-tumour immune responses and overcome resistance to current immunotherapies. Here we show that, in mouse melanoma models, depletion of KDM5B-an H3K4 demethylase that is critical for melanoma maintenance and drug resistance13-15-induces robust adaptive immune responses and enhances responses to immune checkpoint blockade. Mechanistically, KDM5B recruits the H3K9 methyltransferase SETDB1 to repress endogenous retroelements such as MMVL30 in a demethylase-independent manner. Derepression of these retroelements activates cytosolic RNA-sensing and DNA-sensing pathways and the subsequent type-I interferon response, leading to tumour rejection and induction of immune memory. Our results demonstrate that KDM5B suppresses anti-tumour immunity by epigenetic silencing of retroelements. We therefore reveal roles of KDM5B in heterochromatin regulation and immune evasion in melanoma, opening new paths for the development of KDM5B-targeting and SETDB1-targeting therapies to enhance tumour immunogenicity and overcome immunotherapy resistance.

KSHV vIL-6 promotes SIRT3-induced deacetylation of SERBP1 to inhibit ferroptosis and enhance cellular transformation by inducing lipoyltransferase 2 mRNA degradation.

Ferroptosis, a defensive strategy commonly employed by the host cells to restrict pathogenic infections, has been implicated in the development and therapeutic responses of various types of cancer. However, the role of ferroptosis in oncogenic Kaposi's sarcoma-associated herpesvirus (KSHV)-induced cancers remains elusive. While a growing number of non-histone proteins have been identified as acetylation targets, the functions of these modifications have yet to be revealed. Here, we show KSHV reprogramming of host acetylation proteomics following cellular transformation of rat primary mesenchymal precursor. Among them, SERPINE1 mRNA binding protein 1 (SERBP1) deacetylation is increased and required for KSHV-induced cellular transformation. Mechanistically, KSHV-encoded viral interleukin-6 (vIL-6) promotes SIRT3 deacetylation of SERBP1, preventing its binding to and protection of lipoyltransferase 2 (Lipt2) mRNA from mRNA degradation resulting in ferroptosis. Consequently, a SIRT3-specific inhibitor, 3-TYP, suppresses KSHV-induced cellular transformation by inducing ferroptosis. Our findings unveil novel roles of vIL-6 and SERBP1 deacetylation in regulating ferroptosis and KSHV-induced cellular transformation, and establish the vIL-6-SIRT3-SERBP1-ferroptosis pathways as a potential new therapeutic target for KSHV-associated cancers.

DYRK1A promotes viral entry of highly pathogenic human coronaviruses in a kinase-independent manner.

Identifying host genes essential for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has the potential to reveal novel drug targets and further our understanding of Coronavirus Disease 2019 (COVID-19). We previously performed a genome-wide CRISPR/Cas9 screen to identify proviral host factors for highly pathogenic human coronaviruses. Few host factors were required by diverse coronaviruses across multiple cell types, but DYRK1A was one such exception. Although its role in coronavirus infection was previously undescribed, DYRK1A encodes Dual Specificity Tyrosine Phosphorylation Regulated Kinase 1A and is known to regulate cell proliferation and neuronal development. Here, we demonstrate that DYRK1A regulates ACE2 and DPP4 transcription independent of its catalytic kinase function to support SARS-CoV, SARS-CoV-2, and Middle East Respiratory Syndrome Coronavirus (MERS-CoV) entry. We show that DYRK1A promotes DNA accessibility at the ACE2 promoter and a putative distal enhancer, facilitating transcription and gene expression. Finally, we validate that the proviral activity of DYRK1A is conserved across species using cells of nonhuman primate and human origin. In summary, we report that DYRK1A is a novel regulator of ACE2 and DPP4 expression that may dictate susceptibility to multiple highly pathogenic human coronaviruses.

5-Fluorouracil efficacy requires anti-tumor immunity triggered by cancer-cell-intrinsic STING.

5-Fluorouracil (5-FU) is a widely used chemotherapeutic drug, but the mechanisms underlying 5-FU efficacy in immunocompetent hosts in vivo remain largely elusive. Through modeling 5-FU response of murine colon and melanoma tumors, we report that effective reduction of tumor burden by 5-FU is dependent on anti-tumor immunity triggered by the activation of cancer-cell-intrinsic STING. While the loss of STING does not induce 5-FU resistance in vitro, effective 5-FU responsiveness in vivo requires cancer-cell-intrinsic cGAS, STING, and subsequent type I interferon (IFN) production, as well as IFN-sensing by bone-marrow-derived cells. In the absence of cancer-cell-intrinsic STING, a much higher dose of 5-FU is needed to reduce tumor burden. 5-FU treatment leads to increased intratumoral T cells, and T-cell depletion significantly reduces the efficacy of 5-FU in vivo. In human colorectal specimens, higher STING expression is associated with better survival and responsiveness to chemotherapy. Our results support a model in which 5-FU triggers cancer-cell-initiated anti-tumor immunity to reduce tumor burden, and our findings could be harnessed to improve therapeutic effectiveness and toxicity for colon and other cancers.

Immune evasion strategy involving propionylation by the KSHV interferon regulatory factor 1 (vIRF1).

Post-translational modifications (PTMs) are essential for host antiviral immune response and viral immune evasion. Among a set of novel acylations, lysine propionylation (Kpr) has been detected in both histone and non-histone proteins. However, whether protein propionylation occurs in any viral proteins and whether such modifications regulate viral immune evasion remain elusive. Here, we show that Kaposi's sarcoma-associated herpesvirus (KSHV)-encoded viral interferon regulatory factor 1 (vIRF1) can be propionylated in lysine residues, which is required for effective inhibition of IFN-ß production and antiviral signaling. Mechanistically, vIRF1 promotes its own propionylation by blocking SIRT6's interaction with ubiquitin-specific peptidase 10 (USP10) leading to its degradation via a ubiquitin-proteasome pathway. Furthermore, vIRF1 propionylation is required for its function to block IRF3-CBP/p300 recruitment and repress the STING DNA sensing pathway. A SIRT6-specific activator, UBCS039, rescues propionylated vIRF1-mediated repression of IFN-ß signaling. These results reveal a novel mechanism of viral evasion of innate immunity through propionylation of a viral protein. The findings suggest that enzymes involved in viral propionylation could be potential targets for preventing viral infections.

The KDM6A-KMT2D-p300 axis regulates susceptibility to diverse coronaviruses by mediating viral receptor expression.

Identification of host determinants of coronavirus infection informs mechanisms of pathogenesis and may provide novel therapeutic targets. Here, we demonstrate that the histone demethylase KDM6A promotes infection of diverse coronaviruses, including SARS-CoV, SARS-CoV-2, MERS-CoV and mouse hepatitis virus (MHV) in a demethylase activity-independent manner. Mechanistic studies reveal that KDM6A promotes viral entry by regulating expression of multiple coronavirus receptors, including ACE2, DPP4 and Ceacam1. Importantly, the TPR domain of KDM6A is required for recruitment of the histone methyltransferase KMT2D and histone deacetylase p300. Together this KDM6A-KMT2D-p300 complex localizes to the proximal and distal enhancers of ACE2 and regulates receptor expression. Notably, small molecule inhibition of p300 catalytic activity abrogates ACE2 and DPP4 expression and confers resistance to all major SARS-CoV-2 variants and MERS-CoV in primary human airway and intestinal epithelial cells. These data highlight the role for KDM6A-KMT2D-p300 complex activities in conferring diverse coronaviruses susceptibility and reveal a potential pan-coronavirus therapeutic target to combat current and emerging coronaviruses. One Sentence Summary: The KDM6A/KMT2D/EP300 axis promotes expression of multiple viral receptors and represents a potential drug target for diverse coronaviruses.

Systematic evaluation of multiple NGS platforms for structural variants detection.

Structural variations (SV) are critical genome changes affecting human diseases. Although many hybridization-based methods exist, evaluating SVs through next-generation sequencing (NGS) data is still necessary for broader research exploration. Here, we comprehensively compared the performance of 16 SV callers and multiple NGS platforms using NA12878 whole genome sequencing (WGS) datasets. The results indicated that several SV callers performed well relatively, such as Manta, GRIDSS, LUMPY, TARDIS, FermiKit, and Wham. Meanwhile, all NGS platforms have a similar performance using a single software. Additionally, we found that the source of undetected SVs was mostly from long reads datasets, therefore, the more appropriate strategy for accurate SV detection will be an integration of long and shorter reads in the future. At present, in the period of NGS as a mainstream method in bioinformatics, our study would provide helpful and comprehensive guidelines for specific categories of SV research.

Epigenetic markers and therapeutic targets for metastasis.

The last few years have seen an increasing number of discoveries which collectively demonstrate that histone and DNA modifying enzyme modulate different stages of metastasis. Moreover, epigenomic alterations can now be measured at multiple scales of analysis and are detectable in human tumors or liquid biopsies. Malignant cell clones with a proclivity for relapse in certain organs may arise in the primary tumor as a consequence of epigenomic alterations which cause a loss in lineage integrity. These alterations may occur due to genetic aberrations acquired during tumor progression or concomitant to therapeutic response. Moreover, evolution of the stroma can also alter the epigenome of cancer cells. In this review, we highlight current knowledge with a particular emphasis on leveraging chromatin and DNA modifying mechanisms as biomarkers of disseminated disease and as therapeutic targets to treat metastatic cancers.

The C4 photosynthesis bifunctional enzymes, PDRPs, of maize are co-opted to cytoplasmic viral replication complexes to promote infection of a prevalent potyvirus sugarcane mosaic virus.

In maize, two pyruvate orthophosphate dikinase (PPDK) regulatory proteins, ZmPDRP1 and ZmPDRP2, are respectively specific to the chloroplast of mesophyll cells (MCs) and bundle sheath cells (BSCs). Functionally, ZmPDRP1/2 catalyse both phosphorylation/inactivation and dephosphorylation/activation of ZmPPDK, which is implicated as a major rate-limiting enzyme in C4 photosynthesis of maize. Our study here showed that maize plants lacking ZmPDRP1 or silencing of ZmPDRP1/2 confer resistance to a prevalent potyvirus sugarcane mosaic virus (SCMV). We verified that the C-terminal domain (CTD) of ZmPDRP1 plays a key role in promoting viral infection while independent of enzyme activity. Intriguingly, ZmPDRP1 and ZmPDRP2 re-localize to cytoplasmic viral replication complexes (VRCs) following SCMV infection. We identified that SCMV-encoded cytoplasmic inclusions protein CI targets directly ZmPDRP1 or ZmPDRP2 or their CTDs, leading to their re-localization to cytoplasmic VRCs. Moreover, we found that CI could be degraded by the 26S proteasome system, while ZmPDRP1 and ZmPDRP2 could up-regulate the accumulation level of CI through their CTDs by a yet unknown mechanism. Most importantly, with genetic, cell biological and biochemical approaches, we provide evidence that BSCs-specific ZmPDRP2 could accumulate in MCs of Zmpdrp1 knockout (KO) lines, revealing a unique regulatory mechanism crossing different cell types to maintain balanced ZmPPDK phosphorylation, thereby to keep maize normal growth. Together, our findings uncover the genetic link of the two cell-specific maize PDRPs, both of which are co-opted to VRCs to promote viral protein accumulation for robust virus infection.

Correction: Sperm associated antigen 9 promotes oncogenic KSHV-encoded interferon regulatory factor-induced cellular transformation and angiogenesis by activating the JNK/VEGFA pathway.

[This corrects the article DOI: 10.1371/journal.ppat.1008730.].

Alginate oligosaccharide assimilation by gut microorganisms and the potential role in gut inflammation alleviation.

Dietary fiber metabolism by gut microorganisms plays important roles in host physiology and health. Alginate, the major dietary fiber of daily diet seaweeds, is drawing more attention because of multiple biological activities. To advance the understanding of alginate assimilation mechanism in the gut, we show the presence of unsaturated alginate oligosaccharides (uAOS)-specific alginate utilization loci (AUL) in human gut microbiome. As a representative example, a working model of the AUL from the gut microorganism Bacteroides clarus was reconstructed from biochemistry and transcriptome data. The fermentation of resulting monosaccharides through Entner-Doudoroff pathway tunes the metabolism of short-chain fatty acids and amino acids. Furthermore, we show that uAOS feeding protects the mice against dextran sulfate sodium-induced acute colitis probably by remodeling gut microbiota and metabolome. IMPORTANCE: Alginate has been included in traditional Chinese medicine and daily diet for centuries. Recently discovered biological activities suggested that alginate-derived alginate oligosaccharides (AOS) might be an active ingredient in traditional Chinese medicine, but how these AOS are metabolized in the gut and how it affects health need more information. The study on the working mechanism of alginate utilization loci (AUL) by the gut microorganism uncovers the role of unsaturated alginate oligosaccharides (uAOS) assimilation in tuning short-chain fatty acids and amino acids metabolism and demonstrates that uAOS metabolism by gut microorganisms results in a variation of cell metabolites, which potentially contributes to the physiology and health of gut.

Production of recombinant HPV11/16 E6/E7-MBP-His6 fusion proteins and their potential to induce cytokine secretion by immune cells in peripheral blood.

Human papillomavirus (HPV) infection poses a significant threat to public health worldwide. Targeting the function of HPV E6 and E7 proteins and activating the host immune response against these proteins represent promising therapeutic strategies for combating HPV-related diseases. Consequently, the efficient production of soluble, high-purity E6 and E7 proteins is crucial for function and host immune response studies. In this context, we selected the pMCSG19 protein expression vector for Escherichia coli to produce soluble MBP-His6 tagged HPV11/16 E6/E7 proteins, achieving relatively high purity and yield. Notably, these proteins exhibited low toxicity to peripheral blood mononuclear cells (PBMCs) and did not compromise their viability. Additionally, the recombinant proteins were capable of inducing the secretion of multiple cytokines by immune cells in peripheral blood, indicating their potential to elicit immune responses. In conclusion, our study offers a novel approach for the production of HPV11/16 E6/E7 fusion proteins with relatively high purity and yield. The fusing HPV11/16 E6/E7 proteins to MBP-His6 tag may serve as a valuable method for large-scale protein production in future research endeavors.

Purinergic P2X7 receptor involves in anti-retinal photodamage effects of berberine.

Berberine (BBR) is a Chinese herb with antioxidant and anti-inflammatory properties. In a previous study, we found that BBR had a protective effect against light-induced retinal degeneration in BALB/c mice. The purinergic P2X7 receptor (P2X7R) plays a key role in retinal degeneration via inducing oxidative stress, inflammatory changes, and cell death. The aim of this study was to investigate whether BBR can induce protective effects in light damage experiments and whether P2X7R can get involved in these effects. C57BL/6 J mice and P2X7 knockout (KO) mice on the C57BL/6 J background were used. We found that BBR preserved the outer nuclear layer (ONL) thickness and retinal ganglion cells following light stimulation. Furthermore, BBR significantly suppressed photoreceptor apoptosis, pro-apoptotic c-fos expression, pro-inflammatory responses of Mϋller cells, and inflammatory factors (TNF-α, IL-1ß). In addition, protein levels of P2X7R were downregulated in BBR-treated mice. Double immunofluorescence showed that BBR reduced overexpression of P2X7R in retinal ganglion cells and Mϋller cells. Furthermore, BBR combined with the P2X7R agonist BzATP blocked the effects of BBR on retinal morphology and photoreceptor apoptosis. However, in P2X7 KO mice, BBR had an additive effect resulting in thicker ONL and more photoreceptors. The data suggest that the P2X7 receptor is involved in retinal light damage, and BBR inhibits this process by reducing histological impairment, cell death, and inflammatory responses.

CuCo2O4/CuO Heterostructure with Oxygen Vacancies Induced by Plasma for Electrocatalytic Nitrate Reduction to Ammonia.

Electrochemical nitrate reduction (NO3RR) to ammonia production is regarded as one of the potential alternatives for replacing the Haber-Bosch technology for realizing artificial ammonia synthesis. In this study, a CuCo2O4/CuO-Ar heterostructure in the shape of dandelion nanospheres formed by nanoarrays has been successfully constructed, demonstrating excellent NO3RR performance. Experimental results indicate that Ar plasma etching of CuCo2O4/CuO-Ar significantly increases the content of oxygen vacancies compared to the sample of CuCo2O4/CuO-Air etched by air plasma, resulting in improved NO3RR performance. Density functional theory calculations further confirm that the existence of more oxygen vacancies effectively decreases the energy barrier of nitrate adsorption, which is due to the generation of more oxygen vacancies facilitating nitrate adsorption and weakening the N-O bonds of nitrate after plasma treatment. As a result, CuCo2O4/CuO-Ar exhibits a high NH3 yield of 0.55 mmol h-1 cm-2 and a Faraday efficiency of 95.07% at the optimal potential of -0.9 V (vs RHE) in a neutral medium. Importantly, CuCo2O4/CuO-Ar also showcases excellent electrocatalytic stability. This study presents new views on the design and structure regulation of NO3RR electrocatalysts and their potential applications in the future.

Role of Platelet/Lymphocyte, Neutrophil/Lymphocyte, and Interleukin-37/Interleukin-17 Ratios in the Occurrence and Treatment of Rheumatoid Arthritis.

This study was designed to investigate the correlation of neutrophil/lymphocyte ratio (NLR), platelet/lymphocyte ratio (PLR), and interleukin (IL)-37/IL-17 ratio with the incidence/treatment of rheumatoid arthritis (RA). Firstly, fifty-eight patients with RA treated at the first affiliated hospital of Xinjiang Medical University from January 2018 to January 2019 were selected as the RA group; forty-nine healthy volunteers were enrolled in the control group. RA patients were treated with disease-modifying anti-rheumatic drugs (DMARDs). Next, the NLR, PLR, IL-37, IL-17 and 28-joint disease activity score using erythrocyte sedimentation rate (DAS28-ESR) were deleted in two groups. Subsequently, Spearman correlation analysis was adopted for the correlations of various indicators before and after treatment in two groups. According to the analysis results, the levels of NLR, PLR, IL-37, and IL-17 before treatment in the RA group were higher than those in the control group (P < .05), but the difference in the IL-37/IL-17 level between the two groups was not significant (P > .05). After treatment, NLR, PLR, and IL-37/IL-17 levels were significantly reduced in RA patients (P < .05). NLR and PLR were significantly positively correlated with DAS28-ESR, ESR and C-reactive protein (CRP), of which represented the disease activity of RA. NLP was strongly correlated with IL-37/IL-17. Collectively, NLR, PLR, IL-37, and IL-17 are closely related to the occurrence of RA. In addition, NLR and IL-37/IL-17 are more suitable than PLR in reflecting the therapeutic effect. Therefore, IL-37/IL-17 can be considered as a new indicator for reflecting the treatment effectiveness of RA.

The fungal protease BbAorsin contributes to growth, conidiation, germination, virulence, and antiphytopathogenic activities in Beauveria bassiana (Hypocreales: Cordycipitaceae).

The fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae), is one of the most destructive agricultural pests. The entomopathogenic fungus Beauveria bassiana (Hypocreales: Clavicipitaceae) is a biopesticide widely used for biocontrol of various pests. Secreted fungal proteases are critical for insect cuticle destruction and successful infection. We have previously shown that the serine protease BbAorsin in B. bassiana has entomopathogenic and antiphytopathogenic activities. However, the contribution of BbAorsin to fungal growth, conidiation, germination, virulence and antiphytopathogenic activities remains unclear. In this study, the deletion (ΔBbAorsin), complementation (Comp), and overexpression (BbAorsinOE) strains of B. bassiana were generated for comparative studies. The results showed that ΔBbAorsin exhibited slower growth, reduced conidiation, lower germination rate, and longer germination time compared to WT and Comp. In contrast, BbAorsinOE showed higher growth rate, increased conidiation, higher germination rate and shorter germination time. Injection of BbAorsinOE showed the highest virulence against S. frugiperda larvae, while injection of ΔBbAorsin showed the lowest virulence. Feeding BbAorsinOE resulted in lower pupation and adult eclosion rates and malformed adults. 16S rRNA sequencing revealed no changes in the gut microbiota after feeding either WT or BbAorsinOE. However, BbAorsinOE caused a disrupted midgut, leakage of gut microbiota into the hemolymph, and upregulation of apoptosis and immunity-related genes. BbAorsin can disrupt the cell wall of the phytopathogen Fusarium graminearum and alleviate symptoms in wheat seedlings and cherry tomatoes infected with F. graminearum. These results highlight the importance of BbAorsin for B. bassiana and its potential as a multifunctional biopesticide.

Mechanism of Mulberry Leaves and Black Sesame in Alleviating Slow Transit Constipation Revealed by Multi-Omics Analysis.

Traditional Chinese medicine (TCM) possesses the potential of providing good curative effects with no side effects for the effective management of slow transit constipation (STC), an intestinal disease characterized by colonic dyskinesia. Mulberry leaves (Morus alba L.) and black sesame (Sesamum indicum L.), referred to as SH, are processed and conditioned as per standardized protocols. SH has applications as food and medicine. Accordingly, we investigated the therapeutic potential of SH in alleviating STC. The analysis of SH composition identified a total of 504 compounds. The intervention with SH significantly improved intestinal motility, reduced the time for the first black stool, increased antioxidant activity, and enhanced water content, thereby effectively alleviating colon damage caused by STC. Transcriptome analysis revealed the SH in the treatment of STC related to SOD1, MUC2, and AQP1. The analysis of 16S rRNA gene sequences indicated notable differences in the abundance of 10 bacteria between the SH and model. Metabolomic analysis further revealed that SH supplementation increased the levels of nine metabolites associated with STC. Integrative analysis revealed that SH modulated amino acid metabolism, balanced intestinal flora, and targeted key genes (i.e., SOD1, MUC2, AQP1) to exert its effects. SH also inhibited the AQP1 expression and promoted SOD1 and MUC2 expression.

Systematic and benchmarking studies of pipelines for mammal WGBS data in the novel NGS platform.

BACKGROUND: Whole genome bisulfite sequencing (WGBS), possesses the aptitude to dissect methylation status at the nucleotide-level resolution of 5-methylcytosine (5-mC) on a genome-wide scale. It is a powerful technique for epigenome in various cell types, and tissues. As a recently established next-generation sequencing (NGS) platform, GenoLab M is a promising alternative platform. However, its comprehensive evaluation for WGBS has not been reported. We sequenced two bisulfite-converted mammal DNA in this research using our GenoLab M and NovaSeq 6000, respectively. Then, we systematically compared those data via four widely used WGBS tools (BSMAP, Bismark, BatMeth2, BS-Seeker2) and a new bisulfite-seq tool (BSBolt). We interrogated their computational time, genome depth and coverage, and evaluated their percentage of methylated Cs. RESULT: Here, benchmarking a combination of pre- and post-processing methods, we found that trimming improved the performance of mapping efficiency in eight datasets. The data from two platforms uncovered ~ 80% of CpG sites genome-wide in the human cell line. Those data sequenced by GenoLab M achieved a far lower proportion of duplicates (~ 5.5%). Among pipelines, BSMAP provided an intriguing representation of 5-mC distribution at CpG sites with 5-mC levels > ~ 78% in datasets from human cell lines, especially in the GenoLab M. BSMAP performed more advantages in running time, uniquely mapped reads percentages, genomic coverage, and quantitative accuracy. Finally, compared with the previous methylation pattern of human cell line and mouse tissue, we confirmed that the data from GenoLab M performed similar consistency and accuracy in methylation levels of CpG sites with that from NovaSeq 6000. CONCLUSION: Together we confirmed that GenoLab M was a qualified NGS platform for WGBS with high performance. Our results showed that BSMAP was the suitable pipeline that allowed for WGBS studies on the GenoLab M platform.