Single-cell transcriptomics of the naked mole-rat reveals unexpected features of mammalian immunity.

The immune system comprises a complex network of specialized cells that protects against infection, eliminates cancerous cells, and regulates tissue repair, thus serving a critical role in homeostasis, health span, and life span. The subterranean-dwelling naked mole-rat (NM-R; Heterocephalus glaber) exhibits prolonged life span relative to its body size, is unusually cancer resistant, and manifests few physiological or molecular changes with advancing age. We therefore hypothesized that the immune system of NM-Rs evolved unique features that confer enhanced cancer immunosurveillance and prevent the age-associated decline in homeostasis. Using single-cell RNA-sequencing (scRNA-seq) we mapped the immune system of the NM-R and compared it to that of the short-lived, cancer-prone mouse. In contrast to the mouse, we find that the NM-R immune system is characterized by a high myeloid-to-lymphoid cell ratio that includes a novel, lipopolysaccharide (LPS)-responsive, granulocyte cell subset. Surprisingly, we also find that NM-Rs lack canonical natural killer (NK) cells. Our comparative genomics analyses support this finding, showing that the NM-R genome lacks an expanded gene family that controls NK cell function in several other species. Furthermore, we reconstructed the evolutionary history that likely led to this genomic state. The NM-R thus challenges our current understanding of mammalian immunity, favoring an atypical, myeloid-biased mode of innate immunosurveillance, which may contribute to its remarkable health span.

Unraveling the message: insights into comparative genomics of the naked mole-rat.

Animals have evolved to survive, and even thrive, in different environments. Genetic adaptations may have indirectly created phenotypes that also resulted in a longer lifespan. One example of this phenomenon is the preternaturally long-lived naked mole-rat. This strictly subterranean rodent tolerates hypoxia, hypercapnia, and soil-based toxins. Naked mole-rats also exhibit pronounced resistance to cancer and an attenuated decline of many physiological characteristics that often decline as mammals age. Elucidating mechanisms that give rise to their unique phenotypes will lead to better understanding of subterranean ecophysiology and biology of aging. Comparative genomics could be a useful tool in this regard. Since the publication of a naked mole-rat genome assembly in 2011, analyses of genomic and transcriptomic data have enabled a clearer understanding of mole-rat evolutionary history and suggested molecular pathways (e.g., NRF2-signaling activation and DNA damage repair mechanisms) that may explain the extraordinarily longevity and unique health traits of this species. However, careful scrutiny and re-analysis suggest that some identified features result from incorrect or imprecise annotation and assembly of the naked mole-rat genome: in addition, some of these conclusions (e.g., genes involved in cancer resistance and hairlessness) are rejected when the analysis includes additional, more closely related species. We describe how the combination of better study design, improved genomic sequencing techniques, and new bioinformatic and data analytical tools will improve comparative genomics and ultimately bridge the gap between traditional model and nonmodel organisms.

Evolution of T cells in the cancer-resistant naked mole-rat.

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αßT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

Novel insights from a multiomics dissection of the Hayflick limit.

The process wherein dividing cells exhaust proliferative capacity and enter into replicative senescence has become a prominent model for cellular aging in vitro. Despite decades of study, this cellular state is not fully understood in culture and even much less so during aging. Here, we revisit Leonard Hayflick's original observation of replicative senescence in WI-38 human lung fibroblasts equipped with a battery of modern techniques including RNA-seq, single-cell RNA-seq, proteomics, metabolomics, and ATAC-seq. We find evidence that the transition to a senescent state manifests early, increases gradually, and corresponds to a concomitant global increase in DNA accessibility in nucleolar and lamin associated domains. Furthermore, we demonstrate that senescent WI-38 cells acquire a striking resemblance to myofibroblasts in a process similar to the epithelial to mesenchymal transition (EMT) that is regulated by t YAP1/TEAD1 and TGF-ß2. Lastly, we show that verteporfin inhibition of YAP1/TEAD1 activity in aged WI-38 cells robustly attenuates this gene expression program.

BTI1, an azoreductase with pH-dependent substrate specificity.

The group II azoreductase BTI1 utilizes NADPH to directly cleave azo bonds in water-soluble azo dyes, including quenchers of fluorescence. Unexpectedly, optimal reduction was dye specific, ranging from a pH of <5.5 for Janus green B, to pH 6.0 for methyl red, methyl orange, and BHQ-10, to pH >8.3 for flame orange.

Differentiation reveals latent features of aging and an energy barrier in murine myogenesis.

Skeletal muscle experiences a decline in lean mass and regenerative potential with age, in part due to intrinsic changes in progenitor cells. However, it remains unclear how age-related changes in progenitors manifest across a differentiation trajectory. Here, we perform single-cell RNA sequencing (RNA-seq) on muscle mononuclear cells from young and aged mice and profile muscle stem cells (MuSCs) and fibro-adipose progenitors (FAPs) after differentiation. Differentiation increases the magnitude of age-related change in MuSCs and FAPs, but it also masks a subset of age-related changes present in progenitors. Using a dynamical systems approach and RNA velocity, we find that aged MuSCs follow the same differentiation trajectory as young cells but stall in differentiation near a commitment decision. Our results suggest that differentiation reveals latent features of aging and that fate commitment decisions are delayed in aged myogenic cells in vitro.

Solo: Doublet Identification in Single-Cell RNA-Seq via Semi-Supervised Deep Learning.

Single-cell RNA sequencing (scRNA-seq) measurements of gene expression enable an unprecedented high-resolution view into cellular state. However, current methods often result in two or more cells that share the same cell-identifying barcode; these "doublets" violate the fundamental premise of single-cell technology and can lead to incorrect inferences. Here, we describe Solo, a semi-supervised deep learning approach that identifies doublets with greater accuracy than existing methods. Solo embeds cells unsupervised using a variational autoencoder and then appends a feed-forward neural network layer to the encoder to form a supervised classifier. We train this classifier to distinguish simulated doublets from the observed data. Solo can be applied in combination with experimental doublet detection methods to further purify scRNA-seq data to true single cells. It is freely available from https://github.com/calico/solo. A record of this paper's transparent peer review process is included in the Supplemental Information.

Fully Phased Sequence of a Diploid Human Genome Determined de Novo from the DNA of a Single Individual.

In recent years, improved sequencing technology and computational tools have made de novo genome assembly more accessible. Many approaches, however, generate either an unphased or only partially resolved representation of a diploid genome, in which polymorphisms are detected but not assigned to one or the other of the homologous chromosomes. Yet chromosomal phase information is invaluable for the understanding of phenotypic trait inheritance in the cases of compound heterozygosity, allele-specific expression or cis-acting variants. Here we use a combination of tools and sequencing technologies to generate a de novo diploid assembly of the human primary cell line WI-38. First, data from PacBio single molecule sequencing and Bionano Genomics optical mapping were combined to generate an unphased assembly. Next, 10x Genomics linked reads were combined with the hybrid assembly to generate a partially phased assembly. Lastly, we developed and optimized methods to use short-read (Illumina) sequencing of flow cytometry-sorted metaphase chromosomes to provide phase information. The final genome assembly was almost fully (94%) phased with the addition of approximately 2.5-fold coverage of Illumina data from the sequenced metaphase chromosomes. The diploid nature of the final de novo genome assembly improved the resolution of structural variants between the WI-38 genome and the human reference genome. The phased WI-38 sequence data are available for browsing and download at wi38.research.calicolabs.com. Our work shows that assembling a completely phased diploid genome de novo from the DNA of a single individual is now readily achievable.

eIF2B activator prevents neurological defects caused by a chronic integrated stress response.

The integrated stress response (ISR) attenuates the rate of protein synthesis while inducing expression of stress proteins in cells. Various insults activate kinases that phosphorylate the GTPase eIF2 leading to inhibition of its exchange factor eIF2B. Vanishing White Matter (VWM) is a neurological disease caused by eIF2B mutations that, like phosphorylated eIF2, reduce its activity. We show that introduction of a human VWM mutation into mice leads to persistent ISR induction in the central nervous system. ISR activation precedes myelin loss and development of motor deficits. Remarkably, long-term treatment with a small molecule eIF2B activator, 2BAct, prevents all measures of pathology and normalizes the transcriptome and proteome of VWM mice. 2BAct stimulates the remaining activity of mutant eIF2B complex in vivo, abrogating the maladaptive stress response. Thus, 2BAct-like molecules may provide a promising therapeutic approach for VWM and provide relief from chronic ISR induction in a variety of disease contexts.

Detection of minute virus of mice using real time quantitative PCR in assessment of virus clearance during the purification of Mammalian cell substrate derived biotherapeutics.

A real time quantitative PCR assay has been developed for detecting minute virus of mice (MVM). This assay directly quantifies PCR product by monitoring the increase of fluorescence intensity emitted during enzymatic hydrolysis of an oligonucleotide probe labelled covalently with fluorescent reporting and quenching dyes via Taq polymerase 5'-->3' exonuclease activity. The quantity of MVM DNA molecules in the samples was determined using a known amount of MVM standard control DNA fragment cloned into a plasmid (pCR-MVM). We have demonstrated that MVM TaqMan PCR assay is approximately 1000-fold more sensitive than the microplate infectivity assay with the lowest detection limit of approximately one particle per reaction. The reliable detection range is within 100 to 10(9) molecules per reaction with high reproducibility. The intra assay variation is <2.5%, and the inter assays variation is <6.5% when samples contain >100 particles/assay. When we applied the TaqMan PCR to MVM clearance studies done by column chromatography or normal flow viral filtration, we found that the virus removal factors were similar to that of virus infectivity assay. It takes about a day to complete entire assay processes, thus, the TaqMan PCR assay is at least 10-fold faster than the infectivity assay. Therefore, we concluded that this fast, specific, sensitive, and robust assay could replace the infectivity assay for virus clearance evaluation.