Variation in histone configurations correlates with gene expression across nine inbred strains of mice.

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Systems genetics uncover new loci containing functional gene candidates in Mycobacterium tuberculosis-infected Diversity Outbred mice.

Mycobacterium tuberculosis infects two billion people across the globe, and results in 8-9 million new tuberculosis (TB) cases and 1-1.5 million deaths each year. Most patients have no known genetic basis that predisposes them to disease. Here, we investigate the complex genetic basis of pulmonary TB by modelling human genetic diversity with the Diversity Outbred mouse population. When infected with M. tuberculosis, one-third develop early onset, rapidly progressive, necrotizing granulomas and succumb within 60 days. The remaining develop non-necrotizing granulomas and survive longer than 60 days. Genetic mapping using immune and inflammatory mediators; and clinical, microbiological, and granuloma correlates of disease identified five new loci on mouse chromosomes 1, 2, 4, 16; and three known loci on chromosomes 3 and 17. Further, multiple positively correlated traits shared loci on chromosomes 1, 16, and 17 and had similar patterns of allele effects, suggesting these loci contain critical genetic regulators of inflammatory responses to M. tuberculosis. To narrow the list of candidate genes, we used a machine learning strategy that integrated gene expression signatures from lungs of M. tuberculosis-infected Diversity Outbred mice with gene interaction networks to generate scores representing functional relationships. The scores were used to rank candidates for each mapped trait, resulting in 11 candidate genes: Ncf2, Fam20b, S100a8, S100a9, Itgb5, Fstl1, Zbtb20, Ddr1, Ier3, Vegfa, and Zfp318. Although all candidates have roles in infection, inflammation, cell migration, extracellular matrix remodeling, or intracellular signaling, and all contain single nucleotide polymorphisms (SNPs), SNPs in only four genes (S100a8, Itgb5, Fstl1, Zfp318) are predicted to have deleterious effects on protein functions. We performed methodological and candidate validations to (i) assess biological relevance of predicted allele effects by showing that Diversity Outbred mice carrying PWK/PhJ alleles at the H-2 locus on chromosome 17 QTL have shorter survival; (ii) confirm accuracy of predicted allele effects by quantifying S100A8 protein in inbred founder strains; and (iii) infection of C57BL/6 mice deficient for the S100a8 gene. Overall, this body of work demonstrates that systems genetics using Diversity Outbred mice can identify new (and known) QTLs and functionally relevant gene candidates that may be major regulators of complex host-pathogens interactions contributing to granuloma necrosis and acute inflammation in pulmonary TB.

Differential regulation of gene expression pathways with dexamethasone and ACTH after early life seizures.

Early-life seizures (ELS) are associated with persistent cognitive deficits such as ADHD and memory impairment. These co-morbidities have a dramatic negative impact on the quality of life of patients. Therapies that improve cognitive outcomes have enormous potential to improve patients' quality of life. Our previous work in a rat flurothyl-induction model showed that administration of adrenocorticotropic hormone (ACTH) at time of seizure induction led to improved learning and memory in the animals despite no effect on seizure latency or duration. Administration of dexamethasone (Dex), a corticosteroid, did not have the same positive effect on learning and memory and has even been shown to exacerbate injury in a rat model of temporal lobe epilepsy. We hypothesized that ACTH exerted positive effects on cognitive outcomes through beneficial changes to gene expression and proposed that administration of ACTH at seizure induction would return gene-expression in the brain towards the normal pattern of expression in the Control animals whereas Dex would not. Twenty-six Sprague-Dawley rats were randomized into vehicle- Control, and ACTH-, Dex-, and vehicle- ELS. Rat pups were subjected to 60 flurothyl seizures from P5 to P14. After seizure induction, brains were removed and the hippocampus and PFC were dissected, RNA was extracted and sequenced, and differential expression analysis was performed using generalized estimating equations. Differential expression analysis showed that ACTH pushes gene expression in the brain back to a more normal state of expression through enrichment of pathways involved in supporting homeostatic balance and down-regulating pathways that might contribute to excitotoxic cell-damage post-ELS.

Fine Spike Timing in Hippocampal-Prefrontal Ensembles Predicts Poor Encoding and Underlies Behavioral Performance in Healthy and Malformed Brains.

Spatial working memory (SWM) is a central cognitive process during which the hippocampus and prefrontal cortex (PFC) encode and maintain spatial information for subsequent decision-making. This occurs in the context of ongoing computations relating to spatial position, recall of long-term memory, attention, among many others. To establish how intermittently presented information is integrated with ongoing computations we recorded single units, simultaneously in hippocampus and PFC, in control rats and those with a brain malformation during performance of an SWM task. Neurons that encode intermittent task parameters are also well modulated in time and incorporated into a functional network across regions. Neurons from animals with cortical malformation are poorly modulated in time, less likely to encode task parameters, and less likely to be integrated into a functional network. Our results implicate a model in which ongoing oscillatory coordination among neurons in the hippocampal-PFC network describes a functional network that is poised to receive sensory inputs that are then integrated and multiplexed as working memory. The background temporal modulation is systematically altered in disease, but the relationship between these dynamics and behaviorally relevant firing is maintained, thereby providing potential targets for stimulation-based therapies.

Identification of novel loci controlling inflammatory bowel disease susceptibility utilizing the genetic diversity of wild-derived mice.

Inflammatory bowel disease (IBD) is a complex disorder that imposes a growing health burden. Multiple genetic associations have been identified in IBD, but the mechanisms underlying many of these associations are poorly understood. Animal models are needed to bridge this gap, but conventional laboratory mouse strains lack the genetic diversity of human populations. To more accurately model human genetic diversity, we utilized a panel of chromosome (Chr) substitution strains, carrying chromosomes from the wild-derived and genetically divergent PWD/PhJ (PWD) strain on the commonly used C57BL/6J (B6) background, as well as their parental B6 and PWD strains. Two models of IBD were used, TNBS- and DSS-induced colitis. Compared with B6 mice, PWD mice were highly susceptible to TNBS-induced colitis, but resistant to DSS-induced colitis. Using consomic mice, we identified several PWD-derived loci that exhibited profound effects on IBD susceptibility. The most pronounced of these were loci on Chr1 and Chr2, which yielded high susceptibility in both IBD models, each acting at distinct phases of the disease. Leveraging transcriptomic data from B6 and PWD immune cells, together with a machine learning approach incorporating human IBD genetic associations, we identified lead candidate genes, including Itga4, Pip4k2a, Lcn10, Lgmn, and Gpr65.

Detection and Classification of Novel Renal Histologic Phenotypes Using Deep Neural Networks.

With the advent and increased accessibility of deep neural networks (DNNs), complex properties of histologic images can be rigorously and reproducibly quantified. We used DNN-based transfer learning to analyze histologic images of periodic acid-Schiff-stained renal sections from a cohort of mice with different genotypes. We demonstrate that DNN-based machine learning has strong generalization performance on multiple histologic image processing tasks. The neural network extracted quantitative image features and used them as classifiers to look for differences between mice of different genotypes. Excellent performance was observed at segmenting glomeruli from non-glomerular structure and subsequently predicting the genotype of the animal on the basis of glomerular quantitative image features. The DNN-based genotype classifications highly correlate with mesangial matrix expansion scored by a pathologist (R.E.C.), which differed in these animals. In addition, by analyzing non-glomeruli images, the neural network identified novel histologic features that differed by genotype, including the presence of vacuoles, nuclear count, and proximal tubule brush border integrity, which was validated with immunohistologic staining. These features were not identified in systematic pathologic examination. Our study demonstrates the power of DNNs to extract biologically relevant phenotypes and serve as a platform for discovering novel phenotypes. These results highlight the synergistic possibilities for pathologists and DNNs to radically scale up our ability to generate novel mechanistic hypotheses in disease.

Global development and adaptive behaviour in children with early-onset epilepsy: a population-based case-control study.

AIM: There are limited population-based data on global development and adaptive behaviour in children with early-onset epilepsy. The aims of this study were: (1) to identify the prevalence of deficits in global development and adaptive behaviour experienced by children with early-onset epilepsy; (2) to identify factors associated with such deficits; and (3) to compare the relationship between measures of neurodevelopment in the group with epilepsy to a group without epilepsy who had other neurological or neurodevelopmental difficulties. METHOD: The Sussex Early Epilepsy and Neurobehaviour study is a prospective, community-based study involving children (1-7y) with epilepsy. We undertook comprehensive psychological assessment with participants, including measures of global development and adaptive behaviour. We compared the children with epilepsy with a sex, age, and developmentally-matched group of children without epilepsy who had neurodevelopmental or neurological difficulties using correlation matrices. RESULTS: Forty-eight children (91% of the eligible population) with epilepsy underwent assessment. Seventy-one per cent of children displayed delayed global development (<2SD) and 56% showed significant deficits (<2SD) in adaptive behaviour. Our analysis revealed that non-white ethnicity and use of polytherapy were independently associated with decreased scores on measures of global development and adaptive behaviour. The correlations between measures of developmental functioning were higher in children with epilepsy than in those without. INTERPRETATION: Children with early-onset epilepsy frequently have difficulties with global development and adaptive behaviour. The higher correlations between neurodevelopmental measures in children with epilepsy suggest that the profile in children with epilepsy is different. This may have significant implications for both neuropathology and interventions. WHAT THIS PAPER ADDS: Children with early-onset epilepsy are at significant risk of intellectual disability. Developmental impairment is associated with use of polytherapy but not with any seizure parameters. Developmental profiles in young children with epilepsy differ from other conditions.

WONOEP APPRAISAL: The many facets of epilepsy networks.

The brain is a complex system composed of networks of interacting elements, from genes to circuits, whose function (and dysfunction) is not derivable from the superposition of individual components. Epilepsy is frequently described as a network disease, but to date, there is no standardized framework within which network concepts applicable to all levels from genes to whole brain can be used to generate deeper insights into the pathogenesis of seizures or the associated morbidities. To address this shortcoming, the Neurobiology Commission of the International League Against Epilepsy dedicated a Workshop on Neurobiology of Epilepsy (XIV WONOEP 2017) with the aim of formalizing network concepts as they apply to epilepsy and to critically discuss whether and how such concepts could augment current research endeavors. Here, we review concepts and strategies derived by considering epilepsy as a disease of different network hierarchies that range from genes to clinical phenotypes. We propose that the concept of networks is important for understanding epilepsy and is critical for developing new study designs. These approaches could ultimately facilitate the development of novel diagnostic and therapeutic strategies.

The Menopause Visual Analogue Scale: a new tool for measuring the severity and response to treatment of symptoms throughout the menopausal transition.

OBJECTIVES: This study aimed to assess the reliability and validity of the Menopause Visual Analogue Scale (MVAS) in measuring symptoms throughout the menopausal transition. METHODS: Two independent samples of women undergoing the menopausal transition completed both the MVAS and the Greene Climacteric Scale (GCS) at a women's mental health clinic between 2008 and 2016. Data for the first sample were obtained using a retrospective chart review of patients seen between 2008 and 2012 (N1 = 75) and data for the second sample came from a prospective study conducted between 2013 and 2016 (N2 = 86). Internal consistency was assessed using Cronbach's α and Pearson's correlation coefficient was used to evaluate concurrent validity. Bland-Altman plots were developed to assess the degree of agreement between the scales. RESULTS: Internal consistency for the physical and psychological domains of the MVAS was 0.80-0.81 and 0.92-0.94, respectively. Pearson's correlations between the MVAS and the GCS were high for both physical (rphys = 0.74-0.76, p < 0.01) and psychological (rpsych = 0.70-0.72, p < 0.01) components in both samples. Changes in MVAS physical and psychological scores in response to treatment were correlated with changes in GCS physical and psychological scores (rphys = 0.69, p < 0.01; rpsych = 0.49, p < 0.01) in the second sample. Bland-Altman plots indicate low to moderate levels of agreement between most portions of the MVAS and the GCS. CONCLUSIONS: These findings suggest that the MVAS has potential for assessing both severity and change in symptoms throughout the menopausal transition, subject to exploring limitations identified in the analysis and application to other populations.

A multiprotein occupancy map of the mRNP on the 3' end of histone mRNAs.

The animal replication-dependent (RD) histone mRNAs are coordinately regulated with chromosome replication. The RD-histone mRNAs are the only known cellular mRNAs that are not polyadenylated. Instead, the mature transcripts end in a conserved stem-loop (SL) structure. This SL structure interacts with the stem-loop binding protein (SLBP), which is involved in all aspects of RD-histone mRNA metabolism. We used several genomic methods, including high-throughput sequencing of cross-linked immunoprecipitate (HITS-CLIP) to analyze the RNA-binding landscape of SLBP. SLBP was not bound to any RNAs other than histone mRNAs. We performed bioinformatic analyses of the HITS-CLIP data that included (i) clustering genes by sequencing read coverage using CVCA, (ii) mapping the bound RNA fragment termini, and (iii) mapping cross-linking induced mutation sites (CIMS) using CLIP-PyL software. These analyses allowed us to identify specific sites of molecular contact between SLBP and its RD-histone mRNA ligands. We performed in vitro crosslinking assays to refine the CIMS mapping and found that uracils one and three in the loop of the histone mRNA SL preferentially crosslink to SLBP, whereas uracil two in the loop preferentially crosslinks to a separate component, likely the 3'hExo. We also performed a secondary analysis of an iCLIP data set to map UPF1 occupancy across the RD-histone mRNAs and found that UPF1 is bound adjacent to the SLBP-binding site. Multiple proteins likely bind the 3' end of RD-histone mRNAs together with SLBP.

WORMHOLE: Novel Least Diverged Ortholog Prediction through Machine Learning.

The rapid advancement of technology in genomics and targeted genetic manipulation has made comparative biology an increasingly prominent strategy to model human disease processes. Predicting orthology relationships between species is a vital component of comparative biology. Dozens of strategies for predicting orthologs have been developed using combinations of gene and protein sequence, phylogenetic history, and functional interaction with progressively increasing accuracy. A relatively new class of orthology prediction strategies combines aspects of multiple methods into meta-tools, resulting in improved prediction performance. Here we present WORMHOLE, a novel ortholog prediction meta-tool that applies machine learning to integrate 17 distinct ortholog prediction algorithms to identify novel least diverged orthologs (LDOs) between 6 eukaryotic species-humans, mice, zebrafish, fruit flies, nematodes, and budding yeast. Machine learning allows WORMHOLE to intelligently incorporate predictions from a wide-spectrum of strategies in order to form aggregate predictions of LDOs with high confidence. In this study we demonstrate the performance of WORMHOLE across each combination of query and target species. We show that WORMHOLE is particularly adept at improving LDO prediction performance between distantly related species, expanding the pool of LDOs while maintaining low evolutionary distance and a high level of functional relatedness between genes in LDO pairs. We present extensive validation, including cross-validated prediction of PANTHER LDOs and evaluation of evolutionary divergence and functional similarity, and discuss future applications of machine learning in ortholog prediction. A WORMHOLE web tool has been developed and is available at http://wormhole.jax.org/.

Essential Resources for Implementation and Sustainability of Evidence-Based Health Promotion Programs: A Mixed Methods Multi-Site Case Study.

As of October 2016, use of federal Older Americans Act funds for health promotion and disease prevention will be restricted to the Administration on Aging's criteria for high-level evidence-based health promotion programs. Dissemination of these programs to rural communities remains limited. Therefore a strong need exists to identify strategies that facilitate program implementation and sustainability. The objective of this study was to compare organizational readiness and implementation strategies used by rural communities that achieved varying levels of success in sustaining evidence-based health promotion programs for older adults. We utilized a qualitative multi-site case study design to analyze the longitudinal experiences of eight rural sites working to implement evidence-based health promotion program over 3 years (8/2012-7/2015). Multiple sources of data (interviews, documents, reports, surveys) from each site informed the analysis. We used conventional content analysis to conduct a cross-case comparison to identify common features of rural counties that successfully implemented and sustained their target evidence-based health promotion program. Readiness to implement evidence-based programs as low at baseline as all site leaders described needing to secure additional resources for program implementation. Sites that successfully utilized six essential resources implemented and sustained greater numbers of workshops: (1) External Partnerships, (2) Agency Leadership Commitment, (3) Ongoing Source of Workshop Leaders, (4) Health Promotion Coordination Tasks Assigned to Specific Staff, (5) Organizational Stability, and (6) Change Team Engagement. The six essential resources described in this study can help rural communities assess their readiness to implement health promotion programs and work secure the resources necessary for successful implementation.

Systems level analysis of systemic sclerosis shows a network of immune and profibrotic pathways connected with genetic polymorphisms.

Systemic sclerosis (SSc) is a rare systemic autoimmune disease characterized by skin and organ fibrosis. The pathogenesis of SSc and its progression are poorly understood. The SSc intrinsic gene expression subsets (inflammatory, fibroproliferative, normal-like, and limited) are observed in multiple clinical cohorts of patients with SSc. Analysis of longitudinal skin biopsies suggests that a patient's subset assignment is stable over 6-12 months. Genetically, SSc is multi-factorial with many genetic risk loci for SSc generally and for specific clinical manifestations. Here we identify the genes consistently associated with the intrinsic subsets across three independent cohorts, show the relationship between these genes using a gene-gene interaction network, and place the genetic risk loci in the context of the intrinsic subsets. To identify gene expression modules common to three independent datasets from three different clinical centers, we developed a consensus clustering procedure based on mutual information of partitions, an information theory concept, and performed a meta-analysis of these genome-wide gene expression datasets. We created a gene-gene interaction network of the conserved molecular features across the intrinsic subsets and analyzed their connections with SSc-associated genetic polymorphisms. The network is composed of distinct, but interconnected, components related to interferon activation, M2 macrophages, adaptive immunity, extracellular matrix remodeling, and cell proliferation. The network shows extensive connections between the inflammatory- and fibroproliferative-specific genes. The network also shows connections between these subset-specific genes and 30 SSc-associated polymorphic genes including STAT4, BLK, IRF7, NOTCH4, PLAUR, CSK, IRAK1, and several human leukocyte antigen (HLA) genes. Our analyses suggest that the gene expression changes underlying the SSc subsets may be long-lived, but mechanistically interconnected and related to a patients underlying genetic risk.

Cognitive impairment in temporal lobe epilepsy: role of online and offline processing of single cell information.

Cognitive impairment is a common comorbidity in temporal lobe epilepsy (TLE) and is often considered more detrimental to quality of life than seizures. While it has been previously shown that the encoding of memory during behavior is impaired in the pilocarpine model of TLE in rats, how this information is consolidated during the subsequent sleep period remains unknown. In this study, we first report marked deficits in spatial memory performance and severe cell loss in the CA1 layer of the hippocampus lower spatial coherence of firing in TLE rats. We then present the first evidence that the reactivation of behavior-driven patterns of activity of CA1 place cells in the hippocampus is intact in TLE rats. Using a template-matching method, we discovered that real-time (3-5 s) reactivation structure was intact in TLE rats. Furthermore, we estimated the entropy rate of short time scale (∼250 ms) bursting activity using block entropies and found that significant, extended temporal correlations exist in both TLE and control rats. Fitting a first-order Markov Chain model to these bursting time series, we found that long sequences derived from behavior were significantly enriched in the Markov model over corresponding models fit on randomized data confirming the presence of replay in shorter time scales. We propose that the persistent consolidation of poor spatial information in both real time and during bursting activity may contribute to memory impairments in TLE rats.

Adjuvant chemotherapy for bladder cancer-why does level 1 evidence not support it?

Neoadjuvant cisplatin-based combination chemotherapy provides a 5% increase in cure rate, an increase in median survival of about 3 years, and statistically significant and clinically relevant increments in overall survival for patients with invasive bladder cancer. Despite compelling level 1 data, it has become quite clear that facts that are similar to those that changed the paradigm of treatment of breast cancer in the 1970s have not had a similar influence on patterns of practice in bladder cancer care. Instead of using this proven approach, cystectomy alone or surgery followed by adjuvant chemotherapy is often used as a functional alternative for patients with deeply invasive and/or node-metastatic disease discovered at radical cystectomy. However, there is no well-powered level 1 evidence to support routine adjuvant chemotherapy for invasive bladder cancer, and some randomized trials have shown inferior outcomes. There is a clear need for a well-designed, randomized trial that tests the utility of adjuvant chemotherapy for invasive bladder cancer, but until that has been completed, neoadjuvant chemotherapy followed by definitive local treatment should be the standard of care for invasive bladder cancer.

Plasma brain derived neurotrophic factor (BDNF) and response to ketamine in treatment-resistant depression.

Ketamine produces rapid antidepressant effects in treatment-resistant depression (TRD), but the magnitude of response varies considerably between individual patients. Brain-derived neurotrophic factor (BDNF) has been investigated as a biomarker of treatment response in depression and has been implicated in the mechanism of action of ketamine. We evaluated plasma BDNF and associations with symptoms in 22 patients with TRD enrolled in a randomized controlled trial of ketamine compared to an anaesthetic control (midazolam). Ketamine significantly increased plasma BDNF levels in responders compared to non-responders 240 min post-infusion, and Montgomery-Åsberg Depression Rating Scale (MADRS) scores were negatively correlated with BDNF (r=-0.701, p = 0.008). Plasma BDNF levels at 240 min post-infusion were highly negatively associated with MADRS scores at 240 min (r = -0.897, p=.002), 24 h (r = -0.791, p = 0.038), 48 h (r = -0.944, p = 0.001) and 72 h (r = -0.977, p = 0.010). No associations with BDNF were found for patients receiving midazolam. These data support plasma BDNF as a peripheral biomarker relevant to ketamine antidepressant response.

Causal network perturbation analysis identifies known and novel type-2 diabetes driver genes.

The molecular pathogenesis of diabetes is multifactorial, involving genetic predisposition and environmental factors that are not yet fully understood. However, pancreatic ß-cell failure remains among the primary reasons underlying the progression of type-2 diabetes (T2D) making targeting ß-cell dysfunction an attractive pathway for diabetes treatment. To identify genetic contributors to ß-cell dysfunction, we investigated single-cell gene expression changes in ß-cells from healthy (C57BL/6J) and diabetic (NZO/HlLtJ) mice fed with normal or high-fat, high-sugar diet (HFHS). Our study presents an innovative integration of the causal network perturbation assessment (ssNPA) framework with meta-cell transcriptome analysis to explore the genetic underpinnings of type-2 diabetes (T2D). By generating a reference causal network and in silico perturbation, we identified novel genes implicated in T2D and validated our candidates using the Knockout Mouse Phenotyping (KOMP) Project database.

Visual detection of seizures in mice using supervised machine learning.

Seizures are caused by abnormally synchronous brain activity that can result in changes in muscle tone, such as twitching, stiffness, limpness, or rhythmic jerking. These behavioral manifestations are clear on visual inspection and the most widely used seizure scoring systems in preclinical models, such as the Racine scale in rodents, use these behavioral patterns in semiquantitative seizure intensity scores. However, visual inspection is time-consuming, low-throughput, and partially subjective, and there is a need for rigorously quantitative approaches that are scalable. In this study, we used supervised machine learning approaches to develop automated classifiers to predict seizure severity directly from noninvasive video data. Using the PTZ-induced seizure model in mice, we trained video-only classifiers to predict ictal events, combined these events to predict an univariate seizure intensity for a recording session, as well as time-varying seizure intensity scores. Our results show, for the first time, that seizure events and overall intensity can be rigorously quantified directly from overhead video of mice in a standard open field using supervised approaches. These results enable high-throughput, noninvasive, and standardized seizure scoring for downstream applications such as neurogenetics and therapeutic discovery.

A machine learning approach for quantifying age-related histological changes in the mouse kidney.

The ability to quantify aging-related changes in histological samples is important, as it allows for evaluation of interventions intended to effect health span. We used a machine learning architecture that can be trained to detect and quantify these changes in the mouse kidney. Using additional held out data, we show validation of our model, correlation with scores given by pathologists using the Geropathology Research Network aging grading scheme, and its application in providing reproducible and quantifiable age scores for histological samples. Aging quantification also provides the insights into possible changes in image appearance that are independent of specific geropathology-specified lesions. Furthermore, we provide trained classifiers for H&E-stained slides, as well as tutorials on how to use these and how to create additional classifiers for other histological stains and tissues using our architecture. This architecture and combined resources allow for the high throughput quantification of mouse aging studies in general and specifically applicable to kidney tissues.

Inferring single-cell spatial gene expression with tissue morphology via explainable deep learning.

The spatial arrangement of cells is vital in developmental processes and organogenesis in multicellular life forms. Deep learning models trained with spatial omics data uncover complex patterns and relationships among cells, genes, and proteins in a high-dimensional space, providing new insights into biological processes and diseases. State-of-the-art in silico spatial multi-cell gene expression methods using histological images of tissue stained with hematoxylin and eosin (H&E) to characterize cellular heterogeneity. These computational techniques offer the advantage of analyzing vast amounts of spatial data in a scalable and automated manner, thereby accelerating scientific discovery and enabling more precise medical diagnostics and treatments. In this work, we developed a vision transformer (ViT) framework to map histological signatures to spatial single-cell transcriptomic signatures, named SPiRiT ( S patial Omics P rediction and R eproducibility integrated T ransformer). Our framework was enhanced by integrating cross validation with model interpretation during hyper-parameter tuning. SPiRiT predicts single-cell spatial gene expression using the matched histopathological image tiles of human breast cancer and whole mouse pup, evaluated by Xenium (10x Genomics) datasets. Furthermore, ViT model interpretation reveals the high-resolution, high attention area (HAR) that the ViT model uses to predict the gene expression, including marker genes for invasive cancer cells ( FASN ), stromal cells ( POSTN ), and lymphocytes ( IL7R ). In an apple-to-apple comparison with the ST-Net Convolutional Neural Network algorithm, SPiRiT improved predictive accuracy by 40% using human breast cancer Visium (10x Genomics) dataset. Cancer biomarker gene prediction and expression level are highly consistent with the tumor region annotation. In summary, our work highlights the feasibility to infer spatial single-cell gene expression using tissue morphology in multiple-species, i.e., human and mouse, and multi-organs, i.e., mouse whole body morphology. Importantly, incorporating model interpretation and vision transformer is expected to serve as a general-purpose framework for spatial transcriptomics.