Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia

IMPORTANCE

Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Oligotrophic waters of the Northwest Atlantic support taxonomically diverse diatom communities that are distinct from coastal waters

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems.

Microbial metabolites in the marine carbon cycle

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Aureococcus anophagefferens (Pelagophyceae) genomes improve evaluation of nutrient acquisition strategies involved in brown tide dynamics

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment.

Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre

The physical and biological dynamics that influence phytoplankton communities in the oligotrophic ocean are complex, changing across broad temporal and spatial scales. Eukaryotic phytoplankton (e.g., diatoms), despite their relatively low abundance in oligotrophic waters, are responsible for a large component of the organic matter flux to the ocean interior. Mesoscale eddies can impact both microbial community structure and function, enhancing primary production and carbon export, but the mechanisms that underpin these dynamics are still poorly understood. Here, mesoscale eddy influences on the taxonomic diversity and expressed functional profiles of surface communities of microeukaryotes and particle-associated heterotrophic bacteria from the North Pacific Subtropical Gyre were assessed over 2 years (spring 2016 and summer 2017). The taxonomic diversity of the microeukaryotes significantly differed by eddy polarity (cyclonic versus anticyclonic) and between sampling seasons/years and was significantly correlated with the taxonomic diversity of particle-associated heterotrophic bacteria. The expressed functional profile of these taxonomically distinct microeukaryotes varied consistently as a function of eddy polarity, with cyclones having a different expression pattern than anticyclones, and between sampling seasons/years. These data suggest that mesoscale forcing, and associated changes in biogeochemistry, could drive specific physiological responses in the resident microeukaryote community, independent of species composition.

Metatranscriptome Library Preparation Influences Analyses of Viral Community Activity During a Brown Tide Bloom

There is growing interest in the use of metatranscriptomics to study virus community dynamics. We used RNA samples collected from harmful brown tides caused by the eukaryotic alga Aureococcus anophagefferens within New York (United States) estuaries and in the process observed how preprocessing of libraries by either selection for polyadenylation or reduction in ribosomal RNA (rRNA) influenced virus community analyses. As expected, more reads mapped to the A. anophagefferens genome in polyadenylation-selected libraries compared to the rRNA-reduced libraries, with reads mapped in each sample correlating to one another regardless of preprocessing of libraries. Yet, this trend was not seen for reads mapping to the Aureococcus anophagefferens Virus (AaV), where significantly more reads (approximately two orders of magnitude) were mapped to the AaV genome in the rRNA-reduced libraries. In the rRNA-reduced libraries, there was a strong and significant correlation between reads mappings to AaV and A. anophagefferens. Overall, polyadenylation-selected libraries produced fewer viral contigs, fewer reads mapped to viral contigs, and different proportions across viral realms and families, compared to their rRNA-reduced pairs. This study provides evidence that libraries generated by rRNA reduction and not selected for polyadenylation are more appropriate for quantitative characterization of viral communities in aquatic ecosystems by metatranscriptomics.

Progress and promise of omics for predicting the impacts of climate change on harmful algal blooms

Climate change is predicted to increase the severity and prevalence of harmful algal blooms (HABs). In the past twenty years, omics techniques such as genomics, transcriptomics, proteomics and metabolomics have transformed that data landscape of many fields including the study of HABs. Advances in technology have facilitated the creation of many publicly available omics datasets that are complementary and shed new light on the mechanisms of HAB formation and toxin production. Genomics have been used to reveal differences in toxicity and nutritional requirements, while transcriptomics and proteomics have been used to explore HAB species responses to environmental stressors, and metabolomics can reveal mechanisms of allelopathy and toxicity. In this review, we explore how omics data may be leveraged to improve predictions of how climate change will impact HAB dynamics. We also highlight important gaps in our knowledge of HAB prediction, which include swimming behaviors, microbial interactions and evolution that can be addressed by future studies with omics tools. Lastly, we discuss approaches to incorporate current omics datasets into predictive numerical models that may enhance HAB prediction in a changing world. With the ever-increasing omics databases, leveraging these data for understanding climate-driven HAB dynamics will be increasingly powerful.

Transcriptional and Proteomic Choreography Under Phosphorus Deficiency and Re-supply in the N2 Fixing Cyanobacterium Trichodesmium erythraeum

The N₂ fixing cyanobacterium Trichodesmium is a critically important organism in oligotrophic marine ecosystems, supplying “new” nitrogen (N) to the otherwise N-poor tropical and subtropical regions where it occurs. Low concentrations of phosphorus (P) in these regions can constrain Trichodesmium distribution and N₂ fixation rates. Physiological characterization of a single species in a mixed community can be challenging, and ‘omic approaches are increasingly important tools for tracking nutritional physiology in a taxon-specific manner. As such, studies examining the dynamics of gene and protein markers of physiology (e.g., nutrient stress) are critical for the application and interpretation of such ‘omic data in situ. Here we leveraged combined transcriptomics, proteomics, and enzyme activity assays to track the physiological response of Trichodesmium erythraeum IMS101 to P deficiency and subsequent P re-supply over 72 h of sampling. P deficiency resulted in differential gene expression, protein abundance, and enzyme activity that highlighted a synchronous shift in P physiology with increases in the transcripts and corresponding proteins for hydrolyzing organic phosphorus, taking up phosphate with higher affinity, and modulating intracellular P demand. After P deficiency was alleviated, gene expression of these biomarkers was reduced to replete levels within 4 h of P amendment. A number of these gene biomarkers were adjacent to putative pho boxes and their expression patterns were similar to a sphR response regulator. Protein products of the P deficiency biomarkers were slow to decline, with 84% of the original P deficient protein set still significantly differentially expressed after 72 h. Alkaline phosphatase activity tracked with proteins for this enzyme. With the rapid turnover time of transcripts, they appear to be good biomarkers of a P stress phenotype, whereas proteins, with a slower turnover time, may better reflect cellular activities. These results highlight the importance of validating and pairing transcriptome and proteome data that can be applied to physiological studies of key species in situ.

Transcriptional Shifts Highlight the Role of Nutrients in Harmful Brown Tide Dynamics

Harmful algal blooms (HABs) threaten ecosystems and human health worldwide. Controlling nitrogen inputs to coastal waters is a common HAB management strategy, as nutrient concentrations often suggest coastal blooms are nitrogen-limited. However, defining best nutrient management practices is a long-standing challenge: in part, because of difficulties in directly tracking the nutritional physiology of harmful species in mixed communities. Using metatranscriptome sequencing and incubation experiments, we addressed this challenge by assaying the in situ physiological ecology of the ecosystem destructive alga, Aureococcus anophagefferens. Here we show that gene markers of phosphorus deficiency were expressed in situ, and modulated by the enrichment of phosphorus, which was consistent with the observed growth rate responses. These data demonstrate the importance of phosphorus in controlling brown-tide dynamics, suggesting that phosphorus, in addition to nitrogen, should be evaluated in the management and mitigation of these blooms. Given that nutrient concentrations alone were suggestive of a nitrogen-limited ecosystem, this study underscores the value of directly assaying harmful algae in situ for the development of management strategies.

Periodic and coordinated gene expression between a diazotroph and its diatom host

In the surface ocean, light fuels photosynthetic carbon fixation of phytoplankton, playing a critical role in ecosystem processes including carbon export to the deep sea. In oligotrophic oceans, diatom-diazotroph associations (DDAs) play a keystone role in ecosystem function because diazotrophs can provide otherwise scarce biologically available nitrogen to the diatom host, fueling growth and subsequent carbon sequestration. Despite their importance, relatively little is known about the nature of these associations in situ. Here we used metatranscriptomic sequencing of surface samples from the North Pacific Subtropical Gyre (NPSG) to reconstruct patterns of gene expression for the diazotrophic symbiont Richelia and we examined how these patterns were integrated with those of the diatom host over day-night transitions. Richelia exhibited significant diel signals for genes related to photosynthesis, N₂ fixation, and resource acquisition, among other processes. N₂ fixation genes were significantly co-expressed with host nitrogen uptake and metabolism, as well as potential genes involved in carbon transport, which may underpin the exchange of nitrogen and carbon within this association. Patterns of expression suggested cell division was integrated between the host and symbiont across the diel cycle. Collectively these data suggest that symbiont-host physiological ecology is strongly interconnected in the NPSG.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

The impact of elevated CO2 on Prochlorococcus and microbial interactions with 'helper' bacterium Alteromonas

Prochlorococcus is a globally important marine cyanobacterium that lacks the gene catalase and relies on 'helper' bacteria such as Alteromonas to remove reactive oxygen species. Increasing atmospheric CO2 decreases the need for carbon concentrating mechanisms and photorespiration in phytoplankton, potentially altering their metabolism and microbial interactions even when carbon is not limiting growth. Here, Prochlorococcus (VOL4, MIT9312) was co-cultured with Alteromonas (strain EZ55) under ambient (400 p.p.m.) and elevated CO2 (800 p.p.m.). Under elevated CO2, Prochlorococcus had a significantly longer lag phase and greater apparent die-offs after transfers suggesting an increase in oxidative stress. Whole-transcriptome analysis of Prochlorococcus revealed decreased expression of the carbon fixation operon, including carboxysome subunits, corresponding with significantly fewer carboxysome structures observed by electron microscopy. Prochlorococcus co-culture responsive gene 1 had significantly increased expression in elevated CO2, potentially indicating a shift in the microbial interaction. Transcriptome analysis of Alteromonas in co-culture with Prochlorococcus revealed decreased expression of the catalase gene, known to be critical in relieving oxidative stress in Prochlorococcus by removing hydrogen peroxide. The decrease in catalase gene expression was corroborated by a significant ~6-fold decrease in removal rates of hydrogen peroxide from co-cultures. These data suggest Prochlorococcus may be more vulnerable to oxidative stress under elevated CO2 in part from a decrease in ecosystem services provided by heterotrophs like Alteromonas. This work highlights the importance of considering microbial interactions in the context of a changing ocean.

Coordinated gene expression between Trichodesmium and its microbiome over day-night cycles in the North Pacific Subtropical Gyre

Trichodesmium is a widespread, N₂ fixing marine cyanobacterium that drives inputs of newly fixed nitrogen and carbon into the oligotrophic ecosystems where it occurs. Colonies of Trichodesmium ubiquitously occur with heterotrophic bacteria that make up a diverse microbiome, and interactions within this Trichodesmium holobiont could influence the fate of fixed carbon and nitrogen. Metatranscriptome sequencing was performed on Trichodesmium colonies collected during high-frequency Lagrangian sampling in the North Pacific Subtropical Gyre (NPSG) to identify possible interactions between the Trichodesmium host and microbiome over day-night cycles. Here we show significantly coordinated patterns of gene expression between host and microbiome, many of which had significant day-night periodicity. The functions of the co-expressed genes suggested a suite of interactions within the holobiont linked to key resources including nitrogen, carbon, and iron. Evidence of microbiome reliance on Trichodesmium-derived vitamin B12 was also detected in co-expression patterns, highlighting a dependency that could shape holobiont community structure. Collectively, these patterns of expression suggest that biotic interactions could influence colony cycling of resources like nitrogen and vitamin B12, and decouple activities, like N₂ fixation, from typical abiotic drivers of Trichodesmium physiological ecology.

Diverse CO2-Induced Responses in Physiology and Gene Expression among Eukaryotic Phytoplankton

With rising atmospheric CO₂, phytoplankton face shifts in ocean chemistry including increased dissolved CO₂ and acidification that will likely influence the relative competitive fitness of different phytoplankton taxa. Here we compared the physiological and gene expression responses of six species of phytoplankton including a diatom, a raphidophyte, two haptophytes, and two dinoflagellates to ambient (~400 ppm) and elevated (~800 ppm) CO₂. Dinoflagellates had significantly slower growth rates and higher, yet variable, chlorophyll a per cell under elevated CO₂. The other phytoplankton tended to have increased growth rates and/or decreased chlorophyll a per cell. Carbon and nitrogen partitioning of cells shifted under elevated CO₂ in some species, indicating potential changes in energy fluxes due to changes in carbon concentrating mechanisms (CCM) or photorespiration. Consistent with these phenotypic changes, gene set enrichment analyses revealed shifts in energy, carbon and nitrogen metabolic pathways, though with limited overlap between species in the genes and pathways involved. Similarly, gene expression responses across species revealed few conserved CO₂-responsive genes within CCM and photorespiration categories, and a survey of available transcriptomes found high diversity in biophysical CCM and photorespiration expressed gene complements between and within the four phyla represented by these species. The few genes that displayed similar responses to CO₂ across phyla were from understudied gene families, making them targets for further research to uncover the mechanisms of phytoplankton acclimation to elevated CO₂. These results underscore that eukaryotic phytoplankton have diverse gene complements and gene expression responses to CO₂ perturbations and highlight the value of cross-phyla comparisons for identifying gene families that respond to environmental change.

Epibionts dominate metabolic functional potential of Trichodesmium colonies from the oligotrophic ocean

Trichodesmium is a genus of marine diazotrophic colonial cyanobacteria that exerts a profound influence on global biogeochemistry, by injecting 'new' nitrogen into the low nutrient systems where it occurs. Colonies of Trichodesmium ubiquitously contain a diverse assemblage of epibiotic microorganisms, constituting a microbiome on the Trichodesmium host. Metagenome sequences from Trichodesmium colonies were analyzed along a resource gradient in the western North Atlantic to examine microbiome community structure, functional diversity and metabolic contributions to the holobiont. Here we demonstrate the presence of a core Trichodesmium microbiome that is modulated to suit different ocean regions, and contributes over 10 times the metabolic potential of Trichodesmium to the holobiont. Given the ubiquitous nature of epibionts on colonies, the substantial functional diversity within the microbiome is likely an integral facet of Trichodesmium physiological ecology across the oligotrophic oceans where this biogeochemically significant diazotroph thrives.

Transcriptional response of the harmful raphidophyte Heterosigma akashiwo to nitrate and phosphate stress

The marine eukaryotic alga Heterosigma akashiwo (Raphidophyceae) is known for forming ichthyotoxic harmful algal blooms (HABs). In the past 50 years, H. akashiwo blooms have increased, occurring globally in highly eutrophic coastal and estuarine systems. These systems often incur dramatic physicochemical changes, including macronutrient (nitrogen and phosphorus) enrichment and depletion, on short timescales. Here, H. akashiwo cultures grown under nutrient replete, low N and low P growth conditions were examined for changes in biochemical and physiological characteristics in concert with transcriptome sequencing to provide a mechanistic perspective on the metabolic processes involved in responding to N and P stress. There was a marked difference in the overall transcriptional pattern between low N and low P transcriptomes. Both nutrient stresses led to significant changes in the abundance of thousands of contigs related to a wide diversity of metabolic pathways, with limited overlap between the transcriptomic responses to low N and low P. Enriched contigs under low N included many related to nitrogen metabolism, acquisition, and transport. In addition, metabolic modules like photosynthesis and carbohydrate metabolism changed significantly under low N, coincident with treatment-specific changes in photosynthetic efficiency and particulate carbohydrate content. P-specific contigs responsible for P transport and organic P use were more enriched in the low P treatment than in the replete control and low N treatment. These results provide new insight into the genetic mechanisms that distinguish how this HAB species responds to these two common nutrient stresses, and the results can inform future field studies, linking transcriptional patterns to the physiological ecology of H. akashiwo in situ.

Conserved Transcriptional Responses to Nutrient Stress in Bloom-Forming Algae

The concentration and composition of bioavailable nitrogen (N) and phosphorus (P) in the upper ocean shape eukaryotic phytoplankton communities and influence their physiological responses. Phytoplankton are known to exhibit similar physiological responses to limiting N and P conditions such as decreased growth rates, chlorosis, and increased assimilation of N and P. Are these responses similar at the molecular level across multiple species? To interrogate this question, five species from biogeochemically important, bloom-forming taxa (Bacillariophyta, Dinophyta, and Haptophyta) were grown under similar low N, low P, and replete nutrient conditions to identify transcriptional patterns and associated changes in biochemical pools related to N and P stress. Metabolic profiles, revealed through the transcriptomes of these taxa, clustered together based on species rather than nutrient stressor, suggesting that the global metabolic response to nutrient stresses was largely, but not exclusively, species-specific. Nutrient stress led to few transcriptional changes in the two dinoflagellates, consistent with other research. An orthologous group analysis examined functionally conserved (i.e., similarly changed) responses to nutrient stress and therefore focused on the diatom and haptophytes. Most conserved ortholog changes were specific to a single nutrient treatment, but a small number of orthologs were similarly changed under both N and P stress in 2 or more species. Many of these orthologs were related to photosynthesis and may represent generalized stress responses. A greater number of orthologs were conserved across more than one species under low P compared to low N. Screening the conserved orthologs for functions related to N and P metabolism revealed increased relative abundance of orthologs for nitrate, nitrite, ammonium, and amino acid transporters under N stress, and increased relative abundance of orthologs related to acquisition of inorganic and organic P substrates under P stress. Although the global transcriptional responses were dominated by species-specific changes, the analysis of conserved responses revealed functional similarities in resource acquisition pathways among different phytoplankton taxa. This overlap in nutrient stress responses observed among species may be useful for tracking the physiological ecology of phytoplankton field populations.

Virus-host relationships of marine single-celled eukaryotes resolved from metatranscriptomics

Establishing virus–host relationships has historically relied on culture-dependent approaches. Here we report on the use of marine metatranscriptomics to probe virus–host relationships. Statistical co-occurrence analyses of dsDNA, ssRNA and dsRNA viral markers of polyadenylation-selected RNA sequences from microbial communities dominated by Aureococcus anophagefferens (Quantuck Bay, NY), and diatoms (Narragansett Bay, RI) show active infections by diverse giant viruses (NCLDVs) associated with algal and nonalgal hosts. Ongoing infections of A. anophagefferens by a known Mimiviridae (AaV) occur during bloom peak and decline. Bloom decline is also accompanied by increased activity of viruses other than AaV, including (+) ssRNA viruses. In Narragansett Bay, increased temporal resolution reveals active NCLDVs with both ‘boom-and-bust’ and ‘steady-state infection’-like ecologies that include known as well as novel virus–host interactions. Our approach offers a method for screening active viral infections and develops links between viruses and their potential hosts in situ. Our observations further demonstrate that previously unknown virus–host relationships in marine systems are abundant.

Viruses are partners in ecosystem ecology, yet their study has been primarily limited to laboratory models virus-host or derived from metagenomics. Here, Moniruzzaman et al. use metatranscriptomics to resolve interactions between giant viruses and single-celled eukaryotic hosts.

The Genetic Diversity of Mesodinium and Associated Cryptophytes

Ciliates from the genus Mesodinium are globally distributed in marine and freshwater ecosystems and may possess either heterotrophic or mixotrophic nutritional modes. Members of the Mesodinium major/rubrum species complex photosynthesize by sequestering and maintaining organelles from cryptophyte prey, and under certain conditions form periodic or recurrent blooms (= red tides). Here, we present an analysis of the genetic diversity of Mesodinium and cryptophyte populations from 10 environmental samples (eight globally dispersed habitats including five Mesodinium blooms), using group-specific primers for Mesodinium partial 18S, ITS, and partial 28S rRNA genes as well as cryptophyte large subunit RuBisCO genes (rbcL). In addition, 22 new cryptophyte and four new M. rubrum cultures were used to extract DNA and sequence rbcL and 18S-ITS-28S genes, respectively, in order to provide a stronger phylogenetic context for our environmental sequences. Bloom samples were analyzed from coastal Brazil, Chile, two Northeastern locations in the United States, and the Pribilof Islands within the Bering Sea. Additionally, samples were also analyzed from the Baltic and Barents Seas and coastal California under non-bloom conditions. Most blooms were dominated by a single Mesodinium genotype, with coastal Brazil and Chile blooms composed of M. major and the Eastern USA blooms dominated by M. rubrum variant B. Sequences from all four blooms were dominated by Teleaulax amphioxeia-like cryptophytes. Non-bloom communities revealed more diverse assemblages of Mesodinium spp., including heterotrophic species and the mixotrophic Mesodinium chamaeleon. Similarly, cryptophyte diversity was also higher in non-bloom samples. Our results confirm that Mesodinium blooms may be caused by M. major, as well as multiple variants of M. rubrum, and further implicate T. amphioxeia as the key cryptophyte species linked to these phenomena in temperate and subtropical regions.

Variable depth distribution of Trichodesmium clades in the North Pacific Ocean

Populations of nitrogen-fixing cyanobacteria in the genus Trichodesmium are critical to ocean ecosystems, yet predicting patterns of Trichodesmium distribution and their role in ocean biogeochemistry is an ongoing challenge. This may, in part, be due to differences in the physiological ecology of Trichodesmium species, which are not typically considered independently in field studies. In this study, the abundance of the two dominant Trichodesmium clades (Clade I and Clade III) was investigated during a survey at Station ALOHA in the North Pacific Subtropical Gyre (NPSG) using a clade-specific qPCR approach. While Clade I dominated the Trichodesmium community, Clade III abundance was >50% in some NPSG samples, in contrast to the western North Atlantic where Clade III abundance was always <10%. Clade I populations were distributed down to depths >80 m, while Clade III populations were only observed in the mixed layer and found to be significantly correlated with depth and temperature. These data suggest active niche partitioning of Trichodesmium species from different clades, as has been observed in other cyanobacteria. Tracking the distribution and physiology of Trichodesmium spp. would contribute to better predictions of the physiological ecology of this biogeochemically important genus in the present and future ocean.

Mesoscale eddies and Trichodesmium spp. distributions in the southwestern North Atlantic

Trichodesmium abundance was elevated in certain cyclonic and anticyclonic eddiesEnhancement in cyclonic eddies could be driven by Ekman convergenceAnticyclonic eddies with elevated abundance were anomalously fresh.

Expression of a xanthine permease and phosphate transporter in cultures and field populations of the harmful alga Aureococcus anophagefferens: tracking nutritional deficiency during brown tides

Targeted gene expression using quantitative reverse transcription polymerase chain reaction (qRT-PCR) was employed to track patterns in the expression of genes indicative of nitrogen or phosphorus deficiency in the brown tide-forming alga Aureococcus anophagefferens. During culture experiments, a xanthine/uracil/vitamin C permease (XUV) was upregulated ∼20-fold under nitrogen-deficient conditions relative to a nitrogen-replete control and rapidly returned to nitrogen-replete levels after nitrogen-deficient cells were resupplied with nitrate or ammonium. It was not responsive to phosphorus deficiency. Expression of an inorganic phosphate transporter (PTA3) was enriched ∼10-fold under phosphorus-deficient conditions relative to a phosphorus-replete control, and this signal was rapidly lost upon phosphate resupply. PTA3 was not upregulated by nitrogen deficiency. Natural A. anophagefferens populations from a dense brown tide that occurred in Long Island, NY, in 2009 were assayed for XUV and PTA3 expression and compared with nutrient concentrations over the peak of a bloom. Patterns in XUV expression were consistent with nitrogen-replete growth, never reaching the values observed in N-deficient cultures. PTA3 expression was highest prior to peak bloom stages, reaching expression levels within the range of P-deficient cultures. These data highlight the value of molecular-level assessments of nutrient deficiency and suggest that phosphorus deficiency could play a role in the dynamics of destructive A. anophagefferens blooms.

De novo assembly of Aureococcus anophagefferens transcriptomes reveals diverse responses to the low nutrient and low light conditions present during blooms

Transcriptome profiling was performed on the harmful algal bloom-forming pelagophyte Aureococcus anophagefferens strain CCMP 1850 to assess responses to common stressors for dense phytoplankton blooms: low inorganic nitrogen concentrations, low inorganic phosphorus concentrations, low light levels, and a replete control. The de novo assemblies of pooled reads from all treatments reconstructed ~54,000 transcripts using Trinity, and ~31,000 transcripts using ABySS. Comparison to the strain CCMP 1984 genome showed that the majority of the gene models were present in both de novo assemblies and that roughly 95% of contigs from both assemblies mapped to the genome, with Trinity capturing slightly more genome content. Sequence reads were mapped back to the de novo assemblies as well as the gene models and differential expression was analyzed using a Bayesian approach called Analysis of Sequence Counts (ASC). On average, 93% of significantly upregulated transcripts recovered by genome mapping were present in the significantly upregulated pool from both de novo assembly methods. Transcripts related to the transport and metabolism of nitrogen were upregulated in the low nitrogen treatment, transcripts encoding enzymes that hydrolyze organic phosphorus or relieve arsenic toxicity were upregulated in the low phosphorus treatment, and transcripts for enzymes that catabolize organic compounds, restructure lipid membranes, or are involved in sulfolipid biosynthesis were upregulated in the low light treatment. A comparison of this transcriptome to the nutrient regulated transcriptional response of CCMP 1984 identified conserved responses between these two strains. These analyses reveal the transcriptional underpinnings of physiological shifts that could contribute to the ecological success of this species in situ: organic matter processing, metal detoxification, lipid restructuring, and photosynthetic apparatus turnover.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

Identifying reference genes with stable expression from high throughput sequence data

Genes that are constitutively expressed across multiple environmental stimuli are crucial to quantifying differentially expressed genes, particularly when employing quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) assays. However, the identification of these potential reference genes in non-model organisms is challenging and is often guided by expression patterns in distantly related organisms. Here, transcriptome datasets from the diatom Thalassiosira pseudonana grown under replete, phosphorus-limited, iron-limited, and phosphorus and iron co-limited nutrient regimes were analyzed through literature-based searches for homologous reference genes, k-means clustering, and analysis of sequence counts (ASC) to identify putative reference genes. A total of 9759 genes were identified and screened for stable expression. Literature-based searches surveyed 18 generally accepted reference genes, revealing 101 homologs in T. pseudonana with variable expression and a wide range of mean tags per million. k-means analysis parsed the whole transcriptome into 15 clusters. The two most stable clusters contained 709 genes, but still had distinct patterns in expression. ASC analyses identified 179 genes that were stably expressed (posterior probability < 0.1 for 1.25 fold change). Genes known to have a stable expression pattern across the test treatments, like actin, were identified in this pool of 179 candidate genes. ASC can be employed on data without biological replicates and was more robust than the k-means approach in isolating genes with stable expression. The intersection of the genes identified through ASC with commonly used reference genes from the literature suggests that actin and ubiquitin ligase may be useful reference genes for T. pseudonana and potentially other diatoms. With the wealth of transcriptome sequence data becoming available, ASC can be easily applied to transcriptome datasets from other phytoplankton to identify reference genes.

Dissolved Organic Phosphorus Production during Simulated Phytoplankton Blooms in a Coastal Upwelling System

Dissolved organic phosphorus (DOP) is increasingly recognized as an important phosphorus source to marine primary producers. Despite its importance, the production rate and fate of DOP is poorly understood. In this study, patterns of DOP production were evaluated by tracking the evolution of DOP during simulated phytoplankton blooms initiated with nutrient amended surface waters, relative to controls, from the Oregon (USA) coastal upwelling system. Nitrogen (N) and phosphorus (P) additions were used to decouple DOP production and hydrolysis by inducing or repressing, respectively, community alkaline phosphatase activity. In order to examine the progression of nutrient uptake and DOP production under upwelling versus relaxation conditions, two experiments were initiated with waters collected during upwelling events, and two with waters collected during relaxation events. Maximum [under (+P) conditions] and minimum [under (+N) conditions] DOP production rates were calculated and applied to in situ DOP levels to evaluate which end-member rate most closely approximates the in situ DOP production rate at the four study sites in this coastal system. Increases in DOP concentration occurred by day-5 in control treatments in all experiments. N treatments displayed increased chlorophyll a, increased alkaline phosphatase activity, and yielded lower net DOP production rates relative to controls, suggesting that DOP levels were depressed as a consequence of increased hydrolysis of bioavailable DOP substrates. Phosphorus additions resulted in a significant net production of DOP at all stations, but no increase in chlorophyll a relative to control treatments. The contrasting patterns in DOP production between treatments suggests that changes in the ambient dissolved inorganic nitrogen:dissolved inorganic phosphorus (DIN:DIP) ratio could exert profound control over DOP production rates in this system. Patterns of DOP production across the different experiments also suggest that bathymetry-driven differences in water residence times can influence DOP cycling. Taken together, these factors may impact the potential export of DOP to offshore ecosystems.

The transcriptome and proteome of the diatom Thalassiosira pseudonana reveal a diverse phosphorus stress response

Phosphorus (P) is a critical driver of phytoplankton growth and ecosystem function in the ocean. Diatoms are an abundant class of marine phytoplankton that are responsible for significant amounts of primary production. With the control they exert on the oceanic carbon cycle, there have been a number of studies focused on how diatoms respond to limiting macro and micronutrients such as iron and nitrogen. However, diatom physiological responses to P deficiency are poorly understood. Here, we couple deep sequencing of transcript tags and quantitative proteomics to analyze the diatom Thalassiosira pseudonana grown under P-replete and P-deficient conditions. A total of 318 transcripts were differentially regulated with a false discovery rate of <0.05, and a total of 136 proteins were differentially abundant (p<0.05). Significant changes in the abundance of transcripts and proteins were observed and coordinated for multiple biochemical pathways, including glycolysis and translation. Patterns in transcript and protein abundance were also linked to physiological changes in cellular P distributions, and enzyme activities. These data demonstrate that diatom P deficiency results in changes in cellular P allocation through polyphosphate production, increased P transport, a switch to utilization of dissolved organic P through increased production of metalloenzymes, and a remodeling of the cell surface through production of sulfolipids. Together, these findings reveal that T. pseudonana has evolved a sophisticated response to P deficiency involving multiple biochemical strategies that are likely critical to its ability to respond to variations in environmental P availability.

Proteome changes driven by phosphorus deficiency and recovery in the brown tide-forming alga Aureococcus anophagefferens

Shotgun mass spectrometry was used to detect proteins in the harmful alga, Aureococcus anophagefferens, and monitor their relative abundance across nutrient replete (control), phosphate-deficient (-P) and -P refed with phosphate (P-refed) conditions. Spectral counting techniques identified differentially abundant proteins and demonstrated that under phosphate deficiency, A. anophagefferens increases proteins involved in both inorganic and organic phosphorus (P) scavenging, including a phosphate transporter, 5'-nucleotidase, and alkaline phosphatase. Additionally, an increase in abundance of a sulfolipid biosynthesis protein was detected in -P and P-refed conditions. Analysis of the polar membrane lipids showed that cellular concentrations of the sulfolipid sulphoquinovosyldiacylglycerol (SQDG) were nearly two-fold greater in the -P condition versus the control condition, while cellular phospholipids were approximately 8-fold less. Transcript and protein abundances were more tightly coupled for gene products involved in P metabolism compared to those involved in a range of other metabolic functions. Comparison of protein abundances between the -P and P-refed conditions identified differences in the timing of protein degradation and turnover. This suggests that culture studies examining nutrient starvation responses will be valuable in interpreting protein abundance patterns for cellular nutritional status and history in metaproteomic datasets.

Arsenate Resistance in the Unicellular Marine Diazotroph Crocosphaera watsonii

The toxic arsenate ion can behave as a phosphate analog, and this can result in arsenate toxicity especially in areas with elevated arsenate to phosphate ratios like the surface waters of the ocean gyres. In these systems, cellular arsenate resistance strategies would allow phytoplankton to ameliorate the effects of arsenate transport into the cell. Despite the potential coupling between arsenate and phosphate cycling in oligotrophic marine waters, relatively little is known about arsenate resistance in the nitrogen-fixing marine cyanobacteria that are key components of the microbial community in low nutrient systems. The unicellular diazotroph, Crocosphaera watsonii WH8501, was able to grow at reduced rates with arsenate additions up to 30 nM, and estimated arsenate to phosphate ratios of 6:1. The genome of strain WH8501 contains homologs for arsA, arsH, arsB, and arsC, allowing for the reduction of arsenate to arsenite and the pumping of arsenite out of the cell. The short-term addition of arsenate to the growth medium had no effect on nitrogen fixation. However, arsenate addition did result in the up-regulation of the arsB gene with increasing arsenate concentrations, indicating the induction of the arsenate detoxification response. The arsB gene was also up-regulated by phosphorus stress in concert with a gene encoding the high-affinity phosphate binding protein pstS. Both genes were down-regulated when phosphate was re-fed to phosphorus-stressed cells. A field survey of surface water from the low phosphate western North Atlantic detected expression of C. watsoniiarsB, suggestive of the potential importance of arsenate resistance strategies in this and perhaps other systems.

Quorum sensing control of phosphorus acquisition in Trichodesmium consortia

Colonies of the cyanobacterium Trichodesmium are abundant in the oligotrophic ocean, and through their ability to fix both CO(2) and N(2), have pivotal roles in the cycling of carbon and nitrogen in these highly nutrient-depleted environments. Trichodesmium colonies host complex consortia of epibiotic heterotrophic bacteria, and yet, the regulation of nutrient acquisition by these epibionts is poorly understood. We present evidence that epibiotic bacteria in Trichodesmium consortia use quorum sensing (QS) to regulate the activity of alkaline phosphatases (APases), enzymes used by epibionts in the acquisition of phosphate from dissolved-organic phosphorus molecules. A class of QS molecules, acylated homoserine lactones (AHLs), were produced by cultivated epibionts, and adding these AHLs to wild Trichodesmium colonies collected at sea led to a consistent doubling of APase activity. By contrast, amendments of (S)-4,5-dihydroxy-2,3-pentanedione (DPD)-the precursor to the autoinducer-2 (AI-2) family of universal interspecies signaling molecules-led to the attenuation of APase activity. In addition, colonies collected at sea were found by high performance liquid chromatography/mass spectrometry to contain both AHLs and AI-2. Both types of molecules turned over rapidly, an observation we ascribe to quorum quenching. Our results reveal a complex chemical interplay among epibionts using AHLs and AI-2 to control access to phosphate in dissolved-organic phosphorus.

Nutrient-regulated transcriptional responses in the brown tide-forming alga Aureococcus anophagefferens

Long-SAGE (serial analysis of gene expression) was used to profile the transcriptome of the brown tide-forming alga, Aureococcus anophagefferens, under nutrient replete (control), and nitrogen (N) and phosphorus (P) deficiency to understand how this organism responds at the transcriptional level to varying nutrient conditions. This approach has aided A. anophagefferens genome annotation efforts and identified a suite of genes upregulated by N and P deficiency, some of which have known roles in nutrient metabolism. Genes upregulated under N deficiency include an ammonium transporter, an acetamidase/formamidase and two peptidases. This suggests an ability to utilize reduced N compounds and dissolved organic nitrogen, supporting the hypothesized importance of these N sources in A. anophagefferens bloom formation. There are also a broad suite of P-regulated genes, including an alkaline phosphatase, and two 5'-nucleotidases, suggesting A. anophagefferens may use dissolved organic phosphorus under low phosphate conditions. These N- and P-regulated genes may be important targets for exploring nutrient controls on bloom formation in field populations.

Phosphorus supply drives rapid turnover of membrane phospholipids in the diatom Thalassiosira pseudonana

In low-phosphorus (P) marine systems, phytoplankton replace membrane phospholipids with non-phosphorus lipids, but it is not known how rapidly this substitution occurs. Here, when cells of the model diatom Thalassiosira pseudonana were transferred from P-replete medium to P-free medium, the phospholipid content of the cells rapidly declined within 48 h from 45±0.9 to 21±4.5% of the total membrane lipids; the difference was made up by non-phosphorus lipids. Conversely, when P-limited T. pseudonana were resupplied with P, cells reduced the percentage of their total membrane lipids contributed by a non-phosphorus lipid from 43±1.5 to 7.3±0.9% within 24 h, whereas the contribution by phospholipids rose from 2.2±0.1 to 44±3%. This dynamic phospholipid reservoir contained sufficient P to synthesize multiple haploid genomes, suggesting that phospholipid turnover could be an important P source for cells. Field observations of phytoplankton lipid content may thus reflect short-term changes in P supply and cellular physiology, rather than simply long-term adjustment to the environment.

Empirical bayes analysis of sequencing-based transcriptional profiling without replicates

Background

Recent technological advancements have made high throughput sequencing an increasingly popular approach for transcriptome analysis. Advantages of sequencing-based transcriptional profiling over microarrays have been reported, including lower technical variability. However, advances in technology do not remove biological variation between replicates and this variation is often neglected in many analyses.

Results

We propose an empirical Bayes method, titled Analysis of Sequence Counts (ASC), to detect differential expression based on sequencing technology. ASC borrows information across sequences to establish prior distribution of sample variation, so that biological variation can be accounted for even when replicates are not available. Compared to current approaches that simply tests for equality of proportions in two samples, ASC is less biased towards highly expressed sequences and can identify more genes with a greater log fold change at lower overall abundance.

Conclusions

ASC unifies the biological and statistical significance of differential expression by estimating the posterior mean of log fold change and estimating false discovery rates based on the posterior mean. The implementation in R is available at http://www.stat.brown.edu/Zwu/research.aspx.

Parallel analyses of Alexandrium catenella cell concentrations and shellfish toxicity in the Puget Sound

Alexandrium catenella is widespread in western North America and produces a suite of potent neurotoxins that cause paralytic shellfish poisoning (PSP) in humans and have deleterious impacts on public health and economic resources. There are seasonal PSP-related closures of recreational and commercial shellfisheries in the Puget Sound, but the factors that influence cell distribution, abundance, and relationship to paralytic shellfish toxins (PSTs) in this system are poorly described. Here, a quantitative PCR assay was used to detect A. catenella cells in parallel with state shellfish toxicity testing during the 2006 bloom season at 41 sites from April through October. Over 500,000 A. catenella cells liter(-1) were detected at several stations, with two main pulses of cells driving cell distribution, one in June and the other in August. PSTs over the closure limit of 80 mug of PST 100 per g of shellfish tissue were detected at 26 of the 41 sites. Comparison of cell numbers and PST data shows that shellfish toxicity is preceded by an increase in A. catenella cells in 71% of cases. However, cells were also observed in the absence of PSTs in shellfish, highlighting the complex relationship between A. catenella and the resulting shellfish toxicity. These data provide important information on the dynamics of A. catenella cells in the Puget Sound and are a first step toward assessing the utility of plankton monitoring to augment shellfish toxicity testing in this system.

Cell-specific beta-N-acetylglucosaminidase activity in cultures and field populations of eukaryotic marine phytoplankton

It is widely appreciated that eukaryotic marine phytoplankton can hydrolyze a variety of compounds within the dissolved organic matter (DOM) pool in marine environments. Herein, cultures and field populations of marine phytoplankton were assayed for beta-N-acetylglucosaminidase activity, a terminal enzyme of chitin degradation. A traditional bulk assay, which can assess hydrolytic rate, but is not cell-specific, was complemented with a cell-specific assay that images the activity associated with single cells using an enzyme labeled fluorescence (ELF) substrate. beta-N-acetylglucosaminidase activity was widespread across various taxa of marine phytoplankton, and activity was observed both under controlled culture conditions and in field populations. The number of cells with enzyme activity varied with the nutritional physiology of the test species in three of the 17 cultures tested. In these three cases the number of cells with activity in the low nutrient medium was higher than in nutrient replete medium. Taken together, these data suggest that a broad group of marine phytoplankton may be a relevant part of chitin-like DOM degradation and should be incorporated into conceptual models of chitin cycling in marine systems.

Phosphorus scavenging in the unicellular marine diazotroph Crocosphaera watsonii

Through the fixation of atmospheric nitrogen and photosynthesis, marine diazotrophs play a critical role in the global cycling of nitrogen and carbon. Crocosphaera watsonii is a recently described unicellular diazotroph that may significantly contribute to marine nitrogen fixation in tropical environments. One of the many factors that can constrain the growth and nitrogen fixation rates of marine diazotrophs is phosphorus bioavailability. Using genomic and physiological approaches, we examined phosphorus scavenging mechanisms in strains of C. watsonii from both the Atlantic and the Pacific. Observations from the C. watsonii WH8501 genome suggest that this organism has the capacity for high-affinity phosphate transport (e.g., homologs of pstSCAB) in low-phosphate, oligotrophic systems. The pstS gene (high-affinity phosphate binding) is present in strains isolated from both the Atlantic and the Pacific, and its expression was regulated by the exogenous phosphate supply in strain WH8501. Genomic observation also indicated a broad capacity for phosphomonoester hydrolysis (e.g., a putative alkaline phosphatase). In contrast, no clear homologs of genes for phosphonate transport and hydrolysis could be identified. Consistent with these genomic observations, C. watsonii WH8501 is able to grow on phosphomonoesters as a sole source of added phosphorus but not on the phosphonates tested to date. Taken together these data suggest that C. watsonii has a robust capacity for scavenging phosphorus in oligotrophic systems, although this capacity differs from that of other marine cyanobacterial genera, such as Synechococcus, Prochlorococcus, and Trichodesmium.

Phosphonate utilization by the globally important marine diazotroph Trichodesmium

The factors that control the growth and nitrogen fixation rates of marine diazotrophs such as Trichodesmium have been intensively studied because of the role that these processes have in the global cycling of carbon and nitrogen, and in the sequestration of carbon to the deep sea. Because the phosphate concentrations of many ocean gyres are low, the bioavailability of the larger, chemically heterogeneous pool of dissolved organic phosphorus could markedly influence Trichodesmium physiology. Here we describe the induction, by phosphorus stress, of genes from the Trichodesmium erythraeum IMS101 genome that are predicted to encode proteins associated with the high-affinity transport and hydrolysis of phosphonate compounds by a carbon-phosphorus lyase pathway. We show the importance of these genes through expression analyses with T. erythraeum from the Sargasso Sea. Phosphonates are known to be present in oligotrophic marine systems, but have not previously been considered to be bioavailable to marine diazotrophs. The apparent absence of genes encoding a carbon-phosphorus lyase pathway in the other marine cyanobacterial genomes suggests that, relative to other phytoplankton, Trichodesmium is uniquely adapted for scavenging phosphorus from organic sources. This adaptation may help to explain the prevalence of Trichodesmium in low phosphate, oligotrophic systems.

Long serial analysis of gene expression for gene discovery and transcriptome profiling in the widespread marine coccolithophore Emiliania huxleyi

The abundant and widespread coccolithophore Emiliania huxleyi plays an important role in mediating CO2 exchange between the ocean and the atmosphere through its impact on marine photosynthesis and calcification. Here, we use long serial analysis of gene expression (SAGE) to identify E. huxleyi genes responsive to nitrogen (N) or phosphorus (P) starvation. Long SAGE is an elegant approach for examining quantitative and comprehensive gene expression patterns without a priori knowledge of gene sequences via the detection of 21-bp nucleotide sequence tags. E. huxleyi appears to have a robust transcriptional-level response to macronutrient deficiency, with 42 tags uniquely present or up-regulated twofold or greater in the N-starved library and 128 tags uniquely present or up-regulated twofold or greater in the P-starved library. The expression patterns of several tags were validated with reverse transcriptase PCR. Roughly 48% of these differentially expressed tags could be mapped to publicly available genomic or expressed sequence tag (EST) sequence data. For example, in the P-starved library a number of the tags mapped to genes with a role in P scavenging, including a putative phosphate-repressible permease and a putative polyphosphate synthetase. In short, the long SAGE analyses have (i) identified many new differentially regulated gene sequences, (ii) assigned regulation data to EST sequences with no database homology and unknown function, and (iii) highlighted previously uncharacterized aspects of E. huxleyi N and P physiology. To this end, our long SAGE libraries provide a new public resource for gene discovery and transcriptional analysis in this biogeochemically important marine organism.

In situ cytokine production by breast cancer tumor-infiltrating lymphocytes

BACKGROUND

Human breast cancers progressively grow despite the presence of extensive lymphocytic infiltration and specific antitumor immune recognition, thereby calling into question the competency of breast tumor-infiltrating lymphocytes (TIL). The function of breast TILs in vivo and their possible role in the suppression of an antitumor immune response are largely unknown.

METHODS

The cytokines produced in situ by lymphocytes in 89 breast carcinomas and 14 benign breast lesions were assessed using immunohistochemistry.

RESULTS

The majority of tumor and benign breast samples contained T-cell infiltrates, which were disclosed using an anti-CD3 antibody stain. The percentage of tumor samples in which > or =3% of the lymphocytes were producing cytokines was as follows: interleukin (IL)-2 45%, IL-4 36%, tumor necrosis factor-alpha (TNF-alpha) 28%, transforming growth factor-beta 1 (TGF-beta 1) 20%, IL-10 11%, interferon-gamma (IFN-gamma) 4%, and granulocyte-macrophage colony-stimulating factor (GM-CSF) 3%. Production of IL-2, IL-4, and TGF-beta 1 by TILs in breast cancers exceeded that detected in benign breast lesions (p < 0.005). Significantly more tumor samples contained lymphocytes producing IL-2, IL-4, TGF-beta 1, and TNF-alpha than IFN-gamma and GM-CSF (p < 0.002 for each comparison). One or more of the potentially immunoinhibitory cytokines-IL-4, IL-10, or TGF-beta 1-were produced by lymphocytes in 44% of the specimens. No significant associations were seen between lymphocyte production of a particular cytokine and disease-free survival (median follow-up 43 months).

CONCLUSIONS

Immunohistochemical techniques can be used to detect cytokine secretion by TILs in preserved tissue. The relative lack of secretion of IFN-gamma and GM-CSF, rather than a deficiency of IL-2, may explain why the antitumor immune response to breast cancer is impaired.