Publisher Correction: The immune cell landscape in kidneys of patients with lupus nephritis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The immune cell landscape in kidneys of patients with lupus nephritis.

Lupus nephritis is a potentially fatal autoimmune disease for which the current treatment is ineffective and often toxic. To develop mechanistic hypotheses of disease, we analyzed kidney samples from patients with lupus nephritis and from healthy control subjects using single-cell RNA sequencing. Our analysis revealed 21 subsets of leukocytes active in disease, including multiple populations of myeloid cells, T cells, natural killer cells and B cells that demonstrated both pro-inflammatory responses and inflammation-resolving responses. We found evidence of local activation of B cells correlated with an age-associated B-cell signature and evidence of progressive stages of monocyte differentiation within the kidney. A clear interferon response was observed in most cells. Two chemokine receptors, CXCR4 and CX3CR1, were broadly expressed, implying a potentially central role in cell trafficking. Gene expression of immune cells in urine and kidney was highly correlated, which would suggest that urine might serve as a surrogate for kidney biopsies.

Spatial maps of T cell receptors and transcriptomes reveal distinct immune niches and interactions in the adaptive immune response.

T cells mediate antigen-specific immune responses to disease through the specificity and diversity of their clonotypic T cell receptors (TCRs). Determining the spatial distributions of T cell clonotypes in tissues is essential to understanding T cell behavior, but spatial sequencing methods remain unable to profile the TCR repertoire. Here, we developed Slide-TCR-seq, a 10-µm-resolution method, to sequence whole transcriptomes and TCRs within intact tissues. We confirmed the ability of Slide-TCR-seq to map the characteristic locations of T cells and their receptors in mouse spleen. In human lymphoid germinal centers, we identified spatially distinct TCR repertoires. Profiling T cells in renal cell carcinoma and melanoma specimens revealed heterogeneous immune responses: T cell states and infiltration differed intra- and inter-clonally, and adjacent tumor and immune cells exhibited distinct gene expression. Altogether, our method yields insights into the spatial relationships between clonality, neighboring cell types, and gene expression that drive T cell responses.

Dictionary of immune responses to cytokines at single-cell resolution.

Cytokines mediate cell-cell communication in the immune system and represent important therapeutic targets1-3. A myriad of studies have highlighted their central role in immune function4-13, yet we lack a global view of the cellular responses of each immune cell type to each cytokine. To address this gap, we created the Immune Dictionary, a compendium of single-cell transcriptomic profiles of more than 17 immune cell types in response to each of 86 cytokines (>1,400 cytokine-cell type combinations) in mouse lymph nodes in vivo. A cytokine-centric view of the dictionary revealed that most cytokines induce highly cell-type-specific responses. For example, the inflammatory cytokine interleukin-1ß induces distinct gene programmes in almost every cell type. A cell-type-centric view of the dictionary identified more than 66 cytokine-driven cellular polarization states across immune cell types, including previously uncharacterized states such as an interleukin-18-induced polyfunctional natural killer cell state. Based on this dictionary, we developed companion software, Immune Response Enrichment Analysis, for assessing cytokine activities and immune cell polarization from gene expression data, and applied it to reveal cytokine networks in tumours following immune checkpoint blockade therapy. Our dictionary generates new hypotheses for cytokine functions, illuminates pleiotropic effects of cytokines, expands our knowledge of activation states of each immune cell type, and provides a framework to deduce the roles of specific cytokines and cell-cell communication networks in any immune response.

Latent human herpesvirus 6 is reactivated in CAR T cells.

Cell therapies have yielded durable clinical benefits for patients with cancer, but the risks associated with the development of therapies from manipulated human cells are understudied. For example, we lack a comprehensive understanding of the mechanisms of toxicities observed in patients receiving T cell therapies, including recent reports of encephalitis caused by reactivation of human herpesvirus 6 (HHV-6)1. Here, through petabase-scale viral genomics mining, we examine the landscape of human latent viral reactivation and demonstrate that HHV-6B can become reactivated in cultures of human CD4+ T cells. Using single-cell sequencing, we identify a rare population of HHV-6 'super-expressors' (about 1 in 300-10,000 cells) that possess high viral transcriptional activity, among research-grade allogeneic chimeric antigen receptor (CAR) T cells. By analysing single-cell sequencing data from patients receiving cell therapy products that are approved by the US Food and Drug Administration2 or are in clinical studies3-5, we identify the presence of HHV-6-super-expressor CAR T cells in patients in vivo. Together, the findings of our study demonstrate the utility of comprehensive genomics analyses in implicating cell therapy products as a potential source contributing to the lytic HHV-6 infection that has been reported in clinical trials1,6-8 and may influence the design and production of autologous and allogeneic cell therapies.

Landscape of helper and regulatory antitumour CD4+ T cells in melanoma.

Within the tumour microenvironment, CD4+ T cells can promote or suppress antitumour responses through the recognition of antigens presented by human leukocyte antigen (HLA) class II molecules1,2, but how cancers co-opt these physiologic processes to achieve immune evasion remains incompletely understood. Here we performed in-depth analysis of the phenotype and tumour specificity of CD4+ T cells infiltrating human melanoma specimens, finding that exhausted cytotoxic CD4+ T cells could be directly induced by melanoma cells through recognition of HLA class II-restricted neoantigens, and also HLA class I-restricted tumour-associated antigens. CD4+ T regulatory (TReg) cells could be indirectly elicited through presentation of tumour antigens via antigen-presenting cells. Notably, numerous tumour-reactive CD4+ TReg clones were stimulated directly by HLA class II-positive melanoma and demonstrated specificity for melanoma neoantigens. This phenomenon was observed in the presence of an extremely high tumour neoantigen load, which we confirmed to be associated with HLA class II positivity through the analysis of 116 melanoma specimens. Our data reveal the landscape of infiltrating CD4+ T cells in melanoma and point to the presentation of HLA class II-restricted neoantigens and direct engagement of immunosuppressive CD4+ TReg cells as a mechanism of immune evasion that is favoured in HLA class II-positive melanoma.

Phenotype, specificity and avidity of antitumour CD8+ T cells in melanoma.

Interactions between T cell receptors (TCRs) and their cognate tumour antigens are central to antitumour immune responses1-3; however, the relationship between phenotypic characteristics and TCR properties is not well elucidated. Here we show, by linking the antigenic specificity of TCRs and the cellular phenotype of melanoma-infiltrating lymphocytes at single-cell resolution, that tumour specificity shapes the expression state of intratumoural CD8+ T cells. Non-tumour-reactive T cells were enriched for viral specificities and exhibited a non-exhausted memory phenotype, whereas melanoma-reactive lymphocytes predominantly displayed an exhausted state that encompassed diverse levels of differentiation but rarely acquired memory properties. These exhausted phenotypes were observed both among clonotypes specific for public overexpressed melanoma antigens (shared across different tumours) or personal neoantigens (specific for each tumour). The recognition of such tumour antigens was provided by TCRs with avidities inversely related to the abundance of cognate targets in melanoma cells and proportional to the binding affinity of peptide-human leukocyte antigen (HLA) complexes. The persistence of TCR clonotypes in peripheral blood was negatively affected by the level of intratumoural exhaustion, and increased in patients with a poor response to immune checkpoint blockade, consistent with chronic stimulation mediated by residual tumour antigens. By revealing how the quality and quantity of tumour antigens drive the features of T cell responses within the tumour microenvironment, we gain insights into the properties of the anti-melanoma TCR repertoire.

Differentiation of BCMA-specific induced pluripotent stem cells into rejuvenated CD8αß+ T cells targeting multiple myeloma.

ABSTRACT: A major hurdle in adoptive T-cell therapy is cell exhaustion and failure to maintain antitumor responses. Here, we introduce an induced pluripotent stem cell (iPSC) strategy for reprogramming and revitalizing precursor exhausted B-cell maturation antigen (BCMA)-specific T cells to effectively target multiple myeloma (MM). Heteroclitic BCMA72-80 (YLMFLLRKI)-specific CD8+ memory cytotoxic T lymphocytes (CTL) were epigenetically reprogrammed to a pluripotent state, developed into hematopoietic progenitor cells (CD34+ CD43+/CD14- CD235a-), differentiated into the T-cell lineage and evaluated for their polyfunctional activities against MM. The final T-cell products demonstrated (1) mature CD8αß+ memory phenotype, (2) high expression of activation or costimulatory molecules (CD38, CD28, and 41BB), (3) no expression of immune checkpoint and senescence markers (CTLA4, PD1, LAG3, and TIM3; CD57), and (4) robust proliferation and polyfunctional immune responses to MM. The BCMA-specific iPSC-T cells possessed a single T-cell receptor clonotype with cognate BCMA peptide recognition and specificity for targeting MM. RNA sequencing analyses revealed distinct genome-wide shifts and a distinctive transcriptional profile in selected iPSC clones, which can develop CD8αß+ memory T cells. This includes a repertoire of gene regulators promoting T-cell lineage development, memory CTL activation, and immune response regulation (LCK, IL7R, 4-1BB, TRAIL, GZMB, FOXF1, and ITGA1). This study highlights the potential application of iPSC technology to an adaptive T-cell therapy protocol and identifies specific transcriptional patterns that could serve as a biomarker for selection of suitable iPSC clones for the successful development of antigen-specific CD8αß+ memory T cells to improve the outcome in patients with MM.

Mechanisms of response and resistance to combined decitabine and ipilimumab for advanced myeloid disease.

The challenge of eradicating leukemia in patients with acute myelogenous leukemia (AML) after initial cytoreduction has motivated modern efforts to combine synergistic active modalities including immunotherapy. Recently, the ETCTN/CTEP 10026 study tested the combination of the DNA methyltransferase inhibitor decitabine together with the immune checkpoint inhibitor ipilimumab for AML/myelodysplastic syndrome (MDS) either after allogeneic hematopoietic stem cell transplantation (HSCT) or in the HSCT-naïve setting. Integrative transcriptome-based analysis of 304 961 individual marrow-infiltrating cells for 18 of 48 subjects treated on study revealed the strong association of response with a high baseline ratio of T to AML cells. Clinical responses were predominantly driven by decitabine-induced cytoreduction. Evidence of immune activation was only apparent after ipilimumab exposure, which altered CD4+ T-cell gene expression, in line with ongoing T-cell differentiation and increased frequency of marrow-infiltrating regulatory T cells. For post-HSCT samples, relapse could be attributed to insufficient clearing of malignant clones in progenitor cell populations. In contrast to AML/MDS bone marrow, the transcriptomes of leukemia cutis samples from patients with durable remission after ipilimumab monotherapy showed evidence of increased infiltration with antigen-experienced resident memory T cells and higher expression of CTLA-4 and FOXP3. Altogether, activity of combined decitabine and ipilimumab is impacted by cellular expression states within the microenvironmental niche of leukemic cells. The inadequate elimination of leukemic progenitors mandates urgent development of novel approaches for targeting these cell populations to generate long-lasting responses. This trial was registered at www.clinicaltrials.gov as #NCT02890329.

Neoantigen vaccine generates intratumoral T cell responses in phase Ib glioblastoma trial.

Neoantigens, which are derived from tumour-specific protein-coding mutations, are exempt from central tolerance, can generate robust immune responses1,2 and can function as bona fide antigens that facilitate tumour rejection3. Here we demonstrate that a strategy that uses multi-epitope, personalized neoantigen vaccination, which has previously been tested in patients with high-risk melanoma4-6, is feasible for tumours such as glioblastoma, which typically have a relatively low mutation load1,7 and an immunologically 'cold' tumour microenvironment8. We used personalized neoantigen-targeting vaccines to immunize patients newly diagnosed with glioblastoma following surgical resection and conventional radiotherapy in a phase I/Ib study. Patients who did not receive dexamethasone-a highly potent corticosteroid that is frequently prescribed to treat cerebral oedema in patients with glioblastoma-generated circulating polyfunctional neoantigen-specific CD4+ and CD8+ T cell responses that were enriched in a memory phenotype and showed an increase in the number of tumour-infiltrating T cells. Using single-cell T cell receptor analysis, we provide evidence that neoantigen-specific T cells from the peripheral blood can migrate into an intracranial glioblastoma tumour. Neoantigen-targeting vaccines thus have the potential to favourably alter the immune milieu of glioblastoma.

Identification of PPAR-related differentially expressed genes liver hepatocellular carcinoma and construction of a prognostic model based on data analysis and molecular docking.

Liver hepatocellular carcinoma (LIHC) is a significant global health issue with limited treatment options. In this study, single-cell RNA sequencing (scRNA-seq) data were used to explore the molecular mechanisms of LIHC development and identify potential targets for therapy. The expression of peroxisome proliferator-activated receptors (PPAR)-related genes was analysed in LIHC samples, and primary cell populations, including natural killer cells, T cells, B cells, myeloid cells, endothelial cells, fibroblasts and hepatocytes, were identified. Analysis of the differentially expressed genes (DEGs) between normal and tumour tissues revealed significant changes in gene expression in various cell populations. PPAR activity was evaluated using the 'AUCell' R software, which indicated higher scores in the normal versus the malignant hepatocytes. Furthermore, the DEGs showed significant enrichment of pathways related to lipid and glucose metabolism, cell development, differentiation and inflammation. A prognostic model was then constructed using 8 PPARs-related genes, including FABP5, LPL, ACAA1, PPARD, FABP4, PLIN1, HMGCS2 and CYP7A1, identified using least absolute shrinkage and selection operator-Cox regression analysis, and validated in the TCGA-LIHC, ICGI-LIRI and GSE14520 datasets. Patients with low-risk scores had better prognosis in all cohorts. Based on the expression of the eight model genes, two clusters of patients were identified by ConsensusCluster analysis. We also predicted small-molecule drugs targeting the model genes, and identified perfluorohexanesulfonic acid, triflumizole and perfluorononanoic acid as potential candidates. Finally, wound healing assay confirmed that PPARD can promote the migration of liver cancer cells. Overall, our study offers novel perspectives on the molecular mechanisms of LIHC and potential areas for therapeutic intervention, which may facilitate the development of more effective treatment regimens.

Is Phenotypic Evolution Affected by Spiders' Construction Behaviors?

What happens when organisms actively modify their environment? Clarifying the role of construction behavior on a macroevolutionary scale is crucial to interpreting phenotypic evolution. Spiders, an extremely successful group of animals exhibiting a wide variety of functional, morphological, and behavioral diversity, are ideal candidates to test whether animal behaviors are related to their phenotypic evolution. Herein, we reconstructed the phylogenomic relationships of 303 spiders covering 105 families with 99 newly developed molecular markers that universally apply across Araneae, and explicitly tested the potential link between construction behavior and somatic evolution based on extensive morphological data from 3,342 extant species and 682 fossil species. Our dated molecular phylogeny provides the backbone for analyses, revealing the behavioral and ecological processes behind these spiders' morphological adaptations. Evolutionary model analyses showed the artifacts constructed by spiders, especially the suspending webs, expand the spider's ability to inhabit different habitats. These artifacts have more substantial impacts on their somatic traits than habitats and promote different trajectories of morphological evolution. Specifically, for spiders, silk-lined burrowing produced larger bodies, relatively shorter legs, and longer patellae, while web-building produced smaller bodies, relatively longer legs, and shorter patellae, and hunting promoted an intermediate morphological size. Molecular convergence analysis showed that genes related to morphogenesis or response to stimulus and stress are enriched in spiders with construction behavior. Our study demonstrated that the construction behavior of an animal plays a crucial role in determining the direction and intensity of the selection pressure imposed on it and provides strong evidence that construction behaviors are associated with phenotypic evolution on macroevolutionary timescales. [Araneae; body size; habitat change; molecular marker; leg length; phylogenomics.].

T cell activation contributes to purifying selection against the MELAS-associated m.3243A>G pathogenic variant in blood.

T cells have been shown to maintain a lower percentage (heteroplasmy) of the pathogenic m.3243A>G variant (MT-TL1, associated with maternally inherited diabetes and deafness [MIDD] and mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes [MELAS]). The mechanism(s) underlying this purifying selection, however, remain unknown. Here we report that purified patient memory CD4+ T cells have lower bulk m.3243A>G heteroplasmy compared to naïve CD4+ T cells. In vitro activation of naïve CD4+ m.3243A>G patient T cells results in lower bulk m.3243A>G heteroplasmy after proliferation. Finally, m.3243A>G patient T cell receptor repertoire sequencing reveals relative oligoclonality compared to controls. These data support a role for T cell activation in peripheral, purifying selection against high m.3243A>G heteroplasmy T cells at the level of the cell, in a likely cell-autonomous fashion.

Corrigendum: An immunogenic personal neoantigen vaccine for patients with melanoma.

This corrects the article DOI: 10.1038/nature22991.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands.

BACKGROUND: Spiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution. RESULTS: We generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands. CONCLUSIONS: This research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

Molecular and cellular features of CTLA-4 blockade for relapsed myeloid malignancies after transplantation.

Relapsed myeloid disease after allogeneic stem cell transplantation (HSCT) remains largely incurable. We previously demonstrated the potent activity of immune checkpoint blockade in this clinical setting with ipilimumab or nivolumab. To define the molecular and cellular pathways by which CTLA-4 blockade with ipilimumab can reinvigorate an effective graft-versus-leukemia (GVL) response, we integrated transcriptomic analysis of leukemic biopsies with immunophenotypic profiling of matched peripheral blood samples collected from patients treated with ipilimumab following HSCT on the Experimental Therapeutics Clinical Trials Network 9204 trial. Response to ipilimumab was associated with transcriptomic evidence of increased local CD8+ T-cell infiltration and activation. Systemically, ipilimumab decreased naïve and increased memory T-cell populations and increased expression of markers of T-cell activation and costimulation such as PD-1, HLA-DR, and ICOS, irrespective of response. However, responding patients were characterized by higher turnover of T-cell receptor sequences in peripheral blood and showed increased expression of proinflammatory chemokines in plasma that was further amplified by ipilimumab. Altogether, these data highlight the compositional T-cell shifts and inflammatory pathways induced by ipilimumab both locally and systemically that associate with successful GVL outcomes. This trial was registered at www.clinicaltrials.gov as #NCT01822509.

An immunogenic personal neoantigen vaccine for patients with melanoma.

Effective anti-tumour immunity in humans has been associated with the presence of T cells directed at cancer neoantigens, a class of HLA-bound peptides that arise from tumour-specific mutations. They are highly immunogenic because they are not present in normal tissues and hence bypass central thymic tolerance. Although neoantigens were long-envisioned as optimal targets for an anti-tumour immune response, their systematic discovery and evaluation only became feasible with the recent availability of massively parallel sequencing for detection of all coding mutations within tumours, and of machine learning approaches to reliably predict those mutated peptides with high-affinity binding of autologous human leukocyte antigen (HLA) molecules. We hypothesized that vaccination with neoantigens can both expand pre-existing neoantigen-specific T-cell populations and induce a broader repertoire of new T-cell specificities in cancer patients, tipping the intra-tumoural balance in favour of enhanced tumour control. Here we demonstrate the feasibility, safety, and immunogenicity of a vaccine that targets up to 20 predicted personal tumour neoantigens. Vaccine-induced polyfunctional CD4+ and CD8+ T cells targeted 58 (60%) and 15 (16%) of the 97 unique neoantigens used across patients, respectively. These T cells discriminated mutated from wild-type antigens, and in some cases directly recognized autologous tumour. Of six vaccinated patients, four had no recurrence at 25 months after vaccination, while two with recurrent disease were subsequently treated with anti-PD-1 (anti-programmed cell death-1) therapy and experienced complete tumour regression, with expansion of the repertoire of neoantigen-specific T cells. These data provide a strong rationale for further development of this approach, alone and in combination with checkpoint blockade or other immunotherapies.

Difference Between Screw Cement Filling and Adequate Osteotomy With Thick Liner for Primary Total Knee Arthroplasty in Patients With Rand IIb Tibial Defects.

BACKGROUND: The effectiveness of 2 treatment options, screw-cement fill, and adequate osteotomy with a thick liner, in treating patients with Rand IIb tibial defects (tibial plateau defects to a depth of 5 to 10 millimeters) in primary total knee arthroplasty (TKA) has not yet been demonstrated. Therefore, we performed a retrospective study to evaluate the differences between these 2 treatments. METHODS: We retrospectively analyzed patients who underwent primary TKA for Rand IIb tibial plateau defects from 2015 to 2020 from a department database. Patients were categorized into the screw-cement and thick liner groups based on the different options used to repair tibial defects. We evaluated Knee Society Score, range of motion (ROM), Insall-Salvati index (ISI), and Forgotten Joint Score (FJS) in both groups. We also compared differences in prosthesis survival, stiffness, myasthenia, and joint clicking between the 2 groups at mean 2 years postoperatively (range, 2 to 2.3). A power analysis was performed on the number of cases in the cohort. RESULTS: Postoperative femur-tibia mechanical axis (FTMA) correction was significantly higher in the screw-cement group than in the thick-liner group: 18.8 (±5.6°) versus 15.4 (±5.9°) (P < .01); At mean 2 years after surgery, the American Knee Society Functional Score improvement values were higher in the thick-liner group than in the screw-cement group: 36.3 (±12.4) versus 42.4 (±16.4) (P = .05). Postoperative ISI scores were 0.95 (±0.12) points in the screw-cement group and 0.89 (±0.13) points in the-thick liner group (P = .03). There were no statistically significant differences in the Knee Society Clinical Score, ROM, FJS, stiffness, myasthenia, joint clicking, and revision rate. CONCLUSION: The results of this study showed no significant difference in clinical outcomes between the 2 reconstruction strategies of the screw cement fill technique and the adequate osteotomy and thick liner technique for Rand IIb tibial plateau defects. However, in patients who have FTMA deformities greater than 20° or in younger patients who need to preserve bone volume, we recommend the screw cement filling technique to ensure stable postoperative results and to prepare these patients for possible later surgery.

Ancient volcanos as species pumps: A case study of freshwater amphipods in Northeast Asia.

Volcano-tectonic processes have been viewed as primary drivers in the formation of present-day diversity. Volcanos associated with mountain uplifts drive allopatric speciation through vicariance and may impact the surrounding areas like species pump or species attractor. However, the application of these hypotheses to aquatic fauna has rarely been tested explicitly. We tested these hypotheses in the Changbai Mountains (Mts), which are one of the most typical, active volcanic ranges in Northeast (NE) Asia with a long and turbulent geological history. The Gammarus nekkensis species complex of amphipod crustaceans, widely distributed throughout NE Asia with poor dispersal abilities and a long evolutionary history, is a suitable model for testing hypotheses of species pump or species attractor. Phylogenetic and ancestral range reconstructions demonstrated that the studied amphipod originated from the Changbai Mts ~27 Ma and diverged into eastern (Clade I) and western (Clade II) clades, which corresponds well with the initial volcanic eruption of the Changbai Mts in the Late Oligocene. The subsequent diversifications of subclades CI-3, CII-1a and CII-2a were probably driven by second and third eruptions of the Changbai Mts during the Miocene. In particular, the Changbai lineages had spread to the Russian Far East multiple times since the Early Miocene, and widely colonized the region during the Pleistocene. Our discoveries suggest that the ancient volcanos of the Changbai Mts act as species pumps in NE Asia, resulted in burst of diversification around the Changbai Mts and subsequent dispersals into adjacent regions.

Past climate cooling promoted global dispersal of amphipods from Tian Shan montane lakes to circumboreal lakes.

Climate changes have substantial impacts on the geographic distribution of montane lakes and evolutionary dynamics of cold-adapted species. Past climate cooling is hypothesized to have promoted the dispersal of cold-adapted species via montane lakes, while future climate warming is thought to constrain their distributions. We test this hypothesis by using phylogeographic analysis and niche modeling of the Holarctic crustacean Gammarus lacustris with global sampling comprised of 567 sequenced individuals and 3180 occurrence records. We found that the species arose in Tian Shan in Central Asia and dispersed into montane lakes along the Alps, Himalayas, Tibet, East Asia, and the North American Rocky Mountain ranges, with accelerated diversification rates outside Tian Shan. Climatically suitable regions for geographic lineages of G. lacustris were larger during cooling periods (LGM), but smaller during warming periods (Mid-Holocene). In the future (2070) scenario, potential distributions in the Himalayas, North Tibet, South Tibet and North America are predicted to expand, whereas ranges in East Asia, Europe and Tian Shan will decline. Our results suggest that Mid-Miocene-to-Pleistocene continuous cooling promoted multiple independent dispersal events out of Tian Shan due to increased availability of montane lakes via "budding" of lineages. Montane lakes are conduits through which cold-adapted amphipods globally dispersed, dominating circumboreal lakes. However, future climate warming is likely to force organisms to shift upward in altitude and northward in latitude, leading to a future change in local populations. These findings highlight the importance of conservation of montane lakes, especially in the context of climate change.