spatialGE: quantification and visualization of the tumor microenvironment heterogeneity using spatial transcriptomics.

SUMMARY: Spatially resolved transcriptomics promises to increase our understanding of the tumor microenvironment and improve cancer prognosis and therapies. Nonetheless, analytical methods to explore associations between the spatial heterogeneity of the tumor and clinical data are not available. Hence, we have developed spatialGE, a software that provides visualizations and quantification of the tumor microenvironment heterogeneity through gene expression surfaces, spatial heterogeneity statistics that can be compared against clinical information, spot-level cell deconvolution and spatially informed clustering, all using a new data object to store data and resulting analyses simultaneously. AVAILABILITY AND IMPLEMENTATION: The R package and tutorial/vignette are available at https://github.com/FridleyLab/spatialGE. A script to reproduce the analyses in this manuscript is available in Supplementary information. The Thrane study data included in spatialGE was made available from the public available from the website https://www.spatialresearch.org/resources-published-datasets/doi-10-1158-0008-5472-can-18-0747/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

spatialTIME and iTIME: R package and Shiny application for visualization and analysis of immunofluorescence data.

SUMMARY: Multiplex immunofluorescence (mIF) staining combined with quantitative digital image analysis is a novel and increasingly used technique that allows for the characterization of the tumor immune microenvironment (TIME). Generally, mIF data is used to examine the abundance of immune cells in the TIME; however, this does not capture spatial patterns of immune cells throughout the TIME, a metric increasingly recognized as important for prognosis. To address this gap, we developed an R package spatialTIME that enables spatial analysis of mIF data, as well as the iTIME web application that provides a robust but simplified user interface for describing both abundance and spatial architecture of the TIME. The spatialTIME package calculates univariate and bivariate spatial statistics (e.g. Ripley's K, Besag's L, Macron's M and G or nearest neighbor distance) and creates publication quality plots for spatial organization of the cells in each tissue sample. The iTIME web application allows users to statistically compare the abundance measures with patient clinical features along with visualization of the TIME for one tissue sample at a time. AVAILABILITY AND IMPLEMENTATION: spatialTIME is implemented in R and can be downloaded from GitHub (https://github.com/FridleyLab/spatialTIME) or CRAN. An extensive vignette for using spatialTIME can also be found at https://cran.r-project.org/web/packages/spatialTIME/index.html. iTIME is implemented within a R Shiny application and can be accessed online (http://itime.moffitt.org/), with code available on GitHub (https://github.com/FridleyLab/iTIME). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Neurogenomic divergence during speciation by reinforcement of mating behaviors in chorus frogs (Pseudacris).

BACKGROUND: Species interactions can promote mating behavior divergence, particularly when these interactions are costly due to maladaptive hybridization. Selection against hybridization can indirectly cause evolution of reproductive isolation within species, a process termed cascade reinforcement. This process can drive incipient speciation by generating divergent selection pressures among populations that interact with different species assemblages. Theoretical and empirical studies indicate that divergent selection on gene expression networks has the potential to increase reproductive isolation among populations. After identifying candidate synaptic transmission genes derived from neurophysiological studies in anurans, we test for divergence of gene expression in a system undergoing cascade reinforcement, the Upland Chorus Frog (Pseudacris feriarum). RESULTS: Our analyses identified seven candidate synaptic transmission genes that have diverged between ancestral and reinforced populations of P. feriarum, including five that encode synaptic vesicle proteins. Our gene correlation network analyses revealed four genetic modules that have diverged between these populations, two possessing a significant concentration of neurotransmission enrichment terms: one for synaptic membrane components and the other for metabolism of the neurotransmitter nitric oxide. We also ascertained that a greater number of genes have diverged in expression by geography than by sex. Moreover, we found that more genes have diverged within females as compared to males between populations. Conversely, we observed no difference in the number of differentially-expressed genes within the ancestral compared to the reinforced population between the sexes. CONCLUSIONS: This work is consistent with the idea that divergent selection on mating behaviors via cascade reinforcement contributed to evolution of gene expression in P. feriarum. Although our study design does not allow us to fully rule out the influence of environment and demography, the fact that more genes diverged in females than males points to a role for cascade reinforcement. Our discoveries of divergent candidate genes and gene networks related to neurotransmission support the idea that neural mechanisms of acoustic mating behaviors have diverged between populations, and agree with previous neurophysiological studies in frogs. Increasing support for this hypothesis, however, will require additional experiments under common garden conditions. Our work points to the importance of future replicated and tissue-specific studies to elucidate the relative contribution of gene expression divergence to the evolution of reproductive isolation during incipient speciation.

Hierarchical Hybrid Enrichment: Multitiered Genomic Data Collection Across Evolutionary Scales, With Application to Chorus Frogs (Pseudacris).

Determining the optimal targets of genomic subsampling for phylogenomics, phylogeography, and population genomics remains a challenge for evolutionary biologists. Of the available methods for subsampling the genome, hybrid enrichment (sequence capture) has become one of the primary means of data collection for systematics, due to the flexibility and cost efficiency of this approach. Despite the utility of this method, information is lacking as to what genomic targets are most appropriate for addressing questions at different evolutionary scales. In this study, first, we compare the benefits of target loci developed for deep- and shallow scales by comparing these loci at each of three taxonomic levels: within a genus (phylogenetics), within a species (phylogeography), and within a hybrid zone (population genomics). Specifically, we target evolutionarily conserved loci that are appropriate for deeper phylogenetic scales and more rapidly evolving loci that are informative for phylogeographic and population genomic scales. Second, we assess the efficacy of targeting multiple-locus sets for different taxonomic levels in the same hybrid enrichment reaction, an approach we term hierarchical hybrid enrichment. Third, we apply this approach to the North American chorus frog genus Pseudacris to answer key evolutionary questions across taxonomic and temporal scales. We demonstrate that in this system the type of genomic target that produces the most resolved gene trees differs depending on the taxonomic level, although the potential for error is substantially lower for the deep-scale loci at all levels. We successfully recover data for the two different locus sets with high efficiency. Using hierarchical data targeting deep and shallow levels: we 1) resolve the phylogeny of the genus Pseudacris and introduce a novel visual and hypothesis testing method that uses nodal heat maps to examine the robustness of branch support values to the removal of sites and loci; 2) estimate the phylogeographic history of Pseudacris feriarum, which reveals up to five independent invasions leading to sympatry with congener Pseudacris nigrita to form replicated reinforcement contact zones with ongoing gene flow into sympatry; and 3) quantify with high confidence the frequency of hybridization in one of these zones between P. feriarum and P. nigrita, which is lower than microsatellite-based estimates. We find that the hierarchical hybrid enrichment approach offers an efficient, multitiered data collection method for simultaneously addressing questions spanning multiple evolutionary scales. [Anchored hybrid enrichment; heat map; hybridization; phylogenetics; phylogeography; population genomics; reinforcement; reproductive character displacement.].

A likelihood inference of historical biogeography in the world's most diverse terrestrial vertebrate genus: diversification of direct-developing frogs (Craugastoridae: Pristimantis) across the Neotropics.

The geology of the northern Andean region has driven the evolutionary history of Neotropical fauna through the creation of barriers and connections that have resulted in speciation and dispersal events, respectively. One of the most conspicuous groups of anuran fauna in the Andes and surrounding areas is the direct-developing species of the genus Pristimantis. We investigated the molecular phylogenetic placement of 12 species from the montane Andes of Colombia in a broader geographical context with a new genus-level phylogeny in order to identify the role of Andean orogeny over the last 40million years and the effect of elevational differences in diversification of Pristimantis. We examined the biogeographic history of the genus using ancestral range reconstruction by biogeographic regions and elevational ranges. We recognized the middle elevational band (between 1000 and 3000m) in the Northwestern Andes region of Colombia and Ecuador as a focal point for the origin and radiation of Pristimantis species. Additionally, we found several Andean migrations toward new habitats in Central Andes and Merida Andes for some species groups. We suggest that the paleogeological changes in the Northwestern Andes were the main promoter of speciation in Pristimantis, and may have served as a corridor for the dispersion of lowland species.

Differential gene expression analysis of spatial transcriptomic experiments using spatial mixed models.

Spatial transcriptomics (ST) assays represent a revolution in how the architecture of tissues is studied by allowing for the exploration of cells in their spatial context. A common element in the analysis is delineating tissue domains or "niches" followed by detecting differentially expressed genes to infer the biological identity of the tissue domains or cell types. However, many studies approach differential expression analysis by using statistical approaches often applied in the analysis of non-spatial scRNA data (e.g., two-sample t-tests, Wilcoxon's rank sum test), hence neglecting the spatial dependency observed in ST data. In this study, we show that applying linear mixed models with spatial correlation structures using spatial random effects effectively accounts for the spatial autocorrelation and reduces inflation of type-I error rate observed in non-spatial based differential expression testing. We also show that spatial linear models with an exponential correlation structure provide a better fit to the ST data as compared to non-spatial models, particularly for spatially resolved technologies that quantify expression at finer scales (i.e., single-cell resolution).

Spatial transcriptomics analysis identifies a tumor-promoting function of the meningeal stroma in melanoma leptomeningeal disease.

Leptomeningeal disease (LMD) remains a rapidly lethal complication for late-stage melanoma patients. Here, we characterize the tumor microenvironment of LMD and patient-matched extra-cranial metastases using spatial transcriptomics in a small number of clinical specimens (nine tissues from two patients) with extensive in vitro and in vivo validation. The spatial landscape of melanoma LMD is characterized by a lack of immune infiltration and instead exhibits a higher level of stromal involvement. The tumor-stroma interactions at the leptomeninges activate tumor-promoting signaling, mediated through upregulation of SERPINA3. The meningeal stroma is required for melanoma cells to survive in the cerebrospinal fluid (CSF) and promotes MAPK inhibitor resistance. Knocking down SERPINA3 or inhibiting the downstream IGR1R/PI3K/AKT axis results in tumor cell death and re-sensitization to MAPK-targeting therapy. Our data provide a spatial atlas of melanoma LMD, identify the tumor-promoting role of meningeal stroma, and demonstrate a mechanism for overcoming microenvironment-mediated drug resistance in LMD.

WITHDRAWN: Impact of spatial clustering of cytotoxic and tumor infiltrating lymphocytes on overall survival in women with high grade serous ovarian cancer.

The authors have withdrawn their manuscript owing to incorrect handling of multiple measures in the survival analyses. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author.

Spatial transcriptomics analysis identifies a unique tumor-promoting function of the meningeal stroma in melanoma leptomeningeal disease.

Leptomeningeal disease (LMD) remains a rapidly lethal complication for late-stage melanoma patients. The inaccessible nature of the disease site and lack of understanding of the biology of this unique metastatic site are major barriers to developing efficacious therapies for patients with melanoma LMD. Here, we characterize the tumor microenvironment of the leptomeningeal tissues and patient-matched extra-cranial metastatic sites using spatial transcriptomic analyses with in vitro and in vivo validation. We show the spatial landscape of melanoma LMD to be characterized by a lack of immune infiltration and instead exhibit a higher level of stromal involvement. We show that the tumor-stroma interactions at the leptomeninges activate pathways implicated in tumor-promoting signaling, mediated through upregulation of SERPINA3 at the tumor-stroma interface. Our functional experiments establish that the meningeal stroma is required for melanoma cells to survive in the CSF environment and that these interactions lead to a lack of MAPK inhibitor sensitivity in the tumor. We show that knocking down SERPINA3 or inhibiting the downstream IGR1R/PI3K/AKT axis results in re-sensitization of the tumor to MAPK-targeting therapy and tumor cell death in the leptomeningeal environment. Our data provides a spatial atlas of melanoma LMD, identifies the tumor-promoting role of meningeal stroma, and demonstrates a mechanism for overcoming microenvironment-mediated drug resistance unique to this metastatic site.

Branched-chain keto acids promote an immune-suppressive and neurodegenerative microenvironment in leptomeningeal disease.

Leptomeningeal disease (LMD) occurs when tumors seed into the leptomeningeal space and cerebrospinal fluid (CSF), leading to severe neurological deterioration and poor survival outcomes. We utilized comprehensive multi-omics analyses of CSF from patients with lymphoma LMD to demonstrate an immunosuppressive cellular microenvironment and identified dysregulations in proteins and lipids indicating neurodegenerative processes. Strikingly, we found a significant accumulation of toxic branched-chain keto acids (BCKA) in the CSF of patients with LMD. The BCKA accumulation was found to be a pan-cancer occurrence, evident in lymphoma, breast cancer, and melanoma LMD patients. Functionally, BCKA disrupted the viability and function of endogenous T lymphocytes, chimeric antigen receptor (CAR) T cells, neurons, and meningeal cells. Treatment of LMD mice with BCKA-reducing sodium phenylbutyrate significantly improved neurological function, survival outcomes, and efficacy of anti-CD19 CAR T cell therapy. This is the first report of BCKA accumulation in LMD and provides preclinical evidence that targeting these toxic metabolites improves outcomes.

A biographical account of John Paul Richard Thomas, the man who leaves no stone unturned.

This biographical account summarizes the professional career and scientific contributions of John Paul Richard Thomas, a contemporary leading figure in the systematics of West Indian amphibians and non-avian reptiles, especially of blind snakes of the families Typhlopidae and Leptotyphlopidae. Since his first expedition to the West Indies in 1957, Richard's vast field experience (including three trips to Peru between 1968 and 1974), impressive collecting skills, and remarkable ability to detect phenotypic variation among natural populations have resulted in the description of more than 70 species of snakes (24 typhlopids, 4 leptotyphlopids), lizards, and frogs in 16 genera and 11 taxonomic families. Richard joined the faculty of the Department of Biology, University of Puerto Rico, Río Piedras, in 1976 and ever since his efforts significantly advanced organismal biology research at the institution. Although primarily a systematist, his desire to understand multiple aspects of an organism's biology and contagious passion for becoming intimately familiar with animals in their natural environments provided his students the opportunity to conduct research in fields such as behavioral and evolutionary ecology. Richard's mentoring fostered the scientific interests of his graduate students, who were exposed first-hand to every aspect of research, an invaluable experience that served as a springboard for the development of their professional careers inside and outside academia. This Commentary is a fitting tribute to an influential, unassuming scientist whose passion for turning over rocks has led to the discovery of many interesting species.

The herpetological contributions of Richard Thomas.

John Paul Richard Thomas is among the living herpetologists to have described the greatest number of new species of amphibians and reptiles, and his contributions to the herpetology of the West Indies, particularly the Greater Antilles, have been exceptional. His academic career followed an unusual path, having established a strong reputation and described 50 new taxa prior to beginning his doctoral studies. His career was strongly influenced by Albert Schwartz and later was characterized by extensive and fruitful collaboration with S. Blair Hedges. Thomas' contributions to the study of blind snakes have been noteworthy. In addition to describing 28 species of scolecophidians he has been a keen observer of blind snake morphology and his 1976 dissertation remains a valuable source of osteological data. We outline some of the highlights of the career of Richard Thomas and provide a bibliography of his scientific works and a listing of the 108 taxa of amphibians and reptiles described by him.

The status of the anomalepidid snake Liotyphlops albirostris and the revalidation of three taxa based on morphology and ecological niche models.

Liotyphlops is a genus of blindsnakes distributed in Central and South America. We reviewed specimens of Liotyphlops albirostris along its current distribution range and, based on morphological data and ecological niche modeling analyses, we restrict the geographical range of L. albirostris and validate three previously described species. In this revision, we describe the morphological variation in the populations from Panamá, Colombia, Ecuador, and Venezuela, and propose a new taxonomic arrangement. We revalidate three previous synonyms of L. albirostris to full species status, while dividing the populations from Colombia in two subspecies-one attributed to a previously recognized species from the Caribbean region, and a new one from the Andean region. The new species differs from L. albirostris from Panamá in cephalic scale arrangements that effectively reduces the previously reported variability of these scales in L. albirostris. We also explore some osteological differences that are congruent with the variation observed. We hope that the recognition of these new species better represents the diversity within Liotyphlops, helping to bring these new species out of their cryptic status so that they will be considered in future conservation efforts.

Challenges and Opportunities in the Statistical Analysis of Multiplex Immunofluorescence Data.

Immune modulation is considered a hallmark of cancer initiation and progression. The recent development of immunotherapies has ushered in a new era of cancer treatment. These therapeutics have led to revolutionary breakthroughs; however, the efficacy of immunotherapy has been modest and is often restricted to a subset of patients. Hence, identification of which cancer patients will benefit from immunotherapy is essential. Multiplex immunofluorescence (mIF) microscopy allows for the assessment and visualization of the tumor immune microenvironment (TIME). The data output following image and machine learning analyses for cell segmenting and phenotyping consists of the following information for each tumor sample: the number of positive cells for each marker and phenotype(s) of interest, number of total cells, percent of positive cells for each marker, and spatial locations for all measured cells. There are many challenges in the analysis of mIF data, including many tissue samples with zero positive cells or "zero-inflated" data, repeated measurements from multiple TMA cores or tissue slides per subject, and spatial analyses to determine the level of clustering and co-localization between the cell types in the TIME. In this review paper, we will discuss the challenges in the statistical analysis of mIF data and opportunities for further research.

Reduced Diversity in the Bacteriome of the Phytophagous Mite Brevipalpus yothersi (Acari: Tenuipalpidae).

Tenuipalpidae comprises mites that transmit viruses to agriculturally important plants. Several tenuipalpid species present parthenogenesis, and in Brevipalpus yothersi, the endosymbiont Cardinium has been associated with female-only colonies. It is unclear what the bacterial composition of B. yothersi is, and how common Cardinium is in those microbiomes. We performed a comparative analysis of the bacteriomes in three populations of B. yothersi and three additional Tetranychoidea species using sequences from V4-fragment of 16S DNA. The bacteriomes were dominated by Bacteroidetes (especially Cardinium) and Proteobacteria, showing a remarkably low alpha diversity. Cardinium was present in about 22% of all sequences; however, it was not present in R. indica and T. evansi. In B. yothersi, the proportion of Cardinium was higher in adults than eggs, suggesting that proliferation of the bacteria could be the result of selective pressures from the host. This hypothesis was further supported because colonies of B. yothersi from different populations showed different bacterial assemblages, and bacteriomes from different mite species showed similar abundances of Cardinium. A phylogenetic analysis of Cardinium revealed that not only specialization but horizontal transmission has been important for this symbiosis. Together, these results represent a glimpse into the evolution of the Tetranychoidea and Cardinium.

Elevated levels of adaption in Helicobacter pylori genomes from Japan; a link to higher incidences of gastric cancer?

Helicobacter pylori is a bacterium that lives in the human stomach and is a major risk factor for gastric cancer and ulcers. H.pylori is host dependent and has been carried with human populations around the world after their departure from Africa. We wished to investigate how H.pylori has coevolved with its host during that time, focusing on strains from Japanese and European populations, given that gastric cancer incidence is high in Japanese populations, while low in European. A positive selection analysis of eight H.pylori genomes was conducted, using maximum likelihood based pairwise comparisons in order to maximize the number of strain-specific genes included in the study. Using the genic Ka/Ks ratio, comparisons of four Japanese H.pylori genomes suggests 25-34 genes under positive selection, while four European H.pylori genomes suggests 16-21 genes; few of the genes identified were in common between lineages. Of the identified genes which were annotated, 38% possessed homologs associated with pathogenicity and / or host adaptation, consistent with their involvement in a coevolutionary 'arms race' with the host. Given the efficacy of identifying host interaction factors de novo, in the absence of functionally annotated homologs our evolutionary approach may have value in identifying novel genes which H.pylori employs to interact with the human gut environment. In addition, the larger number of genes inferred as being under positive selection in Japanese strains compared to European implies a stronger overall adaptive pressure, potentially resulting from an elevated immune response which may be linked to increased inflammation, an initial stage in the development of gastric cancer.