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.

Urocortin-3 neurons in the mouse perifornical area promote infant-directed neglect and aggression.

While recent studies have uncovered dedicated neural pathways mediating the positive control of parenting, the regulation of infant-directed aggression and how it relates to adult-adult aggression is poorly understood. Here we show that urocortin-3 (Ucn3)-expressing neurons in the hypothalamic perifornical area (PeFAUcn3) are activated during infant-directed attacks in males and females, but not other behaviors. Functional manipulations of PeFAUcn3 neurons demonstrate the role of this population in the negative control of parenting in both sexes. PeFAUcn3 neurons receive input from areas associated with vomeronasal sensing, stress, and parenting, and send projections to hypothalamic and limbic areas. Optogenetic activation of PeFAUcn3 axon terminals in these regions triggers various aspects of infant-directed agonistic responses, such as neglect, repulsion, and aggression. Thus, PeFAUcn3 neurons emerge as a dedicated circuit component controlling infant-directed neglect and aggression, providing a new framework to understand the positive and negative regulation of parenting in health and disease.

The Tug1 lncRNA locus is essential for male fertility.

BACKGROUND: Several long noncoding RNAs (lncRNAs) have been shown to function as components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, studying lncRNA function has been a challenge due to their diverse biological roles and because lncRNA loci can contain multiple molecular modes that may exert function. RESULTS: We previously generated and characterized a cohort of 20 lncRNA loci knockout mice. Here, we extend this initial study and provide a more detailed analysis of the highly conserved lncRNA locus, taurine-upregulated gene 1 (Tug1). We report that Tug1-knockout male mice are sterile with underlying defects including a low number of sperm and abnormal sperm morphology. Because lncRNA loci can contain multiple modes of action, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region have any activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors two distinct noncoding regulatory activities, as a cis-DNA repressor that regulates neighboring genes and as a lncRNA that can regulate genes by a trans-based function. We also show that Tug1 contains an evolutionary conserved open reading frame that when overexpressed produces a stable protein which impacts mitochondrial membrane potential, suggesting a potential third coding function. CONCLUSIONS: Our results reveal an essential role for the Tug1 locus in male fertility and uncover evidence for distinct molecular modes in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Transitions in paternal social status predict patterns of offspring growth and metabolic transcription.

One type of parental effect occurs when changes in parental phenotype or environment trigger changes to offspring phenotype. Such nongenetic parental effects can be precisely triggered in response to an environmental cue in time-locked fashion, or in other cases, persist for multiple generations after the cue has been removed, suggesting multiple timescales of action. For parental effects to serve as reliable signals of current environmental conditions, they should be reversible, such that when cues change, offspring phenotypes change in accordance. Social hierarchy is a prevalent feature of the environment, and current parental social status could signal the environment in which offspring will be born. Here, we sought to address parental effects of social status and their timescale of action in mice. We show that territorial competition in seminatural environments affects offspring growth. Although dominant males are not heavier than nondominant or control males, they produce faster growing offspring, particularly sons. The timing, effect-size, and sex-specificity of this association are modulated by maternal social experience. We show that a change in paternal social status is sufficient to modulate offspring weight: from one breeding cycle to the next, status-ascending males produce heavier sons than before, and status-descending males produce lighter sons than before. Current paternal status is also highly predictive of liver transcription in sons, including molecular pathways controlling oxidative phosphorylation and iron metabolism. These results are consistent with a parental effect of social experience, although alternative explanations are considered. In summary, changes in paternal social status are associated with changes in offspring growth and metabolism.

Murine single-cell RNA-seq reveals cell-identity- and tissue-specific trajectories of aging.

Aging is a pleiotropic process affecting many aspects of mammalian physiology. Mammals are composed of distinct cell type identities and tissue environments, but the influence of these cell identities and environments on the trajectory of aging in individual cells remains unclear. Here, we performed single-cell RNA-seq on >50,000 individual cells across three tissues in young and old mice to allow for direct comparison of aging phenotypes across cell types. We found transcriptional features of aging common across many cell types, as well as features of aging unique to each type. Leveraging matrix factorization and optimal transport methods, we found that both cell identities and tissue environments exert influence on the trajectory and magnitude of aging, with cell identity influence predominating. These results suggest that aging manifests with unique directionality and magnitude across the diverse cell identities in mammals.

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.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata.

BACKGROUND: Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. RESULTS: A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. CONCLUSIONS: The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

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.

Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region.

The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei-including the molecular identity, spatial organization, and function of distinct cell types-is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.

Correction to: A window into extreme longevity; the circulating metabolomic signature of the naked mole-rat, a mammal that shows negligible senescence.

The original version of this article unfortunately contained an error.

Functional circuit architecture underlying parental behaviour.

Parenting is essential for the survival and wellbeing of mammalian offspring. However, we lack a circuit-level understanding of how distinct components of this behaviour are coordinated. Here we investigate how galanin-expressing neurons in the medial preoptic area (MPOAGal) of the hypothalamus coordinate motor, motivational, hormonal and social aspects of parenting in mice. These neurons integrate inputs from a large number of brain areas and the activation of these inputs depends on the animal's sex and reproductive state. Subsets of MPOAGal neurons form discrete pools that are defined by their projection sites. While the MPOAGal population is active during all episodes of parental behaviour, individual pools are tuned to characteristic aspects of parenting. Optogenetic manipulation of MPOAGal projections mirrors this specificity, affecting discrete parenting components. This functional organization, reminiscent of the control of motor sequences by pools of spinal cord neurons, provides a new model for how discrete elements of a social behaviour are generated at the circuit level.

A window into extreme longevity; the circulating metabolomic signature of the naked mole-rat, a mammal that shows negligible senescence.

Mouse-sized naked mole-rats (Heterocephalus glaber), unlike other mammals, do not conform to Gompertzian laws of age-related mortality; adults show no age-related change in mortality risk. Moreover, we observe negligible hallmarks of aging with well-maintained physiological and molecular functions, commonly altered with age in other species. We questioned whether naked mole-rats, living an order of magnitude longer than laboratory mice, exhibit different plasma metabolite profiles, which could then highlight novel mechanisms or targets involved in disease and longevity. Using a comprehensive, unbiased metabolomics screen, we observe striking inter-species differences in amino acid, peptide, and lipid metabolites. Low circulating levels of specific amino acids, particularly those linked to the methionine pathway, resemble those observed during the fasting period at late torpor in hibernating ground squirrels and those seen in longer-lived methionine-restricted rats. These data also concur with metabolome reports on long-lived mutant mice, including the Ames dwarf mice and calorically restricted mice, as well as fruit flies, and even show similarities to circulating metabolite differences observed in young human adults when compared to older humans. During evolution, some of these beneficial nutrient/stress response pathways may have been positively selected in the naked mole-rat. These observations suggest that interventions that modify the aging metabolomic profile to a more youthful one may enable people to lead healthier and longer lives.

Benign and tumor parenchyma metabolomic profiles affect compensatory renal growth in renal cell carcinoma surgical patients.

BACKGROUND AND OBJECTIVES: Pre-operative kidney volume is an independent predictor of glomerular filtration rate in renal cell carcinoma patients. Compensatory renal growth (CRG) can ensue prior to nephrectomy in parallel to tumor growth and benign parenchyma loss. We aimed to test whether renal metabolite abundances significantly associate with CRG, suggesting a causative relationship. DESIGN, SETTING, PARTICIPANTS, AND MEASUREMENTS: Tissue metabolomics data from 49 patients, with a median age of 60 years, were previously collected and the pre-operative fold-change of their contra to ipsi-lateral benign kidney volume served as a surrogate for their CRG. Contra-lateral kidney volume fold-change within a 3.3 +/- 2.1 years follow-up interval was used as a surrogate for long-term CRG. Using a multivariable statistical model, we identified metabolites whose abundances significantly associate with CRG. RESULTS: Our analysis found 13 metabolites in the benign (e.g. L-urobilin, Variable Influence in Projection, VIP, score = 3.02, adjusted p = 0.017) and 163 metabolites in the malignant (e.g. 3-indoxyl-sulfate, VIP score = 1.3, adjusted p = 0.044) tissues that significantly associate with CRG. Benign/tumor fold change in metabolite abundances revealed three additional metabolites with that significantly positively associate with CRG (e.g. p-cresol sulfate, VIP score = 2.945, adjusted p = 0.033). At the pathway level, we show that fatty-acid oxidation is highly enriched with metabolites whose benign tissue abundances strongly positively associate with CRG, both pre-operatively and long term, whereas in the tumor tissue significant enrichment of dipeptides and benzoate (positive association), glycolysis/gluconeogenesis, lysolipid and nucleotide sugar pentose (negative associations) sub-pathways, were observed. CONCLUSION: These data suggest that specific biological processes in the benign as well as in the tumor parenchyma strongly influence compensatory renal growth.

The Multipartite Mitochondrial Genome of Enteromyxum leei (Myxozoa): Eight Fast-Evolving Megacircles.

Myxozoans are a large group of poorly characterized cnidarian parasites. To gain further insight into their evolution, we sequenced the mitochondrial (mt) genome of Enteromyxum leei and reevaluate the mt genome structure of Kudoa iwatai. Although the typical animal mt genome is a compact, 13-25 kb, circular chromosome, the mt genome of E. leei was found to be fragmented into eight circular chromosomes of ∼23 kb, making it the largest described animal mt genome. Each chromosome was found to harbor a large noncoding region (∼15 kb), nearly identical between chromosomes. The protein coding genes show an unusually high rate of sequence evolution and possess little similarity to their cnidarian homologs. Only five protein coding genes could be identified and no tRNA genes. Surprisingly, the mt genome of K. iwatai was also found to be composed of two chromosomes. These observations confirm the remarkable plasticity of myxozoan mt genomes.

New Perspectives on Genomic Imprinting, an Essential and Multifaceted Mode of Epigenetic Control in the Developing and Adult Brain.

Mammalian evolution entailed multiple innovations in gene regulation, including the emergence of genomic imprinting, an epigenetic regulation leading to the preferential expression of a gene from its maternal or paternal allele. Genomic imprinting is highly prevalent in the brain, yet, until recently, its central roles in neural processes have not been fully appreciated. Here, we provide a comprehensive survey of adult and developmental brain functions influenced by imprinted genes, from neural development and wiring to synaptic function and plasticity, energy balance, social behaviors, emotions, and cognition. We further review the widespread identification of parental biases alongside monoallelic expression in brain tissues, discuss their potential roles in dosage regulation of key neural pathways, and suggest possible mechanisms underlying the dynamic regulation of imprinting in the brain. This review should help provide a better understanding of the significance of genomic imprinting in the normal and pathological brain of mammals including humans.

Genomic insights into the evolutionary origin of Myxozoa within Cnidaria.

The Myxozoa comprise over 2,000 species of microscopic obligate parasites that use both invertebrate and vertebrate hosts as part of their life cycle. Although the evolutionary origin of myxozoans has been elusive, a close relationship with cnidarians, a group that includes corals, sea anemones, jellyfish, and hydroids, is supported by some phylogenetic studies and the observation that the distinctive myxozoan structure, the polar capsule, is remarkably similar to the stinging structures (nematocysts) in cnidarians. To gain insight into the extreme evolutionary transition from a free-living cnidarian to a microscopic endoparasite, we analyzed genomic and transcriptomic assemblies from two distantly related myxozoan species, Kudoa iwatai and Myxobolus cerebralis, and compared these to the transcriptome and genome of the less reduced cnidarian parasite, Polypodium hydriforme. A phylogenomic analysis, using for the first time to our knowledge, a taxonomic sampling that represents the breadth of myxozoan diversity, including four newly generated myxozoan assemblies, confirms that myxozoans are cnidarians and are a sister taxon to P. hydriforme. Estimations of genome size reveal that myxozoans have one of the smallest reported animal genomes. Gene enrichment analyses show depletion of expressed genes in categories related to development, cell differentiation, and cell-cell communication. In addition, a search for candidate genes indicates that myxozoans lack key elements of signaling pathways and transcriptional factors important for multicellular development. Our results suggest that the degeneration of the myxozoan body plan from a free-living cnidarian to a microscopic parasitic cnidarian was accompanied by extreme reduction in genome size and gene content.

Regulation of alternative splicing at the single-cell level.

Alternative splicing is a key cellular mechanism for generating distinct isoforms, whose relative abundances regulate critical cellular processes. It is therefore essential that inclusion levels of alternative exons be tightly regulated. However, how the precision of inclusion levels among individual cells is governed is poorly understood. Using single-cell gene expression, we show that the precision of inclusion levels of alternative exons is determined by the degree of evolutionary conservation at their flanking intronic regions. Moreover, the inclusion levels of alternative exons, as well as the expression levels of the transcripts harboring them, also contribute to this precision. We further show that alternative exons whose inclusion levels are considerably changed during stem cell differentiation are also subject to this regulation. Our results imply that alternative splicing is coordinately regulated to achieve accuracy in relative isoform abundances and that such accuracy may be important in determining cell fate.

Quantitative and functional interrogation of parent-of-origin allelic expression biases in the brain.

The maternal and paternal genomes play different roles in mammalian brains as a result of genomic imprinting, an epigenetic regulation leading to differential expression of the parental alleles of some genes. Here we investigate genomic imprinting in the cerebellum using a newly developed Bayesian statistical model that provides unprecedented transcript-level resolution. We uncover 160 imprinted transcripts, including 41 novel and independently validated imprinted genes. Strikingly, many genes exhibit parentally biased--rather than monoallelic--expression, with different magnitudes according to age, organ, and brain region. Developmental changes in parental bias and overall gene expression are strongly correlated, suggesting combined roles in regulating gene dosage. Finally, brain-specific deletion of the paternal, but not maternal, allele of the paternally-biased Bcl-x, (Bcl2l1) results in loss of specific neuron types, supporting the functional significance of parental biases. These findings reveal the remarkable complexity of genomic imprinting, with important implications for understanding the normal and diseased brain.

Diversity and evolution of myxozoan minicollagens and nematogalectins.

BACKGROUND: Myxozoa are a diverse group of metazoan parasites with a very simple organization, which has for decades eluded their evolutionary origin. Their most prominent and characteristic feature is the polar capsule: a complex intracellular structure of the myxozoan spore, which plays a role in host infection. Striking morphological similarities have been found between myxozoan polar capsules and nematocysts, the stinging structures of cnidarians (corals, sea anemones and jellyfish) leading to the suggestion that Myxozoa and Cnidaria share a more recent common ancestry. This hypothesis has recently been supported by phylogenomic evidence and by the identification of a nematocyst specific minicollagen gene in the myxozoan Tetracapsuloides bryosalmonae. Here we searched genomes and transcriptomes of several myxozoan taxa for the presence of additional cnidarian specific genes and characterized these genes within a phylogenetic context. RESULTS: Illumina assemblies of transcriptome or genome data of three myxozoan species (Enteromyxum leei, Kudoa iwatai, and Sphaeromyxa zaharoni) and of the enigmatic cnidarian parasite Polypodium hydriforme (Polypodiozoa) were mined using tBlastn searches with nematocyst-specific proteins as queries. Several orthologs of nematogalectins and minicollagens were identified. Our phylogenetic analyses indicate that myxozoans possess three distinct minicollagens. We found that the cnidarian repertoire of nematogalectins is more complex than previously thought and we identified additional members of the nematogalectin family. Cnidarians were found to possess four nematogalectin/ nematogalectin-related genes, while in myxozoans only three genes could be identified. CONCLUSIONS: Our results demonstrate that myxozoans possess a diverse array of genes that are taxonomically restricted to Cnidaria. Characterization of these genes provide compelling evidence that polar capsules and nematocysts are homologous structures and that myxozoans are highly degenerate cnidarians. The diversity of minicollagens was higher than previously thought, with the presence of three minicollagen genes in myxozoans. Our phylogenetic results suggest that the different myxozoan sequences are the results of ancient divergences within Cnidaria and not of recent specializations of the polar capsule. For both minicollagen and nematogalectin, our results show that myxozoans possess less gene copies than their cnidarian counter parts, suggesting that the polar capsule gene repertoire was simplified with their reduced body plan.

Deep sequencing of mixed total DNA without barcodes allows efficient assembly of highly plastic ascidian mitochondrial genomes.

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.