Cross-platform proteomics signatures of extreme old age.

In previous work we used a Somalogic platform targeting approximately 5000 proteins to generate a serum protein signature of centenarians that we validated in independent studies that used the same technology. We set here to validate and possibly expand the results by profiling the serum proteome of a subset of individuals included in the original study using liquid chromatography tandem mass spectrometry (LC-MS/MS). Following pre-processing, the LC-MS/MS data provided quantification of 398 proteins, with only 266 proteins shared by both platforms. At 1% FDR statistical significance threshold, the analysis of LC-MS/MS data detected 44 proteins associated with extreme old age, including 23 of the original analysis. To identify proteins for which associations between expression and extreme-old age were conserved across platforms, we performed inter-study conservation testing of the 266 proteins quantified by both platforms using a method that accounts for the correlation between the results. From these tests, a total of 80 proteins reached 5% FDR statistical significance, and 26 of these proteins had concordant pattern of gene expression in whole blood. This signature of 80 proteins points to blood coagulation, IGF signaling, extracellular matrix (ECM) organization, and complement cascade as important pathways whose protein level changes provide evidence for age-related adjustments that distinguish centenarians from younger individuals.

STING promotes homeostatic maintenance of tissues and confers longevity with aging.

Local immune processes within aging tissues are a significant driver of aging associated dysfunction, but tissue-autonomous pathways and cell types that modulate these responses remain poorly characterized. The cytosolic DNA sensing pathway, acting through cyclic GMP-AMP synthase (cGAS) and Stimulator of Interferon Genes (STING), is broadly expressed in tissues, and is poised to regulate local type I interferon (IFN-I)-dependent and independent inflammatory processes within tissues. Recent studies suggest that the cGAS/STING pathway may drive pathology in various in vitro and in vivo models of accelerated aging. To date, however, the role of the cGAS/STING pathway in physiological aging processes, in the absence of genetic drivers, has remained unexplored. This remains a relevant gap, as STING is ubiquitously expressed, implicated in multitudinous disorders, and loss of function polymorphisms of STING are highly prevalent in the human population (>50%). Here we reveal that, during physiological aging, STING-deficiency leads to a significant shortening of murine lifespan, increased pro-inflammatory serum cytokines and tissue infiltrates, as well as salient changes in histological composition and organization. We note that aging hearts, livers, and kidneys express distinct subsets of inflammatory, interferon-stimulated gene (ISG), and senescence genes, collectively comprising an immune fingerprint for each tissue. These distinctive patterns are largely imprinted by tissue-specific stromal and myeloid cells. Using cellular interaction network analyses, immunofluorescence, and histopathology data, we show that these immune fingerprints shape the tissue architecture and the landscape of cell-cell interactions in aging tissues. These age-associated immune fingerprints are grossly dysregulated with STING-deficiency, with key genes that define aging STING-sufficient tissues greatly diminished in the absence of STING. Changes in immune signatures are concomitant with a restructuring of the stromal and myeloid fractions, whereby cell:cell interactions are grossly altered and resulting in disorganization of tissue architecture in STING-deficient organs. This altered homeostasis in aging STING-deficient tissues is associated with a cross-tissue loss of homeostatic tissue-resident macrophage (TRM) populations in these tissues. Ex vivo analyses reveal that basal STING-signaling limits the susceptibility of TRMs to death-inducing stimuli and determines their in situ localization in tissue niches, thereby promoting tissue homeostasis. Collectively, these data upend the paradigm that cGAS/STING signaling is primarily pathological in aging and instead indicate that basal STING signaling sustains tissue function and supports organismal longevity. Critically, our study urges caution in the indiscriminate targeting of these pathways, which may result in unpredictable and pathological consequences for health during aging.

ß-catenin/CBP activation of mTORC1 signaling promotes partial epithelial-mesenchymal states in head and neck cancer.

Head and neck cancers, which include oral squamous cell carcinoma (OSCC) as a major subsite, exhibit cellular plasticity that includes features of an epithelial-mesenchymal transition (EMT), referred to as partial-EMT (p-EMT). To identify molecular mechanisms contributing to OSCC plasticity, we performed a multiphase analysis of single cell RNA sequencing (scRNAseq) data from human OSCC. This included a multiresolution characterization of cancer cell subgroups to identify pathways and cell states that are heterogeneously represented, followed by casual inference analysis to elucidate activating and inhibitory relationships between these pathways and cell states. This approach revealed signaling networks associated with hierarchical cell state transitions, which notably included an association between ß-catenin-driven CREB-binding protein (CBP) activity and mTORC1 signaling. This network was associated with subpopulations of cancer cells that were enriched for markers of the p-EMT state and poor patient survival. Functional analyses revealed that ß-catenin/CBP induced mTORC1 activity in part through the transcriptional regulation of a raptor-interacting protein, chaperonin containing TCP1 subunit 5 (CCT5). Inhibition of ß-catenin-CBP activity through the use of the orally active small molecule, E7386, reduced the expression of CCT5 and mTORC1 activity in vitro, and inhibited p-EMT-associated markers and tumor development in a murine model of OSCC. Our study highlights the use of multiresolution network analyses of scRNAseq data to identify targetable signals for therapeutic benefit, thus defining an underappreciated association between ß-catenin/CBP and mTORC1 signaling in head and neck cancer plasticity.

A quantitative trait variant in Gabra2 underlies increased methamphetamine stimulant sensitivity.

Psychostimulant (methamphetamine, cocaine) use disorders have a genetic component that remains mostly unknown. We conducted genome-wide quantitative trait locus (QTL) analysis of methamphetamine stimulant sensitivity. To facilitate gene identification, we employed a Reduced Complexity Cross between closely related C57BL/6 mouse substrains and examined maximum speed and distance traveled over 30 min following methamphetamine (2 mg/kg, i.p.). For maximum methamphetamine-induced speed following the second and third administration, we identified a single genome-wide significant QTL on chromosome 11 that peaked near the Cyfip2 locus (LOD = 3.5, 4.2; peak = 21 cM [36 Mb]). For methamphetamine-induced distance traveled following the first and second administration, we identified a genome-wide significant QTL on chromosome 5 that peaked near a functional intronic indel in Gabra2 coding for the alpha-2 subunit of the GABA-A receptor (LOD = 3.6-5.2; peak = 34-35 cM [66-67 Mb]). Striatal cis-expression QTL mapping corroborated Gabra2 as a functional candidate gene underlying methamphetamine-induced distance traveled. CRISPR/Cas9-mediated correction of the mutant intronic deletion on the C57BL/6J background to the wild-type C57BL/6NJ allele was sufficient to reduce methamphetamine-induced locomotor activity toward the wild-type C57BL/6NJ-like level, thus validating the quantitative trait variant (QTV). These studies show the power and efficiency of Reduced Complexity Crosses in identifying causal variants underlying complex traits. Functionally restoring Gabra2 expression decreased methamphetamine stimulant sensitivity and supports preclinical and human genetic studies implicating the GABA-A receptor in psychostimulant addiction-relevant traits. Importantly, our findings have major implications for studying psychostimulants in the C57BL/6J strain-the gold standard strain in biomedical research.

Multi-resolution characterization of molecular taxonomies in bulk and single-cell transcriptomics data.

As high-throughput genomics assays become more efficient and cost effective, their utilization has become standard in large-scale biomedical projects. These studies are often explorative, in that relationships between samples are not explicitly defined a priori, but rather emerge from data-driven discovery and annotation of molecular subtypes, thereby informing hypotheses and independent evaluation. Here, we present K2Taxonomer, a novel unsupervised recursive partitioning algorithm and associated R package that utilize ensemble learning to identify robust subgroups in a 'taxonomy-like' structure. K2Taxonomer was devised to accommodate different data paradigms, and is suitable for the analysis of both bulk and single-cell transcriptomics, and other '-omics', data. For each of these data types, we demonstrate the power of K2Taxonomer to discover known relationships in both simulated and human tissue data. We conclude with a practical application on breast cancer tumor infiltrating lymphocyte (TIL) single-cell profiles, in which we identified co-expression of translational machinery genes as a dominant transcriptional program shared by T cells subtypes, associated with better prognosis in breast cancer tissue bulk expression data.

5' UTR variants in the quantitative trait gene Hnrnph1 support reduced 5' UTR usage and hnRNP H protein as a molecular mechanism underlying reduced methamphetamine sensitivity.

We previously identified a 210 kb region on chromosome 11 (50.37-50.58 Mb, mm10) containing two protein-coding genes (Hnrnph1, Rufy1) that was necessary for reduced methamphetamine-induced locomotor activity in C57BL/6J congenic mice harboring DBA/2J polymorphisms. Gene editing of a small deletion in the first coding exon supported Hnrnph1 as a quantitative trait gene. We have since shown that Hnrnph1 mutants also exhibit reduced methamphetamine-induced reward, reinforcement, and dopamine release. However, the quantitative trait variants (QTVs) that modulate Hnrnph1 function at the molecular level are not known. Nine single nucleotide polymorphisms and seven indels distinguish C57BL/6J from DBA/2J within Hnrnph1, including four variants within the 5' untranslated region (UTR). Here, we show that a 114 kb introgressed region containing Hnrnph1 and Rufy1 was sufficient to cause a decrease in MA-induced locomotor activity. Gene-level transcriptome analysis of striatal tissue from 114 kb congenics vs Hnrnph1 mutants identified a nearly perfect correlation of fold-change in expression for those differentially expressed genes that were common to both mouse lines, indicating functionally similar effects on the transcriptome and behavior. Exon-level analysis (including noncoding exons) revealed decreased 5' UTR usage of Hnrnph1 and immunoblot analysis identified a corresponding decrease in hnRNP H protein in 114 kb congenic mice. Molecular cloning of the Hnrnph1 5' UTR containing all four variants (but none of them individually) upstream of a reporter induced a decrease in reporter signal in both HEK293 and N2a cells, thus, identifying a set of QTVs underlying molecular regulation of Hnrnph1.

C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception.

Sensitivity to different pain modalities has a genetic basis that remains largely unknown. Employing closely related inbred mouse substrains can facilitate gene mapping of nociceptive behaviors in preclinical pain models. We previously reported enhanced sensitivity to acute thermal nociception in C57BL/6J (B6J) versus C57BL/6N (B6N) substrains. Here, we expanded on nociceptive phenotypes and observed an increase in formalin-induced inflammatory nociceptive behaviors and paw diameter in B6J versus B6N mice (Charles River Laboratories). No strain differences were observed in mechanical or thermal hypersensitivity or in edema following the Complete Freund's Adjuvant model of inflammatory pain, indicating specificity in the inflammatory nociceptive stimulus. In the chronic constrictive nerve injury, a model of neuropathic pain, no strain differences were observed in baseline mechanical threshold or in mechanical hypersensitivity up to one month post-chronic constrictive nerve injury. We replicated the enhanced thermal nociception in the 52.5°C hot plate test in B6J versus B6N mice from The Jackson Laboratory. Using a B6J × B6N-F2 cross (N = 164), we mapped a major quantitative trait locus underlying hot plate sensitivity to chromosome 7 that peaked at 26 Mb (log of the odds [LOD] = 3.81, p < 0.01; 8.74 Mb-36.50 Mb) that was more pronounced in males. Genes containing expression quantitative trait loci associated with the peak nociceptive marker that are implicated in pain and inflammation include Ryr1, Cyp2a5, Pou2f2, Clip3, Sirt2, Actn4, and Ltbp4 (false discovery rate < 0.05). Future studies involving positional cloning and gene editing will determine the quantitative trait gene(s) and potential pleiotropy of this locus across pain modalities.

MicroRNAs in CSF as prodromal biomarkers for Huntington disease in the PREDICT-HD study.

OBJECTIVE: To investigate the feasibility of microRNA (miRNA) levels in CSF as biomarkers for prodromal Huntington disease (HD). METHODS: miRNA levels were measured in CSF from 60 PREDICT-HD study participants using the HTG protocol. Using a CAG-Age Product score, 30 prodromal HD participants were selected based on estimated probability of imminent clinical diagnosis of HD (i.e., low, medium, high; n = 10/group). For comparison, participants already diagnosed (n = 15) and healthy controls (n = 15) were also selected. RESULTS: A total of 2,081 miRNAs were detected and 6 were significantly increased in the prodromal HD gene expansion carriers vs controls at false discovery rate q < 0.05 (miR-520f-3p, miR-135b-3p, miR-4317, miR-3928-5p, miR-8082, miR-140-5p). Evaluating the miRNA levels in each of the HD risk categories, all 6 revealed a pattern of increasing abundance from control to low risk, and from low risk to medium risk, which then leveled off from the medium to high risk and HD diagnosed groups. CONCLUSIONS: This study reports miRNAs as CSF biomarkers of prodromal and diagnosed HD. Importantly, miRNAs were detected in the prodromal HD groups furthest from diagnosis where treatments are likely to be most consequential and meaningful. The identification of potential biomarkers in the disease prodrome may prove useful in evaluating treatments that may postpone disease onset. CLINICALTRIALSGOV IDENTIFIER: NCT00051324.

Cytoplasmic FMR1-Interacting Protein 2 Is a Major Genetic Factor Underlying Binge Eating.

BACKGROUND: Eating disorders are lethal and heritable; however, the underlying genetic factors are unknown. Binge eating is a highly heritable trait associated with eating disorders that is comorbid with mood and substance use disorders. Therefore, understanding its genetic basis will inform therapeutic development that could improve several comorbid neuropsychiatric conditions. METHODS: We assessed binge eating in closely related C57BL/6 mouse substrains and in an F2 cross to identify quantitative trait loci associated with binge eating. We used gene targeting to validate candidate genetic factors. Finally, we used transcriptome analysis of the striatum via messenger RNA sequencing to identify the premorbid transcriptome and the binge-induced transcriptome to inform molecular mechanisms mediating binge eating susceptibility and establishment. RESULTS: C57BL/6NJ but not C57BL/6J mice showed rapid and robust escalation in palatable food consumption. We mapped a single genome-wide significant quantitative trait locus on chromosome 11 (logarithm of the odds = 7.4) to a missense mutation in cytoplasmic FMR1-interacting protein 2 (Cyfip2). We validated Cyfip2 as a major genetic factor underlying binge eating in heterozygous knockout mice on a C57BL/6N background that showed reduced binge eating toward a wild-type C57BL/6J-like level. Transcriptome analysis of premorbid genetic risk identified the enrichment terms morphine addiction and retrograde endocannabinoid signaling, whereas binge eating resulted in the downregulation of a gene set enriched for decreased myelination, oligodendrocyte differentiation, and expression. CONCLUSIONS: We identified Cyfip2 as a major significant genetic factor underlying binge eating and provide a behavioral paradigm for future genome-wide association studies in populations with increased genetic complexity.

Hnrnph1 Is A Quantitative Trait Gene for Methamphetamine Sensitivity.

Psychostimulant addiction is a heritable substance use disorder; however its genetic basis is almost entirely unknown. Quantitative trait locus (QTL) mapping in mice offers a complementary approach to human genome-wide association studies and can facilitate environment control, statistical power, novel gene discovery, and neurobiological mechanisms. We used interval-specific congenic mouse lines carrying various segments of chromosome 11 from the DBA/2J strain on an isogenic C57BL/6J background to positionally clone a 206 kb QTL (50,185,512-50,391,845 bp) that was causally associated with a reduction in the locomotor stimulant response to methamphetamine (2 mg/kg, i.p.; DBA/2J < C57BL/6J)-a non-contingent, drug-induced behavior that is associated with stimulation of the dopaminergic reward circuitry. This chromosomal region contained only two protein coding genes-heterogeneous nuclear ribonucleoprotein, H1 (Hnrnph1) and RUN and FYVE domain-containing 1 (Rufy1). Transcriptome analysis via mRNA sequencing in the striatum implicated a neurobiological mechanism involving a reduction in mesolimbic innervation and striatal neurotransmission. For instance, Nr4a2 (nuclear receptor subfamily 4, group A, member 2), a transcription factor crucial for midbrain dopaminergic neuron development, exhibited a 2.1-fold decrease in expression (DBA/2J < C57BL/6J; p 4.2 x 10-15). Transcription activator-like effector nucleases (TALENs)-mediated introduction of frameshift deletions in the first coding exon of Hnrnph1, but not Rufy1, recapitulated the reduced methamphetamine behavioral response, thus identifying Hnrnph1 as a quantitative trait gene for methamphetamine sensitivity. These results define a novel contribution of Hnrnph1 to neurobehavioral dysfunction associated with dopaminergic neurotransmission. These findings could have implications for understanding the genetic basis of methamphetamine addiction in humans and the development of novel therapeutics for prevention and treatment of substance abuse and possibly other psychiatric disorders.