Neutrophil diversity is associated with T-cell immunity and clinical relevance in patients with thyroid cancer.

Neutrophil heterogeneity is involved in autoimmune diseases, sepsis, and several cancers. However, the link between neutrophil heterogeneity and T-cell immunity in thyroid cancer is incompletely understood. We investigated the circulating neutrophil heterogeneity in 3 undifferentiated thyroid cancer (UTC), 14 differentiated thyroid cancer (DTC) (4 Stage IV, 10 Stage I-II), and healthy controls (n = 10) by transcriptomic data and cytometry. Participants with UTC had a significantly higher proportion of immature high-density neutrophils (HDN) and lower proportion of mature HDN in peripheral blood compared to DTC. The proportion of circulating PD-L1+ immature neutrophils were significantly increased in advanced cancer patients. Unsupervised analysis of transcriptomics data from circulating HDN revealed downregulation of innate immune response and T-cell receptor signaling pathway in cancer patients. Moreover, UTC patients revealed the upregulation of glycolytic process and glutamate receptor signaling pathway. Comparative analysis across tumor types and stages revealed the downregulation of various T-cell-related pathways, such as T-cell receptor signaling pathway and T-cell proliferation in advanced cancer patients. Moreover, the proportions of CD8+ and CD4+ T effector memory CD45RA+ (TEMRA) cells from peripheral blood were significantly decreased in UTC patients compared to DTC patients. Finally, we demonstrated that proportions of tumor-infiltrated neutrophils were increased and related with poor prognosis in advanced thyroid cancer using data from our RNA-seq and TCGA (The Cancer Genome Atlas) data. In conclusion, observed prevalence of circulating immature high-density neutrophils and their immunosuppressive features in undifferentiated thyroid cancers underscore the importance of understanding neutrophil dynamics in the context of tumor progression in thyroid cancer.

Transcriptomic characteristics according to tumor size and SUVmax in papillary thyroid cancer patients.

The SUVmax is a measure of FDG uptake and is related with tumor aggressiveness in thyroid cancer, however, its association with molecular pathways is unclear. Here, we investigated the relationship between SUVmax and gene expression profiles in 80 papillary thyroid cancer (PTC) patients. We conducted an analysis of DEGs and enriched pathways in relation to SUVmax and tumor size. SUVmax showed a positive correlation with tumor size and correlated with glucose metabolic process. The genes that indicate thyroid differentiation, such as SLC5A5 and TPO, were negatively correlated with SUVmax. Unsupervised analysis revealed that SUVmax positively correlated with DNA replication(r = 0.29, p = 0.009), pyrimidine metabolism(r = 0.50, p < 0.0001) and purine metabolism (r = 0.42, p = 0.0001). Based on subgroups analysis, we identified that PSG5, TFF3, SOX2, SL5A5, SLC5A7, HOXD10, FER1L6, and IFNA1 genes were found to be significantly associated with tumor aggressiveness. Both high SUVmax PTMC and macro-PTC are enriched in pathways of DNA replication and cell cycle, however, gene sets for purine metabolic pathways are enriched only in high SUVmax macro-PTC but not in high SUVmax PTMC. Our findings demonstrate the molecular characteristics of high SUVmax tumor and metabolism involved in tumor growth in differentiated thyroid cancer.

Proteogenomic Characterization Reveals Estrogen Signaling as a Target for Never-Smoker Lung Adenocarcinoma Patients without EGFR or ALK Alterations.

Never-smoker lung adenocarcinoma (NSLA) is prevalent in Asian populations, particularly in women. EGFR mutations and anaplastic lymphoma kinase (ALK) fusions are major genetic alterations observed in NSLA, and NSLA with these alterations have been well studied and can be treated with targeted therapies. To provide insights into the molecular profile of NSLA without EGFR and ALK alterations (NENA), we selected 141 NSLA tissues and performed proteogenomic characterization, including whole genome sequencing (WGS), transcriptomic, methylation EPIC array, total proteomic, and phosphoproteomic analyses. Forty patients with NSLA harboring EGFR and ALK alterations and seven patients with NENA with microsatellite instability were excluded. Genome analysis revealed that TP53 (25%), KRAS (22%), and SETD2 (11%) mutations and ROS1 fusions (14%) were the most frequent genetic alterations in NENA patients. Proteogenomic impact analysis revealed that STK11 and ERBB2 somatic mutations had broad effects on cancer-associated genes in NENA. DNA copy number alteration analysis identified 22 prognostic proteins that influenced transcriptomic and proteomic changes. Gene set enrichment analysis revealed estrogen signaling as the key pathway activated in NENA. Increased estrogen signaling was associated with proteogenomic alterations, such as copy number deletions in chromosomes 14 and 21, STK11 mutation, and DNA hypomethylation of LLGL2 and ST14. Finally, saracatinib, an Src inhibitor, was identified as a potential drug for targeting activated estrogen signaling in NENA and was experimentally validated in vitro. Collectively, this study enhanced our understanding of NENA NSLA by elucidating the proteogenomic landscape and proposed saracatinib as a potential treatment for this patient population that lacks effective targeted therapies. SIGNIFICANCE: The proteogenomic landscape in never-smoker lung cancer without known driver mutations reveals prognostic proteins and enhanced estrogen signaling that can be targeted as a potential therapeutic strategy to improve patient outcomes.

Korea4K: whole genome sequences of 4,157 Koreans with 107 phenotypes derived from extensive health check-ups.

BACKGROUND: Phenome-wide association studies (PheWASs) have been conducted on Asian populations, including Koreans, but many were based on chip or exome genotyping data. Such studies have limitations regarding whole genome-wide association analysis, making it crucial to have genome-to-phenome association information with the largest possible whole genome and matched phenome data to conduct further population-genome studies and develop health care services based on population genomics. RESULTS: Here, we present 4,157 whole genome sequences (Korea4K) coupled with 107 health check-up parameters as the largest genomic resource of the Korean Genome Project. It encompasses most of the variants with allele frequency >0.001 in Koreans, indicating that it sufficiently covered most of the common and rare genetic variants with commonly measured phenotypes for Koreans. Korea4K provides 45,537,252 variants, and half of them were not present in Korea1K (1,094 samples). We also identified 1,356 new genotype-phenotype associations that were not found by the Korea1K dataset. Phenomics analyses further revealed 24 significant genetic correlations, 14 pleiotropic associations, and 127 causal relationships based on Mendelian randomization among 37 traits. In addition, the Korea4K imputation reference panel, the largest Korean variants reference to date, showed a superior imputation performance to Korea1K across all allele frequency categories. CONCLUSIONS: Collectively, Korea4K provides not only the largest Korean genome data but also corresponding health check-up parameters and novel genome-phenome associations. The large-scale pathological whole genome-wide omics data will become a powerful set for genome-phenome level association studies to discover causal markers for the prediction and diagnosis of health conditions in future studies.

Analysis of CSF3R mutations in atypical chronic myeloid leukemia and other myeloid malignancies.

We report a series of patients with CSF3R-mutant (CSF3Rmut) atypical chronic myeloid leukemia (aCML), chronic neutrophilic leukemia (CNL) or other hematologic malignancies. We included 25 patients: 5 aCML and 4 CNL CSF3Rmut patients; 1 aCML, 2 CNL, and 2 myelodysplastic/myeloproliferative neoplasm, not otherwise specified patients without CSF3R mutation; and 11 CSF3Rmut patients with other diseases [8 acute myeloid leukemia (AML), 1 chronic myelomonocytic leukemia (CMML), 1 myelodysplastic syndrome (MDS), and 1 acute lymphoblastic leukemia (ALL)]. Patients with aCML or CNL were tested by Sanger sequencing and pyrosequencing to identify CSF3R T618I. Twenty-two patients underwent gene panel analysis. CSF3R mutations, mostly T618I (8/9), were found at high frequencies in both aCML and CNL patients [5/6 aCML and 4/6 CNL]. Two aCML patients in early adulthood with CSF3R T618I and biallelic or homozygous CEBPA mutations without other mutations presented with increased blasts and exhibited remission for >6 years after transplantation. The other 7 CSF3Rmut aCML or CNL patients were elderly adults who all had ASXL1 mutations and frequently presented with SEBP1 and SRSF2 mutations. Five AML patients had CSF3R exon 14 or 15 point mutations, and 6 other patients (3 AML, 1 CMML, 1 MDS, and 1 ALL) had truncating mutations, demonstrating differences in leukocyte counts and mutation status. In conclusion, CSF3R mutations were found at a higher frequency in aCML patients than in previous studies, which might reflect ethnic differences. Additional studies are needed to confirm these findings and the relationship between CSF3R and CEBPA mutations.

Assessing the performance and utility of targeted next-generation sequencing for screening and genotyping of human papillomaviruses.

A next-generation sequencing (NGS)-based Ezplex HPV NGS kit (SML Genetree, Seoul, Korea) was used for human papillomavirus (HPV) screening. Of 885 cervical swab samples, HPV was detected in 162 samples. High-risk HPVs were detected in 82 samples, and other types of HPV were detected in 13 samples (HPV86, 71, 102, 91, and 114). At the read depth ≥ 500, NGS results exhibited 100 % agreement among repeated tests. HPV NGS results were compared with those of real-time PCR assays, Anyplex HPV28 (Seegene, Seoul, Korea) (n = 383) and Cobas HPV (Roche, Mannheim, Germany) (n = 64); concordances were 92.4 % and 95.0 %, respectively. Sanger sequencing of discordant results (n = 13) produced compatible results with those of HPV NGS. Pap smear abnormalities were detected in 31 patients (3.5 %), and 19 patients had high-risk HPV. Using HPV NGS for screening, rare HPV subtypes were detected, and quantitative values were obtained as read depth.

Strategies for Overcoming Immune Evasion in Bladder Cancer.

Tumors intricately shape a highly immunosuppressive microenvironment, hampering effective antitumor immune responses through diverse mechanisms. Consequently, achieving optimal efficacy in cancer immunotherapy necessitates the reorganization of the tumor microenvironment and restoration of immune responses. Bladder cancer, ranking as the second most prevalent malignant tumor of the urinary tract, presents a formidable challenge. Immunotherapeutic interventions including intravesical BCG and immune checkpoint inhibitors such as atezolizumab, avelumab, and pembrolizumab have been implemented. However, a substantial unmet need persists as a majority of bladder cancer patients across all stages do not respond adequately to immunotherapy. Bladder cancer establishes a microenvironment that can actively hinder an efficient anti-tumor immune response. A deeper understanding of immune evasion mechanisms in bladder cancer will aid in suppressing recurrence and identifying viable therapeutic targets. This review seeks to elucidate mechanisms of immune evasion specific to bladder cancer and explore novel pathways and molecular targets that might circumvent resistance to immunotherapy.

Anti-Obesity Effect and Signaling Mechanism of Potassium Poly-γ-Glutamate Produced by Bacillus subtilis Chungkookjang in High-Fat Diet-Induced Obese Mice.

The purpose of this work was to examine the effects of potassium poly-γ-glutamate (PGA-K) on mice fed a high-fat diet consisting of 60% of total calories for 12 weeks. PGA-K administration reduced the increase in body weight, epididymal fat, and liver weight caused by a high-fat diet compared to the obese group. The triglyceride, low-density lipoprotein cholesterol and high-density lipoprotein cholesterol levels, which are blood lipid indicators, were significantly increased in the obese group but were significantly decreased in the PGA-K-treated group. The administration of PGA-K resulted in a significant inhibition of pro-inflammatory cytokines, including tumor necrosis factor α and interleukin 6. Moreover, the levels of leptin and insulin, which are insulin resistance indicators, significantly increased in the obese group but were significantly decreased in the PGA-K-treated group. These results suggest that PGA-K exhibits a protective effect against obesity induced by a high-fat diet, underscoring its potential as a candidate for obesity treatment.

Comparison of structural variant callers for massive whole-genome sequence data.

BACKGROUND: Detecting structural variations (SVs) at the population level using next-generation sequencing (NGS) requires substantial computational resources and processing time. Here, we compared the performances of 11 SV callers: Delly, Manta, GridSS, Wham, Sniffles, Lumpy, SvABA, Canvas, CNVnator, MELT, and INSurVeyor. These SV callers have been recently published and have been widely employed for processing massive whole-genome sequencing datasets. We evaluated the accuracy, sequence depth, running time, and memory usage of the SV callers. RESULTS: Notably, several callers exhibited better calling performance for deletions than for duplications, inversions, and insertions. Among the SV callers, Manta identified deletion SVs with better performance and efficient computing resources, and both Manta and MELT demonstrated relatively good precision regarding calling insertions. We confirmed that the copy number variation callers, Canvas and CNVnator, exhibited better performance in identifying long duplications as they employ the read-depth approach. Finally, we also verified the genotypes inferred from each SV caller using a phased long-read assembly dataset, and Manta showed the highest concordance in terms of the deletions and insertions. CONCLUSIONS: Our findings provide a comprehensive understanding of the accuracy and computational efficiency of SV callers, thereby facilitating integrative analysis of SV profiles in diverse large-scale genomic datasets.

Comprehensive RNA-sequencing analysis of colorectal cancer in a Korean cohort.

Considering the recent increase in the number of colorectal cancer (CRC) cases in South Korea, we aimed to clarify the molecular characteristics of CRC unique to the Korean population. To gain insights into the complexities of CRC and promote the exchange of critical data, RNA-sequencing analysis was performed to reveal the molecular mechanisms that drive the development and progression of CRC; this analysis is critical for developing effective treatment strategies. We performed RNA-sequencing analysis of CRC and adjacent normal tissue samples from 214 Korean participants (comprising a total of 381 including 169 normal and 212 tumor samples) to investigate differential gene expression between the groups. We identified 19,575 genes expressed in CRC and normal tissues, with 3,830 differentially expressed genes (DEGs) between the groups. Functional annotation analysis revealed that the upregulated DEGs were significantly enriched in pathways related to the cell cycle, DNA replication, and IL-17, whereas the downregulated DEGs were enriched in metabolic pathways. We also analyzed the relationship between clinical information and subtypes using the Consensus Molecular Subtype (CMS) classification. Furthermore, we compared groups clustered within our dataset to CMS groups and performed additional analysis of the methylation data between DEGs and CMS groups to provide comprehensive biological insights from various perspectives. Our study provides valuable insights into the molecular mechanisms underlying CRC in Korean patients and serves as a platform for identifying potential target genes for this disease. The raw data and processed results have been deposited in a public repository for further analysis and exploration.

Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses.

The role of the serine/glycine metabolic pathway (SGP) has recently been demonstrated in tumors; however, the pathological relevance of the SGP in thyroid cancer remains unexplored. Here, we perform metabolomic profiling of 17 tumor-normal pairs; bulk transcriptomics of 263 normal thyroid, 348 papillary, and 21 undifferentiated thyroid cancer samples; and single-cell transcriptomes from 15 cases, showing the impact of mitochondrial one-carbon metabolism in thyroid tumors. High expression of serine hydroxymethyltransferase-2 (SHMT2) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is associated with low thyroid differentiation scores and poor clinical features. A subpopulation of tumor cells with high mitochondrial one-carbon pathway activity is observed in the single-cell dataset. SHMT2 inhibition significantly compromises mitochondrial respiration and decreases cell proliferation and tumor size in vitro and in vivo. Collectively, our results highlight the importance of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer and suggest that SHMT2 is a potent therapeutic target.

CYB5R3 functions as a tumor suppressor by inducing ER stress-mediated apoptosis in lung cancer cells via the PERK-ATF4 and IRE1α-JNK pathways.

Cytochrome b5 reductase 3 (CYB5R3) is involved in various cellular metabolic processes, including fatty acid synthesis and drug metabolism. However, the role of CYB5R3 in cancer development remains poorly understood. Here, we show that CYB5R3 expression is downregulated in human lung cancer cell lines and tissues. Adenoviral overexpression of CYB5R3 suppresses lung cancer cell growth in vitro and in vivo. However, CYB5R3 deficiency promotes tumorigenesis and metastasis in mouse models. Transcriptome analysis revealed that apoptosis- and endoplasmic reticulum (ER) stress-related genes are upregulated in CYB5R3-overexpressing lung cancer cells. Metabolomic analysis revealed that CYB5R3 overexpression increased the production of nicotinamide adenine dinucleotide (NAD+) and oxidized glutathione (GSSG). Ectopic CYB5R3 is mainly localized in the ER, where CYB5R3-dependent ER stress signaling is induced via activation of protein kinase RNA-like ER kinase (PERK) and inositol-requiring enzyme 1 alpha (IRE1α). Moreover, NAD+ activates poly (ADP-ribose) polymerase16 (PARP16), an ER-resident protein, to promote ADP-ribosylation of PERK and IRE1α and induce ER stress. In addition, CYB5R3 induces the generation of reactive oxygen species and caspase-9-dependent intrinsic cell death. Our findings highlight the importance of CYB5R3 as a tumor suppressor for the development of CYB5R3-based therapeutics for lung cancer.

The moderating role of circadian gene polymorphisms in the relationship between sleep disturbance and circulating lymphocyte subsets in colorectal cancer patients.

AIM: We investigated the impact of sleep disturbance on immune status in colorectal cancer (CRC) patients with consideration of the moderating role of circadian clock gene polymorphisms. METHODS: A prospective longitudinal study design was used to collect information regarding sleep disturbance. Blood samples for immunologic assays were obtained the day before the first (baseline) and last cycles of 5-fluorouracil, leucovorin, and oxaliplatin (FOLFOX) chemotherapy. Clinical sleep disturbance was compared between the two-time points using the Pittsburgh Sleep Quality Index (PSQI) global score. We analysed single-nucleotide polymorphisms in rs2278749, rs3749474, rs2291738, rs17031614, and rs2287161. The dependent variables included changes in the percentages of CD4+, CD8+, CD19+, and CD16/56+ lymphocytes between the two-time points. The results were analysed using moderated regression analysis; the p-values were adjusted using the false discovery rate. RESULTS: Among the 104 patients, no significant dyadic associations were observed between changes in lymphocyte percentages and the PSQI global score. However, the moderated regression analysis revealed five significant associations (rs2287161 with CD8+, rs2278749 and rs2291738 with CD19+, and rs17031614 with CD4+ and CD16/56+ lymphocytes). The inclusion of each interaction resulted in a significant increase (5.7-10.7%) in the variance explained by changes in lymphocyte percentage. CONCLUSION: Patients with specific circadian gene allele types may be more susceptible to immune dysregulation when experiencing sleep disturbances. Considering that sleep disturbance is a modifiable factor that can impact immune regulation, it is essential to prioritise the management of sleep disturbances in CRC patients receiving FOLFOX chemotherapy.

Clinical Manifestations and Genotype of Primary Ciliary Dyskinesia Diagnosed in Korea: Multicenter Study.

PURPOSE: Primary ciliary dyskinesia (PCD) is a genetically heterogeneous disorder that leads to secondary ciliary dysfunction. PCD is a rare disease, and data on it are limited in Korea. This study systematically evaluated the clinical symptoms, diagnostic characteristics, and treatment modalities of pediatric PCD in Korea. METHODS: This Korean nationwide, multicenter study, conducted between January 2000 and August 2022, reviewed the medical records of pediatric patients diagnosed with PCD. Prospective studies have been added to determine whether additional genetic testing is warranted in some patients. RESULTS: Overall, 41 patients were diagnosed with PCD in 15 medical institutions. The mean age at diagnosis was 11.8 ± 5.4 years (range: 0.5 months-18.9 years). Most patients (40/41) were born full term, 15 (36.6%) had neonatal respiratory symptoms, and 12 (29.3%) had a history of admission to the neonatal intensive care unit. The most common complaint (58.5%) was chronic nasal symptoms. Thirty-three patients were diagnosed with transmission electron microscopy (TEM) and 12 patients by genetic studies. TEM mostly identified outer dynein arm defects (alone or combined with inner dynein arm defects, n = 17). The genes with the highest mutation rates were DNAH5 (3 cases) and DNAAF1 (3 cases). Rare genotypes (RPGR, HYDIN, NME5) were found as well. Chest computed tomography revealed bronchiectasis in 33 out of 41 patients. Among them, 15 patients had a PrImary CiliAry DyskinesiA Rule score of over 5 points. CONCLUSIONS: To our knowledge, this is the first multicenter study to report the clinical characteristics, diagnostic methods, and genotypes of PCD in Korea. These results can be used as basic data for further PCD research.

Diagnosis of Crohn's disease and ulcerative colitis using the microbiome.

BACKGROUND: Inflammatory bowel disease (IBD) is a multifactorial chronic inflammatory disease resulting from dysregulation of the mucosal immune response and gut microbiota. Crohn's disease (CD) and ulcerative colitis (UC) are difficult to distinguish, and differential diagnosis is essential for establishing a long-term treatment plan for patients. Furthermore, the abundance of mucosal bacteria is associated with the severity of the disease. This study aimed to differentiate and diagnose these two diseases using the microbiome and identify specific biomarkers associated with disease activity. RESULTS: Differences in the abundance and composition of the microbiome between IBD patients and healthy controls (HC) were observed. Compared to HC, the diversity of the gut microbiome in patients with IBD decreased; the diversity of the gut microbiome in patients with CD was significantly lower. Sixty-eight microbiota members (28 for CD and 40 for UC) associated with these diseases were identified. Additionally, as the disease progressed through different stages, the diversity of the bacteria decreased. The abundances of Alistipes shahii and Pseudodesulfovibrio aespoeensis were negatively correlated with the severity of CD, whereas the abundance of Polynucleobacter wianus was positively correlated. The severity of UC was negatively correlated with the abundance of A. shahii, Porphyromonas asaccharolytica and Akkermansia muciniphilla, while it was positively correlated with the abundance of Pantoea candidatus pantoea carbekii. A regularized logistic regression model was used for the differential diagnosis of the two diseases. The area under the curve (AUC) was used to examine the performance of the model. The model discriminated UC and CD at an AUC of 0.873 (train set), 0.778 (test set), and 0.633 (validation set) and an area under the precision-recall curve (PRAUC) of 0.888 (train set), 0.806 (test set), and 0.474 (validation set). CONCLUSIONS: Based on fecal whole-metagenome shotgun (WMS) sequencing, CD and UC were diagnosed using a machine-learning predictive model. Microbiome biomarkers associated with disease activity (UC and CD) are also proposed.

Hypoxia stabilizes SETDB1 to maintain genome stability.

Von Hippel-Lindau (VHL) is a tumor suppressor that functions as the substrate recognition subunit of the CRL2VHL E3 complex. While substrates of VHL have been identified, its tumor suppressive role remains to be fully understood. For further determination of VHL substrates, we analyzed the physical interactome of VHL and identified the histone H3K9 methyltransferase SETBD1 as a novel target. SETDB1 undergoes oxygen-dependent hydroxylation by prolyl hydroxylase domain proteins and the CRL2VHL complex recognizes hydroxylated SETDB1 for ubiquitin-mediated degradation. Under hypoxic conditions, SETDB1 accumulates by escaping CRL2VHL activity. Loss of SETDB1 in hypoxia compared with that in normoxia escalates the production of transposable element-derived double-stranded RNAs, thereby hyperactivating the immune-inflammatory response. In addition, strong derepression of TEs in hypoxic cells lacking SETDB1 triggers DNA damage-induced death. Our collective results support a molecular mechanism of oxygen-dependent SETDB1 degradation by the CRL2VHL E3 complex and reveal a role of SETDB1 in genome stability under hypoxia.

Enrichment of Activated Fibroblasts as a Potential Biomarker for a Non-Durable Response to Anti-Tumor Necrosis Factor Therapy in Patients with Crohn's Disease.

We investigated whether the response to anti-tumor necrosis factor (anti-TNF) treatment varied according to inflammatory tissue characteristics in Crohn's disease (CD). Bulk RNA sequencing (RNA-seq) data were obtained from inflamed and non-inflamed tissues from 170 patients with CD. The samples were clustered based on gene expression profiles using principal coordinate analysis (PCA). Cellular heterogeneity was inferred using CiberSortx, with bulk RNA-seq data. The PCA results displayed two clusters of CD-inflamed samples: one close to (Inflamed_1) and the other far away (Inflamed_2) from the non-inflamed samples. Inflamed_1 was rich in anti-TNF durable responders (DRs), and Inflamed_2 was enriched in non-durable responders (NDRs). The CiberSortx results showed that the cell fraction of activated fibroblasts was six times higher in Inflamed_2 than in Inflamed_1. Validation with public gene expression datasets (GSE16879) revealed that the activated fibroblasts were enriched in NDRs over Next, we used DRs by 1.9 times pre-treatment and 7.5 times after treatment. Fibroblast activation protein (FAP) was overexpressed in the Inflamed_2 and was also overexpressed in the NDRs in both the RISK and GSE16879 datasets. The activation of fibroblasts may play a role in resistance to anti-TNF therapy. Characterizing fibroblasts in inflamed tissues at diagnosis may help to identify patients who are likely to respond to anti-TNF therapy.

Effect of TRAF6-knockout on gene expression and lncRNA expression in Epithelioma papulosum cyprini (EPC) cells.

TRAF6 is a key immune gene that plays a significant role in toll-like receptor signal transduction and activates downstream immune genes involved in antiviral immunity in fish. To explore the role of TRAF6 in Epithelioma papulosum cyprini (EPC) cells, we knocked out the TRAF6 gene using the Clustered Regularly Interspaced Short Palindromic Repeats-Cas9 (CRISPR-Cas9) technique and then analyzed the transcriptomes of the knockout cells. In this study, we identified that 232 transcripts were differentially expressed in naive cells. Using the pipeline, we identified 381 novel lncRNAs in EPC cells, 23 of which were differentially expressed. Gene Ontology enrichment analysis demonstrated that differentially expressed genes (DEG) are implicated in various immune processes, such as neutrophil chemotaxis and mitogen-activated protein kinase binding. In addition, the KEGG pathway analysis revealed enrichment in immune-related pathways (Interleukin-17 signaling pathway, cytokine-cytokine receptor interaction, and TNF signaling pathway). Furthermore, the target genes of the differentially expressed lncRNAs were implicated in the negative regulation of interleukin-6 and tumor necrosis factor production. These results indicate that lncRNAs and protein-coding genes participate in the regulation of immune and metabolic processes in fish.