COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63+ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

PSMA-targeting TGFβ-insensitive armored CAR T cells in metastatic castration-resistant prostate cancer: a phase 1 trial

Chimeric antigen receptor (CAR) T cells have demonstrated promising efficacy, particularly in hematologic malignancies. One challenge regarding CAR T cells in solid tumors is the immunosuppressive tumor microenvironment (TME), characterized by high levels of multiple inhibitory factors, including transforming growth factor (TGF)-β. We report results from an in-human phase 1 trial of castration-resistant, prostate cancer-directed CAR T cells armored with a dominant-negative TGF-β receptor (NCT03089203). Primary endpoints were safety and feasibility, while secondary objectives included assessment of CAR T cell distribution, bioactivity and disease response. All prespecified endpoints were met. Eighteen patients enrolled, and 13 subjects received therapy across four dose levels. Five of the 13 patients developed grade ≥2 cytokine release syndrome (CRS), including one patient who experienced a marked clonal CAR T cell expansion, >98% reduction in prostate-specific antigen (PSA) and death following grade 4 CRS with concurrent sepsis. Acute increases in inflammatory cytokines correlated with manageable high-grade CRS events. Three additional patients achieved a PSA reduction of ≥30%, with CAR T cell failure accompanied by upregulation of multiple TME-localized inhibitory molecules following adoptive cell transfer. CAR T cell kinetics revealed expansion in blood and tumor trafficking. Thus, clinical application of TGF-β-resistant CAR T cells is feasible and generally safe. Future studies should use superior multipronged approaches against the TME to improve outcomes.

Visualizing in deceased COVID-19 patients how SARS-CoV-2 attacks the respiratory and olfactory mucosae but spares the olfactory bulb

Anosmia, the loss of smell, is a common and often the sole symptom of COVID-19. The onset of the sequence of pathobiological events leading to olfactory dysfunction remains obscure. Here, we have developed a postmortem bedside surgical procedure to harvest endoscopically samples of respiratory and olfactory mucosae and whole olfactory bulbs. Our cohort of 85 cases included COVID-19 patients who died a few days after infection with SARS-CoV-2, enabling us to catch the virus while it was still replicating. We found that sustentacular cells are the major target cell type in the olfactory mucosa. We failed to find evidence for infection of olfactory sensory neurons, and the parenchyma of the olfactory bulb is spared as well. Thus, SARS-CoV-2 does not appear to be a neurotropic virus. We postulate that transient insufficient support from sustentacular cells triggers transient olfactory dysfunction in COVID-19. Olfactory sensory neurons would become affected without getting infected.

The spatial landscape of lung pathology during COVID-19 progression

Recent studies have provided insights into the pathology of and immune response to COVID-191-8. However, a thorough investigation of the interplay between infected cells and the immune system at sites of infection has been lacking. Here we use high-parameter imaging mass cytometry9 that targets the expression of 36 proteins to investigate the cellular composition and spatial architecture of acute lung injury in humans (including injuries derived from SARS-CoV-2 infection) at single-cell resolution. These spatially resolved single-cell data unravel the disordered structure of the infected and injured lung, alongside the distribution of extensive immune infiltration. Neutrophil and macrophage infiltration are hallmarks of bacterial pneumonia and COVID-19, respectively. We provide evidence that SARS-CoV-2 infects predominantly alveolar epithelial cells and induces a localized hyperinflammatory cell state that is associated with lung damage. We leverage the temporal range of fatal outcomes of COVID-19 in relation to the onset of symptoms, which reveals increased macrophage extravasation and increased numbers of mesenchymal cells and fibroblasts concomitant with increased proximity between these cell types as the disease progresses-possibly as a result of attempts to repair the damaged lung tissue. Our data enable us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. We use this landscape to characterize the pathophysiology of the human lung from its macroscopic presentation to the single-cell level, which provides an important basis for understanding COVID-19 and lung pathology in general.

Goldilocks Meets Santa Rosalia: An Ephemeral Speciation Model Explains Patterns of Diversification Across Time Scales

Understanding the rate at which new species form is a key question in studying the evolution of life on earth. Here we review our current understanding of speciation rates, focusing on studies based on the fossil record, phylogenies, and mathematical models. We find that speciation rates estimated from these different studies can be dramatically different: some studies find that new species form quickly and often, while others find that new species form much less frequently. We suggest that instead of being contradictory, differences in speciation rates across different scales can be reconciled by a common model. Under the "ephemeral speciation model", speciation is very common and very rapid but the new species produced almost never persist. Evolutionary studies should therefore focus on not only the formation but also the persistence of new species.

The CD155/TIGIT axis promotes and maintains immune evasion in neoantigen-expressing pancreatic cancer

The CD155/TIGIT axis can be co-opted during immune evasion in chronic viral infections and cancer. Pancreatic adenocarcinoma (PDAC) is a highly lethal malignancy, and immune-based strategies to combat this disease have been largely unsuccessful to date. We corroborate prior reports that a substantial portion of PDAC harbors predicted high-affinity MHC class I-restricted neoepitopes and extend these findings to advanced/metastatic disease. Using multiple preclinical models of neoantigen-expressing PDAC, we demonstrate that intratumoral neoantigen-specific CD8+ T cells adopt multiple states of dysfunction, resembling those in tumor-infiltrating lymphocytes of PDAC patients. Mechanistically, genetic and/or pharmacologic modulation of the CD155/TIGIT axis was sufficient to promote immune evasion in autochthonous neoantigen-expressing PDAC. Finally, we demonstrate that the CD155/TIGIT axis is critical in maintaining immune evasion in PDAC and uncover a combination immunotherapy (TIGIT/PD-1 co-blockade plus CD40 agonism) that elicits profound anti-tumor responses in preclinical models, now poised for clinical evaluation.

Shotgun transcriptome, spatial omics, and isothermal profiling of SARS-CoV-2 infection reveals unique host responses, viral diversification, and drug interactions

In less than nine months, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) killed over a million people, including >25,000 in New York City (NYC) alone. The COVID-19 pandemic caused by SARS-CoV-2 highlights clinical needs to detect infection, track strain evolution, and identify biomarkers of disease course. To address these challenges, we designed a fast (30-minute) colorimetric test (LAMP) for SARS-CoV-2 infection from naso/oropharyngeal swabs and a large-scale shotgun metatranscriptomics platform (total-RNA-seq) for host, viral, and microbial profiling. We applied these methods to clinical specimens gathered from 669 patients in New York City during the first two months of the outbreak, yielding a broad molecular portrait of the emerging COVID-19 disease. We find significant enrichment of a NYC-distinctive clade of the virus (20C), as well as host responses in interferon, ACE, hematological, and olfaction pathways. In addition, we use 50,821 patient records to find that renin-angiotensin-aldosterone system inhibitors have a protective effect for severe COVID-19 outcomes, unlike similar drugs. Finally, spatial transcriptomic data from COVID-19 patient autopsy tissues reveal distinct ACE2 expression loci, with macrophage and neutrophil infiltration in the lungs. These findings can inform public health and may help develop and drive SARS-CoV-2 diagnostic, prevention, and treatment strategies.

SARS-CoV-2 infection produces chronic pulmonary epithelial and immune cell dysfunction with fibrosis in mice

A subset of individuals who recover from coronavirus disease 2019 (COVID-19) develop post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal tissue samples. The mouse-adapted SARS-CoV-2 strain MA10 produces an acute respiratory distress syndrome in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute to clinical recovery phases. At 15 to 120 days after virus clearance, pulmonary histologic findings included subpleural lesions composed of collagen, proliferative fibroblasts, and chronic inflammation, including tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal up-regulation of profibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early antifibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

Emergent EEG: indications and diagnostic yield

The authors reviewed the reports of all emergent EEG (EmEEG) performed in our hospital within 1 hour of the test being ordered over a period of 52 months. Two hundred sixty-one EmEEG (12.8% of all EEG) were performed. The most common reason to order the test was a change in mental status or coma (17.6%). Although EmEEG was ordered to rule out status epilepticus (SE) in 60.2% of cases, this diagnosis was made in only 10.7% of patients. The only independent predictor for SE was a history of cardiac or respiratory arrest (odds [95% CI], 6.8 [2.7 to 16.9]).

Genomics and introgression: Discovery and mapping of thousands of species-diagnostic SNPs using RAD sequencing

Invasive hybridization and introgression pose a serious threat to the persistence of many native species. Understanding the effects of hybridization on native populations (e.g., fitness consequences) requires numerous species-diagnostic loci distributed genome-wide. Here we used RAD sequencing to discover thousands of single-nucleotide polymorphisms (SNPs) that are diagnostic between rainbow trout (RBT, Oncorhynchus mykiss), the world’s most widely introduced fish, and native westslope cutthroat trout (WCT, O. clarkii lewisi) in the northern Rocky Mountains, USA. We advanced previous work that identified 4,914 species-diagnostic loci by using longer sequence reads (100 bp vs. 60 bp) and a larger set of individuals (n = 84). We sequenced RAD libraries for individuals from diverse sampling sources, including native populations of WCT and hatchery broodstocks of WCT and RBT. We also took advantage of a newly released reference genome assembly for RBT to align our RAD loci. In total, we discovered 16,788 putatively diagnostic SNPs, 10,267 of which we mapped to anchored chromosome locations on the RBT genome. A small portion of previously discovered putative diagnostic loci (325 of 4,914) were no longer diagnostic (i.e., fixed between species) based on our wider survey of non-hybridized RBT and WCT individuals. Our study suggests that RAD loci mapped to a draft genome assembly could provide the marker density required to identify genes and chromosomal regions influencing selection in admixed populations of conservation concern and evolutionary interest.

System-wide transcriptome damage and tissue identity loss in COVID-19 patients

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Identifying and Tracking Low-Frequency Virus-Specific TCR Clonotypes Using High-Throughput Sequencing

Tracking antigen-specific T cell responses over time within individuals is difficult because of lack of knowledge of antigen-specific TCR sequences, limitations in sample size, and assay sensitivities. We hypothesized that analyses of high-throughput sequencing of TCR clonotypes could provide functional readouts of individuals' immunological histories. Using high-throughput TCR sequencing, we develop a database of TCRβ sequences from large cohorts of mice before (naive) and after smallpox vaccination. We computationally identify 315 vaccine-associated TCR sequences (VATS) that are used to train a diagnostic classifier that distinguishes naive from vaccinated samples in mice up to 9 months post-vaccination with >99% accuracy. We determine that the VATS library contains virus-responsive TCRs by in vitro expansion assays and virus-specific tetramer sorting. These data outline a platform for advancing our capabilities to identify pathogen-specific TCR sequences, which can be used to identify and quantitate low-frequency pathogen-specific TCR sequences in circulation over time with exceptional sensitivity.

Machine learning identifies specific habitats associated with genetic connectivity in Hyla squirella

The goal of this study was to identify and differentiate the influence of multiple habitat types that span a spectrum of suitability for Hyla squirella, a widespread frog species that occurs in a broad range of habitat types. We collected microsatellite data from 675 samples representing 20 localities from the southeastern USA and used machine-learning methodologies to identify significant habitat features associated with genetic structure. In simulation, we confirm that our machine-learning algorithm can successfully identify landscape features responsible for generating between-population genetic differentiation, suggesting that it can be a useful hypothesis-generating tool for landscape genetics. In our study system, we found that H. squirella were spatially structured and models including specific habitat types (i.e. upland oak forest and urbanization) consistently explained more variation in genetic distance (median pR(2)  = 47.78) than spatial distance alone (median pR(2)  = 23.81). Moreover, we estimate the relative importance that spatial distance, upland oak and urbanized habitat have in explaining genetic structure of H. squirella. We discuss how these habitat types may mechanistically facilitate dispersal in H. squirella. This study provides empirical support for the hypothesis that habitat-use can be an informative correlate of genetic differentiation, even for species that occur in a wide range of habitats.

The tumor immune microenvironment of primary and metastatic HER2- positive breast cancers utilizing gene expression and spatial proteomic profiling

BACKGROUND

The characterization of the immune component of the tumor microenvironment (TME) of human epidermal growth factor receptor 2 positive (HER2+) breast cancer has been limited. Molecular and spatial characterization of HER2+ TME of primary, recurrent, and metastatic breast tumors has the potential to identify immune mediated mechanisms and biomarker targets that could be used to guide selection of therapies.

METHODS

We examined 15 specimens from eight patients with HER2+ breast cancer: 10 primary breast tumors (PBT), two soft tissue, one lung, and two brain metastases (BM). Using molecular profiling by bulk gene expression TME signatures, including the Tumor Inflammation Signature (TIS) and PAM50 subtyping, as well as spatial characterization of immune hot, warm, and cold regions in the stroma and tumor epithelium using 64 protein targets on the GeoMx Digital Spatial Profiler.

RESULTS

PBT had higher infiltration of immune cells relative to metastatic sites and higher protein and gene expression of immune activation markers when compared to metastatic sites. TIS scores were lower in metastases, particularly in BM. BM also had less immune infiltration overall, but in the stromal compartment with the highest density of immune infiltration had similar levels of T cells that were less activated than PBT stromal regions suggesting immune exclusion in the tumor epithelium.

CONCLUSIONS

Our findings show stromal and tumor localized immune cells in the TME are more active in primary versus metastatic disease. This suggests patients with early HER2+ breast cancer could have more benefit from immune-targeting therapies than patients with advanced disease.

Cryptic chytridiomycosis linked to climate and genetic variation in amphibian populations of the southeastern United States

North American amphibians have recently been impacted by two major emerging pathogens, the fungus Batrachochytrium dendrobatidis (Bd) and iridoviruses in the genus Ranavirus (Rv). Environmental factors and host genetics may play important roles in disease dynamics, but few studies incorporate both of these components into their analyses.

Here, we investigated the role of environmental and genetic factors in driving Bd and Rv infection prevalence and severity in a biodiversity hot spot, the southeastern United States.

We used quantitative PCR to characterize Bd and Rv dynamics in natural populations of three amphibian species: Notophthalmus perstriatus, Hyla squirella and Pseudacris ornata. We combined pathogen data, genetic diversity metrics generated from neutral markers, and environmental variables into general linear models to evaluate how these factors impact infectious disease dynamics. Occurrence, prevalence and intensity of Bd and Rv varied across species and populations, but only one species, Pseudacris ornata, harbored high Bd intensities in the majority of sampled populations. Genetic diversity and climate variables both predicted Bd prevalence, whereas climatic variables alone predicted infection intensity.

We conclude that Bd is more abundant in the southeastern United States than previously thought and that genetic and environmental factors are both important for predicting amphibian pathogen dynamics. Incorporating both genetic and environmental information into conservation plans for amphibians is necessary for the development of more effective management strategies to mitigate the impact of emerging infectious diseases.

Genetic regulatory network motifs constrain adaptation through curvature in the landscape of mutational (co)variance

Systems biology is accumulating a wealth of understanding about the structure of genetic regulatory networks, leading to a more complete picture of the complex genotype-phenotype relationship. However, models of multivariate phenotypic evolution based on quantitative genetics have largely not incorporated a network-based view of genetic variation. Here we model a set of two-node, two-phenotype genetic network motifs, covering a full range of regulatory interactions. We find that network interactions result in different patterns of mutational (co)variance at the phenotypic level (the M-matrix), not only across network motifs but also across phenotypic space within single motifs. This effect is due almost entirely to mutational input of additive genetic (co)variance. Variation in M has the effect of stretching and bending phenotypic space with respect to evolvability, analogous to the curvature of space-time under general relativity, and similar mathematical tools may apply in each case. We explored the consequences of curvature in mutational variation by simulating adaptation under divergent selection with gene flow. Both standing genetic variation (the G-matrix) and rate of adaptation are constrained by M, so that G and adaptive trajectories are curved across phenotypic space. Under weak selection the phenotypic mean at migration-selection balance also depends on M.

Systemic Tissue and Cellular Disruption from SARS-CoV-2 Infection revealed in COVID-19 Autopsies and Spatial Omics Tissue Maps

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus has infected over 115 million people and caused over 2.5 million deaths worldwide. Yet, the molecular mechanisms underlying the clinical manifestations of COVID-19, as well as what distinguishes them from common seasonal influenza virus and other lung injury states such as Acute Respiratory Distress Syndrome (ARDS), remains poorly understood. To address these challenges, we combined transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues, matched with spatial protein and expression profiling (GeoMx) across 357 tissue sections. These results define both body-wide and tissue-specific (heart, liver, lung, kidney, and lymph nodes) damage wrought by the SARS-CoV-2 infection, evident as a function of varying viral load (high vs. low) during the course of infection and specific, transcriptional dysregulation in splicing isoforms, T cell receptor expression, and cellular expression states. In particular, cardiac and lung tissues revealed the largest degree of splicing isoform switching and cell expression state loss. Overall, these findings reveal a systemic disruption of cellular and transcriptional pathways from COVID-19 across all tissues, which can inform subsequent studies to combat the mortality of COVID-19, as well to better understand the molecular dynamics of lethal SARS-CoV-2 infection and other viruses.

A model of persistent post SARS-CoV-2 induced lung disease for target identification and testing of therapeutic strategies

COVID-19 survivors develop post-acute sequelae of SARS-CoV-2 (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal samples. Mouse-adapted SARS-CoV-2 MA10 produces an acute respiratory distress syndrome (ARDS) in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute disease through clinical recovery. At 15-120 days post-virus clearance, histologic evaluation identified subpleural lesions containing collagen, proliferative fibroblasts, and chronic inflammation with tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal upregulation of pro-fibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early anti-fibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Fat frogs, mobile genes: unexpected phylogeographic patterns for the ornate chorus frog (Pseudacris ornata)

The southeastern coastal plain of the United States is a region marked by extraordinary phylogeographic congruence that is frequently attributed to the changing sea levels that occurred during the glacial-interglacial cycles of the Pleistocene epoch. A phylogeographic break corresponding to the Apalachicola River has been suggested for many species studied to date that are endemic to this region. Here, we used this pattern of phylogeographic congruence to develop and test explicit hypotheses about the genetic structure in the ornate chorus frog (Pseudacris ornata). Using 1299 bp of mtDNA sequence and seven nuclear microsatellite markers in 13 natural populations of P. ornata, we found three clades corresponding to geographically distinct regions; one spans the Apalachicola River (Southern Clade), one encompasses Georgia and South Carolina (Central Clade) and a third comprises more northerly individuals (Northern Clade). However, it does not appear that typical phylogeographic barriers demarcate these clades. Instead, isolation by distance across the range of the entire species explained the pattern of genetic variation that we observed. We propose that P. ornata was historically widespread in the southeastern United States, and that a balance between genetic drift and migration was the root of the genetic divergence among populations. Additionally, we investigated fine-scale patterns of genetic structure and found the spatial scale at which there was significant genetic structure varied among the regions studied. Furthermore, we discuss our results in light of other phylogeographic studies of southeastern coastal plain organisms and in relation to amphibian conservation and management.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

BACKGROUND

The molecular underpinnings of organ dysfunction in severe COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we perform single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents.

RESULTS

We identify hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells, and a central role in a pro-fibrotic TGFβ signaling cell-cell communications network. Integrated analysis and comparisons with healthy controls reveal extensive changes in the cellular composition and expression states in COVID-19 liver, providing the underpinning of hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis characteristic of COVID-19 cholangiopathy. We also observe Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition is dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells.

CONCLUSIONS

Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Abstract 66: A multi-institution examination of concordance in spatial transcriptomics using the GeoMx Cancer Transcriptome Atlas

Multiplexed spatial profiling can enable biological insights by characterizing gene expression within discrete physical locations of a tissue. However, these advanced techniques can be confounded by variability of sample collection, storage, or profiling protocols, making it difficult to accurately compare data generated by different laboratories. Further, as more spatial profiling datasets become publicly available, having methods to enable meta-analysis of samples collected on different studies will maximize the learnings to support the advancement of treatments or identifying patient segments.

To quantify the variability of spatial profiling data at different laboratories and to advance data normalization methodologies, 4 independent laboratories used the NanoString® GeoMx® Digital Spatial Profiling (DSP) to profile serial 5 µm sections of tissue and cell pellet arrays (CPAs). GeoMx DSP enables high throughput, spatially resolved analysis of gene or protein expression from fresh or archival human tissues. In this study, the GeoMx Cancer Transcriptome Atlas was used to profile >1800 genes simultaneously. We examined the concordance of GeoMx data generated in the different laboratories when controlling for methodical variation (e.g., reagents, tissue source) and experimentally varying region of interest (ROI) size, collection site, and sample preservation methods.

Sections of tonsil, colon, and 2 CPAs were profiled separately at the 4 laboratories. Each analyzed fresh cut (FC) tissues and two sites examined sample stability by analyzing the impact of storing slides at -80°C for 1 month prior to spatial profiling. Concordance analysis was performed using the Horn-Morisita Index on raw data comparing paired and unpaired ROIs across each set of slides.

In CPA samples where each pellet was a different tumor type (e.g., NSCLC, melanoma), we observed strong clustering by cell line. While data initially showed varying degrees of clustering by slide, factoring out this variable removed the association of slide, allowing integration of the data across profiling locations without affecting concordance within slides. In tonsil, ROIs with increasing area were profiled. Comparing expression between pairs of samples for a given area, concordance increased with ROI size (R = -0.40, p<6e-06). Finally, we observed little impact of preservation method (FC vs -80°C) in these data.

In this study, we quantify slide-specific variation observed in high-plex RNA profiling by the DSP platform and detail methods for accounting for this variation. We note that many downstream analyses (e.g., differential expression) already model slide effects during the analysis, but modeling it explicitly allows for direct comparison of concordance with other approaches (e.g., clustering, PCA). These methods support the use of multi-institution studies leveraging the GeoMx platform.

Citation Format: Tyler Hether, Tim Howes, David Scoville, Charlie Glaser, Yanyun Li, Rami Vanguri, Neeman Mohibullah, Wan-Jung Chang, Todd Yoder, Minnal Gupta, Kathy Ton, Yan Liang, Ying Huang, Zach Herbert, Jason Reeves, Elizabeth Mittendorf, Simon Lacey, Travis Hollmann, Sarah Warren, Theresa LaVallee. A multi-institution examination of concordance in spatial transcriptomics using the GeoMx Cancer Transcriptome Atlas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 66.