Cancer Evolution: No Room for Negative Selection.

In this issue of Cell, Martincorena et al. and Campbell et al. interrogated the selection dynamics during tumor evolution using large-scale genomics datasets. They found that somatic mutations in cancer are largely neutral, highlighting a near-complete absence of negative selection. Neutral evolution enables tolerance of hypermutation, which defines a surprisingly large fraction of adult cancer.

Mapping genotypes to chromatin accessibility profiles in single cells.

In somatic tissue differentiation, chromatin accessibility changes govern priming and precursor commitment towards cellular fates1-3. Therefore, somatic mutations are likely to alter chromatin accessibility patterns, as they disrupt differentiation topologies leading to abnormal clonal outgrowth. However, defining the impact of somatic mutations on the epigenome in human samples is challenging due to admixed mutated and wild-type cells. Here, to chart how somatic mutations disrupt epigenetic landscapes in human clonal outgrowths, we developed genotyping of targeted loci with single-cell chromatin accessibility (GoT-ChA). This high-throughput platform links genotypes to chromatin accessibility at single-cell resolution across thousands of cells within a single assay. We applied GoT-ChA to CD34+ cells from patients with myeloproliferative neoplasms with JAK2V617F-mutated haematopoiesis. Differential accessibility analysis between wild-type and JAK2V617F-mutant progenitors revealed both cell-intrinsic and cell-state-specific shifts within mutant haematopoietic precursors, including cell-intrinsic pro-inflammatory signatures in haematopoietic stem cells, and a distinct profibrotic inflammatory chromatin landscape in megakaryocytic progenitors. Integration of mitochondrial genome profiling and cell-surface protein expression measurement allowed expansion of genotyping onto DOGMA-seq through imputation, enabling single-cell capture of genotypes, chromatin accessibility, RNA expression and cell-surface protein expression. Collectively, we show that the JAK2V617F mutation leads to epigenetic rewiring in a cell-intrinsic and cell type-specific manner, influencing inflammation states and differentiation trajectories. We envision that GoT-ChA will empower broad future investigations of the critical link between somatic mutations and epigenetic alterations across clonal populations in malignant and non-malignant contexts.

Targeting histone acetylation dynamics and oncogenic transcription by catalytic P300/CBP inhibition.

To separate causal effects of histone acetylation on chromatin accessibility and transcriptional output, we used integrated epigenomic and transcriptomic analyses following acute inhibition of major cellular lysine acetyltransferases P300 and CBP in hematological malignancies. We found that catalytic P300/CBP inhibition dynamically perturbs steady-state acetylation kinetics and suppresses oncogenic transcriptional networks in the absence of changes to chromatin accessibility. CRISPR-Cas9 screening identified NCOR1 and HDAC3 transcriptional co-repressors as the principal antagonists of P300/CBP by counteracting acetylation turnover kinetics. Finally, deacetylation of H3K27 provides nucleation sites for reciprocal methylation switching, a feature that can be exploited therapeutically by concomitant KDM6A and P300/CBP inhibition. Overall, this study indicates that the steady-state histone acetylation-methylation equilibrium functions as a molecular rheostat governing cellular transcription that is amenable to therapeutic exploitation as an anti-cancer regimen.

Evolution and impact of subclonal mutations in chronic lymphocytic leukemia.

Clonal evolution is a key feature of cancer progression and relapse. We studied intratumoral heterogeneity in 149 chronic lymphocytic leukemia (CLL) cases by integrating whole-exome sequence and copy number to measure the fraction of cancer cells harboring each somatic mutation. We identified driver mutations as predominantly clonal (e.g., MYD88, trisomy 12, and del(13q)) or subclonal (e.g., SF3B1 and TP53), corresponding to earlier and later events in CLL evolution. We sampled leukemia cells from 18 patients at two time points. Ten of twelve CLL cases treated with chemotherapy (but only one of six without treatment) underwent clonal evolution, predominantly involving subclones with driver mutations (e.g., SF3B1 and TP53) that expanded over time. Furthermore, presence of a subclonal driver mutation was an independent risk factor for rapid disease progression. Our study thus uncovers patterns of clonal evolution in CLL, providing insights into its stepwise transformation, and links the presence of subclones with adverse clinical outcomes.

Integrating genetic and non-genetic determinants of cancer evolution by single-cell multi-omics.

Cancer represents an evolutionary process through which growing malignant populations genetically diversify, leading to tumour progression, relapse and resistance to therapy. In addition to genetic diversity, the cell-to-cell variation that fuels evolutionary selection also manifests in cellular states, epigenetic profiles, spatial distributions and interactions with the microenvironment. Therefore, the study of cancer requires the integration of multiple heritable dimensions at the resolution of the single cell - the atomic unit of somatic evolution. In this Review, we discuss emerging analytic and experimental technologies for single-cell multi-omics that enable the capture and integration of multiple data modalities to inform the study of cancer evolution. These data show that cancer results from a complex interplay between genetic and non-genetic determinants of somatic evolution.

Tristetraprolin impairs myc-induced lymphoma and abolishes the malignant state.

Myc oncoproteins directly regulate transcription by binding to target genes, yet this only explains a fraction of the genes affected by Myc. mRNA turnover is controlled via AU-binding proteins (AUBPs) that recognize AU-rich elements (AREs) found within many transcripts. Analyses of precancerous and malignant Myc-expressing B cells revealed that Myc regulates hundreds of ARE-containing (ARED) genes and select AUBPs. Notably, Myc directly suppresses transcription of Tristetraprolin (TTP/ZFP36), an mRNA-destabilizing AUBP, and this circuit is also operational during B lymphopoiesis and IL7 signaling. Importantly, TTP suppression is a hallmark of cancers with MYC involvement, and restoring TTP impairs Myc-induced lymphomagenesis and abolishes maintenance of the malignant state. Further, there is a selection for TTP loss in malignancy; thus, TTP functions as a tumor suppressor. Finally, Myc/TTP-directed control of select cancer-associated ARED genes is disabled during lymphomagenesis. Thus, Myc targets AUBPs to regulate ARED genes that control tumorigenesis.

TARGET-seq Takes Aim at Cancer Evolution through Multi-omics Single-Cell Genotyping and Transcriptomics.

In this issue of Molecular Cell, Rodriguez-Meira et al. (2019) present TARGET-seq, an elegant single-cell method that genotypes somatic mutations and captures whole transcriptomes in the same tumor cells, thus paving the way to directly link somatic mutations with resulting transcriptional phenotypes in clonally diverse cancer populations.

Loss of the stress sensor GADD45A promotes stem cell activity and ferroptosis resistance in LGR4/HOXA9-dependent AML.

ABSTRACT: The overall prognosis of acute myeloid leukemia (AML) remains dismal, largely because of the inability of current therapies to kill leukemia stem cells (LSCs) with intrinsic resistance. Loss of the stress sensor growth arrest and DNA damage-inducible 45 alpha (GADD45A) is implicated in poor clinical outcomes, but its role in LSCs and AML pathogenesis is unknown. Here, we define GADD45A as a key downstream target of G protein-coupled receptor (LGR)4 pathway and discover a regulatory role for GADD45A loss in promoting leukemia-initiating activity and oxidative resistance in LGR4/HOXA9-dependent AML, a poor prognosis subset of leukemia. Knockout of GADD45A enhances AML progression in murine and patient-derived xenograft (PDX) mouse models. Deletion of GADD45A induces substantial mutations, increases LSC self-renewal and stemness in vivo, and reduces levels of reactive oxygen species (ROS), accompanied by a decreased response to ROS-associated genotoxic agents (eg, ferroptosis inducer RSL3) and acquisition of an increasingly aggressive phenotype on serial transplantation in mice. Our single-cell cellular indexing of transcriptomes and epitopes by sequencing analysis on patient-derived LSCs in PDX mice and subsequent functional studies in murine LSCs and primary AML patient cells show that loss of GADD45A is associated with resistance to ferroptosis (an iron-dependent oxidative cell death caused by ROS accumulation) through aberrant activation of antioxidant pathways related to iron and ROS detoxification, such as FTH1 and PRDX1, upregulation of which correlates with unfavorable outcomes in patients with AML. These results reveal a therapy resistance mechanism contributing to poor prognosis and support a role for GADD45A loss as a critical step for leukemia-initiating activity and as a target to overcome resistance in aggressive leukemia.

Clustered versus catastrophic global vertebrate declines.

Recent analyses have reported catastrophic global declines in vertebrate populations1,2. However, the distillation of many trends into a global mean index obscures the variation that can inform conservation measures and can be sensitive to analytical decisions. For example, previous analyses have estimated a mean vertebrate decline of more than 50% since 1970 (Living Planet Index2). Here we show, however, that this estimate is driven by less than 3% of vertebrate populations; if these extremely declining populations are excluded, the global trend switches to an increase. The sensitivity of global mean trends to outliers suggests that more informative indices are needed. We propose an alternative approach, which identifies clusters of extreme decline (or increase) that differ statistically from the majority of population trends. We show that, of taxonomic-geographic systems in the Living Planet Index, 16 systems contain clusters of extreme decline (comprising around 1% of populations; these extreme declines occur disproportionately in larger animals) and 7 contain extreme increases (around 0.4% of populations). The remaining 98.6% of populations across all systems showed no mean global trend. However, when analysed separately, three systems were declining strongly with high certainty (all in the Indo-Pacific region) and seven were declining strongly but with less certainty (mostly reptile and amphibian groups). Accounting for extreme clusters fundamentally alters the interpretation of global vertebrate trends and should be used to help to prioritize conservation efforts.

An Mtb-Human Protein-Protein Interaction Map Identifies a Switch between Host Antiviral and Antibacterial Responses.

Although macrophages are armed with potent antibacterial functions, Mycobacterium tuberculosis (Mtb) replicates inside these innate immune cells. Determinants of macrophage intrinsic bacterial control, and the Mtb strategies to overcome them, are poorly understood. To further study these processes, we used an affinity tag purification mass spectrometry (AP-MS) approach to identify 187 Mtb-human protein-protein interactions (PPIs) involving 34 secreted Mtb proteins. This interaction map revealed two factors involved in Mtb pathogenesis-the secreted Mtb protein, LpqN, and its binding partner, the human ubiquitin ligase CBL. We discovered that an lpqN Mtb mutant is attenuated in macrophages, but growth is restored when CBL is removed. Conversely, Cbl-/- macrophages are resistant to viral infection, indicating that CBL regulates cell-intrinsic polarization between antibacterial and antiviral immunity. Collectively, these findings illustrate the utility of this Mtb-human PPI map for developing a deeper understanding of the intricate interactions between Mtb and its host.

Kelp carbon sink potential decreases with warming due to accelerating decomposition.

Cycling of organic carbon in the ocean has the potential to mitigate or exacerbate global climate change, but major questions remain about the environmental controls on organic carbon flux in the coastal zone. Here, we used a field experiment distributed across 28° of latitude, and the entire range of 2 dominant kelp species in the northern hemisphere, to measure decomposition rates of kelp detritus on the seafloor in relation to local environmental factors. Detritus decomposition in both species were strongly related to ocean temperature and initial carbon content, with higher rates of biomass loss at lower latitudes with warmer temperatures. Our experiment showed slow overall decomposition and turnover of kelp detritus and modeling of coastal residence times at our study sites revealed that a significant portion of this production can remain intact long enough to reach deep marine sinks. The results suggest that decomposition of these kelp species could accelerate with ocean warming and that low-latitude kelp forests could experience the greatest increase in remineralization with a 9% to 42% reduced potential for transport to long-term ocean sinks under short-term (RCP4.5) and long-term (RCP8.5) warming scenarios. However, slow decomposition at high latitudes, where kelp abundance is predicted to expand, indicates potential for increasing kelp-carbon sinks in cooler (northern) regions. Our findings reveal an important latitudinal gradient in coastal ecosystem function that provides an improved capacity to predict the implications of ocean warming on carbon cycling. Broad-scale patterns in organic carbon decomposition revealed here can be used to identify hotspots of carbon sequestration potential and resolve relationships between carbon cycling processes and ocean climate at a global scale.

Somatic mutations and cell identity linked by Genotyping of Transcriptomes.

Defining the transcriptomic identity of malignant cells is challenging in the absence of surface markers that distinguish cancer clones from one another, or from admixed non-neoplastic cells. To address this challenge, here we developed Genotyping of Transcriptomes (GoT), a method to integrate genotyping with high-throughput droplet-based single-cell RNA sequencing. We apply GoT to profile 38,290 CD34+ cells from patients with CALR-mutated myeloproliferative neoplasms to study how somatic mutations corrupt the complex process of human haematopoiesis. High-resolution mapping of malignant versus normal haematopoietic progenitors revealed an increasing fitness advantage with myeloid differentiation of cells with mutated CALR. We identified the unfolded protein response as a predominant outcome of CALR mutations, with a considerable dependency on cell identity, as well as upregulation of the NF-κB pathway specifically in uncommitted stem cells. We further extended the GoT toolkit to genotype multiple targets and loci that are distant from transcript ends. Together, these findings reveal that the transcriptional output of somatic mutations in myeloproliferative neoplasms is dependent on the native cell identity.

Epigenetic evolution and lineage histories of chronic lymphocytic leukaemia.

Genetic and epigenetic intra-tumoral heterogeneity cooperate to shape the evolutionary course of cancer1. Chronic lymphocytic leukaemia (CLL) is a highly informative model for cancer evolution as it undergoes substantial genetic diversification and evolution after therapy2,3. The CLL epigenome is also an important disease-defining feature4,5, and growing populations of cells in CLL diversify by stochastic changes in DNA methylation known as epimutations6. However, previous studies using bulk sequencing methods to analyse the patterns of DNA methylation were unable to determine whether epimutations affect CLL populations homogeneously. Here, to measure the epimutation rate at single-cell resolution, we applied multiplexed single-cell reduced-representation bisulfite sequencing to B cells from healthy donors and patients with CLL. We observed that the common clonal origin of CLL results in a consistently increased epimutation rate, with low variability in the cell-to-cell epimutation rate. By contrast, variable epimutation rates across healthy B cells reflect diverse evolutionary ages across the trajectory of B cell differentiation, consistent with epimutations serving as a molecular clock. Heritable epimutation information allowed us to reconstruct lineages at high-resolution with single-cell data, and to apply this directly to patient samples. The CLL lineage tree shape revealed earlier branching and longer branch lengths than in normal B cells, reflecting rapid drift after the initial malignant transformation and a greater proliferative history. Integration of single-cell bisulfite sequencing analysis with single-cell transcriptomes and genotyping confirmed that genetic subclones mapped to distinct clades, as inferred solely on the basis of epimutation information. Finally, to examine potential lineage biases during therapy, we profiled serial samples during ibrutinib-associated lymphocytosis, and identified clades of cells that were preferentially expelled from the lymph node after treatment, marked by distinct transcriptional profiles. The single-cell integration of genetic, epigenetic and transcriptional information thus charts the lineage history of CLL and its evolution with therapy.

Deciphering the Exact Sequence of Endogenous Soluble B Cell Maturation Antigen and Unbiased Quantitation in Multiple Myeloma Patient Samples by LC-MS.

BACKGROUND: B-cell maturation antigen is a pivotal therapeutic target for multiple myeloma (MM). Membrane-bound BCMA can be cleaved by γ-secretase and shed as soluble BCMA (sBCMA). sBCMA can act as a neutralizing sink to compete with drug, as well as serve as a diagnostic/prognostic biomarker for MM. Antibody-capture based methods, such as enzyme-linked immunosorbent assay (ELISA) and immunoaffinity-liquid chromatography-multiple reaction monitoring (IA-LC-MRM), have been reported and well adopted to measure sBCMA in clinical samples. However, both methods are biased by capturing antibodies. METHODS: We have used various LC-MS workflows to characterize and quantify endogenous sBCMA in MM patient samples, including bottom-up peptide mapping, intact analysis, IA-based, and reagent-free (RF)-LC-MRM quantitation. RESULTS: We have confirmed that sBCMA contains a variable N-terminus and a C-terminus that extends to the transmembrane domain, ending at amino acid 61. Leveraging an in-house synthesized G-1-61 sBCMA recombinant standard, we developed a RF-LC-MRM method for unbiased sBCMA quantitation in MM patient samples. By comparing the results from RF-LC-MRM with ELISA and IA-LC-MRM, we demonstrated that RF-LC-MRM measures a more complete pool of endogenous sBCMA compared to the antibody-based methods. CONCLUSIONS: This work fills the knowledge gap of the exact sequence of endogenous sBCMA for the first time, which differs from the current commercially available standard. Additionally, this work highlights the necessity of identifying the actual sequence of an endogenous soluble target such as sBCMA, both for bioanalytical purposes and to underpin pharmacodynamic measurements.

The ERK2-DBP domain opposes pathogenesis of a mouse JAK2V617F-driven myeloproliferative neoplasm.

Although Ras/mitogen-activated protein kinase (MAPK) signaling is activated in most human cancers, attempts to target this pathway using kinase-active site inhibitors have not typically led to durable clinical benefit. To address this shortcoming, we sought to test the feasibility of an alternative targeting strategy, focused on the ERK2 substrate binding domains, D and DEF binding pocket (DBP). Disabling the ERK2-DBP domain in mice caused baseline erythrocytosis. Consequently, we investigated the role of the ERK2-D and -DBP domains in disease, using a JAK2-dependent model of polycythemia vera (PV). Of note, inactivation of the ERK2-DBP domain promoted the progression of disease from PV to myelofibrosis, suggesting that the ERK2-DBP domain normally opposes progression. ERK2-DBP inactivation also prevented oncogenic JAK2 kinase (JAK2V617F) from promoting oncogene-induced senescence in vitro. The ERK2-DBP mutation attenuated JAK2-mediated oncogene-induced senescence by preventing the physical interaction of ERK2 with the transcription factor Egr1. Because inactivation of the ERK2-DBP created a functional ERK2 kinase limited to binding substrates through its D domain, these data suggested that the D domain substrates were responsible for promoting oncogene-induced progenitor growth and tumor progression and that pharmacologic targeting of the ERK2-D domain may attenuate cancer cell growth. Indeed, pharmacologic agents targeting the ERK2-D domain were effective in attenuating the growth of JAK2-dependent myeloproliferative neoplasm cell lines. Taken together, these data indicate that the ERK-D and -DBP domains can play distinct roles in the progression of neoplasms and that the D domain has the potential to be a potent therapeutic target in Ras/MAPK-dependent cancers.

Correlation of In Vitro Kinetic Stability to Preclinical In Vivo Pharmacokinetics for a Panel of Anti-PD-1 Monoclonal Antibody Interleukin 21 Mutein Immunocytokines.

The development of therapeutic fusion protein drugs is often impeded by the unintended consequences that occur from fusing together domains from independent naturally occurring proteins, consequences such as altered biodistribution, tissue uptake, or rapid clearance and potential immunogenicity. For therapeutic fusion proteins containing globular domains, we hypothesized that aberrant in vivo behavior could be related to low kinetic stability of these domains leading to local unfolding and susceptibility to partial proteolysis and/or salvage and uptake. Herein we describe an assay to measure kinetic stability of therapeutic fusion proteins by way of their sensitivity to the protease thermolysin. The results indicate that in vivo pharmacokinetics of a panel of anti-programmed cell death protein 1 monocolonal antibody:interleukin 21 immunocytokines in both mice and nonhuman primates are highly correlated with their in vitro susceptibility to thermolysin-mediated proteolysis. This assay can be used as a tool to quickly identify in vivo liabilities of globular domains of therapeutic proteins, thus aiding in the optimization and development of new multispecific drug candidates. SIGNIFICANCE STATEMENT: This work describes a novel assay utilizing protein kinetic stability to identify preclinical in vivo pharmacokinetic liabilities of multispecific therapeutic fusion proteins. This provides an efficient, inexpensive method to ascertain inherent protein stability in vitro before conducting in vivo studies, which can rapidly increase the speed of preclinical drug development.

Spatiotemporal variability in population demography and morphology of the habitat-forming macroalga Saccorhiza polyschides in the Western English Channel.

BACKGROUND AND AIMS: Large brown macroalgae serve as foundation organisms along temperate and polar coastlines, providing a range of ecosystem services. Saccorhiza polyschides is a warm-temperate kelp-like species found in the northeast Atlantic, which is suggested to have proliferated in recent decades across the southern UK, possibly in response to increasing temperatures, physical disturbance and reduced competition. However, little is known about S. polyschides with regard to ecological functioning and population dynamics across its geographical range. Here we examined the population demography of S. polyschides populations in southwest UK, located within the species' range centre, to address a regional knowledge gap and to provide a baseline against which to detect future changes. METHODS: Intertidal surveys were conducted during spring low tides at three sites along a gradient of wave exposure in Plymouth Sound (Western English Channel) over a period of 15 months. Density, cover, age, biomass and morphology of S. polyschides were quantified. Additionally, less frequent sampling of shallow subtidal reefs was conducted to compare intertidal and subtidal populations. KEY RESULTS: We recorded pronounced seasonality, with fairly consistent demographic patterns across sites and depths. By late summer, S. polyschides was a dominant habitat-former on both intertidal and subtidal reefs, with maximum standing stock exceeding 13 000 g wet weight m-2. CONCLUSIONS: Saccorhiza polyschides is a conspicuous and abundant member of rocky reef assemblages in the region, providing complex and abundant biogenic habitat for associated organisms and high rates of primary productivity. However, its short-lived pseudo-annual life strategy is in stark contrast to dominant long-lived perennial laminarian kelps. As such, any replacement or reconfiguration of habitat-forming macroalgae due to ocean warming will probably have implications for local biodiversity and community composition. More broadly, our study demonstrates the importance of high-resolution cross-habitat surveys to generate robust baselines of kelp population demography, against which the ecological impacts of climate change and other stressors can be reliably detected.

Marine heatwave duration and intensity interact to reduce physiological tipping points of kelp species with contrasting thermal affinities.

BACKGROUND AND AIMS: Marine heatwaves (MHWs) are widely recognized as pervasive drivers of ecosystem change, yet our understanding of how different MHW properties mediate ecological responses remains largely unexplored. Understanding MHW impacts on foundation species is particularly important, given their structural role in communities and ecosystems. METHODS: We simulated a series of realistic MHWs with different levels of intensity (Control: 14 °C, Moderate: 18 °C, Extreme: 22 °C) and duration (14 or 28 d) and examined responses of two habitat-forming kelp species in the southwest UK. Here, Laminaria digitata reaches its trailing edge and is undergoing a range contraction, whereas Laminaria ochroleuca reaches its leading edge and is undergoing a range expansion. KEY RESULTS: For both species, sub-lethal stress responses induced by moderate-intensity MHWs were exacerbated by longer duration. Extreme-intensity MHWs caused dramatic declines in growth and photosynthetic performance, and elevated bleaching, which were again exacerbated by longer MHW duration. Stress responses were most pronounced in L. ochroleuca, where almost complete tissue necrosis was observed by the end of the long-duration MHW. This was unexpected given the greater thermal safety margins assumed with leading edge populations. It is likely that prolonged exposure to sub-lethal thermal stress exceeded a physiological tipping point for L. ochroleuca, presumably due to depletion of internal reserves. CONCLUSIONS: Overall, our study showed that exposure to MHW profiles projected to occur in the region in the coming decades can have significant deleterious effects on foundation kelp species, regardless of their thermal affinities and location within respective latitudinal ranges, which would probably have consequences for entire communities and ecosystems.

Variation in farming damselfish behaviour creates a competitive landscape of risk on coral reefs.

Interspecific interactions are fundamental drivers of animal space use. Yet while non-consumptive effects of predation risk on prey space use are well-known, the risk of aggressive interactions on space use of competitors is largely unknown. We apply the landscape of risk framework to competition-driven space use for the first time, with the hypothesis that less aggressive competitors may alter their behaviour to avoid areas of high competitor density. Specifically, we test how aggressive risk from territorial algal-farming damselfishes can shape the spatial distribution of herbivore fish competitors. We found that only the most aggressive damselfish had fewer competitors in their surrounding area, demonstrating that individual-level behavioural variation can shape spatial distributions. In contradiction to the landscape of risk framework, abundances of farming damselfish and other fishes were positively associated. Our results suggest that reef fishes do not simply avoid areas of high damselfish abundance, but that spatial variation in aggressive behaviour, rather than of individuals, created a competitive landscape of risk. We emphasize the importance of individual-level behaviour in identifying patterns of space use and propose expanding the landscape of risk framework to non-predatory interactions to explore cascading behavioural responses to aggressive risk.

Can Medication be a Factor that Can Negatively Affect the Effect of Transcranial Magnetic Stimulation in Depression?

BACKGROUND: This study aims to evaluate the efficacy of transcranial magnetic stimulation (TMS) in patients with depression and whether concurrent psychotropic medication use negatively affects the treatment outcome of TMS. Patients' characteristics, predictors of treatment response, the relationship between demographics, and the selection of TMS as a treatment modality were also analyzed. STUDY QUESTION: Can psychotropic medication be a factor that can negatively affect the efficacy of TMS in patients with depression? STUDY DESIGN: This pilot-controlled study included 40 subjects from Romanian clinical practice who were treated with pharmacological treatment and TMS for major depressive disorder. The severity of depression and anxiety symptoms was measured using validated scales at baseline (day 1) and follow-up (day 30). DATA SOURCES: All patients' characteristics and information were collected manually from the clinic's medical records, deidentified, and then introduced into an electronic database for analysis. LIMITATIONS: Conducting the study in a clinical routine practice, it was not possible to include an active and/or sham control group. In addition, because TMS is not used as a monotherapy in this type of practice, we could not evaluate its safety and efficacy without concomitant pharmacological treatment. The study sample is small; therefore, the results cannot be generalized. RESULTS: Sixty percentage of patients (n = 24) included in this study obtained a clinical response, and 30% of patients (n = 12) obtained remission of depression. The group with pharmacological treatment obtained clinical responses in 80% of patients (n = 16) and remission of depression in 45% of patients (n = 9). The group with pharmacological treatment and TMS obtained clinical responses in 40% of patients (n = 8) and remission of depression in 15% (n = 3) of cases. CONCLUSIONS: The study results show a lack of efficacy for TMS as an adjunctive therapy to pharmacological treatment for patients with depression. In addition, a negative impact of psychotropic medication on TMS efficacy is observed in our study sample.