Photothermal Prussian blue nanoparticles generate potent multi-targeted tumor-specific T cells as an adoptive cell therapy.

Prussian blue nanoparticle-based photothermal therapy (PBNP-PTT) is an effective tumor treatment capable of eliciting an antitumor immune response. Motivated by the ability of PBNP-PTT to potentiate endogenous immune responses, we recently demonstrated that PBNP-PTT could be used ex vivo to generate tumor-specific T cells against glioblastoma (GBM) cell lines as an adoptive T cell therapy (ATCT). In this study, we further developed this promising T cell development platform. First, we assessed the phenotype and function of T cells generated using PBNP-PTT. We observed that PBNP-PTT facilitated CD8+ T cell expansion from healthy donor PBMCs that secreted IFNγ and TNFα and upregulated CD107a in response to engagement with target U87 cells, suggesting specific antitumor T cell activation and degranulation. Further, CD8+ effector and effector memory T cell populations significantly expanded after co-culture with U87 cells, consistent with tumor-specific effector responses. In orthotopically implanted U87 GBM tumors in vivo, PBNP-PTT-derived T cells effectively reduced U87 tumor growth and generated long-term survival in >80% of tumor-bearing mice by Day 100, compared to 0% of mice treated with PBS, non-specific T cells, or T cells expanded from lysed U87 cells, demonstrating an enhanced antitumor efficacy of this ATCT platform. Finally, we tested the generalizability of our approach by generating T cells targeting medulloblastoma (D556), breast cancer (MDA-MB-231), neuroblastoma (SH-SY5Y), and acute monocytic leukemia (THP-1) cell lines. The resulting T cells secreted IFNγ and exerted increased tumor-specific cytolytic function relative to controls, demonstrating the versatility of PBNP-PTT in generating tumor-specific T cells for ATCT.

Putative climate adaptation in American pikas (Ochotona princeps) is associated with copy number variation across environmental gradients.

Improved understanding of the genetic basis of adaptation to climate change is necessary for maintaining global biodiversity moving forward. Studies to date have largely focused on sequence variation, yet there is growing evidence that suggests that changes in genome structure may be an even more significant source of adaptive potential. The American pika (Ochotona princeps) is an alpine specialist that shows some evidence of adaptation to climate along elevational gradients, but previous work has been limited to single nucleotide polymorphism based analyses within a fraction of the species range. Here, we investigated the role of copy number variation underlying patterns of local adaptation in the American pika using genome-wide data previously collected across the entire species range. We identified 37-193 putative copy number variants (CNVs) associated with environmental variation (temperature, precipitation, solar radiation) within each of the six major American pika lineages, with patterns of divergence largely following elevational and latitudinal gradients. Genes associated (n = 158) with independent annotations across lineages, variables, and/or CNVs had functions related to mitochondrial structure/function, immune response, hypoxia, olfaction, and DNA repair. Some of these genes have been previously linked to putative high elevation and/or climate adaptation in other species, suggesting they may serve as important targets in future studies.

Phylogenomics of American pika (Ochotona princeps) lineage diversification.

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Evaluating genotyping-in-thousands by sequencing as a genetic monitoring tool for a climate sentinel mammal using non-invasive and archival samples.

Genetic tools for wildlife monitoring can provide valuable information on spatiotemporal population trends and connectivity, particularly in systems experiencing rapid environmental change. Multiplexed targeted amplicon sequencing techniques, such as genotyping-in-thousands by sequencing (GT-seq), can provide cost-effective approaches for collecting genetic data from low-quality and quantity DNA samples, making them potentially useful for long-term wildlife monitoring using non-invasive and archival samples. Here, we developed a GT-seq panel as a potential monitoring tool for the American pika (Ochotona princeps) and evaluated its performance when applied to traditional, non-invasive, and archival samples, respectively. Specifically, we optimized a GT-seq panel (307 single nucleotide polymorphisms (SNPs)) that included neutral, sex-associated, and putatively adaptive SNPs using contemporary tissue samples (n = 77) from the Northern Rocky Mountains lineage of American pikas. The panel demonstrated high genotyping success (94.7%), low genotyping error (0.001%), and excellent performance identifying individuals, sex, relatedness, and population structure. We subsequently applied the GT-seq panel to archival tissue (n = 17) and contemporary fecal pellet samples (n = 129) collected within the Canadian Rocky Mountains to evaluate its effectiveness. Although the panel demonstrated high efficacy with archival tissue samples (90.5% genotyping success, 0.0% genotyping error), this was not the case for the fecal pellet samples (79.7% genotyping success, 28.4% genotyping error) likely due to the exceptionally low quality/quantity of recovered DNA using the approaches implemented. Overall, our study reinforced GT-seq as an effective tool using contemporary and archival tissue samples, providing future opportunities for temporal applications using historical specimens. Our results further highlight the need for additional optimization of sample and genetic data collection techniques prior to broader-scale implementation of a non-invasive genetic monitoring tool for American pikas.

Generation of locus coeruleus norepinephrine neurons from human pluripotent stem cells.

Central norepinephrine (NE) neurons, located mainly in the locus coeruleus (LC), are implicated in diverse psychiatric and neurodegenerative diseases and are an emerging target for drug discovery. To facilitate their study, we developed a method to generate 40-60% human LC-NE neurons from human pluripotent stem cells. The approach depends on our identification of ACTIVIN A in regulating LC-NE transcription factors in dorsal rhombomere 1 (r1) progenitors. In vitro generated human LC-NE neurons display extensive axonal arborization; release and uptake NE; and exhibit pacemaker activity, calcium oscillation and chemoreceptor activity in response to CO2. Single-nucleus RNA sequencing (snRNA-seq) analysis at multiple timepoints confirmed NE cell identity and revealed the differentiation trajectory from hindbrain progenitors to NE neurons via an ASCL1-expressing precursor stage. LC-NE neurons engineered with an NE sensor reliably reported extracellular levels of NE. The availability of functional human LC-NE neurons enables investigation of their roles in psychiatric and neurodegenerative diseases and provides a tool for therapeutics development.

Multimodal analysis reveals genes driving neuronal maturation in the primate prefrontal cortex.

The dorsolateral prefrontal cortex (dlPFC) is an evolutionarily derived cortical region in primates critical for high-level cognitive functions and implicated in various neuropsychiatric disorders. The cells that compose the dlPFC, especially excitatory and inhibitory neurons, undergo extensive maturation throughout midfetal and late-fetal development, during which critical neurodevelopmental events, such as circuit assembly and electrophysiological maturation of neurons, occur. Despite the relevance of neuronal maturation in several neurodevelopmental and psychiatric disorders, the molecular mechanisms underlying this process remain largely unknown. Here, we performed an integrated Patch-seq and single-nucleus multiomic analysis of the rhesus macaque dlPFC to identify genes governing neuronal maturation from midfetal to late-fetal development. Our multimodal analysis identified gene pathways and regulatory networks important for the maturation of distinct neuronal populations, including upper-layer intratelencephalicprojecting neurons. We identified genes underlying the maturation of specific electrophysiological properties of these neurons. Furthermore, gene knockdown in organotypic slices revealed that RAPGEF4 regulates the maturation of resting membrane potential and inward sodium current. Using this strategy, we also found that the knockdown of CHD8, a high-confidence autism spectrum disorder risk gene, in human slices led to deficits in neuronal maturation, via the downstream downregulation of several key genes, including RAPGEF4. Our study revealed novel regulators of neuronal maturation during a critical period of prefrontal development in primates and implicated such regulators in molecular processes underlying autism.

Molecular and cellular evolution of the primate dorsolateral prefrontal cortex.

The granular dorsolateral prefrontal cortex (dlPFC) is an evolutionary specialization of primates that is centrally involved in cognition. We assessed more than 600,000 single-nucleus transcriptomes from adult human, chimpanzee, macaque, and marmoset dlPFC. Although most cell subtypes defined transcriptomically are conserved, we detected several that exist only in a subset of species as well as substantial species-specific molecular differences across homologous neuronal, glial, and non-neural subtypes. The latter are exemplified by human-specific switching between expression of the neuropeptide somatostatin and tyrosine hydroxylase, the rate-limiting enzyme in dopamine production in certain interneurons. The above molecular differences are also illustrated by expression of the neuropsychiatric risk gene FOXP2, which is human-specific in microglia and primate-specific in layer 4 granular neurons. We generated a comprehensive survey of the dlPFC cellular repertoire and its shared and divergent features in anthropoid primates.

Novel Injury Scoring Tool for Assessing Brain Injury following Neonatal Hypoxia-Ischemia in Mice.

The variability of severity in hypoxia-ischemia (HI)-induced brain injury among research subjects is a major challenge in developmental brain injury research. Our laboratory developed a novel injury scoring tool based on our gross pathological observations during hippocampal extraction. The hippocampi received scores of 0-6 with 0 being no injury and 6 being severe injury post-HI. The hippocampi exposed to sham surgery were grouped as having no injury. We have validated the injury scoring tool with T2-weighted MRI analysis of percent hippocampal/hemispheric tissue loss and cell survival/death markers after exposing the neonatal mice to Vannucci's rodent model of neonatal HI. In addition, we have isolated hippocampal nuclei and quantified the percent good quality nuclei to provide an example of utilization of our novel injury scoring tool. Our novel injury scores correlated significantly with percent hippocampal and hemispheric tissue loss, cell survival/death markers, and percent good quality nuclei. Caspase-3 and Poly (ADP-ribose) polymerase-1 (PARP1) have been implicated in different cell death pathways in response to neonatal HI. Another gene, sirtuin1 (SIRT1), has been demonstrated to have neuroprotective and anti-apoptotic properties. To assess the correlation between the severity of injury and genes involved in cell survival/death, we analyzed caspase-3, PARP1, and SIRT1 mRNA expressions in hippocampi 3 days post-HI and sham surgery, using quantitative reverse transcription polymerase chain reaction. The ipsilateral (IL) hippocampal caspase-3 and SIRT1 mRNA expressions post-HI were significantly higher than sham IL hippocampi and positively correlated with the novel injury scores in both males and females. We detected a statistically significant sex difference in IL hippocampal caspase-3 mRNA expression with comparable injury scores between males and females with higher expression in females.

Genome-wide analysis reveals associations between climate and regional patterns of adaptive divergence and dispersal in American pikas.

Understanding the role of adaptation in species' responses to climate change is important for evaluating the evolutionary potential of populations and informing conservation efforts. Population genomics provides a useful approach for identifying putative signatures of selection and the underlying environmental factors or biological processes that may be involved. Here, we employed a population genomic approach within a space-for-time study design to investigate the genetic basis of local adaptation and reconstruct patterns of movement across rapidly changing environments in a thermally sensitive mammal, the American pika (Ochotona princeps). Using genotypic data at 49,074 single-nucleotide polymorphisms (SNPs), we analyzed patterns of genome-wide diversity, structure, and migration along three independent elevational transects located at the northern extent (Tweedsmuir South Provincial Park, British Columbia, Canada) and core (North Cascades National Park, Washington, USA) of the Cascades lineage. We identified 899 robust outlier SNPs within- and among-transects. Of those annotated to genes with known function, many were linked with cellular processes related to climate stress including ATP-binding, ATP citrate synthase activity, ATPase activity, hormone activity, metal ion-binding, and protein-binding. Moreover, we detected evidence for contrasting patterns of directional migration along transects across geographic regions that suggest an increased propensity for American pikas to disperse among lower elevation populations at higher latitudes where environments are generally cooler. Ultimately, our data indicate that fine-scale demographic patterns and adaptive processes may vary among populations of American pikas, providing an important context for evaluating biotic responses to climate change in this species and other alpine-adapted mammals.

Offline World: the Internet as Social Infrastructure among the Unconnected in Quasi-Rural Illinois.

The United States continues to experience a persistent rural-urban digital divide. However, in this area of research, less attention has been paid to the divide in regions between these two demographic and geographic extremes. In this paper, we examine the perceived effects of internet inaccessibility in this in-between space, which we term "quasi-rural." Using quasi-rural Illinois as a case study, semi-structured interview data is used to highlight the experiences of those who are directly affected by the digital divide, as well as those who provide service alternatives and advocate for internet connectivity. With this data, we describe the personal experience of at-home internet inaccessibility. We then focus on how limited, or a lack of access shapes the perceptions of community connectedness and disadvantage among those affected by the divide. Our findings demonstrate the internet's function as social infrastructure; differences in access are then conceptualized as a form of socioeconomic inequality.

Genotyping-in-Thousands by sequencing reveals marked population structure in Western Rattlesnakes to inform conservation status.

Delineation of units below the species level is critical for prioritizing conservation actions for species at-risk. Genetic studies play an important role in characterizing patterns of population connectivity and diversity to inform the designation of conservation units, especially for populations that are geographically isolated. The northernmost range margin of Western Rattlesnakes (Crotalus oreganus) occurs in British Columbia, Canada, where it is federally classified as threatened and restricted to five geographic regions. In these areas, Western Rattlesnakes hibernate (den) communally, raising questions about connectivity within and between den complexes. At present, Western Rattlesnake conservation efforts are hindered by a complete lack of information on genetic structure and degree of isolation at multiple scales, from the den to the regional level. To fill this knowledge gap, we used Genotyping-in-Thousands by sequencing (GT-seq) to genotype an optimized panel of 362 single nucleotide polymorphisms (SNPs) from individual samples (n = 461) collected across the snake's distribution in western Canada and neighboring Washington (USA). Hierarchical STRUCTURE analyses found evidence for population structure within and among the five geographic regions in BC, as well as in Washington. Within these regions, 11 genetically distinct complexes of dens were identified, with some regions having multiple complexes. No significant pattern of isolation-by-distance and generally low levels of migration were detected among den complexes across regions. Additionally, snakes within dens generally were more related than those among den complexes within a region, indicating limited movement. Overall, our results suggest that the single, recognized designatable unit for Western Rattlesnakes in Canada should be re-assessed to proactively focus conservation efforts on preserving total genetic variation detected range-wide. More broadly, our study demonstrates a novel application of GT-seq for investigating patterns of diversity in wild populations at multiple scales to better inform conservation management.

Genotyping-in-Thousands by sequencing (GT-seq) panel development and application to minimally invasive DNA samples to support studies in molecular ecology.

Minimally invasive sampling (MIS) is widespread in wildlife studies; however, its utility for massively parallel DNA sequencing (MPS) is limited. Poor sample quality and contamination by exogenous DNA can make MIS challenging to use with modern genotyping-by-sequencing approaches, which have been traditionally developed for high-quality DNA sources. Given that MIS is often more appropriate in many contexts, there is a need to make such samples practical for harnessing MPS. Here, we test the ability for Genotyping-in-Thousands by sequencing (GT-seq), a multiplex amplicon sequencing approach, to effectively genotype minimally invasive cloacal DNA samples collected from the Western Rattlesnake (Crotalus oreganus), a threatened species in British Columbia, Canada. As there was no previous genetic information for this species, an optimized panel of 362 SNPs was selected for use with GT-seq from a de novo restriction site-associated DNA sequencing (RADseq) assembly. Comparisons of genotypes generated within and among RADseq and GT-seq for the same individuals found low rates of genotyping error (GT-seq: 0.50%; RADseq: 0.80%) and discordance (2.57%), the latter likely due to the different genotype calling models employed. GT-seq mean genotype discordance between blood and cloacal swab samples collected from the same individuals was also minimal (1.37%). Estimates of population diversity parameters were similar across GT-seq and RADseq data sets, as were inferred patterns of population structure. Overall, GT-seq can be effectively applied to low-quality DNA samples, minimizing the inefficiencies presented by exogenous DNA typically found in minimally invasive samples and continuing the expansion of molecular ecology and conservation genetics in the genomics era.

Clinical and serological analysis of patients with positive anticyclic citrullinated Peptide antibodies referred through a Rheumatology Central Triage System.

OBJECTIVE: Anticitrullinated protein antibodies (ACPA) are a highly specific and sensitive biomarker for the diagnosis of rheumatoid arthritis (RA). Some patients who were found to have a positive ACPA test were referred to our Rheumatology Central Triage (CT; Calgary, Alberta, Canada) for assessment by a rheumatologist. The objectives of our study were to determine the clinical accuracy of ACPA in establishing a diagnosis of RA in a real-time clinical setting. METHODS: Cases that met 3 criteria were included in the study: (1) referred to the CT over 3 calendar years (n = 20,389), (2) reason for referral was a positive ACPA test (n = 568), and (3) evaluated by a certified rheumatologist (n = 314). An administrative serological database was used to retrieve specific ACPA results. RESULTS: Of patients referred through our CT for evaluation of a positive ACPA test, 57.6% received a diagnosis of RA; the remainder had a variety of other diagnoses, some of which might be considered early RA (9%). The predictive values of ACPA for the diagnosis of RA were increased when rheumatoid factor (RF) results were included in the analysis. When definite and possible RA were combined and the prevalence of moderate/high ACPA was compared to all other individuals, the positive and negative predictive values for moderate/high ACPA for RA were 74.3% and 68.4%, respectively. CONCLUSION: About 58% of patients with a positive ACPA referred through a triage system for a rheumatologist opinion received a diagnosis of RA at their first visit. RF provides additional useful information to guide the diagnosis and urgency of referral.