Fundamental immune-oncogenicity trade-offs define driver mutation fitness

Missense driver mutations in cancer are concentrated in a few hotspots1. Various mechanisms have been proposed to explain this skew, including biased mutational processes2, phenotypic differences3-6 and immunoediting of neoantigens7,8; however, to our knowledge, no existing model weighs the relative contribution of these features to tumour evolution. We propose a unified theoretical 'free fitness' framework that parsimoniously integrates multimodal genomic, epigenetic, transcriptomic and proteomic data into a biophysical model of the rate-limiting processes underlying the fitness advantage conferred on cancer cells by driver gene mutations. Focusing on TP53, the most mutated gene in cancer1, we present an inference of mutant p53 concentration and demonstrate that TP53 hotspot mutations optimally solve an evolutionary trade-off between oncogenic potential and neoantigen immunogenicity. Our model anticipates patient survival in The Cancer Genome Atlas and patients with lung cancer treated with immunotherapy as well as the age of tumour onset in germline carriers of TP53 variants. The predicted differential immunogenicity between hotspot mutations was validated experimentally in patients with cancer and in a unique large dataset of healthy individuals. Our data indicate that immune selective pressure on TP53 mutations has a smaller role in non-cancerous lesions than in tumours, suggesting that targeted immunotherapy may offer an early prophylactic opportunity for the former. Determining the relative contribution of immunogenicity and oncogenic function to the selective advantage of hotspot mutations thus has important implications for both precision immunotherapies and our understanding of tumour evolution.

126. Magnitude and Dynamics of the T-Cell Response to SARS-CoV-2 Infection and Vaccination

Background

T cells are central to the early identification and clearance of viral infections and support antibody generation by B cells, making them desirable for assessing the immune response to SARS-CoV-2 infection and vaccines. We combined 2 high-throughput immune profiling methods to create a quantitative picture of the SARS-CoV-2 T-cell response that is highly sensitive, durable, diagnostic, and discriminatory between natural infection and vaccination.

Methods

We deeply characterized 116 convalescent COVID-19 subjects by experimentally mapping CD8 and CD4 T-cell responses via antigen stimulation to 545 Human Leukocyte Antigen (HLA) class I and 284 class II viral peptides. We also performed T-cell receptor (TCR) repertoire sequencing on 1815 samples from 1521 PCR-confirmed SARS-CoV-2 cases and 3500 controls to identify shared public TCRs from SARS-CoV-2-associated CD8 and CD4 T cells. Combining these approaches with additional samples from vaccinated individuals, we characterized the response to natural infection as well as vaccination by separating responses to spike protein from other viral targets.

Results

We find that T-cell responses are often driven by a few immunodominant, HLA-restricted epitopes. As expected, the SARS-CoV-2 T-cell response peaks about 1-2 weeks after infection and is detectable at least several months after recovery. Applying these data, we trained a classifier to diagnose past SARS-CoV-2 infection based solely on TCR sequencing from blood samples and observed, at 99.8% specificity, high sensitivity soon after diagnosis (Day 3–7 = 85.1%; Day 8–14 = 94.8%) that persists after recovery (Day 29+/convalescent = 95.4%). Finally, by evaluating TCRs binding epitopes targeting all non-spike SARS-CoV-2 proteins, we were able to separate natural infection from vaccination with > 99% specificity.

Conclusion

TCR repertoire sequencing from whole blood reliably measures the adaptive immune response to SARS-CoV-2 soon after viral antigenic exposure (before antibodies are typically detectable) as well as at later time points, and distinguishes post-infection vs. vaccine immune responses with high specificity. This approach to characterizing the cellular immune response has applications in clinical diagnostics as well as vaccine development and monitoring.

Disclosures

Thomas M. Snyder, PhD, Adaptive Biotechnologies (Employee, Shareholder) Rachel M. Gittelman, PhD, Adaptive Biotechnologies (Employee, Shareholder) Mark Klinger, PhD, Adaptive Biotechnologies (Employee, Shareholder) Damon H. May, PhD, Adaptive Biotechnologies (Employee, Shareholder) Edward J. Osborne, PhD, Adaptive Biotechnologies (Employee, Shareholder) Ruth Taniguchi, PhD, Adaptive Biotechnologies (Employee, Shareholder) H. Jabran Zahid, PhD, Microsoft Research (Employee, Shareholder) Rebecca Elyanow, PhD, Adaptive Biotechnologies (Employee, Shareholder) Sudeb C. Dalai, MD, PhD, Adaptive Biotechnologies (Employee, Shareholder) Ian M. Kaplan, PhD, Adaptive Biotechnologies (Employee, Shareholder) Jennifer N. Dines, MD, Adaptive Biotechnologies (Employee, Shareholder) Matthew T. Noakes, PhD, Adaptive Biotechnologies (Employee, Shareholder) Ravi Pandya, PhD, Microsoft Research (Employee, Shareholder) Lance Baldo, MD, Adaptive Biotechnologies (Employee, Shareholder, Leadership Interest) James R. Heath, PhD, Merck (Research Grant or Support, Funding (from BARDA) for the ISB INCOV project, but had no role in planning the research or in writing the paper.) Joaquin Martinez-Lopez, MD, PhD, Adaptive Biotechnologies (Consultant) Jonathan M. Carlson, PhD, Microsoft Research (Employee, Shareholder) Harlan S. Robins, PhD, Adaptive Biotechnologies (Board Member, Employee, Shareholder)

Magnitude and Dynamics of the T-Cell Response to SARS-CoV-2 Infection at Both Individual and Population Levels

T cells are involved in the early identification and clearance of viral infections and also support the development of antibodies by B cells. This central role for T cells makes them a desirable target for assessing the immune response to SARS-CoV-2 infection. Here, we combined two high-throughput immune profiling methods to create a quantitative picture of the T-cell response to SARS-CoV-2. First, at the individual level, we deeply characterized 3 acutely infected and 58 recovered COVID-19 subjects by experimentally mapping their CD8 T-cell response through antigen stimulation to 545 Human Leukocyte Antigen (HLA) class I presented viral peptides (class II data in a forthcoming study). Then, at the population level, we performed T-cell repertoire sequencing on 1,815 samples (from 1,521 COVID-19 subjects) as well as 3,500 controls to identify shared "public" T-cell receptors (TCRs) associated with SARS-CoV-2 infection from both CD8 and CD4 T cells. Collectively, our data reveal that CD8 T-cell responses are often driven by a few immunodominant, HLA-restricted epitopes. As expected, the T-cell response to SARS-CoV-2 peaks about one to two weeks after infection and is detectable for at least several months after recovery. As an application of these data, we trained a classifier to diagnose SARS-CoV-2 infection based solely on TCR sequencing from blood samples, and observed, at 99.8% specificity, high early sensitivity soon after diagnosis (Day 3-7 = 85.1% [95% CI = 79.9-89.7]; Day 8-14 = 94.8% [90.7-98.4]) as well as lasting sensitivity after recovery (Day 29+/convalescent = 95.4% [92.1-98.3]). These results demonstrate an approach to reliably assess the adaptive immune response both soon after viral antigenic exposure (before antibodies are typically detectable) as well as at later time points. This blood-based molecular approach to characterizing the cellular immune response has applications in clinical diagnostics as well as in vaccine development and monitoring.

A copy number variant is associated with a spectrum of pigmentation patterns in the rock pigeon (Columba livia)

Rock pigeons (Columba livia) display an extraordinary array of pigment pattern variation. One such pattern, Almond, is characterized by a variegated patchwork of plumage colors that are distributed in an apparently random manner. Almond is a sex-linked, semi-dominant trait controlled by the classical Stipper (St) locus. Heterozygous males (Z^StZ⁺ sex chromosomes) and hemizygous Almond females (Z^StW) are favored by breeders for their attractive plumage. In contrast, homozygous Almond males (Z^StZ^St) develop severe eye defects and often lack plumage pigmentation, suggesting that higher dosage of the mutant allele is deleterious. To determine the molecular basis of Almond, we compared the genomes of Almond pigeons to non-Almond pigeons and identified a candidate St locus on the Z chromosome. We found a copy number variant (CNV) within the differentiated region that captures complete or partial coding sequences of four genes, including the melanosome maturation gene Mlana. We did not find fixed coding changes in genes within the CNV, but all genes are misexpressed in regenerating feather bud collar cells of Almond birds. Notably, six other alleles at the St locus are associated with depigmentation phenotypes, and all exhibit expansion of the same CNV. Structural variation at St is linked to diversity in plumage pigmentation and gene expression, and thus provides a potential mode of rapid phenotypic evolution in pigeons.

The genetic changes responsible for different animal color patterns are poorly understood, due in part to a paucity of research organisms that are both genetically tractable and phenotypically diverse. Domestic pigeons (Columba livia) have been artificially selected for many traits, including an enormous variety of color patterns that are variable both within and among different breeds of this single species. We investigated the genetic basis of a sex-linked color pattern in pigeons called Almond that is characterized by a sprinkled pattern of plumage pigmentation. Pigeons with one copy of the Almond allele have desirable color pattern; however, male pigeons with two copies of the Almond mutation have severely depleted pigmentation and congenital eye defects. By comparing the genomes of Almond and non-Almond pigeons, we discovered that Almond pigeons have extra copies of a chromosome region that contains a gene that is critical for the formation of pigment granules. We also found that different numbers of copies of this region are associated with varying degrees of pigment reduction. The Almond phenotype in pigeons bears a remarkable resemblance to Merle coat color mutants in dogs, and our new results from pigeons suggest that similar genetic mechanisms underlie these traits in both species. Our work highlights the role of gene copy number variation as a potential driver of rapid phenotypic evolution.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Plant organ evolution revealed by phylotranscriptomics in Arabidopsis thaliana

The evolution of phenotypes occurs through changes both in protein sequence and gene expression levels. Though much of plant morphological evolution can be explained by changes in gene expression, examining its evolution has challenges. To gain a new perspective on organ evolution in plants, we applied a phylotranscriptomics approach. We combined a phylostratigraphic approach with gene expression based on the strand-specific RNA-seq data from seedling, floral bud, and root of 19 Arabidopsis thaliana accessions to examine the age and sequence divergence of transcriptomes from these organs and how they adapted over time. Our results indicate that, among the sense and antisense transcriptomes of these organs, the sense transcriptomes of seedlings are the evolutionarily oldest across all accessions and are the most conserved in amino acid sequence for most accessions. In contrast, among the sense transcriptomes from these same organs, those from floral bud are evolutionarily youngest and least conserved in sequence for most accessions. Different organs have adaptive peaks at different stages in their evolutionary history; however, all three show a common adaptive signal from the Magnoliophyta to Brassicale stage. Our research highlights how phylotranscriptomic analyses can be used to trace organ evolution in the deep history of plant species.

Wham: Identifying Structural Variants of Biological Consequence

Existing methods for identifying structural variants (SVs) from short read datasets are inaccurate. This complicates disease-gene identification and efforts to understand the consequences of genetic variation. In response, we have created Wham (Whole-genome Alignment Metrics) to provide a single, integrated framework for both structural variant calling and association testing, thereby bypassing many of the difficulties that currently frustrate attempts to employ SVs in association testing. Here we describe Wham, benchmark it against three other widely used SV identification tools–Lumpy, Delly and SoftSearch–and demonstrate Wham’s ability to identify and associate SVs with phenotypes using data from humans, domestic pigeons, and vaccinia virus. Wham and all associated software are covered under the MIT License and can be freely downloaded from github (https://github.com/zeeev/wham), with documentation on a wiki (http://zeeev.github.io/wham/). For community support please post questions to https://www.biostars.org/.

DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation

Epigenome modulation potentially provides a mechanism for organisms to adapt, within and between generations. However, neither the extent to which this occurs, nor the mechanisms involved are known. Here we investigate DNA methylation variation in Swedish Arabidopsis thaliana accessions grown at two different temperatures. Environmental effects were limited to transposons, where CHH methylation was found to increase with temperature. Genome-wide association studies (GWAS) revealed that the extensive CHH methylation variation was strongly associated with genetic variants in both cis and trans, including a major trans-association close to the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body methylation (GBM) was not affected by growth temperature, but was instead correlated with the latitude of origin. Accessions from colder regions had higher levels of GBM for a significant fraction of the genome, and this was associated with increased transcription for the genes affected. GWAS revealed that this effect was largely due to trans-acting loci, many of which showed evidence of local adaptation.

DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation

Epigenome modulation potentially provides a mechanism for organisms to adapt, within and between generations. However, neither the extent to which this occurs, nor the mechanisms involved are known. Here we investigate DNA methylation variation in Swedish Arabidopsis thaliana accessions grown at two different temperatures. Environmental effects were limited to transposons, where CHH methylation was found to increase with temperature. Genome-wide association studies (GWAS) revealed that the extensive CHH methylation variation was strongly associated with genetic variants in both cis and trans, including a major trans-association close to the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body methylation (GBM) was not affected by growth temperature, but was instead correlated with the latitude of origin. Accessions from colder regions had higher levels of GBM for a significant fraction of the genome, and this was associated with increased transcription for the genes affected. GWAS revealed that this effect was largely due to trans-acting loci, many of which showed evidence of local adaptation.

DOI:http://dx.doi.org/10.7554/eLife.05255.001

Organisms need to adapt quickly to changes in their environment. Mutations in the DNA sequence of genes can lead to new adaptations, but this can take many generations. Instead, altering how genes are switched on by changing how the DNA is packaged in cells can allow organisms to adapt within and between generations. One way that genes are controlled in organisms is by a process known as DNA methylation, where ‘methyl’ tags are added to DNA and act as markers for other proteins involved in activating genes.

DNA is made of four different molecules called ‘nucleotides’ that are arranged in different orders to produce a vast variety of DNA sequences. One type of DNA methylation can happen at sites where a nucleotide called cytosine is followed by two other non-cytosine nucleotides. Another type of methylation can take place at sites where a cytosine is followed by a guanine nucleotide. However, it is not clear how big a role DNA methylation plays in allowing organisms to adapt to their changing environment.

Here, Dubin, Zhang, Meng, Remigereau et al. studied DNA methylation in a plant called Arabidopsis thaliana. Several different varieties of A. thaliana plants from Sweden were grown at two different temperatures. The experiments showed that the A. thaliana plants grown at higher temperatures were more likely to have methyl tags attached to sections of DNA called transposons, which are able to move around the genome. There was a lot of variety in the levels of this DNA methylation in the different plants, and some of it was shown to be associated with variation in a gene that is involved in DNA methylation.

However, not all of the DNA methylation in these plants was sensitive to the temperature the plants were grown in. Dubin, Zhang, Meng, Remigereau et al. show that the pattern of a type of DNA methylation that is found within genes depends on how far north in Sweden the plants' ancestors came from rather than the temperature the plants were grown in. Plants that originated from colder regions, farther north, had more DNA methylation within many genes and these genes were more active.

These findings suggest that genetic differences in these plants strongly influence the levels of DNA methylation, and they provide the first direct link between DNA methylation and adaption to the environment. Future studies should reveal how DNA methylation is regulated in these plants, and whether it plays a key role in adaptation, or merely reflects other changes in the genome.

DOI:http://dx.doi.org/10.7554/eLife.05255.002

DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation

Epigenome modulation potentially provides a mechanism for organisms to adapt, within and between generations. However, neither the extent to which this occurs, nor the mechanisms involved are known. Here we investigate DNA methylation variation in Swedish Arabidopsis thaliana accessions grown at two different temperatures. Environmental effects were limited to transposons, where CHH methylation was found to increase with temperature. Genome-wide association studies (GWAS) revealed that the extensive CHH methylation variation was strongly associated with genetic variants in both cis and trans, including a major trans-association close to the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body methylation (GBM) was not affected by growth temperature, but was instead correlated with the latitude of origin. Accessions from colder regions had higher levels of GBM for a significant fraction of the genome, and this was associated with increased transcription for the genes affected. GWAS revealed that this effect was largely due to trans-acting loci, many of which showed evidence of local adaptation.

High resolution genetic mapping uncovers chitin synthase-1 as the target-site of the structurally diverse mite growth inhibitors clofentezine, hexythiazox and etoxazole in Tetranychus urticae

The acaricides clofentezine, hexythiazox and etoxazole are commonly referred to as 'mite growth inhibitors', and clofentezine and hexythiazox have been used successfully for the integrated control of plant mite pests for decades. Although they are still important today, their mode of action has remained elusive. Recently, a mutation in chitin synthase 1 (CHS1) was linked to etoxazole resistance. In this study, we identified and investigated a Tetranychus urticae strain (HexR) harboring recessive, monogenic resistance to each of hexythiazox, clofentezine, and etoxazole. To elucidate if there is a common genetic basis for the observed cross-resistance, we adapted a previously developed bulk segregant analysis method to map with high resolution a single, shared resistance locus for all three compounds. This finding indicates that the underlying molecular basis for resistance to all three compounds is identical. This locus is centered on the CHS1 gene, and as supported by additional genetic and biochemical studies, a non-synonymous variant (I1017F) in CHS1 associates with resistance to each of the tested acaricides in HexR. Our findings thus demonstrate a shared molecular mode of action for the chemically diverse mite growth inhibitors clofentezine, hexythiazox and etoxazole as inhibitors of an essential, non-catalytic activity of CHS1. Given the previously documented cross-resistance between clofentezine, hexythiazox and the benzyolphenylurea (BPU) compounds flufenoxuron and cycloxuron, CHS1 should be also considered as a potential target-site of insecticidal BPUs.

Reciprocal responses in the interaction between Arabidopsis and the cell-content-feeding chelicerate herbivore spider mite

Most molecular-genetic studies of plant defense responses to arthropod herbivores have focused on insects. However, plant-feeding mites are also pests of diverse plants, and mites induce different patterns of damage to plant tissues than do well-studied insects (e.g. lepidopteran larvae or aphids). The two-spotted spider mite (Tetranychus urticae) is among the most significant mite pests in agriculture, feeding on a staggering number of plant hosts. To understand the interactions between spider mite and a plant at the molecular level, we examined reciprocal genome-wide responses of mites and its host Arabidopsis (Arabidopsis thaliana). Despite differences in feeding guilds, we found that transcriptional responses of Arabidopsis to mite herbivory resembled those observed for lepidopteran herbivores. Mutant analysis of induced plant defense pathways showed functionally that only a subset of induced programs, including jasmonic acid signaling and biosynthesis of indole glucosinolates, are central to Arabidopsis's defense to mite herbivory. On the herbivore side, indole glucosinolates dramatically increased mite mortality and development times. We identified an indole glucosinolate dose-dependent increase in the number of differentially expressed mite genes belonging to pathways associated with detoxification of xenobiotics. This demonstrates that spider mite is sensitive to Arabidopsis defenses that have also been associated with the deterrence of insect herbivores that are very distantly related to chelicerates. Our findings provide molecular insights into the nature of, and response to, herbivory for a representative of a major class of arthropod herbivores.

Reciprocal responses in the interaction between Arabidopsis and the cell-content-feeding chelicerate herbivore spider mite

Most molecular-genetic studies of plant defense responses to arthropod herbivores have focused on insects. However, plant-feeding mites are also pests of diverse plants, and mites induce different patterns of damage to plant tissues than do well-studied insects (e.g. lepidopteran larvae or aphids). The two-spotted spider mite (Tetranychus urticae) is among the most significant mite pests in agriculture, feeding on a staggering number of plant hosts. To understand the interactions between spider mite and a plant at the molecular level, we examined reciprocal genome-wide responses of mites and its host Arabidopsis (Arabidopsis thaliana). Despite differences in feeding guilds, we found that transcriptional responses of Arabidopsis to mite herbivory resembled those observed for lepidopteran herbivores. Mutant analysis of induced plant defense pathways showed functionally that only a subset of induced programs, including jasmonic acid signaling and biosynthesis of indole glucosinolates, are central to Arabidopsis's defense to mite herbivory. On the herbivore side, indole glucosinolates dramatically increased mite mortality and development times. We identified an indole glucosinolate dose-dependent increase in the number of differentially expressed mite genes belonging to pathways associated with detoxification of xenobiotics. This demonstrates that spider mite is sensitive to Arabidopsis defenses that have also been associated with the deterrence of insect herbivores that are very distantly related to chelicerates. Our findings provide molecular insights into the nature of, and response to, herbivory for a representative of a major class of arthropod herbivores.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana

Genetic differences between Arabidopsis thaliana accessions underlie the plant's extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

Heat shock response of killifish (Fundulus heteroclitus): candidate gene and heterologous microarray approaches

Northern and southern subspecies of the Atlantic killifish, Fundulus heteroclitus, differ in maximal thermal tolerance. To determine whether these subspecies also differ in their heat shock response (HSR), we exposed 20°C acclimated killifish to a 2 h heat shock at 34°C and examined gene expression in fish from both subspecies during heat shock and recovery using real-time quantitative PCR and a heterologous cDNA microarray designed for salmonid fishes. The heat shock proteins Hsp70-1, hsp27, and hsp30 were upregulated to a greater extent in the high temperature-tolerant southern subspecies than in the less tolerant northern subspecies, whereas hsp70-2 (which showed the largest upregulation of all the heat shock proteins) in both gill and muscle and hsp90α in muscle was upregulated to a greater extent in northern than in southern fish. These data demonstrate that differences in the HSR between subspecies cannot be due to changes in a single global regulator but must occur via gene-specific mechanisms. They also suggest that the role, if any, of hsps in establishing thermal tolerance is complex and varies from gene to gene. Heterologous microarray hybridization provided interpretable gene expression signatures, detecting differential regulation of genes known to be involved in the heat shock response in other species. Under control conditions, a variety of genes were differentially expressed in muscle between subspecies that suggest differences in muscle fiber type and could relate to previously observed differences between subspecies in the thermal sensitivity of swimming performance and metabolism.

Hypersensitivity of human testis-tumour cell lines to chemotherapeutic drugs

Metastatic testis tumours, in contrast to most other types of cancer, can be cured by drugs. To investigate which classes of chemotherapeutic drug are differentially toxic to testis-tumour cells, we compared the in vitro dose-response curves of 5 human testis and 5 bladder-cancer cell lines to 12 compounds. The testis cells were hypersensitive to drugs that interact directly with DNA (m-amsa, bleomycin, cisplatin, doxorubicin, methylnitrosourea, mitozolomide, etoposide, mitomycin-C), but little or no difference between the 2 cell types was seen following exposure to drugs whose mechanisms of action do not involve direct interaction with DNA (methotrexate, 5-fluorouracil, colchicine, vinblastine). We conclude that testis tumour cells are either less tolerant of, or have a reduced capacity to repair, DNA damage.

The effect of exposure to elevated temperatures on membrane permeability to adriamycin in Chinese hamster ovary cells in vitro

The effect of environmental temperature on plasma membrane permeability to Adriamycin was studied in Chinese hamster ovary (CHO) cells in vitro using flow cytometry. Initial rates of uptake increased steadily between 23 degrees C and 47 degrees C and there was a greater than 2-fold increase in permeability between 37 degrees C and 47 degrees C. The increase in permeability with increasing temperature was greater than expected based on the model of passive drug influx. Adriamycin uptake was also measured at 37 degrees C following previous exposure of the cells to elevated temperatures. Twenty-minute preexposures to temperatures above 41 degrees C caused a significant decrease in membrane permeability, which fell to approximately 60% of control levels after exposure to 45 degrees C. Longer periods of pre-exposure to temperatures as low as 40 degrees C were also shown to decrease membrane permeability to Adriamycin subsequently measured at 37 degrees C. The state of decreased permeability to Adriamycin induced by hyperthermia persisted for at least 2 h. The thermally induced decrease in membrane permeability to Adriamycin is of potential importance to the design of optimal schedules for thermochemotherapy.