Cellular development and evolution of the mammalian cerebellum.

The expansion of the neocortex, a hallmark of mammalian evolution1,2, was accompanied by an increase in cerebellar neuron numbers3. However, little is known about the evolution of the cellular programmes underlying the development of the cerebellum in mammals. In this study we generated single-nucleus RNA-sequencing data for around 400,000 cells to trace the development of the cerebellum from early neurogenesis to adulthood in human, mouse and the marsupial opossum. We established a consensus classification of the cellular diversity in the developing mammalian cerebellum and validated it by spatial mapping in the fetal human cerebellum. Our cross-species analyses revealed largely conserved developmental dynamics of cell-type generation, except for Purkinje cells, for which we observed an expansion of early-born subtypes in the human lineage. Global transcriptome profiles, conserved cell-state markers and gene-expression trajectories across neuronal differentiation show that cerebellar cell-type-defining programmes have been overall preserved for at least 160 million years. However, we also identified many orthologous genes that gained or lost expression in cerebellar neural cell types in one of the species or evolved new expression trajectories during neuronal differentiation, indicating widespread gene repurposing at the cell-type level. In sum, our study unveils shared and lineage-specific gene-expression programmes governing the development of cerebellar cells and expands our understanding of mammalian brain evolution.

The molecular evolution of spermatogenesis across mammals.

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals1-6, probably owing to the evolutionary pressure on males to be reproductively successful7. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome.

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Organoid single-cell genomic atlas uncovers human-specific features of brain development.

The human brain has undergone substantial change since humans diverged from chimpanzees and the other great apes1,2. However, the genetic and developmental programs that underlie this divergence are not fully understood. Here we have analysed stem cell-derived cerebral organoids using single-cell transcriptomics and accessible chromatin profiling to investigate gene-regulatory changes that are specific to humans. We first analysed cell composition and reconstructed differentiation trajectories over the entire course of human cerebral organoid development from pluripotency, through neuroectoderm and neuroepithelial stages, followed by divergence into neuronal fates within the dorsal and ventral forebrain, midbrain and hindbrain regions. Brain-region composition varied in organoids from different iPSC lines, but regional gene-expression patterns remained largely reproducible across individuals. We analysed chimpanzee and macaque cerebral organoids and found that human neuronal development occurs at a slower pace relative to the other two primates. Using pseudotemporal alignment of differentiation paths, we found that human-specific gene expression resolved to distinct cell states along progenitor-to-neuron lineages in the cortex. Chromatin accessibility was dynamic during cortex development, and we identified divergence in accessibility between human and chimpanzee that correlated with human-specific gene expression and genetic change. Finally, we mapped human-specific expression in adult prefrontal cortex using single-nucleus RNA sequencing analysis and identified developmental differences that persist into adulthood, as well as cell-state-specific changes that occur exclusively in the adult brain. Our data provide a temporal cell atlas of great ape forebrain development, and illuminate dynamic gene-regulatory features that are unique to humans.

Gene expression across mammalian organ development.

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, rhesus macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified correspondences of developmental stages across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex.

BACKGROUND: Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques. RESULTS: We found a discrete functional association, in which H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while H3K4me3HP loss played a vital role in synaptic activity. Moreover, H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by H3K4me3HP and H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes. CONCLUSIONS: Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

Single-cell-resolution transcriptome map of human, chimpanzee, bonobo, and macaque brains.

Identification of gene expression traits unique to the human brain sheds light on the molecular mechanisms underlying human evolution. Here, we searched for uniquely human gene expression traits by analyzing 422 brain samples from humans, chimpanzees, bonobos, and macaques representing 33 anatomical regions, as well as 88,047 cell nuclei composing three of these regions. Among 33 regions, cerebral cortex areas, hypothalamus, and cerebellar gray and white matter evolved rapidly in humans. At the cellular level, astrocytes and oligodendrocyte progenitors displayed more differences in the human evolutionary lineage than the neurons. Comparison of the bulk tissue and single-nuclei sequencing revealed that conventional RNA sequencing did not detect up to two-thirds of cell-type-specific evolutionary differences.

Untargeted Lipidomics after D2O Administration Reveals the Turnover Rate of Individual Lipids in Various Organs of Living Organisms.

The administration of low doses of D2O to living organisms was used for decades for the investigation of metabolic pathways and for the measurement of the turnover rate for specific compounds. Usually, the investigation of the deuterium uptake in lipids is performed by measuring the deuteration level of the palmitic acid residue using GC-MS instruments, and to our knowledge, the application of the modern untargeted LC-MS/MS lipidomics approaches was only reported a few times. Here, we investigated the deuterium uptake for >500 lipids for 13 organs and body liquids of mice (brain, lung, heart, liver, kidney, spleen, plasma, urine, etc.) after 4 days of 100% D2O administration. The maximum deuteration level was observed in the liver, plasma, and lung, while in the brain and heart, the deuteration level was lower. Using MS/MS, we demonstrated the incorporation of deuterium in palmitic and stearic fragments in lipids (PC, PE, TAG, PG, etc.) but not in the corresponding free forms. Our results were analyzed based on the metabolic pathways of lipids.

LINC00116 codes for a mitochondrial peptide linking respiration and lipid metabolism.

Genes coding for small peptides have been frequently misannotated as long noncoding RNA (lncRNA) genes. Here we have demonstrated that one such transcript is translated into a 56-amino-acid-long peptide conserved in chordates, corroborating the work published while this manuscript was under review. The Mtln peptide could be detected in mitochondria of mouse cell lines and tissues. In line with its mitochondrial localization, lack of the Mtln decreases the activity of mitochondrial respiratory chain complex I. Unlike the integral components and assembly factors of NADH:ubiquinone oxidoreductase, Mtln does not alter its enzymatic activity directly. Interaction of Mtln with NADH-dependent cytochrome b5 reductase stimulates complex I functioning most likely by providing a favorable lipid composition of the membrane. Study of Mtln illuminates the importance of small peptides, whose genes might frequently be misannotated as lncRNAs, for the control of vitally important cellular processes.

Time-Dependent Effect of Sciatic Nerve Injury on Rat Plasma Lipidome.

Neuropathic pain is a condition affecting the quality of life of a substantial part of the population, but biomarkers and treatment options are still limited. While this type of pain is caused by nerve damage, in which lipids play key roles, lipidome alterations related to nerve injury remain poorly studied. Here, we assessed blood lipidome alterations in a common animal model, the rat sciatic nerve crush injury. We analyzed alterations in blood lipid abundances between seven rats with nerve injury (NI) and eight control (CL) rats in a time-course experiment. For these rats, abundances of 377 blood lipid species were assessed at three distinct time points: immediately after, two weeks, and five weeks post injury. Although we did not detect significant differences between NI and CL at the first two time points, 106 lipids were significantly altered in NI five weeks post injury. At this time point, we found increased levels of triglycerides (TGs) and lipids containing esterified palmitic acid (16:0) in the blood plasma of NI animals. Lipids containing arachidonic acid (20:4), by contrast, were significantly decreased after injury, aligning with the crucial role of arachidonic acid reported for NI. Taken together, these results indicate delayed systematic alterations in fatty acid metabolism after nerve injury, potentially reflecting nerve tissue restoration dynamics.

Variation of microRNA expression in the human placenta driven by population identity and sex of the newborn.

BACKGROUND: Analysis of lymphocyte cell lines revealed substantial differences in the expression of mRNA and microRNA (miRNA) among human populations. The extent of such population-associated differences in actual human tissues remains largely unexplored. The placenta is one of the few solid human tissues that can be collected in substantial numbers in a controlled manner, enabling quantitative analysis of transient biomolecules such as RNA transcripts. Here, we analyzed microRNA (miRNA) expression in human placental samples derived from 36 individuals representing four genetically distinct human populations: African Americans, European Americans, South Asians, and East Asians. All samples were collected at the same hospital following a unified protocol, thus minimizing potential biases that might influence the results. RESULTS: Sequence analysis of the miRNA fraction yielded 938 annotated and 70 novel miRNA transcripts expressed in the placenta. Of them, 82 (9%) of annotated and 11 (16%) of novel miRNAs displayed quantitative expression differences among populations, generally reflecting reported genetic and mRNA-expression-based distances. Several co-expressed miRNA clusters stood out from the rest of the population-associated differences in terms of miRNA evolutionary age, tissue-specificity, and disease-association characteristics. Among three non-environmental influenced demographic parameters, the second largest contributor to miRNA expression variation after population was the sex of the newborn, with 32 miRNAs (3% of detected) exhibiting significant expression differences depending on whether the newborn was male or female. Male-associated miRNAs were evolutionarily younger and correlated inversely with the expression of target mRNA involved in neuron-related functions. In contrast, both male and female-associated miRNAs appeared to mediate different types of hormonal responses. Demographic factors further affected reported imprinted expression of 66 placental miRNAs: the imprinting strength correlated with the mother's weight, but not height. CONCLUSIONS: Our results showed that among 12 assessed demographic variables, population affiliation and fetal sex had a substantial influence on miRNA expression variation among human placental samples. The effect of newborn-sex-associated miRNA differences further led to expression inhibition of the target genes clustering in specific functional pathways. By contrast, population-driven miRNA differences might mainly represent neutral changes with minimal functional impacts.

Genotyping and lipid profiling of 601 cultivated sunflower lines reveals novel genetic determinants of oil fatty acid content.

BACKGROUND: Sunflower is an important oilseed crop domesticated in North America approximately 4000 years ago. During the last century, oil content in sunflower was under strong selection. Further improvement of oil properties achieved by modulating its fatty acid composition is one of the main directions in modern oilseed crop breeding. RESULTS: We searched for the genetic basis of fatty acid content variation by genotyping 601 inbred sunflower lines and assessing their lipid and fatty acid composition. Our genome-wide association analysis based on the genotypes for 15,483 SNPs and the concentrations of 23 fatty acids, including minor fatty acids, revealed significant genetic associations for eleven of them. Identified genomic regions included the loci involved in rare fatty acids variation on chromosomes 3 and 14, explaining up to 34.5% of the total variation of docosanoic acid (22:0) in sunflower oil. CONCLUSIONS: This is the first large scale implementation of high-throughput lipidomic profiling to sunflower germplasm characterization. This study contributes to the genetic characterization of Russian sunflower collections, which made a substantial contribution to the development of sunflower as the oilseed crop worldwide, and provides new insights into the genetic control of oil composition that can be implemented in future studies.

Structure-Preserving and Perceptually Consistent Approach for Visualization of Mass Spectrometry Imaging Datasets.

Mass spectrometry imaging (MSI) has become an important tool for 2D profiling of biological tissues, allowing for the visualization of individual compound distributions in the sample. Based on this information, it is possible to investigate the molecular organization within any particular tissue and detect abnormal regions (such as tumor regions) and many other biologically relevant phenomena. However, the large number of compounds present in the spectra hinders the productive analysis of large MSI datasets when utilizing standard tools. The heterogeneity of samples makes exploratory visualization (a presentation of the general idea of the molecular and structural organization of the inspected tissues) challenging. Here, we explore the application of various dimensionality reduction techniques that have been used extensively in the visualization of hyperspectral images and the MSI data specifically, such as principal component analysis, independent component analysis, non-negative matrix factorization, t-distributed stochastic neighbor embedding, and uniform manifold approximation and projection. Further, we propose a new approach based on a combination of structure preserving visualization with nonlinear manifold embedding of normalized spectral data. This way, we aim to preserve as much spatially overlapping signals as possible while augmenting them with information on compositional (spectral) variation. The proposed approach can be used for exploratory visualization of MSI datasets without prior deep chemical or histological knowledge of the sample. Thus, different datasets can be visually compared employing the proposed method. The proposed approach allowed for the clear visualization of the molecular layer, granular layer, and white matter in chimpanzee and macaque cerebellum slices.

Human-specific features of spatial gene expression and regulation in eight brain regions.

Molecular maps of the human brain alone do not inform us of the features unique to humans. Yet, the identification of these features is important for understanding both the evolution and nature of human cognition. Here, we approached this question by analyzing gene expression and H3K27ac chromatin modification data collected in eight brain regions of humans, chimpanzees, gorillas, a gibbon, and macaques. An analysis of spatial transcriptome trajectories across eight brain regions in four primate species revealed 1851 genes showing human-specific transcriptome differences in one or multiple brain regions, in contrast to 240 chimpanzee-specific differences. More than half of these human-specific differences represented elevated expression of genes enriched in neuronal and astrocytic markers in the human hippocampus, whereas the rest were enriched in microglial markers and displayed human-specific expression in several frontal cortical regions and the cerebellum. An analysis of the predicted regulatory interactions driving these differences revealed the role of transcription factors in species-specific transcriptome changes, and epigenetic modifications were linked to spatial expression differences conserved across species.

Lipidome alterations in human prefrontal cortex during development, aging, and cognitive disorders.

Lipids are essential to brain functions, yet they remain largely unexplored. Here we investigated the lipidome composition of prefrontal cortex gray matter in 396 cognitively healthy individuals with ages spanning 100 years, as well as 67 adult individuals diagnosed with autism (ASD), schizophrenia (SZ), and Down syndrome (DS). Of the 5024 detected lipids, 95% showed significant age-dependent concentration differences clustering into four temporal stages, and resulting in a gradual increase in membrane fluidity in individuals ranging from newborn to nonagenarian. Aging affects 14% of the brain lipidome with late-life changes starting predominantly at 50-55 years of age-a period of general metabolic transition. All three diseases alter the brain lipidome composition, leading-among other things-to a concentration decrease in glycerophospholipid metabolism and endocannabinoid signaling pathways. Lipid concentration decreases in SZ were further linked to genetic variants associated with disease, indicating the relevance of the lipidome changes to disease progression.

Long-Term Fluoxetine Administration Causes Substantial Lipidome Alteration of the Juvenile Macaque Brain.

Fluoxetine is an antidepressant commonly prescribed not only to adults but also to children for the treatment of depression, obsessive-compulsive disorder, and neurodevelopmental disorders. The adverse effects of the long-term treatment reported in some patients, especially in younger individuals, call for a detailed investigation of molecular alterations induced by fluoxetine treatment. Two-year fluoxetine administration to juvenile macaques revealed effects on impulsivity, sleep, social interaction, and peripheral metabolites. Here, we built upon this work by assessing residual effects of fluoxetine administration on the expression of genes and abundance of lipids and polar metabolites in the prelimbic cortex of 10 treated and 11 control macaques representing two monoamine oxidase A (MAOA) genotypes. Analysis of 8871 mRNA transcripts, 3608 lipids, and 1829 polar metabolites revealed substantial alterations of the brain lipid content, including significant abundance changes of 106 lipid features, accompanied by subtle changes in gene expression. Lipid alterations in the drug-treated animals were most evident for polyunsaturated fatty acids (PUFAs). A decrease in PUFAs levels was observed in all quantified lipid classes excluding sphingolipids, which do not usually contain PUFAs, suggesting systemic changes in fatty acid metabolism. Furthermore, the residual effect of the drug on lipid abundances was more pronounced in macaques carrying the MAOA-L genotype, mirroring reported behavioral effects of the treatment. We speculate that a decrease in PUFAs may be associated with adverse effects in depressive patients and could potentially account for the variation in individual response to fluoxetine in young people.

Lipidome analysis of milk composition in humans, monkeys, bovids, and pigs.

BACKGROUND: Lipids contained in milk are an essential source of energy and structural materials for a growing neonate. Furthermore, lipids' long-chain unsaturated fatty acid residues can directly participate in neonatal tissue formation. Here, we used untargeted mass spectrometric measurements to assess milk lipid composition in seven mammalian species: humans, two macaque species, cows, goats, yaks, and pigs. RESULTS: Analysis of the main milk lipid class, triacylglycerides (TAGs), revealed species-specific quantitative differences in the composition of fatty acid residues for each of seven species. Overall, differences in milk lipid composition reflect evolutionary distances among species, with each species group demonstrating specific lipidome features. Among them, human milk contained more medium and long-chain unsaturated fatty acids compared to other species, while pig milk was the most distinct, featuring the highest proportion of long-chain polyunsaturated fatty acids. CONCLUSIONS: We show that milk lipidome composition is dynamic across mammalian species, changed extensively in pigs, and contains features particular to humans.

Predominant patterns of splicing evolution on human, chimpanzee and macaque evolutionary lineages.

Although splicing is widespread and evolves rapidly among species, the mechanisms driving this evolution, as well as its functional implications, are not yet fully understood. We analyzed the evolution of splicing patterns based on transcriptome data from five tissues of humans, chimpanzees, rhesus macaques and mice. In total, 1526 exons and exon sets from 1236 genes showed significant splicing differences among primates. More than 60% of these differences represent constitutive-to-alternative exon transitions while an additional 25% represent changes in exon inclusion frequency. These two dominant evolutionary patterns have contrasting conservation, regulation and functional features. The sum of these features indicates that, despite their prevalence, constitutive-to-alternative exon transitions do not substantially contribute to long-term functional transcriptome changes. Conversely, changes in exon inclusion frequency appear to be functionally relevant, especially for changes taking place in the brain on the human evolutionary lineage.

Conservation, evolution, and regulation of splicing during prefrontal cortex development in humans, chimpanzees, and macaques.

Changes in splicing are known to affect the function and regulation of genes. We analyzed splicing events that take place during the postnatal development of the prefrontal cortex in humans, chimpanzees, and rhesus macaques based on data obtained from 168 individuals. Our study revealed that among the 38,822 quantified alternative exons, 15% are differentially spliced among species, and more than 6% splice differently at different ages. Mutations in splicing acceptor and/or donor sites might explain more than 14% of all splicing differences among species and up to 64% of high-amplitude differences. A reconstructed trans-regulatory network containing 21 RNA-binding proteins explains a further 4% of splicing variations within species. While most age-dependent splicing patterns are conserved among the three species, developmental changes in intron retention are substantially more pronounced in humans.

Lipidome Evolution in Mammalian Tissues.

Lipids are essential structural and functional components of cells. Little is known, however, about the evolution of lipid composition in different tissues. Here, we report a large-scale analysis of the lipidome evolution in six tissues of 32 species representing primates, rodents, and bats. While changes in genes' sequence and expression accumulate proportionally to the phylogenetic distances, <2% of the lipidome evolves this way. Yet, lipids constituting this 2% cluster in specific functions shared among all tissues. Among species, human show the largest amount of species-specific lipidome differences. Many of the uniquely human lipidome features localize in the brain cortex and cluster in specific pathways implicated in cognitive disorders.