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.

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.

Pheromone of grouped female mice impairs genome stability in male mice through stress-mediated pathways

Population density is known to affect the health and survival of many species, and is especially important for social animals. In mice, living in crowded conditions results in the disruption of social interactions, chronic stress, and immune and reproductive suppression; however, the underlying mechanisms remain unclear. Here, we investigated the role of chemosignals in the regulation of mouse physiology and behavior in response to social crowding. The pheromone 2,5-dimethylpyrazine (2,5-DMP), which is released by female mice in crowded conditions, induced aversion, glucocorticoid elevation and, when chronic, resulted in reproductive and immune suppression. 2,5-DMP olfaction induced genome destabilization in bone marrow cells in a stress-dependent manner, providing a plausible mechanism for crowding-induced immune dysfunction. Interestingly, the genome-destabilizing effect of 2,5-DMP was comparable to a potent mouse stressor (immobilization), and both stressors led to correlated expression changes in genes regulating cellular stress response. Thus, our findings demonstrate that, in mice, the health effects of crowding may be explained at least in part by chemosignals and also propose a significant role of stress and genome destabilization in the emergence of crowding effects.

QTL mapping of oleic acid content in modern VNIIMK sunflower (Helianthus annuus L.) lines by using GBS-based SNP map

Oleic acid is a monounsaturated fatty acid increasing oil oxidative stability. High content of oleic acid is thus a valuable trait in oilseed crops. Sunflower (Helianthus annuus L.) normally accumulates linoleic acid as a major fatty acid, but a mutant expressing a high oleic phenotype form was previously obtained by chemical mutagenesis and mapped on the sunflower genome. Several studies suggest the presence of additional genes involved in the control of the high content of oleic acid, with their expression possibly depending on the genetic background. To test this hypothesis, we performed a QTL mapping of the high oleic acid trait within two independent F2 crosses involving lines with contrasting oleic acid content from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We applied genotyping-by-sequencing (GBS) to construct single nucleotide polymorphism-based genetic maps and performed QTL mapping using quantitative and qualitative encoding for oleic acid content. Our results support the finding that the oleic acid content in the assessed crosses is controlled by one major effect locus. However, different dominant/recessive effects of the major locus were reported for both crosses. Additionally, a possible translocation between chromosome 7 and 14 was reported in one assessed cross. We defined a set of single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

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

The administration of low doses of D₂O 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% D₂O 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.

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.

Lipid Alteration Signature in the Blood Plasma of Individuals With Schizophrenia, Depression, and Bipolar Disorder

Importance

No clinically applicable diagnostic test exists for severe mental disorders. Lipids harbor potential as disease markers.

Objective

To define a reproducible profile of lipid alterations in the blood plasma of patients with schizophrenia (SCZ) independent of demographic and environmental variables and to investigate its specificity in association with other psychiatric disorders, ie, major depressive disorder (MDD) and bipolar disorder (BPD).

Design, Setting, and Participants

This was a multicohort case-control diagnostic analysis involving plasma samples from psychiatric patients and control individuals collected between July 17, 2009, and May 18, 2018. Study participants were recruited as consecutive and volunteer samples at multiple inpatient and outpatient mental health hospitals in Western Europe (Germany and Austria [DE-AT]), China (CN), and Russia (RU). Individuals with DSM-IV or International Statistical Classification of Diseases and Related Health Problems, Tenth Revision diagnoses of SCZ, MDD, BPD, or a first psychotic episode, as well as age- and sex-matched healthy controls without a mental health-related diagnosis were included in the study. Samples and data were analyzed from January 2018 to September 2020.

Main Outcomes and Measures

Plasma lipidome composition was assessed using liquid chromatography coupled with untargeted mass spectrometry.

Results

Blood lipid levels were assessed in 980 individuals (mean [SD] age, 36 [13] years; 510 male individuals [52%]) diagnosed with SCZ, BPD, MDD, or those with a first psychotic episode and in 572 controls (mean [SD] age, 34 [13] years; 323 male individuals [56%]). A total of 77 lipids were found to be significantly altered between those with SCZ (n = 436) and controls (n = 478) in all 3 sample cohorts. Alterations were consistent between cohorts (CN and RU: [Pearson correlation] r = 0.75; DE-AT and CN: r = 0.78; DE-AT and RU: r = 0.82; P < 10-38). A lipid-based predictive model separated patients with SCZ from controls with high diagnostic ability (area under the receiver operating characteristic curve = 0.86-0.95). Lipidome alterations in BPD and MDD, assessed in 184 and 256 individuals, respectively, were found to be similar to those of SCZ (BPD: r = 0.89; MDD: r = 0.92; P < 10-79). Assessment of detected alterations in individuals with a first psychotic episode, as well as patients with SCZ not receiving medication, demonstrated only limited association with medication restricted to particular lipids.

Conclusions and Relevance

In this study, SCZ was accompanied by a reproducible profile of plasma lipidome alterations, not associated with symptom severity, medication, and demographic and environmental variables, and largely shared with BPD and MDD. This lipid alteration signature may represent a trait marker of severe psychiatric disorders, indicating its potential to be transformed into a clinically applicable testing procedure.

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals^1-6, probably owing to the evolutionary pressure on males to be reproductively successful⁷. 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.

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.

Genetic mapping of loci involved in oil tocopherol composition control in Russian sunflower (Helianthus annuus L.) lines

Tocopherols are antioxidants that preserve oil lipids against oxidation and serve as a natural source of vitamin E in the human diet. Compared with other major oilseeds like rapeseed and soybean, sunflower (Helianthus annuus L.) exhibits low phenotypic diversity of tocopherol composition, both in wild and cultivated accessions from germplasm collections. Two major mutations that alter tocopherol composition were identified in genetic collections, and several studies suggested additional loci controlling tocopherol composition, with their expression possibly depending on the genetic background. In the present study, we performed QTL mapping of tocopherol composition in two independent F2 crosses between lines with contrasting tocopherol composition from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We used genotyping-bysequencing (GBS) to construct single nucleotide polymorphism-based genetic maps, and performed QTL mapping using quantitative and qualitative encoding for phenotypic traits. Our results support the notion that the tocopherol composition in the assessed crosses is controlled by two loci. We additionally selected and validated two single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity

Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence, or DiCo. We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell type composition and also by cell-autonomous expression changes within particular cell types.

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.

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.

Alternative splicing during mammalian organ development

Alternative splicing (AS) is pervasive in mammalian genomes, yet cross-species comparisons have been largely restricted to adult tissues and the functionality of most AS events remains unclear. We assessed AS patterns across pre- and postnatal development of seven organs in six mammals and a bird. Our analyses revealed that developmentally dynamic AS events, which are especially prevalent in the brain, are substantially more conserved than nondynamic ones. Cassette exons with increasing inclusion frequencies during development show the strongest signals of conserved and regulated AS. Newly emerged cassette exons are typically incorporated late in testis development, but those retained during evolution are predominantly brain specific. Our work suggests that an intricate interplay of programs controlling gene expression levels and AS is fundamental to organ development, especially for the brain and heart. In these regulatory networks, AS affords substantial functional diversification of genes through the generation of tissue- and time-specific isoforms from broadly expressed genes.

Reduced purine biosynthesis in humans after their divergence from Neandertals

We analyze the metabolomes of humans, chimpanzees, and macaques in muscle, kidney and three different regions of the brain. Although several compounds in amino acid metabolism occur at either higher or lower concentrations in humans than in the other primates, metabolites downstream of adenylosuccinate lyase, which catalyzes two reactions in purine synthesis, occur at lower concentrations in humans. This enzyme carries an amino acid substitution that is present in all humans today but absent in Neandertals. By introducing the modern human substitution into the genomes of mice, as well as the ancestral, Neandertal-like substitution into the genomes of human cells, we show that this amino acid substitution contributes to much or all of the reduction of de novo synthesis of purines in humans.

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.

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.

Genetic Characterization of Russian Rapeseed Collection and Association Mapping of Novel Loci Affecting Glucosinolate Content

Rapeseed is the second most common oilseed crop worldwide. While the start of rapeseed breeding in Russia dates back to the middle of the 20th century, its widespread cultivation began only recently. In contrast to the world’s rapeseed genetic variation, the genetic composition of Russian rapeseed lines remained unexplored. We have addressed this question by performing genome-wide genotyping of 90 advanced rapeseed accessions provided by the All-Russian Research Institute of Oil Crops (VNIIMK). Genome-wide genetic analysis demonstrated a clear difference between Russian rapeseed varieties and the rapeseed varieties from the rest of the world, including the European ones, indicating that rapeseed breeding in Russia proceeded in its own independent direction. Hence, genetic determinants of agronomical traits might also be different in Russian rapeseed lines. To assess it, we collected the glucosinolate content data for the same 90 genotyped accessions obtained during three years and performed an association mapping of this trait. We indeed found that the loci significantly associated with glucosinolate content variation in the Russian rapeseed collection differ from those previously reported for the non-Russian rapeseed lines.

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.

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.

Accelerated evolution of an Lhx2 enhancer shapes mammalian social hierarchies

Social hierarchies emerged during evolution, and social rank influences behavior and health of individuals. However, the evolutionary mechanisms of social hierarchy are still unknown in amniotes. Here we developed a new method and performed a genome-wide screening for identifying regions with accelerated evolution in the ancestral lineage of placental mammals, where mammalian social hierarchies might have initially evolved. Then functional analyses were conducted for the most accelerated region designated as placental-accelerated sequence 1 (PAS1, P = 3.15 × 10-18). Multiple pieces of evidence show that PAS1 is an enhancer of the transcription factor gene Lhx2 involved in brain development. PAS1s isolated from various amniotes showed different cis-regulatory activity in vitro, and affected the expression of Lhx2 differently in the nervous system of mouse embryos. PAS1 knock-out mice lack social stratification. PAS1 knock-in mouse models demonstrate that PAS1s determine the social dominance and subordinate of adult mice, and that social ranks could even be turned over by mutated PAS1. All homozygous mutant mice had normal huddled sleeping behavior, motor coordination and strength. Therefore, PAS1-Lhx2 modulates social hierarchies and is essential for establishing social stratification in amniotes, and positive Darwinian selection on PAS1 plays pivotal roles in the occurrence of mammalian social hierarchies.

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.

Differences in lipidome and metabolome organization of prefrontal cortex among human populations

Human populations, despite their overwhelming similarity, contain some distinct phenotypic, genetic, epigenetic, and gene expression features. In this study, we explore population differences at yet another level of molecular phenotype: the abundance of non-polar and polar low molecular weight compounds, lipids and metabolites in the prefrontal cortical region of the brain. We assessed the abundance of 1,670 lipids and 258 metabolites in 146 Han Chinese, 97 Western European, and 60 African American individuals of varying ages, covering most of the lifespan. The statistical analysis and logistic regression models both demonstrated extensive lipid and metabolic divergence of the Han Chinese individuals from the other two populations. This divergence was age-dependent, peaking in young adults, and involved metabolites and lipids clustering in specific metabolic pathways.

Hydrogen/Deuterium Exchange Aiding Compound Identification for LC-MS and MALDI Imaging Lipidomics

We present a novel approach for the increasing reliability of compound identification for LC-MS and MALDI imaging lipidomics. Our approach is based on the characterization of compounds not only by the elution time, accurate mass, and fragmentation spectra but also by the number of labile hydrogens that can be measured using the hydrogen/deuterium (H/D) exchange approach. The number of labile hydrogens (those from -OH and -NH groups) serves as an additional structural descriptor used when performing a database search. For LC-MS experiment, the H/D exchange was performed in the heating capillary of the modified electrospray ionization (ESI) source, while for MALDI imaging, the exchange was performed in the ion funnel at 10 Torr pressure. It was observed that such an approach allowed one to achieve a considerable degree of deuteration, enough to unambiguously distinguish between different classes of lipids. The proposed analytical approach may be successfully used for the identification not only of lipids but also of peptides and metabolites. A special software for the automatic filtration of molecules based on the number of functional groups was also developed.

Molecular footprint of Medawar's mutation accumulation process in mammalian aging

Medawar's mutation accumulation hypothesis explains aging by the declining force of natural selection with age: Slightly deleterious germline mutations expressed in old age can drift to fixation and thereby lead to aging‐related phenotypes. Although widely cited, empirical evidence for this hypothesis has remained limited. Here, we test one of its predictions that genes relatively highly expressed in old adults should be under weaker purifying selection than genes relatively highly expressed in young adults. Combining 66 transcriptome datasets (including 16 tissues from five mammalian species) with sequence conservation estimates across mammals, here we report that the overall conservation level of expressed genes is lower at old age compared to young adulthood. This age‐related decrease in transcriptome conservation (ADICT) is systematically observed in diverse mammalian tissues, including the brain, liver, lung, and artery, but not in others, most notably in the muscle and heart. Where observed, ADICT is driven partly by poorly conserved genes being up‐regulated during aging. In general, the more often a gene is found up‐regulated with age among tissues and species, the lower its evolutionary conservation. Poorly conserved and up‐regulated genes have overlapping functional properties that include responses to age‐associated tissue damage, such as apoptosis and inflammation. Meanwhile, these genes do not appear to be under positive selection. Hence, genes contributing to old age phenotypes are found to harbor an excess of slightly deleterious alleles, at least in certain tissues. This supports the notion that genetic drift shapes aging in multicellular organisms, consistent with Medawar's mutation accumulation hypothesis.

"Redox lipidomics technology: Looking for a needle in a haystack"

Aerobic life is based on numerous metabolic oxidation reactions as well as biosynthesis of oxygenated signaling compounds. Among the latter are the myriads of oxygenated lipids including a well-studied group of polyunsaturated fatty acids (PUFA) - octadecanoids, eicosanoids, and docosanoids. During the last two decades, remarkable progress in liquid-chromatography-mass spectrometry has led to significant progress in the characterization of oxygenated PUFA-containing phospholipids, thus designating the emergence of a new field of lipidomics, redox lipidomics. Although non-enzymatic free radical reactions of lipid peroxidation have been mostly associated with the aberrant metabolism typical of acute injury or chronic degenerative processes, newly accumulated evidence suggests that enzymatically catalyzed (phospho)lipid oxygenation reactions are essential mechanisms of many physiological pathways. In this review, we discuss a variety of contemporary protocols applicable for identification and quantitative characterization of different classes of peroxidized (phospho)lipids. We describe applications of different types of LCMS for analysis of peroxidized (phospho)lipids, particularly cardiolipins and phosphatidylethanolalmines, in two important types of programmed cell death - apoptosis and ferroptosis. We discuss the role of peroxidized phosphatidylserines in phagocytotic signaling. We exemplify the participation of peroxidized neutral lipids, particularly tri-acylglycerides, in immuno-suppressive signaling in cancer. We also consider new approaches to exploring the spatial distribution of phospholipids in the context of their oxidizability by MS imaging, including the latest achievements in high resolution imaging techniques. We present innovative approaches to the interpretation of LC-MS data, including audio-representation analysis. Overall, we emphasize the role of redox lipidomics as a communication language, unprecedented in diversity and richness, through the analysis of peroxidized (phospho)lipids.

Metabolome signature of autism in the human prefrontal cortex

Autism spectrum disorder (ASD) is a common neurodevelopmental disorder with yet incompletely uncovered molecular determinants. Alterations in the abundance of low molecular weight compounds (metabolites) in ASD could add to our understanding of the disease. Indeed, such alterations take place in the urine, plasma and cerebellum of ASD individuals. In this work, we investigated mass-spectrometric signal intensities of 1,366 metabolites in the prefrontal cortex grey matter of 32 ASD and 40 control individuals. 15% of these metabolites showed significantly different intensities in ASD and clustered in 16 metabolic pathways. Of them, ten pathways were altered in urine and blood of ASD individuals (Fisher test, p < 0.05), opening an opportunity for the design of new diagnostic instruments. Furthermore, metabolic measurements conducted in 40 chimpanzees and 40 macaques showed an excess of metabolite intensity differences unique to humans, supporting the hypothesized disruption of evolutionary novel cortical mechanisms in ASD.

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.

Identification of Immediate Early Genes in the Nervous System of Snail Helix lucorum

Immediate early genes (IEGs) are useful markers of neuronal activation and essential components of neuronal response. While studies of gastropods have provided many insights into the basic learning and memory mechanisms, the genome-wide assessment of IEGs has been mainly restricted to vertebrates. In this study, we identified IEGs in the terrestrial snail Helix lucorum In the absence of the genome, we conducted de novo transcriptome assembly using reads with short and intermediate lengths cumulatively covering more than 98 billion nucleotides. Based on this assembly, we identified 37 proteins corresponding to contigs differentially expressed (DE) in either the parietal ganglia (PaG) or two giant interneurons located within the PaG of the snail in response to the neuronal stimulation. These proteins included homologues of well-known mammalian IEGs, such as c-jun/jund, C/EBP, c-fos/fosl2, and Egr1, as well as homologues of genes not yet implicated in the neuronal response.

Recently Evolved Tumor Suppressor Transcript TP73-AS1 Functions as Sponge of Human-Specific miR-941

MicroRNA (miRNA) sponges are vital components of posttranscriptional gene regulation. Yet, only a limited number of miRNA sponges have been identified. Here, we show that the recently evolved noncoding tumor suppressor transcript, antisense RNA to TP73 gene (TP73-AS1), functions as a natural sponge of human-specific miRNA miR-941. We find unusually nine high-affinity miR-941 binding sites clustering within 1 kb region on TP73-AS1, which forms miR-941 sponge region. This sponge region displays increased sequence constraint only in humans, and its formation can be traced to the tandem expansion of a 71-nt-long sequence containing a single miR-941 binding site in old world monkeys. We further confirm TP73-AS1 functions as an efficient miR-941 sponge based on massive transcriptome data analyses, wound-healing assay, and Argonaute protein immunoprecipitation experiments conducted in cell lines. The expression of miR-941 and its sponge correlate inversely across multiple healthy and cancerous tissues, with miR-941 being highly expressed in tumors and preferentially repressing tumor suppressors. Thus, the TP73-AS1 and miR-941 duo represents an unusual case of the extremely rapid evolution of noncoding regulators controlling cell migration, proliferation, and tumorigenesis.

Ultra-performance liquid chromatography-mass spectrometry for precise fatty acid profiling of oilseed crops

Oilseed crops are one of the most important sources of vegetable oils for food and industry. Nutritional and technical properties of vegetable oil are primarily determined by its fatty acid (FA) composition. The content and composition of FAs in plants are commonly determined using gas chromatography-mass spectrometry (GS-MS) or gas chromatography-flame ionization detection (GC-FID) techniques. In the present work, we applied ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) technique to FA profiling of sunflower and rapeseed seeds and compared this method with the GC-FID technique. GC-FID detected 11 FAs in sunflower and 13 FAs in rapeseed, while UPLC-MS appeared to be more sensitive, detecting about 2.5 times higher numbers of FAs in both plants. In addition to even-chain FAs, UPLC-MS was able to detect odd-chain FAs. The longest FA detected using GC-FID was an FA with 24 carbon atoms, whereas UPLC-MS could reveal the presence of longer FAs with the tails of up to 28 carbon atoms. Based on our results, we may conclude that UPLC-MS has great potential to be used for the assessment of FA profiles of oil crops.

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.

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.

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.

Changes in Lipidome Composition during Brain Development in Humans, Chimpanzees, and Macaque Monkeys

Lipids are essential components of the brain. Here, we conducted a comprehensive mass spectrometry-based analysis of lipidome composition in the prefrontal cortex of 40 humans, 40 chimpanzees, and 40 rhesus monkeys over postnatal development and adulthood. Of the 11,772 quantified lipid peaks, 7,589 change significantly along the lifespan. More than 60% of these changes occur prior to adulthood, with less than a quarter associated with myelination progression. Evolutionarily, 36% of the age-dependent lipids exhibit concentration profiles distinct to one of the three species; 488 (18%) of them were unique to humans. In both humans and chimpanzees, the greatest extent of species-specific differences occurs in early development. Human-specific lipidome differences, however, persist over most of the lifespan and reach their peak from 20 to 35 years of age, when compared with chimpanzee-specific ones.

Gene expression reversal toward pre-adult levels in the aging human brain and age-related loss of cellular identity

It was previously reported that mRNA expression levels in the prefrontal cortex at old age start to resemble pre-adult levels. Such expression reversals could imply loss of cellular identity in the aging brain, and provide a link between aging-related molecular changes and functional decline. Here we analyzed 19 brain transcriptome age-series datasets, comprising 17 diverse brain regions, to investigate the ubiquity and functional properties of expression reversal in the human brain. Across all 19 datasets, 25 genes were consistently up-regulated during postnatal development and down-regulated in aging, displaying an "up-down" pattern that was significant as determined by random permutations. In addition, 113 biological processes, including neuronal and synaptic functions, were consistently associated with genes showing an up-down tendency among all datasets. Genes up-regulated during in vitro neuronal differentiation also displayed a tendency for up-down reversal, although at levels comparable to other genes. We argue that reversals may not represent aging-related neuronal loss. Instead, expression reversals may be associated with aging-related accumulation of stochastic effects that lead to loss of functional and structural identity in neurons.

Comprehensive transcriptome analysis of neocortical layers in humans, chimpanzees and macaques

While human cognitive abilities are clearly unique, underlying changes in brain organization and function remain unresolved. Here we characterized the transcriptome of the cortical layers and adjacent white matter in the prefrontal cortexes of humans, chimpanzees and rhesus macaques using unsupervised sectioning followed by RNA sequencing. More than 20% of detected genes were expressed predominantly in one layer, yielding 2,320 human layer markers. While the bulk of the layer markers were conserved among species, 376 switched their expression to another layer in humans. By contrast, only 133 of such changes were detected in the chimpanzee brain, suggesting acceleration of cortical reorganization on the human evolutionary lineage. Immunohistochemistry experiments further showed that human-specific expression changes were not limited to neurons but affected a broad spectrum of cortical cell types. Thus, despite apparent histological conservation, human neocortical organization has undergone substantial changes affecting more than 5% of its transcriptome.

Disruption of an Evolutionarily Novel Synaptic Expression Pattern in Autism

Cognitive defects in autism spectrum disorder (ASD) include socialization and communication: key behavioral capacities that separate humans from other species. Here, we analyze gene expression in the prefrontal cortex of 63 autism patients and control individuals, as well as 62 chimpanzees and macaques, from natal to adult age. We show that among all aberrant expression changes seen in ASD brains, a single aberrant expression pattern overrepresented in genes involved synaptic-related pathways is enriched in nucleotide variants linked to autism. Furthermore, only this pattern contains an excess of developmental expression features unique to humans, thus resulting in the disruption of human-specific developmental programs in autism. Several members of the early growth response (EGR) transcription factor family can be implicated in regulation of this aberrant developmental change. Our study draws a connection between the genetic risk architecture of autism and molecular features of cortical development unique to humans.

Mice carrying a human GLUD2 gene recapitulate aspects of human transcriptome and metabolome development

Whereas all mammals have one glutamate dehydrogenase gene (GLUD1), humans and apes carry an additional gene (GLUD2), which encodes an enzyme with distinct biochemical properties. We inserted a bacterial artificial chromosome containing the human GLUD2 gene into mice and analyzed the resulting changes in the transcriptome and metabolome during postnatal brain development. Effects were most pronounced early postnatally, and predominantly genes involved in neuronal development were affected. Remarkably, the effects in the transgenic mice partially parallel the transcriptome and metabolome differences seen between humans and macaques analyzed. Notably, the introduction of GLUD2 did not affect glutamate levels in mice, consistent with observations in the primates. Instead, the metabolic effects of GLUD2 center on the tricarboxylic acid cycle, suggesting that GLUD2 affects carbon flux during early brain development, possibly supporting lipid biosynthesis.

Fluoxetine treatment prevents the inflammatory response in a mouse model of posttraumatic stress disorder

Despite intense research efforts the molecular mechanisms affecting stress-vulnerable brain regions in posttraumatic stress disorder (PTSD) remain elusive. In the current study we have applied global transcriptomic profiling to a PTSD mouse model induced by foot shock fear conditioning. We compared the transcriptomes of prelimbic cortex, anterior cingulate cortex (ACC), basolateral amygdala, central nucleus of amygdala, nucleus accumbens (NAc) and CA1 of the dorsal hippocampus between shocked and non-shocked (control) mice, with and without fluoxetine treatment by RNA sequencing. Differentially expressed (DE) genes were identified and clustered for in silico pathway analysis. Findings in relevant brain regions were further validated with immunohistochemistry. DE genes belonging to 11 clusters were identified including increased inflammatory response in ACC in shocked mice. In line with this finding, we noted higher microglial activation in ACC of shocked mice. Chronic fluoxetine treatment initiated in the aftermath of the trauma prevented inflammatory gene expression alterations in ACC and ameliorated PTSD-like symptoms, implying an important role of the immune response in PTSD pathobiology. Our results provide novel insights into molecular mechanisms affected in PTSD and suggest therapeutic applications with anti-inflammatory agents.

Fluoxetine Treatment Rescues Energy Metabolism Pathway Alterations in a Posttraumatic Stress Disorder Mouse Model

Posttraumatic stress disorder (PTSD) is a prevalent psychiatric disorder. Several studies have attempted to characterize molecular alterations associated with PTSD, but most findings were limited to the investigation of specific cellular markers in the periphery or defined brain regions. In the current study, we aimed to unravel affected molecular pathways/mechanisms in the fear circuitry associated with PTSD. We interrogated a foot shock-induced PTSD mouse model by integrating proteomics and metabolomics profiling data. Alterations at the proteome level were analyzed using in vivo (15)N metabolic labeling combined with mass spectrometry in the prelimbic cortex (PrL), anterior cingulate cortex (ACC), basolateral amygdala, central nucleus of the amygdala and CA1 of the hippocampus between shocked and nonshocked (control) mice, with and without fluoxetine treatment. In silico pathway analyses revealed an upregulation of the citric acid cycle pathway in PrL, and downregulation in ACC and nucleus accumbens (NAc). Chronic fluoxetine treatment prevented decreased citric acid cycle activity in NAc and ACC and ameliorated conditioned fear response in shocked mice. Our results shed light on the role of energy metabolism in PTSD pathogenesis and suggest potential therapy through mitochondrial targeting.

Changes in snoRNA and snRNA Abundance in the Human, Chimpanzee, Macaque, and Mouse Brain

Small nuclear and nucleolar RNAs (snRNAs and snoRNAs) are known to be functionally and evolutionarily conserved elements of transcript processing machinery. Here, we investigated the expression evolution of snRNAs and snoRNAs by measuring their abundance in the frontal cortex of humans, chimpanzees, rhesus monkeys, and mice. Although snRNA expression is largely conserved, 44% of the 185 measured snoRNA and 40% of the 134 snoRNA families showed significant expression divergence among species. The snRNA and snoRNA expression divergence included drastic changes unique to humans: A 10-fold elevated expression of U1 snRNA and a 1,000-fold drop in expression of SNORA29. The decreased expression of SNORA29 might be due to two mutations that affect secondary structure stability. Using in situ hybridization, we further localized SNORA29 expression to nucleolar regions of neuronal cells. Our study presents the first observation of snoRNA abundance changes specific to the human lineage and suggests a possible mechanism underlying these changes.

Transcript and protein expression decoupling reveals RNA binding proteins and miRNAs as potential modulators of human aging

BACKGROUND

In studies of development and aging, the expression of many genes has been shown to undergo drastic changes at mRNA and protein levels. The connection between mRNA and protein expression level changes, as well as the role of posttranscriptional regulation in controlling expression level changes in postnatal development and aging, remains largely unexplored.

RESULTS

Here, we survey mRNA and protein expression changes in the prefrontal cortex of humans and rhesus macaques over developmental and aging intervals of both species' lifespans. We find substantial decoupling of mRNA and protein expression levels in aging, but not in development. Genes showing increased mRNA/protein disparity in primate brain aging form expression patterns conserved between humans and macaques and are enriched in specific functions involving mammalian target of rapamycin (mTOR) signaling, mitochondrial function and neurodegeneration. Mechanistically, aging-dependent mRNA/protein expression decoupling could be linked to a specific set of RNA binding proteins and, to a lesser extent, to specific microRNAs.

CONCLUSIONS

Increased decoupling of mRNA and protein expression profiles observed in human and macaque brain aging results in specific co-expression profiles composed of genes with shared functions and shared regulatory signals linked to specific posttranscriptional regulators. Genes targeted and predicted to be targeted by the aging-dependent posttranscriptional regulation are associated with biological processes known to play important roles in aging and lifespan extension. These results indicate the potential importance of posttranscriptional regulation in modulating aging-dependent changes in humans and other species.

Organization and evolution of brain lipidome revealed by large-scale analysis of human, chimpanzee, macaque, and mouse tissues

Lipids are prominent components of the nervous system. Here we performed a large-scale mass spectrometry-based analysis of the lipid composition of three brain regions as well as kidney and skeletal muscle of humans, chimpanzees, rhesus macaques, and mice. The human brain shows the most distinct lipid composition: 76% of 5,713 lipid compounds examined in our study are either enriched or depleted in the human brain. Concentration levels of lipids enriched in the brain evolve approximately four times faster among primates compared with lipids characteristic of non-neural tissues and show further acceleration of change in human neocortical regions but not in the cerebellum. Human-specific concentration changes are supported by human-specific expression changes for corresponding enzymes. These results provide the first insights into the role of lipids in human brain evolution.

De novo identification and quantification of single amino-acid variants in human brain

The detection of single amino-acid variants (SAVs) usually depends on single-nucleotide polymorphisms (SNPs) database. Here, we describe a novel method that discovers SAVs at proteome level independent of SNPs data. Using mass spectrometry-based de novo sequencing algorithm, peptide-candidates are identified and compared with theoretical protein database to generate SAVs under pairing strategy, which is followed by database re-searching to control false discovery rate. In human brain tissues, we can confidently identify known and novel protein variants with diverse origins. Combined with DNA/RNA sequencing, we verify SAVs derived from DNA mutations, RNA alternative splicing, and unknown post-transcriptional mechanisms. Furthermore, quantitative analysis in human brain tissues reveals several tissue-specific differential expressions of SAVs. This approach provides a novel access to high-throughput detection of protein variants, which may offer the potential for clinical biomarker discovery and mechanistic research.

Conserved expression of lincRNA during human and macaque prefrontal cortex development and maturation

The current annotation of the human genome includes more than 12,000 long intergenic noncoding RNAs (lincRNA). While a handful of lincRNA have been shown to play important regulatory roles, the functionality of most remains unclear. Here, we examined the expression conservation and putative functionality of lincRNA in human and macaque prefrontal cortex (PFC) development and maturation. We analyzed transcriptome sequence (RNA-seq) data from 38 human and 40 macaque individuals covering the entire postnatal development interval. Using the human data set, we detected the expression of 5835 lincRNA annotated in GENCODE and further identified 1888 novel lincRNA. Most of these lincRNA show low DNA sequence conservation, as well as low expression levels. Remarkably, developmental expression patterns of these lincRNA were as conserved between humans and macaques as those of protein-coding genes. Transfection of development-associated lincRNA into human SH-SY5Y cells affected gene expression, indicating their regulatory potential. In brain, expression of these putative target genes correlated with the expression of the corresponding lincRNA during human and macaque PFC development. These results support the potential functionality of lincRNA in primate PFC development.

Deep sequencing reveals a novel class of bidirectional promoters associated with neuronal genes

Background

Comprehensive annotation of transcripts expressed in a given tissue is a critical step towards the understanding of regulatory and functional pathways that shape the transcriptome.

Results

Here, we reconstructed a cumulative transcriptome of the human prefrontal cortex (PFC) based on approximately 300 million strand-specific RNA sequence (RNA-seq) reads collected at different stages of postnatal development. We find that more than 50% of reconstructed transcripts represent novel transcriptome elements, including 8,343 novel exons and exon extensions of annotated coding genes, 11,217 novel antisense transcripts and 29,541 novel intergenic transcripts or their fragments showing canonical features of long non-coding RNAs (lncRNAs). Our analysis further led to a surprising discovery of a novel class of bidirectional promoters (NBiPs) driving divergent transcription of mRNA and novel lncRNA pairs and displaying a distinct set of sequence and epigenetic features. In contrast to known bidirectional and unidirectional promoters, NBiPs are strongly associated with genes involved in neuronal functions and regulated by neuron-associated transcription factors.

Conclusions

Taken together, our results demonstrate that large portions of the human transcriptome remain uncharacterized. The distinct sequence and epigenetic features of NBiPs, as well as their specific association with neuronal genes, further suggest existence of regulatory pathways specific to the human brain.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-457) contains supplementary material, which is available to authorized users.