Altered thalamic volumes and functional connectivity in the recovered patients with psychosis.

BACKGROUND: Investigating neural correlates in recovered patients with psychosis is important in terms of identifying biological markers associated with recovery status or predicting a possible future relapse. We sought to examine thalamic nuclei volumes and thalamus-centered functional connectivity (FC) in recovered patients with psychosis who discontinued their medication. METHODS: Thirty patients with psychosis who satisfied the criteria for full recovery and 50 healthy controls (HC) matched for age, sex, and education underwent magnetic resonance imaging and clinical evaluation. The recovered patients were divided into the maintained and relapsed subjects according to their clinical status on the follow-ups. Thalamic nuclei volumes and thalamus-centered FC were measured between the recovered patients and HC. Correlations between the thalamic nuclei or altered FC, and clinical symptoms and cognitive functioning were explored. RESULTS: Modest cognitive impairments and reduced thalamic nuclei volumes were evident in the recovered patients. Moreover, we found altered thalamo-cortical connectivity and its associations with negative symptoms and cognitive functioning in the recovered patients compared with HC. CONCLUSION: These findings suggest that there are still cognitive impairments, and aberrant neuronal changes in the recovered patients. The implication of differential FC patterns between the maintained and the relapsed patients remain to be further explored.

Deep Learning-based Brain Age Prediction in Patients With Schizophrenia Spectrum Disorders.

BACKGROUND AND HYPOTHESIS: The brain-predicted age difference (brain-PAD) may serve as a biomarker for neurodegeneration. We investigated the brain-PAD in patients with schizophrenia (SCZ), first-episode schizophrenia spectrum disorders (FE-SSDs), and treatment-resistant schizophrenia (TRS) using structural magnetic resonance imaging (sMRI). STUDY DESIGN: We employed a convolutional network-based regression (SFCNR), and compared its performance with models based on three machine learning (ML) algorithms. We pretrained the SFCNR with sMRI data of 7590 healthy controls (HCs) selected from the UK Biobank. The parameters of the pretrained model were transferred to the next training phase with a new set of HCs (n = 541). The brain-PAD was analyzed in independent HCs (n = 209) and patients (n = 233). Correlations between the brain-PAD and clinical measures were investigated. STUDY RESULTS: The SFCNR model outperformed three commonly used ML models. Advanced brain aging was observed in patients with SCZ, FE-SSDs, and TRS compared to HCs. A significant difference in brain-PAD was observed between FE-SSDs and TRS with ridge regression but not with the SFCNR model. Chlorpromazine equivalent dose and cognitive function were correlated with the brain-PAD in SCZ and FE-SSDs. CONCLUSIONS: Our findings indicate that there is advanced brain aging in patients with SCZ and higher brain-PAD in SCZ can be used as a surrogate marker for cognitive dysfunction. These findings warrant further investigations on the causes of advanced brain age in SCZ. In addition, possible psychosocial and pharmacological interventions targeting brain health should be considered in early-stage SCZ patients with advanced brain age.

Network biomarkers in recovered psychosis patients who discontinued antipsychotics.

There are no studies investigating topological properties of resting-state fMRI (rs-fMRI) in patients who have recovered from psychosis and discontinued medication (hereafter, recovered patients [RP]). This study aimed to explore topological organization of the functional brain connectome in the RP using graph theory approach. We recruited 30 RP and 50 age and sex-matched healthy controls (HC). The RP were further divided into the subjects who were relapsed after discontinuation of antipsychotics (RP-R) and who maintained recovered state without relapse (RP-M). Using graph-based network analysis of rs-fMRI signals, global and local metrics and hub information were obtained. The robustness of the network was tested with random failure and targeted attack. As an ancillary analysis, Network-Based Statistic (NBS) was performed. Association of significant findings with psychopathology and cognitive functioning was also explored. The RP showed intact network properties in terms of global and local metrics. However, higher global functional connectivity strength and hyperconnectivity in the interconnected component were observed in the RP compared to HC. In the subgroup analysis, the RP-R were found to have lower global efficiency, longer characteristic path length and lower robustness whereas no such abnormalities were identified in the RP-M. Associations of the degree centrality of some hubs with cognitive functioning were identified in the RP-M. Even though network properties of the RP were intact, subgroup analysis revealed more altered topological organizations in the RP-R. The findings in the RP-R and RP-M may serve as network biomarkers for predicting relapse or maintained recovery after the discontinuation of antipsychotics.

Isolation and characterization of single domain antibodies from banded houndshark (Triakis scyllium) targeting SARS-CoV-2 spike RBD protein.

The COVID-19 pandemic has significantly impacted human health for three years. To mitigate the spread of SARS-CoV-2, the development of neutralizing antibodies has been accelerated, including the exploration of alternative antibody formats such as single-domain antibodies. In this study, we identified variable new antigen receptors (VNARs) specific for the receptor binding domain (RBD) of SARS-CoV-2 by immunizing a banded houndshark (Triakis scyllium) with recombinant wild-type RBD. Notably, the CoV2NAR-1 clone showed high binding affinities in the nanomolar range to various RBDs and demonstrated neutralizing activity against SARS-CoV-2 pseudoviruses. These results highlight the potential of the banded houndshark as an animal model for the development of VNAR-based therapeutics or diagnostics against future pandemics.

Neuronal signatures of anger and fear in patients with psychosis.

The present study investigated the functional neuroanatomy in response to sentence stimuli related to anger-provoking situations and fear of negative evaluation in patients with psychosis. The tasks consisted of four active conditions, Self-Anger (SA), Self-Fear, Other-Anger (OA), and Other-Fear (OF), and two neutral conditions, Neutral-Anger (NA) and Neutral-Fear (NF). Several relevant clinical measures were obtained. Under all contrasts, significantly higher activation in the left inferior parietal gyrus or superior parietal gyrus and the left middle occipital gyrus or superior occipital gyrus was observed in patients compared to healthy controls (HCs). However, we observed significantly lower activation in the left angular gyrus (AG) and left middle temporal gyrus (MTG) under the OA vs. NA contrast, as well as in the left precuneus and left posterior cingulate gyrus (PCG) under the OF vs. NF contrast in patients. The mean beta values for the significant regions under the SA vs. NA and OF vs. NF contrasts were significantly associated with the total PI and PANSS scores, respectively. These findings indicate that patients with psychosis exhibit hypoactivation in the AG, MTG, precuneus, and PCG compared to HCs. The findings suggest that patients with psychosis are less efficient at recruiting neural responses in those regions for semantic processing and social evaluation.

Differential effects of parental socioeconomic status on cortical thickness in patients with schizophrenia spectrum disorders and healthy controls.

OBJECTIVES: Widespread changes in cortical thickness (CT) have been repeatedly reported in schizophrenia (SZ). The nature of the pathophysiologic process underlying such changes remains to be elucidated. The aims of the present study were to measure the CT; evaluate parent socioeconomic status (pSES), childhood trauma (ChT) and premorbid adjustment (PA) in patients with schizophrenia spectrum disorders (SSDs); and investigate group differences in CT (i.e., SSD vs. healthy controls (HCs)), pSES, PA, and/or ChT, as well as the interactions among these factors. METHODS: 164 patients with SSD and 245 age-, sex- and education-matched healthy controls have participated. The pSES, ChT and PA were evaluated using Korean version of Polyenvironmental Risk Score, Early Trauma Inventory Self Report-Short Form and Premorbid Adjustment Scale, respectively. Vertex-wise measure of CT was estimated using the FreeSurfer. To investigate the main effects and interactions, multilevel regression was employed. RESULTS: We found widespread cortical thinning in patients with SSDs compared to HCs. The cortical thinning was associated with ChT, symptom severity and chlorpromazine equivalent dose and duration of illness in patients. In multilevel regression, main effects of group and pSES and interaction between group and pSES were found whereas a significant interaction between ChT and CPZ equivalent was found in patients. CONCLUSION: Our findings indicate that compared to HCs, patients with SSDs have cortical structural abnormalities, and that group and pSES interaction determines CT. Further studies are needed to explore the effects of psychosocial factors on brain structural and functional abnormalities in SZ.

Predictors of full recovery in patients with early stage schizophrenia spectrum disorders.

To promote recovery in psychosis, targeting modifiable factors related to recovery is critical. Using more strict definition of full recovery, we examined predictors for recovery in patients with early stage schizophrenia spectrum disorders (SSD) followed up to 6.5 years. The target subjects were 375 patients with early stage SSD who had been over at least 1-year after registration and evaluated. The criteria for full recovery were having the score of the Positive and Negative Syndrome Scale (PANSS) 8-item ≤ 2 and adequate functional recovery for at least 1-year. We performed univariate Cox and stepwise Cox regression in both total and acute patients. In stepwise Cox regression, several independent predictors for recovery, i.e., negative symptoms of the PANSS, duration of untreated psychosis (DUP) and non-professional job were identified in patients with early stage SSD. In acute patients, other factors such as professional job and subjective well-being under neuroleptics were more important. The present study identified independent predictors for recovery modifiable by various psychosocial intervention and early intervention services. Moreover, it highlights the need of providing different treatment strategies depending on clinical status.

Stage-Specific Brain Aging in First-Episode Schizophrenia and Treatment-Resistant Schizophrenia.

BACKGROUND: Brain age is a popular brain-based biomarker that offers a powerful strategy for using neuroscience in clinical practice. We investigated the brain-predicted age difference (PAD) in patients with schizophrenia (SCZ), first-episode schizophrenia spectrum disorders (FE-SSDs), and treatment-resistant schizophrenia (TRS) using structural magnetic resonance imaging data. The association between brain-PAD and clinical parameters was also assessed. METHODS: We developed brain age prediction models for the association between 77 average structural brain measures and age in a training sample of controls (HCs) using ridge regression, support vector regression, and relevance vector regression. The trained models in the controls were applied to the test samples of the controls and 3 patient groups to obtain brain-based age estimates. The correlations were tested between the brain PAD and clinical measures in the patient groups. RESULTS: Model performance indicated that, regardless of the type of regression metric, the best model was support vector regression and the worst model was relevance vector regression for the training HCs. Accelerated brain aging was identified in patients with SCZ, FE-SSDs, and TRS compared with the HCs. A significant difference in brain PAD was observed between FE-SSDs and TRS using the ridge regression algorithm. Symptom severity, the Social and Occupational Functioning Assessment Scale, chlorpromazine equivalents, and cognitive function were correlated with the brain PAD in the patient groups. CONCLUSIONS: These findings suggest additional progressive neuronal changes in the brain after SCZ onset. Therefore, pharmacological or psychosocial interventions targeting brain health should be developed and provided during the early course of SCZ.

Brain ageing in schizophrenia: evidence from 26 international cohorts via the ENIGMA Schizophrenia consortium.

Schizophrenia (SZ) is associated with an increased risk of life-long cognitive impairments, age-related chronic disease, and premature mortality. We investigated evidence for advanced brain ageing in adult SZ patients, and whether this was associated with clinical characteristics in a prospective meta-analytic study conducted by the ENIGMA Schizophrenia Working Group. The study included data from 26 cohorts worldwide, with a total of 2803 SZ patients (mean age 34.2 years; range 18-72 years; 67% male) and 2598 healthy controls (mean age 33.8 years, range 18-73 years, 55% male). Brain-predicted age was individually estimated using a model trained on independent data based on 68 measures of cortical thickness and surface area, 7 subcortical volumes, lateral ventricular volumes and total intracranial volume, all derived from T1-weighted brain magnetic resonance imaging (MRI) scans. Deviations from a healthy brain ageing trajectory were assessed by the difference between brain-predicted age and chronological age (brain-predicted age difference [brain-PAD]). On average, SZ patients showed a higher brain-PAD of +3.55 years (95% CI: 2.91, 4.19; I2 = 57.53%) compared to controls, after adjusting for age, sex and site (Cohen's d = 0.48). Among SZ patients, brain-PAD was not associated with specific clinical characteristics (age of onset, duration of illness, symptom severity, or antipsychotic use and dose). This large-scale collaborative study suggests advanced structural brain ageing in SZ. Longitudinal studies of SZ and a range of mental and somatic health outcomes will help to further evaluate the clinical implications of increased brain-PAD and its ability to be influenced by interventions.

Functional Expression of the Recombinant Spike Receptor Binding Domain of SARS-CoV-2 Omicron in the Periplasm of Escherichia coli.

A new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant known as Omicron has caused a rapid increase in recent global patients with coronavirus infectious disease 2019 (COVID-19). To overcome the COVID-19 Omicron variant, production of a recombinant spike receptor binding domain (RBD) is vital for developing a subunit vaccine or a neutralizing antibody. Although bacterial expression has many advantages in the production of recombinant proteins, the spike RBD expressed in a bacterial system experiences a folding problem related to disulfide bond formation. In this study, the soluble Omicron RBD was obtained by a disulfide isomerase-assisted periplasmic expression system in Escherichia coli. The Omicron RBD purified from E. coli was very well recognized by anti-SARS-CoV-2 antibodies, sotrovimab (S309), and CR3022, which were previously reported to bind to various SARS-CoV-2 variants. In addition, the kinetic parameters of the purified Omicron RBD upon binding to the human angiotensin-converting enzyme 2 (ACE2) were similar to those of the Omicron RBD produced in the mammalian expression system. These results suggest that an E. coli expression system would be suitable to produce functional and correctly folded spike RBDs of the next emerging SARS-CoV-2 variants quickly and inexpensively.

Altered functional connectivity in psychotic disorder not otherwise specified.

BACKGROUND: Few studies have investigated functional connectivity (FC) in patients with psychotic disorder not otherwise specified (PNOS). We sought to identify distinct FC differentiating PNOS from schizophrenia (SZ). METHODS: In total, 49 patients with PNOS, 42 with SZ, and 55 healthy controls (HC) matched for age, sex, and education underwent functional magnetic resonance imaging (fMRI) brain scans and clinical evaluation. Using six functional networks consisting of 40 regions of interest (ROIs), we conducted ROI to ROI and intra- and inter-network FC analyses using resting-state fMRI (rs-fMRI) data. Correlations of altered FC with symptomatology were explored. RESULTS: We found common brain connectomics in PNOS and SZ including thalamo-cortical (especially superior temporal gyrus) hyperconnectivity, thalamo-cerebellar hypoconnectivity, and reduced within-thalamic connectivity compared to HC. Additionally, features differentiating the two patient groups included hyperconnectivity between the thalamic subregion and anterior cingulate cortex in PNOS compared to SZ and hyperconnectivity of the thalamic subregions with the posterior cingulate cortex and precentral gyrus in SZ compared to PNOS. CONCLUSIONS: These findings suggest that PNOS and SZ exhibit both common and differentiating changes in neuronal connectivity. Furthermore, they may support the hypothesis that PNOS should be treated as a separate clinical syndrome with distinct neural connectomics.

Three-year outcomes and predictors for full recovery in patients with early-stage psychosis.

In the present study, various outcomes over 3-year period in patients with early stage psychosis including remission, recovery, relapse and medication adherence were investigated. Predictor for full recovery at year 3 was also examined. Three-year follow-up data in 534 patients with schizophrenia spectrum disorders (SSD) and psychotic disorder not otherwise specified (PNOS) were examined for overall outcome trajectories. The data of completers at year 3 (n = 157) were used to identify predictors for recovery using logistic regression. The rates of symptomatic remission and full recovery at 6-, 12-, 24-, and 36-month follow-up were 76.10, 69.20, 79.50, and 79.10%, and 22.80, 26.40, 28.60, and 39.60%, respectively. The rates of drop-out and relapse at 6-, 12-, 24-, and 36-month follow-up were 25.4, 29.5, 38.6, and 51.1%, and 3.7, 8.9, 19.0, and 38.9%, respectively. The rates of good adherence and prescription of Long-Acting Injectable Antipsychotics (LAIA) at 6-, 12-, 24- and 36-month follow-up were 87.8, 88.0, 91.9, and 93.9%, and 18.3, 21.7, 22.0, and 25.5%, respectively. Significant predictors for full recovery were duration of untreated psychosis (DUP), family intimacy and physical activity. We observed similar or better results on remission, recovery, and relapse rates compared to other previous studies. Effective psychosocial intervention should be provided to shorten the gap between remission and recovery rates and to address DUP, family issues, and exercise to enhance recovery.

Altered thalamic subregion functional networks in patients with treatment-resistant schizophrenia.

BACKGROUND: The thalamus plays a key role in filtering information and has extensive interconnectivity with other brain regions. A large body of evidence points to impaired functional connectivity (FC) of the thalamocortical pathway in schizophrenia. However, the functional network of the thalamic subregions has not been investigated in patients with treatment-resistant schizophrenia (TRS). AIM: To identify the neural mechanisms underlying TRS, we investigated FC of thalamic sub-regions with cortical networks and voxels, and the associations of this FC with clinical symptoms. We hypothesized that the FC of thalamic sub-regions with cortical networks and voxels would differ between TRS patients and HCs. METHODS: In total, 50 patients with TRS and 61 healthy controls (HCs) matched for age, sex, and education underwent resting-state functional magnetic resonance imaging (rs-fMRI) and clinical evaluation. Based on the rs-fMRI data, we conducted a FC analysis between thalamic subregions and cortical functional networks and voxels, and within thalamic subregions and cortical functional networks, in the patients with TRS. A functional parcellation atlas was used to segment the thalamus into nine subregions. Correlations between altered FC and TRS symptoms were explored. RESULTS: We found differences in FC within thalamic subregions and cortical functional networks between patients with TRS and HCs. In addition, increased FC was observed between thalamic subregions and the sensorimotor cortex, frontal medial cortex, and lingual gyrus. These abnormalities were associated with the pathophysiology of TRS. CONCLUSION: Our findings suggest that disrupted FC within thalamic subregions and cortical functional networks, and within the thalamocortical pathway, has potential as a marker for TRS. Our findings also improve our understanding of the relationship between the thalamocortical pathway and TRS symptoms.

Altered Neuroanatomical Signatures of Patients With Treatment-Resistant Schizophrenia Compared to Patients With Early-Stage Schizophrenia and Healthy Controls.

Background: The relationship between brain structural changes and cognitive dysfunction in schizophrenia is strong. However, few studies have investigated both neuroanatomical abnormalities and cognitive dysfunction in treatment-resistant schizophrenia (TRS). We examined neuroanatomical markers and cognitive function between patients with TRS or early-stage schizophrenia (ES-S) and healthy controls (HCs). Relationships between neuroanatomical markers and cognitive function in the patient groups were also investigated. Methods: A total of 46 and 45 patients with TRS and ES-S and 61 HCs underwent structural magnetic resonance imaging (MRI) brain scanning and comprehensive cognitive tests. MRI scans were analyzed using the FreeSurfer to investigate differences in cortical surface area (CSA), cortical thickness (CT), cortical volume (CV), and subcortical volume (SCV) among the groups. Four cognitive domains (attention, verbal memory, executive function, and language) were assessed. Comparisons of neuroanatomical and cognitive function results among the three groups were performed. Results: A widespread reduction in CT was observed in patients with TRS compared to HCs, but differences in cortical thinning between TRS and ES-S patients were mainly limited to the inferior frontal gyrus and insula. Several subcortical structures (accumbens, amygdala, hippocampus, putamen, thalamus and ventricles) were significantly altered in TRS patients compared to both ES-S patients and HCs. Performance in the verbal memory domain was significantly worse in TRS patients compared to ES-S patients. A positive relationship between the thickness of the left middle temporal gyrus and the composite score for language was identified in patients with ES-S. Conclusions: Our findings suggest significant cognitive impairment and reductions in CT and SCV in individuals with TRS compared to those with ES-S and HCs. These abnormalities could act as biomarkers for earlier identification of TRS.

Comparison of clinical features and 1-year outcomes between patients with psychotic disorder not otherwise specified and those with schizophrenia.

AIM: Research on psychotic disorder not otherwise specified (PNOS) that clearly mentions its subgroups is very rare. This study was conducted to identify the demographic and clinical features, cognitive function, and 1-year outcomes of patients with early stage PNOS compared with those with early stage schizophrenia (SZ). METHODS: The study subjects were 54 and 321 patients with PNOS and SZ, respectively, who were registered at least more than 1 year ago. Due to drop out, only 37 and 210 patients with PNOS and SZ were evaluated at the 1-year follow-up. We compared clinical variables (duration of untreated psychosis, symptom severity, self-rating scales, and so on), cognitive function, and short-term outcomes (treatment response, remission, compliance, drop out, relapse) between the two groups. RESULTS: The patients with PNOS were associated with higher diagnostic stability (53.7%) compared with those in previous studies. They had lower symptom severity, better treatment response at 2 months and higher remission rates at 12 months, but poorer compliance at 6 months compared with patients with SZ. Level of cognitive impairment in PNOS was intermediate between those of SZ patients and healthy controls. CONCLUSIONS: These findings indicate that PNOS has unique clinical features, suggesting that it should be treated as a distinct clinical syndrome. At the same time, however, prevention of its possible progression to other psychotic disorders in some patients with PNOS is also important.

Diagnosis of schizophrenia with functional connectome data: a graph-based convolutional neural network approach.

Previous deep learning methods have not captured graph or network representations of brain structural or functional connectome data. To address this, we developed the BrainNet-Global Covariance Pooling-Attention Convolutional Neural Network (BrainNet-GA CNN) by incorporating BrainNetCNN and global covariance pooling into the self-attention mechanism. Resting-state functional magnetic resonance imaging data were obtained from 171 patients with schizophrenia spectrum disorders (SSDs) and 161 healthy controls (HCs). We conducted an ablation analysis of the proposed BrainNet-GA CNN and quantitative performance comparisons with competing methods using the nested tenfold cross validation strategy. The performance of our model was compared with competing methods. Discriminative connections were visualized using the gradient-based explanation method and compared with the results obtained using functional connectivity analysis. The BrainNet-GA CNN showed an accuracy of 83.13%, outperforming other competing methods. Among the top 10 discriminative connections, some were associated with the default mode network and auditory network. Interestingly, these regions were also significant in the functional connectivity analysis. Our findings suggest that the proposed BrainNet-GA CNN can classify patients with SSDs and HCs with higher accuracy than other models. Visualization of salient regions provides important clinical information. These results highlight the potential use of the BrainNet-GA CNN in the diagnosis of schizophrenia.

Longitudinal symptom network structure in first-episode psychosis: a possible marker for remission.

BACKGROUND: Network approach has been applied to a wide variety of psychiatric disorders. The aim of the present study was to identify network structures of remitters and non-remitters in patients with first-episode psychosis (FEP) at baseline and the 6-month follow-up. METHODS: Participants (n = 252) from the Korean Early Psychosis Study (KEPS) were enrolled. They were classified as remitters or non-remitters using Andreasen's criteria. We estimated network structure with 10 symptoms (three symptoms from the Positive and Negative Syndrome Scale, one depressive symptom, and six symptoms related to schema and rumination) as nodes using a Gaussian graphical model. Global and local network metrics were compared within and between the networks over time. RESULTS: Global network metrics did not differ between the remitters and non-remitters at baseline or 6 months. However, the network structure and nodal strengths associated with positive-self and positive-others scores changed significantly in the remitters over time. Unique central symptoms for remitters and non-remitters were cognitive brooding and negative-self, respectively. The correlation stability coefficients for nodal strength were within the acceptable range. CONCLUSION: Our findings indicate that network structure and some nodal strengths were more flexible in remitters. Negative-self could be an important target for therapeutic intervention.

Network analysis of trauma in patients with early-stage psychosis.

Childhood trauma (ChT) is a risk factor for psychosis. Negative lifestyle factors such as rumination, negative schemas, and poor diet and exercise are common in psychosis. The present study aimed to perform a network analysis of interactions between ChT and negative lifestyle in patients and controls. We used data of patients with early-stage psychosis (n = 500) and healthy controls (n = 202). Networks were constructed using 12 nodes from five scales: the Brief Core Schema Scale (BCSS), Brooding Scale (BS), Dietary Habits Questionnaire, Physical Activity Rating, and Early Trauma Inventory Self Report-Short Form (ETI). Graph metrics were calculated. The nodes with the highest predictability and expected influence in both patients and controls were cognitive and emotional components of the BS and emotional abuse of the ETI. The emotional abuse was a mediator in the shortest pathway connecting the ETI and negative lifestyle for both groups. The negative others and negative self of the BCSS mediated emotional abuse to other BCSS or BS for patients and controls, respectively. Our findings suggest that rumination and emotional abuse were central symptoms in both groups and that negative others and negative self played important mediating roles for patients and controls, respectively.Trial Registration: ClinicalTrials.gov identifier: CUH201411002.

Different prognostic implications of hepatic metastasis according to front-line treatment in non-small cell lung cancer: a real-world retrospective study.

BACKGROUND: Although liver metastasis occurs in approximately 15% of metastatic non-small cell lung cancer (NSCLC) patients with poor prognosis, its prognostic effect in patients who receive immunotherapy is unclear. This study aimed to verify the effects of liver metastasis on the prognosis of metastatic NSCLC patients according to their first-line treatment. METHODS: Patients who were initially diagnosed with stage 4 NSCLC from January 2015 to December 2019 were analyzed in this retrospective real-world data-based study. The patients were divided into three groups according to the type of first-line chemotherapy they received: cytotoxic, targeted, and immunotherapy. Prognosis was then compared depending on the presence of liver metastasis in each treatment group. RESULTS: Among the 1,470 patients, 723 (49.2%) received cytotoxic chemotherapy, 678 (46.1%) received targeted therapy, and 69 (4.7%) received immunotherapy as their first-line chemotherapy. A total of 234 (15.9%) patients had liver metastasis at the initial diagnosis. The mean patient age was 63.7 years, and 59.1% were male. There was no difference in overall survival (OS) in the immunotherapy group in patients with or without liver metastasis (11.7 vs. 13.0 months, P=0.968); however, patients with liver metastasis had worse outcomes in the cytotoxic and targeted therapy groups compared to patients without liver metastasis. Furthermore, in patients with liver metastasis, the immunotherapy group had a longer OS than the cytotoxic chemotherapy group (11.7 vs. 4.4 months, P<0.001). Liver metastasis was associated with poor outcomes (hazard ratio of 1.438), as were age, male sex, bone, adrenal gland, or soft tissue metastasis, and three or more metastatic sites; however, lymph node, brain, collateral lung, and pleura metastasis did not affect prognosis. CONCLUSIONS: Although liver metastasis was associated with poor outcomes, it did not affect prognosis in patients who received immunotherapy.