A phase 1 study of biweekly nab-paclitaxel/oxaliplatin/S-1/LV for advanced upper gastrointestinal cancers: TCOG T1216 study.

BACKGROUND: Oxaliplatin- and fluoropyrimidine-based triplet regimens have demonstrated feasibility and efficacy in the treatment of upper gastrointestinal (UGI) cancers. Herein, we evaluate the feasibility and preliminary efficacy of biweekly nab-paclitaxel plus oxaliplatin and S-1/leucovorin (SOLAR) in chemonaïve UGI cancers. METHODS: A 3 + 3 phase 1 study was conducted to determine the maximal tolerated dose (MTD) of oxaliplatin in SOLAR (nab-paclitaxel [150 mg/m2 in D1], oxaliplatin [60, 75, or 85 mg/m2 in D1], and oral S-1/leucovorin [35 mg/m2 and 30 mg bid from D1 to D7]). The secondary endpoints were overall response rate (ORR), progression-free survival (PFS), overall survival (OS), and safety. RESULTS: Thirteen and 6 accruals were in the dose-escalation and MTD expansion cohorts, respectively. One of 6 patients at level III experienced dose-limiting toxicity (grade 3 diarrhea), which revealed that the MTD of oxaliplatin was 85 mg/m2. After a mean of 15.9 cycles of treatment, the most common treatment-related grade 3/4 toxicities were neutropenia (57.9%) and diarrhea (21.1%). The ORR was 63.2%. The median PFS and OS were 12.5 and 24.7 months, respectively. CONCLUSION: The current study revealed the MTD of oxaliplatin and demonstrated the preliminary efficacy of SOLAR in UGI cancers, which deserves further investigation. CLINICALTRIALS.GOV IDENTIFIER: NCT03162510.

Survival of advanced/recurrent gastrointestinal stromal tumors treated with tyrosine kinase inhibitors in Taiwan: a nationwide registry study.

BACKGROUND: Most gastrointestinal stromal tumors (GISTs) harbor c-KIT or PDGFRA mutations. Administration of tyrosine kinase inhibitors (TKIs) has significantly improved the survival of patients with GISTs. We aimed to evaluate the clinical outcome of advanced or recurrent GIST patients in Taiwan. METHODS: Patients diagnosed between 2010 and 2020 were enrolled. The collected data included baseline characteristics, treatment pattern, treatment outcome, genetic aberrations and survival status. Progression-free survival (PFS) and overall survival (OS) were analyzed and plotted with the Kaplan-Meier method. Cox regression analysis was used to analyze the prognostic factors of survival. RESULTS: A total of 224 patients with advanced or recurrent GISTs treated with TKIs were enrolled. All patients received imatinib treatment. Ninety-three and 42 patients received sunitinib and regorafenib treatment, respectively. The 48-month PFS and OS rates for patients treated with imatinib were 50.5% and 79.5%, respectively. c-KIT exon 9 and PDGFRA mutations were prognostic factors for a poor PFS and PDGFRA mutation was a prognostic factor for a poor OS in patients treated with imatinib in multivariate Cox regression analysis. The median PFS of patients who received sunitinib treatment was 12.76 months (95% confidence interval (CI), 11.01-14.52). Patients with c-KIT exon 9 mutations had a longer PFS than those with other genetic aberrations. The median PFS of patients treated with regorafenib was 7.14 months (95% CI, 3.39-10.89). CONCLUSIONS: We present real-world clinical outcomes for advanced GIST patients treated with TKIs and identify mutational status as an independent prognostic factor for patient survival.

Bayesian hierarchical models for adaptive basket trial designs.

Basket trials evaluate a single drug targeting a single genetic variant in multiple cancer cohorts. Empirical findings suggest that treatment efficacy across baskets may be heterogeneous. Most modern basket trial designs use Bayesian methods. These methods require the prior specification of at least one parameter that permits information sharing across baskets. In this study, we provide recommendations for selecting a prior for scale parameters for adaptive basket trials by using Bayesian hierarchical modeling. Heterogeneity among baskets attracts much attention in basket trial research, and substantial heterogeneity challenges the basic assumption of exchangeability of Bayesian hierarchical approach. Thus, we also allowed each stratum-specific parameter to be exchangeable or nonexchangeable with similar strata by using data observed in an interim analysis. Through a simulation study, we evaluated the overall performance of our design based on statistical power and type I error rates. Our research contributes to the understanding of the properties of Bayesian basket trial designs.

Use of likelihood estimates for variances for the design and evaluation of multiregional clinical trials with heterogeneous variances.

Globalized drug development studies, such as multiregional clinical trials (MRCTs), have attracted much attention due to their ability to expedite drug development and shorten the time lag of drug release. While observing the overall effect of a new drug, the region-specific effects to support drug registration in constituent regions can also be evaluated. Several challenges arise in conducting MRCTs, such as the heterogeneity in the variability of the primary endpoint across regions. However, most of the existing statistical methods assume a common variability, which may not be valid in practice due to differences across regions (eg, diversities in ethnicity or disparities in medical culture/practice). We present a statistical method for the design and evaluation of MRCTs to consider the heterogeneous variability across regions. We assessed the overall sample size requirement and addressed the region-specific sample size determination to establish the consistency of treatment effects between the specific region and the entire group. We demonstrate the proposed approach with numerical examples.

Use of a two-sided tolerance interval in the design and evaluation of biosimilarity in clinical studies.

In assessing biosimilarity between two products, the question to ask is always "How similar is similar?" Traditionally, the equivalence of the means between products is the primary consideration in a clinical trial. This study suggests an alternative assessment for testing a certain percentage of the population of differences lying within a prespecified interval. In doing so, the accuracy and precision are assessed simultaneously by judging whether a two-sided tolerance interval falls within a prespecified acceptance range. We further derive an asymptotic distribution of the tolerance limits to determine the sample size for achieving a targeted level of power. Our numerical study shows that the proposed two-sided tolerance interval test controls the type I error rate and provides sufficient power. A real example is presented to illustrate our proposed approach.

An Open-Label, Single-Arm, Two-Stage, Multicenter, Phase II Study to Evaluate the Efficacy of TLC388 and Genomic Analysis for Poorly Differentiated Neuroendocrine Carcinomas.

BACKGROUND: The discovery of effective therapeutic options for treating metastatic poorly differentiated neuroendocrine carcinoma (NEC) after prior platinum-based chemotherapy remains elusive. This study analyzed the efficacy of TLC388 (Lipotecan) Hydrochloride, a novel camptothecin analog, for pretreated patients with metastatic NEC. METHODS: This single-arm, two-stage, phase II clinical trial was conducted at four community and academic centers in Taiwan. Patients aged 20 years or older with confirmed metastatic NEC and who had received prior systemic therapy with etoposide plus cisplatin were enrolled between July 2015 and May 2018. Patients received 40 mg/m2 of TLC388 intravenously on days 1, 8, and 15 of a 28-day cycle until disease progression or unacceptable toxic effects. Gene mutations were analyzed by next-generation sequencing. RESULTS: Twenty-three patients with a median age of 61 (range, 44-73) years, 18 of whom were men (78%), were enrolled. Patients received a median of 2 (range, 0-6) treatment cycles. Among 20 evaluable patients, 3 patients exhibited stable disease and no patient experienced a complete or partial remission, resulting in a disease control rate of 15%. Median progression-free survival was 1.8 (95% confidence interval [CI], 0.4-15) months, and the median overall survival was 4.3 (95% CI, 1.7-15) months. The most common treatment-related hematologic adverse events at grade 3 or higher were leukopenia (22.7%), anemia (31.8%), and thrombocytopenia (18.2%). The most frequent mutated genes in 35 patients with NEC were ARSA, DPYD, HEXB, BRCA1, HPD, MYBPC3, BBS2, IL7R, HSD17B4, and PRODH. CONCLUSION: TLC388 demonstrates limited antitumor activity in metastatic NEC. ClinicalTrials.gov identifier: NCT02457273. IMPLICATIONS FOR PRACTICE: Poorly differentiated neuroendocrine carcinomas (NECs) are rare and aggressive. Currently, effective therapeutic options for treating metastatic poorly differentiated NECs beyond platinum-based chemotherapy remain elusive. In this single-arm, multicenter, phase II study, 23 patients with NEC were enrolled and received TLC388 (Lipotecan) Hydrochloride, which is a novel camptothecin analog. The results demonstrated the disease control rate of 15%, the median progression-free survival of 1.8 (95% confidence interval [CI], 0.4-15) months, and the median overall survival of 4.3 (95% CI, 1.7-15) months. Most importantly, several novel genetic mutations and pathways were identified. These results offer the opportunity to develop future treatment strategies in this rare cancer.

Use of tolerance intervals for assessing biosimilarity.

A biosimilar is a biological product that is highly similar to an existing approved reference drug and has no clinically meaningful difference from it. Biosimilars are composed of or derived from living cells or organisms. Therefore, they are often sensitive to slight variations in the manufacturing process. Consequently, in demonstrating biosimilarity, it might be inappropriate to focus solely on the mean difference, or ratio of means, while ignoring the variabilities associated with the test and reference products. It is important to account for the entire population of clinical outcomes. Thus, we propose using the concept of tolerance intervals and related hypothesis testing for assessing biosimilarity. Our approach has the advantage of considering entire populations associated with both groups. A real example is used to illustrate our proposed method, and our approach is more stringent than those that employ confidence intervals. This is specifically the case when the mean difference of two drugs is not sufficiently large, but the biosimilar has a higher variability than that in the reference drug.

Use of a tolerance interval approach as a statistical quality control tool for traditional Chinese medicine.

Raw materials for traditional Chinese medicine (TCM) are often from different resources and its final product may also be made by different sites. Therefore, variabilities from different resources such as site-to-site or within site component-to-component may be expected. Consequently, test for consistency in raw materials, in-process materials, and/or final product has become an important issue in the quality control (QC) process in TCM development. In this paper, a statistical QC process for raw materials and/or the final product of TCM is proposed based on a two sided [Formula: see text]-content, [Formula: see text]-confidence tolerance interval. More specifically, we construct the tolerance interval for a random-effects model to assess the QC of TCM products from different regions and possibly different product batches. The products can be claimed to be consistency when the constructed tolerance interval is within the permitted range. Given the region and batch effects, sample sizes can also be calculated to ensure the desired measure of goodness. An example is presented to illustrate the proposed approach.

Cost-effectiveness of preventing hepatitis B virus reactivation in patients with lymphoma and resolved HBV infection.

BACKGROUND/PURPOSE: Hepatitis B virus (HBV) reactivation may occur in >10% of patients with lymphoma and resolved HBV infection who undergo rituximab-containing chemotherapy. Preventive strategies may have marked impact on resource allocation in HBV endemic areas. This study aims to compare the cost-effectiveness between prophylactic antiviral therapy and HBV DNA monitoring for the prevention of HBV-related complications. METHODS: Data sources are studies of HBV-related events and survival for patients with lymphoma and resolved HBV infection published since 2006. Decision tree analysis was used to compare the incremental cost-effectiveness ratio (ICER) of preventing HBV-related death or liver decompensation for patients who undergo first-line rituximab-containing chemotherapy. Sensitivity analysis was performed to examine the impact of the preventive efficacy, the duration of prophylactic antiviral therapy, and the cost of different interventions. The direct medical cost was derived from the database of the NHI Administration, Taiwan. The time frame of our analysis was set to 3 years after the completion of chemotherapy. RESULTS: The median ICER of prophylactic antiviral therapy, according to current practice guidelines, ranged between USD 150,000 and 250,000 if we apply the guidelines generally. When a long-course (12 months after completion of chemotherapy according to clinical guidelines) prophylactic therapy was assumed, Option A was cheaper and more effective only in the anti-HBs-negative subgroup (median ICER US$149,932 vs. US$161,526, p = 0.013). CONCLUSION: Identification of anti-HBs-negative subgroups is critical to improve the cost-effectiveness of prophylactic antiviral therapy in lymphoma patients with resolved HBV infection.

Genetic Modifiers of Progression-Free Survival in Never-Smoking Lung Adenocarcinoma Patients Treated with First-Line Tyrosine Kinase Inhibitors.

RATIONALE: Patients with non-small cell lung cancer (NSCLC) with mutated epidermal growth factor receptor (EGFR) are relatively sensitive to EGFR-tyrosine kinase inhibitor (TKI) treatment and have longer progression-free survival (PFS) when treated with EGFR-TKI compared with platinum-based chemotherapy. However, many patients with advanced NSCLC who have mutated EGFR do not respond to first-line EGFR-TKI treatment and still have shorter PFS. OBJECTIVES: The aim of this study was to identify genetic variants associated with PFS among patients with lung adenocarcinoma who were treated with first-line EGFR-TKIs. METHODS: A genome-wide association study on PFS was performed in never-smoking women diagnosed with lung adenocarcinoma and who were treated with first-line EGFR-TKIs (n = 128). Significant single-nucleotide polymorphisms (SNPs) were selected for follow-up association analysis (n = 198) and for replication assay in another independent cohort (n = 153). MEASUREMENTS AND MAIN RESULTS: We identified SNPs at 4q12 associated with PFS at genome-wide significance (P < 10-8) and with an estimated hazard ratio of more than 4. This association was also replicated in a larger but similar cohort and in an independent NSCLC cohort. Follow-up functional analyses showed that these SNPs were associated with the expression of EGFR, which encodes the TKI target, and with a nearby gene neuromedin-U, which encodes a G protein-coupled receptor ligand known to be involved in the progression of NSCLC. Considering these as possible prognostic biomarkers for the treatment of patients with late-stage lung cancer, we found that these SNPs were not associated with EGFR mutation status or with polymorphism of the Bcl2-interacting mediator of cell death gene. CONCLUSIONS: Genetic variants in 4q12 merit further investigation to assess their potential as pharmacogenomic predictors for and to understand the biology underlying its influence on PFS in patients treated with TKI therapy.

An approach for sample size determination of average bioequivalence based on interval estimation.

In 1992, the US Food and Drug Administration declared that two drugs demonstrate average bioequivalence (ABE) if the log-transformed mean difference of pharmacokinetic responses lies in (-0.223, 0.223). The most widely used approach for assessing ABE is the two one-sided tests procedure. More specifically, ABE is concluded when a 100(1 - 2α) % confidence interval for mean difference falls within (-0.223, 0.223). As known, bioequivalent studies are usually conducted by crossover design. However, in the case that the half-life of a drug is long, a parallel design for the bioequivalent study may be preferred. In this study, a two-sided interval estimation - such as Satterthwaite's, Cochran-Cox's, or Howe's approximations - is used for assessing parallel ABE. We show that the asymptotic joint distribution of the lower and upper confidence limits is bivariate normal, and thus the sample size can be calculated based on the asymptotic power so that the confidence interval falls within (-0.223, 0.223). Simulation studies also show that the proposed method achieves sufficient empirical power. A real example is provided to illustrate the proposed method. Copyright © 2017 John Wiley & Sons, Ltd.

Assessing the consistency of the treatment effect under the discrete random effects model in multiregional clinical trials.

In recent years, developing pharmaceutical products via multiregional clinical trials (MRCTs) has become standard. Traditionally, an MRCT would assume that a treatment effect is uniform across regions. However, heterogeneity among regions may have impact upon the evaluation of a medicine's effect. In this study, we consider a random effects model using discrete distribution (DREM) to account for heterogeneous treatment effects across regions for the design and evaluation of MRCTs. We derive an power function for a treatment that is beneficial under DREM and illustrate determination of the overall sample size in an MRCT. We use the concept of consistency based on Method 2 of the Japanese Ministry of Health, Labour, and Welfare's guidance to evaluate the probability for treatment benefit and consistency under DREM. We further derive an optimal sample size allocation over regions to maximize the power for consistency. Moreover, we provide three algorithms for deriving sample size at the desired level of power for benefit and consistency. In practice, regional treatment effects are unknown. Thus, we provide some guidelines on the design of MRCTs with consistency when the regional treatment effect are assumed to fall into a specified interval. Numerical examples are given to illustrate applications of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.

Genetic variants associated with longer telomere length are associated with increased lung cancer risk among never-smoking women in Asia: a report from the female lung cancer consortium in Asia.

Recent evidence from several relatively small nested case-control studies in prospective cohorts shows an association between longer telomere length measured phenotypically in peripheral white blood cell (WBC) DNA and increased lung cancer risk. We sought to further explore this relationship by examining a panel of seven telomere-length associated genetic variants in a large study of 5,457 never-smoking female Asian lung cancer cases and 4,493 never-smoking female Asian controls using data from a previously reported genome-wide association study. Using a group of 1,536 individuals with phenotypically measured telomere length in WBCs in the prospective Shanghai Women's Health study, we demonstrated the utility of a genetic risk score (GRS) of seven telomere-length associated variants to predict telomere length in an Asian population. We then found that GRSs used as instrumental variables to predict longer telomere length were associated with increased lung cancer risk (OR = 1.51 (95% CI = 1.34-1.69) for upper vs. lower quartile of the weighted GRS, p value = 4.54 × 10(-14) ) even after removing rs2736100 (p value = 4.81 × 10(-3) ), a SNP in the TERT locus robustly associated with lung cancer risk in prior association studies. Stratified analyses suggested the effect of the telomere-associated GRS is strongest among younger individuals. We found no difference in GRS effect between adenocarcinoma and squamous cell subtypes. Our results indicate that a genetic background that favors longer telomere length may increase lung cancer risk, which is consistent with earlier prospective studies relating longer telomere length with increased lung cancer risk.

Interactions between household air pollution and GWAS-identified lung cancer susceptibility markers in the Female Lung Cancer Consortium in Asia (FLCCA).

We previously carried out a multi-stage genome-wide association study (GWAS) on lung cancer among never smokers in the Female Lung Cancer Consortium in Asia (FLCCA) (6,609 cases, 7,457 controls) that identified novel susceptibility loci at 10q25.2, 6q22.2, and 6p21.32, and confirmed two previously identified loci at 5p15.33 and 3q28. Household air pollution (HAP) attributed to solid fuel burning for heating and cooking, is the leading cause of the overall disease burden in Southeast Asia, and is known to contain lung carcinogens. To evaluate the gene-HAP interactions associated with lung cancer in loci independent of smoking, we analyzed data from studies participating in FLCCA with fuel use information available (n = 3; 1,731 cases; 1,349 controls). Coal use was associated with a 30% increased risk of lung cancer (OR 1.3, 95% CI 1.0-1.6). Among the five a priori SNPs identified by our GWAS, two showed a significant interaction with coal use (HLA Class II rs2395185, p = 0.02; TP63 rs4488809 (rs4600802), p = 0.04). The risk of lung cancer associated with coal exposure varied with the respective alleles for these two SNPs. Our observations provide evidence that genetic variation in HLA Class II and TP63 may modify the association between HAP and lung cancer risk. The roles played in the cell cycle and inflammation pathways by the proteins encoded by these two genes provide biological plausibility for these interactions; however, additional replication studies are needed in other non-smoking populations.

Quality of life changes in patients undergoing treatment for hepatocellular carcinoma.

PURPOSE: Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related death worldwide. One of the primary treatment goals for incurable advanced cases is to prolong quality of life (QoL). Thus, to determine which HCC therapies may be linked to a more favorable QoL, we assessed the association between QoL changes and different treatments in HCC patients. METHODS: We analyzed a non-randomized multicenter longitudinal study, which included 171 patients treated with surgery (n = 53), ablation (n = 53) or embolization (n = 65) from seven centers: four Asian and three European sites. All participants completed the EORTC QLQ-C30 and QLQ-HCC18 questionnaires before and after treatment. Propensity scores were calculated and used in addition to race for adjustment in the logistic regression model to account for the confounding effects of patient characteristics including age, gender, race, employment, living with family, at least one comorbid condition, years since diagnosis, prior treatment history, BCLC stage, Child-Pugh grade, cirrhosis, bilirubin levels and QoL score before treatment. RESULTS: After adjustment for confounders, patients tended to have higher odds of QoL deterioration when treated with ablation versus embolization (dyspnea: p = 0.019; appetite loss: p = 0.018; body image: p = 0.035) or ablation versus surgery (dyspnea: p = 0.099; appetite loss: p = 0.100; body image: p = 0.038). CONCLUSIONS: There were significant differences in QoL deterioration across different treatment groups. This information may assist patients and providers when selecting patient-centered treatment approaches for HCC.

The association of genetic polymorphisms in the κ-opioid receptor 1 gene with body weight, alcohol use, and withdrawal symptoms in patients with methadone maintenance.

Methadone is a synthetic opioid that binds to the κ-opioid receptor with a low affinity. This study tested the hypotheses that the genetic polymorphisms in the κ-opioid receptor 1 (OPRK1) gene region are associated with methadone treatment responses in a Taiwan methadone maintenance treatment (MMT) cohort. Seventeen single nucleotide polymorphisms (SNPs) in OPRK1 were selected and genotyped on DNA of 366 MMT patients. Six SNPs from rs7843965 to rs1051660 (intron 2 to exon 2) were significantly associated with body weight (P < 0.007). A haplotype of 4 SNPs rs7832417-rs16918853-rs702764-rs7817710 (exon 4 to intron 3) was associated with bone or joint aches (P ≤ 0.004) and with the amount of alcohol use (standard drinks per day; global P < 0.0001). The haplotype rs10958350-rs7016778-rs12675595 was associated with gooseflesh skin (global P < 0.0001), yawning (global P = 0.0001), and restlessness (global P < 0.0001) withdrawal symptoms. The findings suggest that genetic polymorphisms in OPRK1 were associated with the body weight, alcohol use, and opioid withdrawal symptoms in MMT patients.

Celecoxib plus chemoradiotherapy for locally advanced rectal cancer: a phase II TCOG study.

BACKGROUND: To report the results of a phase II trial combining celecoxib and preoperative chemoradiotherapy (CRT) for locally advanced rectal cancer. PATIENTS AND METHODS: Patients with clinical stage II or III rectal cancer were treated with radiotherapy of 44 Gy in 22 fractions. Concurrent chemotherapy consisted of oral tegafur-uracil and folinate on days 1-30 and 38-65. Celecoxib (400 mg/day) given from days 1 to 65. Surgery was done on day 70. The expression of cyclooxygenase 2 (COX-2) in tumor tissues was evaluated microscopically as a prognostic factor. RESULTS: From 2008 to 2011, 53 patients completed CRT+ celecoxib therapy and 47 received radical surgery. Grade 3 diarrhea developed in 5 (9%). Grade 4 anemia was seen in 2 (4%). Pathological complete response (pCR) was seen in 6 (13%). T or N downstaging found in 38 (81%). Sphincter preservation was achieved in 77% of low-positioned tumors. Patients with tumors expressing high-level COX-2 after CRT + celecoxib treatment had inferior pelvic control (P = 0.01), disease-free survival (P = 0.04), and overall survival (P = 0.03) than those with low-level expression. CONCLUSIONS: Celecoxib can be safely combined with preoperative CRT for rectal cancer. More intensified adjuvant therapy may be considered for tumors expressing high-level COX-2 after CRT and surgery.

Sample size estimation of multiregional clinical trials with heterogeneous variability across regions.

In recent years, global collaboration has become a conventional strategy for new drug development. To accelerate the development process and shorten approval time, the design of multiregional clinical trials (MRCTs) incorporates subjects from many countries around the world under the same protocol. After showing the overall efficacy of a drug in all global regions, one can also simultaneously evaluate the possibility of applying the overall trial results to all regions and subsequently support drug registration in each of them. Several statistical methods have been proposed for the design and evaluation of MRCTs. Most of these approaches, however, assume a common variability of the primary endpoint across regions. In practice, this assumption may not be true, due to differences across regions (e.g., differences in ethnic factors and/or medical culture/practice). In this article, we use a random-effect model for modeling heterogeneous variability across regions for the design and evaluation of MRCTs. We also address consideration on the determination of the number of subjects in a specific region to establish the consistency of treatment effects between the specific region and the entire group.

The Clinical Characteristics and Manifestation of Anxious Depression Among Patients With Major Depressive Disorders-Results From a Taiwan Multicenter Study.

OBJECTIVE: Anxious depression is a prevalent characteristic observed in Asian psychiatric patients diagnosed with major depressive disorder (MDD). This study aims to investigate the prevalence and clinical presentation of anxious depression in Taiwanese individuals diagnosed with MDD. METHODS: We recruited psychiatric outpatients aged over 18 who had been diagnosed with MDD through clinical interviews. This recruitment took place at five hospitals located in northern Taiwan. We gathered baseline clinical and demographic information from the participants. Anxious depression was identified using a threshold of an anxiety/somatization factor score ≥7 on the 21-item Hamilton Rating Scale for Depression (HAM-D). RESULTS: In our study of 399 patients (84.21% female), 64.16% met the criteria for anxious depression. They tended to be older, married, less educated, with more children, and an older age of onset. Anxious depression patients had higher HAM-D and Clinical Global Impression-Severity scale score, more panic disorder (without agoraphobia), and exhibited symptoms like agitation, irritability, concentration difficulties, psychological and somatic anxiety, somatic complaints, hypochondriasis, weight loss, and increased insight. Surprisingly, their suicide rates did not significantly differ from non-anxious depression patients. This highlights the importance of recognizing and addressing these unique characteristics. CONCLUSION: Our study findings unveiled that the prevalence of anxious depression among Taiwanese outpatients diagnosed with MDD was lower compared to inpatients but substantially higher than the reported rates in European countries and the United States. Furthermore, patients with anxious depression exhibited a greater occurrence of somatic symptoms.