Archaeogenomic distinctiveness of the Isthmo-Colombian area.

The recently enriched genomic history of Indigenous groups in the Americas is still meager concerning continental Central America. Here, we report ten pre-Hispanic (plus two early colonial) genomes and 84 genome-wide profiles from seven groups presently living in Panama. Our analyses reveal that pre-Hispanic demographic events contributed to the extensive genetic structure currently seen in the area, which is also characterized by a distinctive Isthmo-Colombian Indigenous component. This component drives these populations on a specific variability axis and derives from the local admixture of different ancestries of northern North American origin(s). Two of these ancestries were differentially associated to Pleistocene Indigenous groups that also moved into South America, leaving heterogenous genetic footprints. An additional Pleistocene ancestry was brought by a still unsampled population of the Isthmus (UPopI) that remained restricted to the Isthmian area, expanded locally during the early Holocene, and left genomic traces up to the present day.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci.

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

The Genetic Origin of Daunians and the Pan-Mediterranean Southern Italian Iron Age Context.

The geographical location and shape of Apulia, a narrow land stretching out in the sea at the South of Italy, made this region a Mediterranean crossroads connecting Western Europe and the Balkans. Such movements culminated at the beginning of the Iron Age with the Iapygian civilization which consisted of three cultures: Peucetians, Messapians, and Daunians. Among them, the Daunians left a peculiar cultural heritage, with one-of-a-kind stelae and pottery, but, despite the extensive archaeological literature, their origin has been lost to time. In order to shed light on this and to provide a genetic picture of Iron Age Southern Italy, we collected and sequenced human remains from three archaeological sites geographically located in Northern Apulia (the area historically inhabited by Daunians) and radiocarbon dated between 1157 and 275 calBCE. We find that Iron Age Apulian samples are still distant from the genetic variability of modern-day Apulians, they show a degree of genetic heterogeneity comparable with the cosmopolitan Republican and Imperial Roman civilization, even though a few kilometers and centuries separate them, and they are well inserted into the Iron Age Pan-Mediterranean genetic landscape. Our study provides for the first time a window on the genetic make-up of pre-Roman Apulia, whose increasing connectivity within the Mediterranean landscape, would have contributed to laying the foundation for modern genetic variability. In this light, the genetic profile of Daunians may be compatible with an at least partial autochthonous origin, with plausible contributions from the Balkan peninsula.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy.

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

Analysis of the human Y-chromosome haplogroup Q characterizes ancient population movements in Eurasia and the Americas.

BACKGROUND: Recent genome studies of modern and ancient samples have proposed that Native Americans derive from a subset of the Eurasian gene pool carried to America by an ancestral Beringian population, from which two well-differentiated components originated and subsequently mixed in different proportion during their spread in the Americas. To assess the timing, places of origin and extent of admixture between these components, we performed an analysis of the Y-chromosome haplogroup Q, which is the only Pan-American haplogroup and accounts for virtually all Native American Y chromosomes in Mesoamerica and South America. RESULTS: Our analyses of 1.5 Mb of 152 Y chromosomes, 34 re-sequenced in this work, support a "coastal and inland routes scenario" for the first entrance of modern humans in North America. We show a major phase of male population growth in the Americas after 15 thousand years ago (kya), followed by a period of constant population size from 8 to 3 kya, after which a secondary sign of growth was registered. The estimated dates of the first expansion in Mesoamerica and the Isthmo-Colombian Area, mainly revealed by haplogroup Q-Z780, suggest an entrance in South America prior to 15 kya. During the global constant population size phase, local South American hints of growth were registered by different Q-M848 sub-clades. These expansion events, which started during the Holocene with the improvement of climatic conditions, can be ascribed to multiple cultural changes rather than a steady population growth and a single cohesive culture diffusion as it occurred in Europe. CONCLUSIONS: We established and dated a detailed haplogroup Q phylogeny that provides new insights into the geographic distribution of its Eurasian and American branches in modern and ancient samples.

Y-chromosome and Surname Analyses for Reconstructing Past Population Structures: The Sardinian Population as a Test Case.

Many anthropological, linguistic, genetic and genomic analyses have been carried out to evaluate the potential impact that evolutionary forces had in shaping the present-day Sardinian gene pool, the main outlier in the genetic landscape of Europe. However, due to the homogenizing effect of internal movements, which have intensified over the past fifty years, only partial information has been obtained about the main demographic events. To overcome this limitation, we analyzed the male-specific region of the Y chromosome in three population samples obtained by reallocating a large number of Sardinian subjects to the place of origin of their monophyletic surnames, which are paternally transmitted through generations in most of the populations, much like the Y chromosome. Three Y-chromosome founding lineages, G2-L91, I2-M26 and R1b-V88, were identified as strongly contributing to the definition of the outlying position of Sardinians in the European genetic context and marking a significant differentiation within the island. The present distribution of these lineages does not always mirror that detected in ancient DNAs. Our results show that the analysis of the Y-chromosome gene pool coupled with a sampling method based on the origin of the family name, is an efficient approach to unravelling past heterogeneity, often hidden by recent movements, in the gene pool of modern populations. Furthermore, the reconstruction and comparison of past genetic isolates represent a starting point to better assess the genetic information deriving from the increasing number of available ancient DNA samples.

Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective.

BACKGROUND: Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then. AIM: This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation. SUBJECTS AND METHODS: High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals. RESULTS: Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones. CONCLUSIONS: A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

Phylogeography of human Y-chromosome haplogroup Q3-L275 from an academic/citizen science collaboration.

BACKGROUND: The Y-chromosome haplogroup Q has three major branches: Q1, Q2, and Q3. Q1 is found in both Asia and the Americas where it accounts for about 90% of indigenous Native American Y-chromosomes; Q2 is found in North and Central Asia; but little is known about the third branch, Q3, also named Q1b-L275. Here, we combined the efforts of population geneticists and genetic genealogists to use the potential of full Y-chromosome sequencing for reconstructing haplogroup Q3 phylogeography and suggest possible linkages to events in population history. RESULTS: We analyzed 47 fully sequenced Y-chromosomes and reconstructed the haplogroup Q3 phylogenetic tree in detail. Haplogroup Q3-L275, derived from the oldest known split within Eurasian/American haplogroup Q, most likely occurred in West or Central Asia in the Upper Paleolithic period. During the Mesolithic and Neolithic epochs, Q3 remained a minor component of the West Asian Y-chromosome pool and gave rise to five branches (Q3a to Q3e), which spread across West, Central and parts of South Asia. Around 3-4 millennia ago (Bronze Age), the Q3a branch underwent a rapid expansion, splitting into seven branches, some of which entered Europe. One of these branches, Q3a1, was acquired by a population ancestral to Ashkenazi Jews and grew within this population during the 1st millennium AD, reaching up to 5% in present day Ashkenazi. CONCLUSIONS: This study dataset was generated by a massive Y-chromosome genotyping effort in the genetic genealogy community, and phylogeographic patterns were revealed by a collaboration of population geneticists and genetic genealogists. This positive experience of collaboration between academic and citizen science provides a model for further joint projects. Merging data and skills of academic and citizen science promises to combine, respectively, quality and quantity, generalization and specialization, and achieve a well-balanced and careful interpretation of the paternal-side history of human populations.

Global diversity, recurrent evolution, and recent selection on amylase structural haplotypes in humans.

The adoption of agriculture, first documented ~12,000 years ago in the Fertile Crescent, triggered a rapid shift toward starch-rich diets in human populations. Amylase genes facilitate starch digestion and increased salivary amylase copy number has been observed in some modern human populations with high starch intake, though evidence of recent selection is lacking. Here, using 52 long-read diploid assemblies and short read data from ~5,600 contemporary and ancient humans, we resolve the diversity, evolutionary history, and selective impact of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in populations with agricultural subsistence compared to fishing, hunting, and pastoral groups. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each exhibit multiple duplications/deletions with mutation rates >10,000-fold the SNP mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome graph-based approach to infer structural haplotypes across thousands of humans, we identify extensively duplicated haplotypes present at higher frequencies in modern day populations with traditionally agricultural diets. Leveraging 533 ancient human genomes we find that duplication-containing haplotypes (i.e. haplotypes with more amylase gene copies than the ancestral haplotype) have increased in frequency more than seven-fold over the last 12,000 years providing evidence for recent selection in West Eurasians. Together, our study highlights the potential impacts of the agricultural revolution on human genomes and the importance of long-read sequencing in identifying signatures of selection at structurally complex loci.

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

The multifaceted genomic history of Ashaninka from Amazonian Peru.

Despite its crucial location, the western side of Amazonia between the Andes and the source(s) of the Amazon River is still understudied from a genomic and archaeogenomic point of view, albeit possibly harboring essential information to clarify the complex genetic history of local Indigenous groups and their interactions with nearby regions,1,2,3,4,5,6,7,8 including central America and the Caribbean.9,10,11,12 Focusing on this key region, we analyzed the genome-wide profiles of 51 Ashaninka individuals from Amazonian Peru, observing an unexpected extent of genomic variation. We identified at least two Ashaninka subgroups with distinctive genomic makeups, which were differentially shaped by the degree and timing of external admixtures, especially with the Indigenous groups from the Andes and the Pacific coast. On a continental scale, Ashaninka ancestors probably derived from a south-north migration of Indigenous groups moving into the Amazonian rainforest from a southeastern area with contributions from the Southern Cone and the Atlantic coast. These ancestral populations diversified in the variegated geographic regions of interior South America, on the eastern side of the Andes, differentially interacting with surrounding coastal groups. In this complex scenario, we also revealed strict connections between the ancestors of present-day Ashaninka, who belong to the Arawakan language family,13 and those Indigenous groups that moved further north into the Caribbean, contributing to the early Ceramic (Saladoid) tradition in the islands.14,15.

Identification and Characterization of SOG1 (Suppressor of Gamma Response 1) Homologues in Plants Using Data Mining Resources and Gene Expression Profiling.

SOG1 (Suppressor of the Gamma response 1) is the master-regulator of plant DNA damage response (DDR), a highly coordinated network of DNA damage sensors, transducers, mediators, and effectors, with highly coordinated activities. SOG1 transcription factor belongs to the NAC/NAM protein family, containing the well-conserved NAC domain and five serine-glutamine (SQ) motifs, preferential targets for phosphorylation by ATM and ATR. So far, the information gathered for the SOG1 function comes from studies on the model plant Arabidopsis thaliana. To expand the knowledge on plant-specific DDR, it is opportune to gather information on other SOG1 orthologues. The current study identified plants where multiple SOG1 homologues are present and evaluated their functions by leveraging the information contained in publicly available transcriptomics databases. This analysis revealed the presence of multiple SOG1 sequences in thirteen plant species, and four (Medicago truncatula, Glycine max, Kalankoe fedtschenkoi, Populus trichocarpa) were selected for gene expression data mining based on database availability. Additionally, M. truncatula seeds and seedlings exposed to treatments known to activate DDR pathways were used to evaluate the expression profiles of MtSOG1a and MtSOG1b. The experimental workflow confirmed the data retrieved from transcriptomics datasets, suggesting that the SOG1 homologues have redundant functions in different plant species.

A single-cell transcriptomic landscape of innate and adaptive intratumoral immunity in triple negative breast cancer during chemo- and immunotherapies.

Breast cancer (BC) constitutes a major health problem worldwide, making it the most common malignancy in women. Current treatment options for BC depend primarily on histological type, molecular markers, clinical aggressiveness and stage of disease. Immunotherapy, such as αPD-1, have shown combinatorial clinical activity with chemotherapy in triple negative breast cancer (TNBC) delineating some therapeutic combinations as more effective than others. However, a clear overview of the main immune cell populations involved in these treatments has never been provided.Here, an assessment of the immune landscape in the tumor microenvironment (TME) of two TNBC mouse models has been performed using single-cell RNA sequencing technology. Specifically, immune cells were evaluated in untreated conditions and after treatments with chemotherapy or immunotherapy used as single agents or in combination. A decrease of Treg was found in treatments with in vivo efficacy as well as γδ T cells, which have a pro-tumoral activity in mice. Focusing on Cd8 T cells, across all the conditions, a general increase of exhausted-like Cd8 T cells was confirmed in pre-clinical treatments with low efficacy and an opposite trend was found for the proliferative Cd8 T cells. Regarding macrophages, M2-like cells were enriched in treatments with low efficacy while M1-like macrophages followed an opposite trend. For both models, similar proportions of B cells were detected with an increase of proliferative B cells in treatments involving cisplatin in combination with αPD-1. The fine-scale characterization of the immune TME in this work can lead to new insights on the diagnosis and treatment of TNBC.

Overview of the Americas' First Peopling from a Patrilineal Perspective: New Evidence from the Southern Continent.

Uniparental genetic systems are unique sex indicators and complement the study of autosomal diversity by providing landmarks of human migrations that repeatedly shaped the structure of extant populations. Our knowledge of the variation of the male-specific region of the Y chromosome in Native Americans is still rather scarce and scattered, but by merging sequence information from modern and ancient individuals, we here provide a comprehensive and updated phylogeny of the distinctive Native American branches of haplogroups C and Q. Our analyses confirm C-MPB373, C-P39, Q-Z780, Q-M848, and Q-Y4276 as the main founding haplogroups and identify traces of unsuccessful (pre-Q-F1096) or extinct (C-L1373*, Q-YP4010*) Y-chromosome lineages, indicating that haplogroup diversity of the founder populations that first entered the Americas was greater than that observed in the Indigenous component of modern populations. In addition, through a diachronic and phylogeographic dissection of newly identified Q-M848 branches, we provide the first Y-chromosome insights into the early peopling of the South American hinterland (Q-BY104773 and Q-BY15730) and on overlying inland migrations (Q-BY139813).

Fasting renders immunotherapy effective against low-immunogenic breast cancer while reducing side effects.

Immunotherapy is improving the prognosis and survival of cancer patients, but despite encouraging outcomes in different cancers, the majority of tumors are resistant to it, and the immunotherapy combinations are often accompanied by severe side effects. Here, we show that a periodic fasting-mimicking diet (FMD) can act on the tumor microenvironment and increase the efficacy of immunotherapy (anti-PD-L1 and anti-OX40) against the poorly immunogenic triple-negative breast tumors (TNBCs) by expanding early exhausted effector T cells, switching the cancer metabolism from glycolytic to respiratory, and reducing collagen deposition. Furthermore, FMD reduces the occurrence of immune-related adverse events (irAEs) by preventing the hyperactivation of the immune response. These results indicate that FMD cycles have the potential to enhance the efficacy of anti-cancer immune responses, expand the portion of tumors sensitive to immunotherapy, and reduce its side effects.

Genomic variation in baboons from central Mozambique unveils complex evolutionary relationships with other Papio species.

BACKGROUND: Gorongosa National Park in Mozambique hosts a large population of baboons, numbering over 200 troops. Gorongosa baboons have been tentatively identified as part of Papio ursinus on the basis of previous limited morphological analysis and a handful of mitochondrial DNA sequences. However, a recent morphological and morphometric analysis of Gorongosa baboons pinpointed the occurrence of several traits intermediate between P. ursinus and P. cynocephalus, leaving open the possibility of past and/or ongoing gene flow in the baboon population of Gorongosa National Park. In order to investigate the evolutionary history of baboons in Gorongosa, we generated high and low coverage whole genome sequence data of Gorongosa baboons and compared it to available Papio genomes. RESULTS: We confirmed that P. ursinus is the species closest to Gorongosa baboons. However, the Gorongosa baboon genomes share more derived alleles with P. cynocephalus than P. ursinus does, but no recent gene flow between P. ursinus and P. cynocephalus was detected when available Papio genomes were analyzed. Our results, based on the analysis of autosomal, mitochondrial and Y chromosome data, suggest complex, possibly male-biased, gene flow between Gorongosa baboons and P. cynocephalus, hinting to direct or indirect contributions from baboons belonging to the "northern" Papio clade, and signal the presence of population structure within P. ursinus. CONCLUSIONS: The analysis of genome data generated from baboon samples collected in central Mozambique highlighted a complex set of evolutionary relationships with other baboons. Our results provided new insights in the population dynamics that have shaped baboon diversity.

The Dual Role of Innate Lymphoid and Natural Killer Cells in Cancer. from Phenotype to Single-Cell Transcriptomics, Functions and Clinical Uses.

The role of innate lymphoid cells (ILCs), including natural killer (NK) cells, is pivotal in inflammatory modulation and cancer. Natural killer cell activity and count have been demonstrated to be regulated by the expression of activating and inhibitory receptors together with and as a consequence of different stimuli. The great majority of NK cell populations have an anti-tumor activity due to their cytotoxicity, and for this reason have been used for cellular therapies in cancer patients. On the other hand, the recently classified helper ILCs are fundamentally involved in inflammation and they can be either helpful or harmful in cancer development and progression. Tissue niche seems to play an important role in modulating ILC function and conversion, as observed at the transcriptional level. In the past, these cell populations have been classified by the presence of specific cellular receptor markers; more recently, due to the advent of single-cell RNA sequencing (scRNA-seq), it has been possible to also explore them at the transcriptomic level. In this article we review studies on ILC (and NK cell) classification, function and their involvement in cancer. We also summarize the potential application of NK cells in cancer therapy and give an overview of the most recent studies involving ILCs and NKs at scRNA-seq, focusing on cancer. Finally, we provide a resource for those who wish to start single-cell transcriptomic analysis on the context of these innate lymphoid cell populations.

Drug Repurposing in Oncology, an Attractive Opportunity for Novel Combinatorial Regimens.

The costs of developing, validating and buying new drugs are dramatically increasing. On the other hand, sobering economies have difficulties in sustaining their healthcare systems, particularly in countries with an elderly population requiring increasing welfare. This conundrum requires immediate action, and a possible option is to study the large, already present arsenal of drugs approved and to use them for innovative therapies. This possibility is particularly interesting in oncology, where the complexity of the cancer genome dictates in most patients a multistep therapeutic approach. In this review, we discuss a) Computational approaches; b) preclinical models; c) currently ongoing or already published clinical trials in the drug repurposing field in oncology; and d) drug repurposing to overcome resistance to previous therapies.

Projecting Ancient Ancestry in Modern-Day Arabians and Iranians: A Key Role of the Past Exposed Arabo-Persian Gulf on Human Migrations.

The Arabian Peninsula is strategic for investigations centered on the early structuring of modern humans in the wake of the out-of-Africa migration. Despite its poor climatic conditions for the recovery of ancient human DNA evidence, the availability of both genomic data from neighboring ancient specimens and informative statistical tools allow modeling the ancestry of local modern populations. We applied this approach to a data set of 741,000 variants screened in 291 Arabians and 78 Iranians, and obtained insightful evidence. The west-east axis was a strong forcer of population structure in the Peninsula, and, more importantly, there were clear continuums throughout time linking western Arabia with the Levant, and eastern Arabia with Iran and the Caucasus. Eastern Arabians also displayed the highest levels of the basal Eurasian lineage of all tested modern-day populations, a signal that was maintained even after correcting for a possible bias due to a recent sub-Saharan African input in their genomes. Not surprisingly, eastern Arabians were also the ones with highest similarity with Iberomaurusians, who were, so far, the best proxy for the basal Eurasians amongst the known ancient specimens. The basal Eurasian lineage is the signature of ancient non-Africans who diverged from the common European-eastern Asian pool before 50,000 years ago, prior to the later interbred with Neanderthals. Our results appear to indicate that the exposed basin of the Arabo-Persian Gulf was the possible home of basal Eurasians, a scenario to be further investigated by searching ancient Arabian human specimens.

Cyclophosphamide and Vinorelbine Activate Stem-Like CD8+ T Cells and Improve Anti-PD-1 Efficacy in Triple-Negative Breast Cancer.

Checkpoint inhibitors (CI) instigate anticancer immunity in many neoplastic diseases, albeit only in a fraction of patients. The clinical success of cyclophosphamide (C)-based haploidentical stem-cell transplants indicates that this drug may re-orchestrate the immune system. Using models of triple-negative breast cancer (TNBC) with different intratumoral immune contexture, we demonstrate that a combinatorial therapy of intermittent C, CI, and vinorelbine activates antigen-presenting cells (APC), and abrogates local and metastatic tumor growth by a T-cell-related effect. Single-cell transcriptome analysis of >50,000 intratumoral immune cells after therapy treatment showed a gene signature suggestive of a change resulting from exposure to a mitogen, ligand, or antigen for which it is specific, as well as APC-to-T-cell adhesion. This transcriptional program also increased intratumoral Tcf1+ stem-like CD8+ T cells and altered the balance between terminally and progenitor-exhausted T cells favoring the latter. Overall, our data support the clinical investigation of this therapy in TNBC. SIGNIFICANCE: A combinatorial therapy in mouse models of breast cancer increases checkpoint inhibition by activating antigen-presenting cells, enhancing intratumoral Tcf1+ stem-like CD8+ T cells, and increasing progenitor exhausted CD8+ T cells.