RPANDA: an R package for macroevolutionary analyses on phylogenetic trees

A number of approaches for studying macroevolution using phylogenetic trees have been developed in the last few years. Here, we present RPANDA, an R package that implements model‐free and model‐based phylogenetic comparative methods for macroevolutionary analyses. The model‐free approaches implemented in RPANDA are recently developed approaches stemming from graph theory that allow summarizing the information contained in phylogenetic trees, computing distances between trees, and clustering them accordingly. They also allow identifying distinct branching patterns within single trees. RPANDA also implements likelihood‐based models for fitting various diversification models to phylogenetic trees. It includes birth–death models with i) constant, ii) time‐dependent and iii) environmental‐dependent speciation and extinction rates. It also includes models with equilibrium diversity derived from the coalescent process, as well as a likelihood‐based inference framework to fit the individual‐based model of Speciation by Genetic Differentiation, which is an extension of Hubbell's neutral theory of biodiversity. RPANDA can be used to (i) characterize trees by plotting their spectral density profiles (ii) compare trees and cluster them according to their similarities, (iii) identify and plot distinct branching patterns within trees, (iv) compare the fit of alternative diversification models to phylogenetic trees, (v) estimate rates of speciation and extinction, (vi) estimate and plot how these rates have varied with time and environmental variables and (vii) deduce and plot estimates of species richness through geological time. RPANDA provides investigators with a set of tools for exploring patterns in phylogenetic trees and fitting various models to these trees, thereby contributing to the ongoing development of phylogenetics in the life sciences

Cancer-Associated Fibroblasts Neutralize the Anti-tumor Effect of CSF1 Receptor Blockade by Inducing PMN-MDSC Infiltration of Tumors.

Tumor-associated macrophages (TAM) contribute to all aspects of tumor progression. Use of CSF1R inhibitors to target TAM is therapeutically appealing, but has had very limited anti-tumor effects. Here, we have identified the mechanism that limited the effect of CSF1R targeted therapy. We demonstrated that carcinoma-associated fibroblasts (CAF) are major sources of chemokines that recruit granulocytes to tumors. CSF1 produced by tumor cells caused HDAC2-mediated downregulation of granulocyte-specific chemokine expression in CAF, which limited migration of these cells to tumors. Treatment with CSF1R inhibitors disrupted this crosstalk and triggered a profound increase in granulocyte recruitment to tumors. Combining CSF1R inhibitor with a CXCR2 antagonist blocked granulocyte infiltration of tumors and showed strong anti-tumor effects

Toward a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa.

Rapidly growing biological data-including molecular sequences and fossils-hold an unprecedented potential to reveal how evolutionary processes generate and maintain biodiversity. However, researchers often have to develop their own idiosyncratic workflows to integrate and analyze these data for reconstructing time-calibrated phylogenies. In addition, divergence times estimated under different methods and assumptions, and based on data of various quality and reliability, should not be combined without proper correction. Here we introduce a modular framework termed SUPERSMART (Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa), and provide a proof of concept for dealing with the moving targets of evolutionary and biogeographical research. This framework assembles comprehensive data sets of molecular and fossil data for any taxa and infers dated phylogenies using robust species tree methods, also allowing for the inclusion of genomic data produced through next-generation sequencing techniques. We exemplify the application of our method by presenting phylogenetic and dating analyses for the mammal order Primates and for the plant family Arecaceae (palms). We believe that this framework will provide a valuable tool for a wide range of hypothesis-driven research questions in systematics, biogeography, and evolution. SUPERSMART will also accelerate the inference of a "Dated Tree of Life" where all node ages are directly comparable. [Bayesian phylogenetics; data mining; divide-and-conquer methods; GenBank; multilocus multispecies coalescent; next-generation sequencing; palms; primates; tree calibration.]

ER stress regulates myeloid-derived suppressor cell fate through TRAIL-R–mediated apoptosis

Myeloid-derived suppressor cells (MDSCs) dampen the immune response thorough inhibition of T cell activation and proliferation and often are expanded in pathological conditions. Here, we studied the fate of MDSCs in cancer. Unexpectedly, MDSCs had lower viability and a shorter half-life in tumor-bearing mice compared with neutrophils and monocytes. The reduction of MDSC viability was due to increased apoptosis, which was mediated by increased expression of TNF-related apoptosis–induced ligand receptors (TRAIL-Rs) in these cells. Targeting TRAIL-Rs in naive mice did not affect myeloid cell populations, but it dramatically reduced the presence of MDSCs and improved immune responses in tumor-bearing mice. Treatment of myeloid cells with proinflammatory cytokines did not affect TRAIL-R expression; however, induction of ER stress in myeloid cells recapitulated changes in TRAIL-R expression observed in tumor-bearing hosts. The ER stress response was detected in MDSCs isolated from cancer patients and tumor-bearing mice, but not in control neutrophils or monocytes, and blockade of ER stress abrogated tumor-associated changes in TRAIL-Rs. Together, these data indicate that MDSC pathophysiology is linked to ER stress, which shortens the lifespan of these cells in the periphery and promotes expansion in BM. Furthermore, TRAIL-Rs can be considered as potential targets for selectively inhibiting MDSCs

Phylogeny and Historical Biogeography of Asian Pterourus Butterflies (Lepidoptera: Papilionidae): A Case of Intercontinental Dispersal from North America to East Asia

The phylogenetic status of the well-known Asian butterflies often known as Agehana (a species group, often treated as a genus or a subgenus, within Papilio sensu lato) has long remained unresolved. Only two species are included, and one of them especially, Papilio maraho, is not only rare but near-threatened, being monophagous on its vulnerable hostplant, Sassafras randaiense (Lauraceae). Although the natural history and population conservation of “Agehana” has received much attention, the biogeographic origin of this group still remains enigmatic. To clarify these two questions, a total of 86 species representatives within Papilionidae were sampled, and four genes (concatenated length 3842 bp) were used to reconstruct their phylogenetic relationships and historical scenarios. Surprisingly, “Agehana” fell within the American Papilio subgenus Pterourus and not as previously suggested, phylogenetically close to the Asian Papilio subgenus Chilasa. We therefore formally synonymize Agehana with Pterourus. Dating and biogeographic analysis allow us to infer an intercontinental dispersal of an American ancestor of Asian Pterourus in the early Miocene, which was coincident with historical paleo-land bridge connections, resulting in the present “East Asia-America” disjunction distribution. We emphasize that species exchange between East Asia and America seems to be a quite frequent occurrence in butterflies during the Oligocene to Miocene climatic optima.© 2015 Wu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage

Rice yield is most sensitive to salinity stress imposed during the panicle initiation (PI) stage. In this study, we have focused on physiological and transcriptional responses of four rice genotypes exposed to salinity stress during PI. The genotypes selected included a pair of indicas (IR63731 and IR29) and a pair of japonica (Agami and M103) rice subspecies with contrasting salt tolerance. Physiological characterization showed that tolerant genotypes maintained a much lower shoot Na(+) concentration relative to sensitive genotypes under salinity stress. Global gene expression analysis revealed a strikingly large number of genes which are induced by salinity stress in sensitive genotypes, IR29 and M103 relative to tolerant lines. We found 19 probe sets to be commonly induced in all four genotypes. We found several salinity modulated, ion homeostasis related genes from our analysis. We also studied the expression of SKC1, a cation transporter reported by others as a major source of variation in salt tolerance in rice. The transcript abundance of SKC1 did not change in response to salinity stress at PI stage in the shoot tissue of all four genotypes. However, we found the transcript abundance of SKC1 to be significantly higher in tolerant japonica Agami relative to sensitive japonica M103 under control and stressed conditions during PI stage. ELECTRONIC SUPPLEMENTARY MATERIAL: Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s11103-006-9112-0 and is accessible for authorized users

The Speed of Sound in Methane under Conditions of the Thermal Boundary Layer of Uranus

We present the first direct observations of acoustic waves in warm dense matter. We analyze wavenumber- and energy-resolved X-ray spectra taken from warm dense methane created by laser-heating a cryogenic liquid jet. X-ray diffraction and inelastic free electron scattering yield sample conditions of 0.3$\pm$0.1 eV and 0.8$\pm$0.1 g/cm$^3$, corresponding to a pressure of $\sim$13 GPa and matching the conditions predicted in the thermal boundary layer between the inner and outer envelope of Uranus. Inelastic X-ray scattering was used to observe the collective oscillations of the ions. With a highly improved energy resolution of $\sim$50 meV, we could clearly distinguish the Brillouin peaks from the quasi-elastic Rayleigh feature. Data at different wavenumbers were used to obtain a sound speed of 5.9$\pm$0.5 km/s, which enabled us to validate the use of Birch's law in this new parameter regime.Comment: 7 pages, 4 figures with supplementary informatio

IRE1α–XBP1 controls T cell function in ovarian cancer by regulating mitochondrial activity

Tumours evade immune control by creating hostile microenvironments that perturb T cell metabolism and effector function 1?4 . However, it remains unclear how intra-tumoral T cells integrate and interpret metabolic stress signals. Here we report that ovarian cancer?an aggressive malignancy that is refractory to standard treatments and current immunotherapies 5?8 ?induces endoplasmic reticulum stress and activates the IRE1α?XBP1 arm of the unfolded protein response 9,10 in T cells to control their mitochondrial respiration and anti-tumour function. In T cells isolated from specimens collected from patients with ovarian cancer, upregulation of XBP1 was associated with decreased infiltration of T cells into tumours and with reduced IFNG mRNA expression. Malignant ascites fluid obtained from patients with ovarian cancer inhibited glucose uptake and caused N-linked protein glycosylation defects in T cells, which triggered IRE1α?XBP1 activation that suppressed mitochondrial activity and IFNγ production. Mechanistically, induction of XBP1 regulated the abundance of glutamine carriers and thus limited the influx of glutamine that is necessary to sustain mitochondrial respiration in T cells under glucose-deprived conditions. Restoring N-linked protein glycosylation, abrogating IRE1α?XBP1 activation or enforcing expression of glutamine transporters enhanced mitochondrial respiration in human T cells exposed to ovarian cancer ascites. XBP1-deficient T cells in the metastatic ovarian cancer milieu exhibited global transcriptional reprogramming and improved effector capacity. Accordingly, mice that bear ovarian cancer and lack XBP1 selectively in T cells demonstrate superior anti-tumour immunity, delayed malignant progression and increased overall survival. Controlling endoplasmic reticulum stress or targeting IRE1α?XBP1 signalling may help to restore the metabolic fitness and anti-tumour capacity of T cells in cancer hosts.Fil: Song, Minkyung. Weill Cornell Medicine; Estados UnidosFil: Sandoval, Tito A.. Weill Cornell Medicine; Estados UnidosFil: Chae, Chang-Suk. Weill Cornell Medicine; Estados UnidosFil: Chopra, Sahil. Weill Cornell Medicine; Estados UnidosFil: Tan, Chen. Weill Cornell Medicine; Estados UnidosFil: Rutkowski, Melanie R.. University of Virginia; Estados UnidosFil: Raundhal, Mahesh. Dana Farber Cancer Institute; Estados Unidos. Harvard Medical School; Estados UnidosFil: Chaurio, Ricardo A.. H. Lee Moffitt Cancer Center & Research Institute; Estados UnidosFil: Payne, Kyle K.. H. Lee Moffitt Cancer Center & Research Institute; Estados UnidosFil: Konrad, Csaba. Weill Cornell Medicine; Estados UnidosFil: Bettigole, Sarah E.. Quentis Therapeutics Inc.; Estados UnidosFil: Shin, Hee Rae. Quentis Therapeutics Inc.; Estados UnidosFil: Crowley, Michael J. P.. Weill Cornell Graduate School of Medical Sciences; Estados UnidosFil: Cerliani, Juan Pablo. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Biología y Medicina Experimental. Fundación de Instituto de Biología y Medicina Experimental. Instituto de Biología y Medicina Experimental; ArgentinaFil: Kossenkov, Andrew V.. The Wistar Institute; Estados UnidosFil: Motorykin, Ievgen. Weill Cornell Medicine,; Estados UnidosFil: Zhang, Sheng. Weill Cornell Medicine,; Estados UnidosFil: Manfredi, Giovanni. Weill Cornell Medicine,; Estados UnidosFil: Zamarin, Dmitriy. Memorial Sloan Kettering Cancer Center; Estados UnidosFil: Holcomb, Kevin. Weill Cornell Medicine,; Estados UnidosFil: Rodriguez, Paulo C.. H. Lee Moffitt Cancer Center & Research Institute; Estados UnidosFil: Rabinovich, Gabriel Adrián. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Biología y Medicina Experimental. Fundación de Instituto de Biología y Medicina Experimental. Instituto de Biología y Medicina Experimental; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Química Biológica; ArgentinaFil: Conejo Garcia, Jose R.. H. Lee Moffitt Cancer Center & Research Institute; Estados UnidosFil: Glimcher, Laurie H.. Dana Farber Cancer Institute; Estados Unidos. Harvard Medical School; Estados UnidosFil: Cubillos-Ruiz, Juan R.. Weill Graduate School Of Medical Sciences; Estados Unidos. Weill Graduate School Of Medical Sciences; Estados Unido

Congruence and diversity of butterfly-host plant associations at higher taxonomic levels

We aggregated data on butterfly-host plant associations from existing sources in order to address the following questions: (1) is there a general correlation between host diversity and butterfly species richness?, (2) has the evolution of host plant use followed consistent patterns across butterfly lineages?, (3) what is the common ancestral host plant for all butterfly lineages? The compilation included 44,148 records from 5,152 butterfly species (28.6% of worldwide species of Papilionoidea) and 1,193 genera (66.3%). The overwhelming majority of butterflies use angiosperms as host plants. Fabales is used by most species (1,007 spp.) from all seven butterfly families and most subfamilies, Poales is the second most frequently used order, but is mostly restricted to two species-rich subfamilies: Hesperiinae (56.5% of all Hesperiidae), and Satyrinae (42.6% of all Nymphalidae). We found a significant and strong correlation between host plant diversity and butterfly species richness. A global test for congruence (Parafit test) was sensitive to uncertainty in the butterfly cladogram, and suggests a mixed system with congruent associations between Papilionidae and magnoliids, Hesperiidae and monocots, and the remaining subfamilies with the eudicots (fabids and malvids), but also numerous random associations. The congruent associations are also recovered as the most probable ancestral states in each node using maximum likelihood methods. The shift from basal groups to eudicots appears to be more likely than the other way around, with the only exception being a Satyrine-clade within the Nymphalidae that feed on monocots. Our analysis contributes to the visualization of the complex pattern of interactions at superfamily level and provides a context to discuss the timing of changes in host plant utilization that might have promoted diversification in some butterfly lineages