Heterogeneity of Pancreatic Cancer Metastases in a Single Patient Revealed by Quantitative Proteomics

Many patients with pancreatic cancer have metastases to distant organs at the time of initial presentation. Recent studies examining the evolution of pancreatic cancer at the genetic level have shown that clonal complexity of metastatic pancreatic cancer is already initiated within primary tumors, and organ-specific metastases are derived from different subclones. However, we do not yet understand to what extent the evolution of pancreatic cancer contributes to proteomic and signaling alterations. We hypothesized that genetic heterogeneity of metastatic pancreatic cancer results in heterogeneity at the proteome level. To address this, we employed a model system in which cells isolated from three sites of metastasis (liver, lung, and peritoneum) from a single patient were compared. We used a SILAC-based accurate quantitative proteomic strategy combined with high-resolution mass spectrometry to analyze the total proteome and tyrosine phosphoproteome of each of the distal metastases. Our data revealed distinct patterns of both overall proteome expression and tyrosine kinase activities across the three different metastatic lesions. This heterogeneity was significant because it led to differential sensitivity of the neoplastic cells to small molecule inhibitors targeting various kinases and other pathways. For example, R428, a tyrosine kinase inhibitor that targets Axl receptor tyrosine kinase, was able to inhibit cells derived from lung and liver metastases much more effectively than cells from the peritoneal metastasis. Finally, we confirmed that administration of R428 in mice bearing xenografts of cells derived from the three different metastatic sites significantly diminished tumors formed from liver- and lung-metastasis-derived cell lines as compared with tumors derived from the peritoneal metastasis cell line. Overall, our data provide proof-of-principle support that personalized therapy of multiple organ metastases in a single patient should involve the administration of a combination of agents, with each agent targeted to the features of different subclones.Approximately half of the patients with pancreatic cancer are initially diagnosed with metastases to distal sites, with the commonest sites being the liver, lung, and peritoneum (1). Therapeutic strategies against metastases could help reduce the high mortality rates associated with this cancer (2). Understanding the nature of metastatic pancreatic cancer at a systems level can enable the discovery of potential targets for the development of targeted therapies.Pancreatic cancer has been shown to be a genetically evolving and heterogeneous disease (3–5). Clonal diversity and evolution of cancer genomes have also been demonstrated based on the isolation of distinct clonal populations purified directly from patient biopsies by means of flow cytometry followed by genomic characterization (6). A number of reports have documented the adoption of a proteomic approach for the discovery of potential biomarkers in pancreatic cancer (7, 8). However, these studies generally assume pancreatic cancers to be homogeneous, and the emphasis is placed on identifying molecules that are common across a broad array of tumors. There is a lack of studies systematically examining the proteomic changes or signaling pathways across pancreatic cancers to dissect the nature of the heterogeneity of each clone. An excellent setting in which the heterogeneity of tumors can be studied systematically is in a patient harboring metastases to several distant sites. To this end, we chose cells isolated from three metastatic pancreatic lesions of a single patient. The exomes of each tumor site were previously sequenced to study the progression of pancreatic cancer, and the results showed that all cell lines were identical for the genetic status of driver mutations (e.g. KRAS, TP53, and SMAD4) (9). Our hypothesis was that a better understanding of the proteomic consequences of the heterogeneity derived from genetic changes, and possibly other types of alterations, might provide additional opportunities to identify therapeutic targets.In order to precisely quantify differences across the proteomes of multiple metastatic pancreatic cancer lesions, we employed a SILAC-based¹ quantitative proteomics strategy combined with high-resolution mass spectrometry (10, 11). Based on changes observed at the whole-proteome level, we found that a class of cell surface receptors showed significant enrichment with the highest alteration of their expression among the three metastatic pancreatic cancer cell lines examined (i.e. peritoneum, lung, and liver). Because the total protein levels provide information about the static levels of proteins and not their activity per se, we decided to examine the activation of phosphorylation-driven pathways, many of which are activated by cell surface receptors. To globally examine tyrosine phosphorylation-based signaling pathways, we carried out mass spectrometric analysis of purified tyrosine phosphorylated peptides enriched using anti-phosphotyrosine antibodies. As a result, we observed differential activation of tyrosine kinases in the three different sites of metastases. For example, Axl receptor tyrosine kinase was found to be hyperphosphorylated in lung and liver metastases relative to peritoneal metastasis. Expression of Axl receptor tyrosine kinase in primary and matched pancreatic cancers on tissue microarrays was validated by immunohistochemistry. Given such unique patterns of activation of pathways, it was possible that tumors derived from different sites could show differences in their sensitivity to pathway inhibitors. To test this, we performed experiments in which we screened cell lines derived from each metastatic site against a panel of small molecule inhibitors. We observed that the three metastatic pancreatic cancers had differential sensitivities to different inhibitors. For example, cells derived from the peritoneal metastasis were highly sensitive to lapatinib, whereas greater sensitivity to the Axl inhibitor R428 was observed in the lung metastasis cell line. Finally, we showed that treatment of mice bearing xenografts from these different pancreatic cancer cell lines with R428, an inhibitor of Axl receptor tyrosine kinase, led to reduction of tumors with evidence of activation of Axl.     

Subcellular Localization of Rat Brain Insulin Binding Sites

Fractions and subcellular structures were prepared from rat brain homogenate and their purity was assessed using enzyme markers, gamma-aminobutyric acid binding, DNA content, and electron microscopy. Insulin binding was highest on the plasma membrane preparations and approximately 50% less so on brain homogenate crude mitochondrial (P2), myelinated axon, and synaptosome preparations. Very low levels of binding were found on mitochondria and nuclei. Differences in binding between fractions were due to numbers of binding sites, and not variable binding affinity. There was a close relationship between insulin binding and the activity of Na/K ATPase (E.C. 3.6.1.4) in all fractions (r = 0.98). Insulin binding to the P2 was compared with plasma membrane fractions in seven brain regions, and the results demonstrated the same close relationship between insulin binding and plasma membrane content in all regions except hypothalamus. Plasma membrane insulin binding was well represented by the binding on P2 membranes in all regions except hypothalamus and brainstem. It was concluded that insulin binding is distributed evenly over the surface of brain cells and is not increased on nerve endings.     

A multi-trait systems approach reveals a response cascade to bleaching in corals

Background

Climate change causes the breakdown of the symbiotic relationships between reef-building corals and their photosynthetic symbionts (genus Symbiodinium), with thermal anomalies in 2015–2016 triggering the most widespread mass coral bleaching on record and unprecedented mortality on the Great Barrier Reef. Targeted studies using specific coral stress indicators have highlighted the complexity of the physiological processes occurring during thermal stress, but have been unable to provide a clear mechanistic understanding of coral bleaching.

Results

Here, we present an extensive multi-trait-based study in which we compare the thermal stress responses of two phylogenetically distinct and widely distributed coral species, Acropora millepora and Stylophora pistillata, integrating 14 individual stress indicators over time across a simulated thermal anomaly. We found that key stress responses were conserved across both taxa, with the loss of symbionts and the activation of antioxidant mechanisms occurring well before collapse of the physiological parameters, including gross oxygen production and chlorophyll a. Our study also revealed species-specific traits, including differences in the timing of antioxidant regulation, as well as drastic differences in the production of the sulfur compound dimethylsulfoniopropionate during bleaching. Indeed, the concentration of this antioxidant increased two-fold in A. millepora after the corals started to bleach, while it decreased 70% in S. pistillata.

Conclusions

We identify a well-defined cascading response to thermal stress, demarking clear pathophysiological reactions conserved across the two species, which might be central to fully understanding the mechanisms triggering thermally induced coral bleaching. These results highlight that bleaching is a conserved mechanism, but specific adaptations linked to the coral’s antioxidant capacity drive differences in the sensitivity and thus tolerance of each coral species to thermal stress.

     

Sexuality and Genetic Identity in the Agaricus Section Arvenses

Twelve wild collections and one commercial strain were used to characterize breeding systems and to develop molecular identities in the Arvenses section of the genus Agaricus, which includes the “horse mushroom” A. arvensis. Two morphotypes were identified based on macro- and micromorphological features. However, not all collections could be delimited by conventional taxonomic characters. Sequencing of the small subunit intergenic spacer (ITS) region (368 to 370 bp) of the rRNA genes clearly resolved the 13 collections into two clusters consistent with the identified morphotypes. Single-spore progenies and mating type testers were established and used to test intra- and interstock compatibility. The two compatibility groups identified were consistent with ITS clusters. Compatibility group I stocks readily interbred within the constraints of a unifactorial heterothallic system with a multiallelic mating type factor. Compatibility group II had a more restricted breeding pattern, and interactions were difficult to predict on the basis of mating type. Morphological data, ITS sequences, and the ability to interbreed suggest that these collections are part of a complex of interrelated species. Single-spore, homokaryotic isolates from both compatibility groups were able to fruit in compost culture, and two of the collections may represent natural homokaryotic fruiting. We conclude that species from the section Arvenses have versatile unifactorial heterothallic life cycles that permit both interbreeding and homokaryotic fruiting.     

Climate change promotes transitions to tall evergreen vegetation in tropical Asia

Vegetation in tropical Asia is highly diverse due to large environmental gradients and heterogeneity of landscapes. This biodiversity is threatened by intense land use and climate change. However, despite the rich biodiversity and the dense human population, tropical Asia is often underrepresented in global biodiversity assessments. Understanding how climate change influences the remaining areas of natural vegetation is therefore highly important for conservation planning. Here, we used the adaptive Dynamic Global Vegetation Model version 2 (aDGVM2) to simulate impacts of climate change and elevated CO₂ on vegetation formations in tropical Asia for an ensemble of climate change scenarios. We used climate forcing from five different climate models for representative concentration pathways RCP4.5 and RCP8.5. We found that vegetation in tropical Asia will remain a carbon sink until 2099, and that vegetation biomass increases of up to 28% by 2099 are associated with transitions from small to tall woody vegetation and from deciduous to evergreen vegetation. Patterns of phenology were less responsive to climate change and elevated CO₂ than biomes and biomass, indicating that the selection of variables and methods used to detect vegetation changes is crucial. Model simulations revealed substantial variation within the ensemble, both in biomass increases and in distributions of different biome types. Our results have important implications for management policy, because they suggest that large ensembles of climate models and scenarios are required to assess a wide range of potential future trajectories of vegetation change and to develop robust management plans. Furthermore, our results highlight open ecosystems with low tree cover as most threatened by climate change, indicating potential conflicts of interest between biodiversity conservation in open ecosystems and active afforestation to enhance carbon sequestration.     

Estrogen and Progestogen Correlates of the Structure of Female Copulation Calls in Semi-Free-Ranging Barbary Macaques (<Emphasis Type="BoldItalic">Macaca sylvanus</Emphasis>)

Females of many Old World primates produce conspicuous vocalizations in combination with copulations. Indirect evidence exists that in Barbary macaques (Macaca sylvanus), the structure of these copulation calls is related to changes in reproductive hormone levels. However, the structure of these calls does not vary significantly around the timing of ovulation when estrogen and progestogen levels show marked changes. We here aimed to clarify this paradox by investigating how the steroid hormones estrogen and progesterone are related to changes in the acoustic structure of copulation calls. We collected data on semi-free-ranging Barbary macaques in Gibraltar and at La Forêt des Singes in Rocamadour, France. We determined estrogen and progestogen concentrations from fecal samples and combined them with a fine-grained structural analysis of female copulation calls (N = 775 calls of 11 females). Our analysis indicates a time lag of 3 d between changes in fecal hormone levels, adjusted for the excretion lag time, and in the acoustic structure of copulation calls. Specifically, we found that estrogen increased the duration and frequency of the calls, whereas progestogen had an antagonistic effect. Importantly, however, variation in acoustic variables did not track short-term changes such as the peak in estrogen occurring around the timing of ovulation. Taken together, our results help to explain why female Barbary macaque copulation calls are related to changes in hormone levels but fail to indicate the fertile phase.     

Trophic cul-de-sac, Pyrazus ebeninus, limits trophic transfer through an estuarine detritus-based food web

The importance to food‐webs of trophic cul‐de‐sacs, species that channel energy flow away from higher trophic levels, is seldom considered outside of the pelagic systems in which they were first identified. On intertidal mudflats, inputs of detritus from saltmarshes, macroalgae or microphytobenthos are generally regarded as a major structuring force underpinning food‐webs and there has been no consideration of trophic cul‐de‐sacs to date. A fully orthogonal three‐factor experiment manipulating the density of the abundant gastropod, Pyrazus ebeninus, detritus and macrobenthic predators on a Sydney mudflat revealed large deleterious effects of the gastropod, irrespective of detrital loading or the presence of predators. Two months after experimental manipulation, the standing‐stock of microphytobenthos in plots with high (44 per m²) densities of P. ebeninus was 20% less than in plots with low (4 per m²) densities. Increasing densities of P. ebeninus from low to high halved the abundance of macroinvertebrates and the average number of species. In contrast, the addition of detritus had differing effects on microphytobenthos (positively affected) and macroinvertebrates (negatively affected). Over the two‐months of our experiment, no predatory mortality of P. ebeninus was observed and high densities of P. ebeninus decreased impacts of predators on macroinvertebrate abundances. Given that the dynamics of southeast Australian mudflats are driven more by disturbance than seasonality in predators and their interactions with prey, it is likely that Pyrazus would be similarly resistant to predation and have negative effects on benthic assemblages at other times of the year, outside of our study period. Thus, in reducing microphytobenthos and the abundance and species richness of macrofauna, high abundances of the detritivore P. ebeninus may severely limit the flow of energy up the food chain to commercially‐important species. This study therefore suggests that trophic cul‐de‐sacs are not limited to the eutrophied pelagic systems in which they were first identified, but may exist in other systems as well.     

Mechanisms of zoonotic severe acute respiratory syndrome coronavirus host range expansion in human airway epithelium

In 2003, severe acute respiratory syndrome coronavirus (SARS-CoV) emerged and caused over 8,000 human cases of infection and more than 700 deaths worldwide. Zoonotic SARS-CoV likely evolved to infect humans by a series of transmission events between humans and animals for sale in China. Using synthetic biology, we engineered the spike protein (S) from a civet strain, SZ16, into our epidemic strain infectious clone, creating the chimeric virus icSZ16-S, which was infectious but yielded progeny viruses incapable of propagating in vitro. After introducing a K479N mutation within the S receptor binding domain (RBD) of SZ16, the recombinant virus (icSZ16-S K479N) replicated in Vero cells but was severely debilitated in growth. The in vitro evolution of icSZ16-S K479N on human airway epithelial (HAE) cells produced two viruses (icSZ16-S K479N D8 and D22) with enhanced growth on HAE cells and on delayed brain tumor cells expressing the SARS-CoV receptor, human angiotensin I converting enzyme 2 (hACE2). The icSZ16-S K479N D8 and D22 virus RBDs contained mutations in ACE2 contact residues, Y442F and L472F, that remodeled S interactions with hACE2. Further, these viruses were neutralized by a human monoclonal antibody (MAb), S230.15, but the parent icSZ16-S K479N strain was eight times more resistant than the mutants. These data suggest that the human adaptation of zoonotic SARS-CoV strains may select for some variants that are highly susceptible to select MAbs that bind to RBDs. The epidemic, icSZ16-S K479N, and icSZ16-S K479N D22 viruses replicate similarly in the BALB/c mouse lung, highlighting the potential use of these zoonotic spike SARS-CoVs to assess vaccine or serotherapy efficacy in vivo.