Search

Abstract Lipids are the most abundant but poorly explored components of the human brain. Here, we present a lipidome map of the human brain comprising 75 regions, including 52 neocortical ones. The lipidome composition varies greatly among the brain regions, affecting 93% of the 419 analyzed lipids. These differences reflect the brain’s structural characteristics, such as myelin content (345 lipids) and cell type composition (353 lipids), but also functional traits: functional connectivity (76 lipids) and information processing hierarchy (60 lipids). Combining lipid composition and mRNA expression data further enhances functional connectivity association. Biochemically, lipids linked with structural and functional brain features display distinct lipid class distribution, unsaturation extent, and prevalence of omega-3 and omega-6 fatty acid residues. We verified our conclusions by parallel analysis of three adult macaque brains, targeted analysis of 216 lipids, mass spectrometry imaging, and lipidome assessment of sorted murine neurons.

Abstract 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.

Abstract 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.

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.

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.

Abstract 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.

Abstract 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.

AbstractTocopherols 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

Abstract 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.

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 (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.

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.

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.

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.

Abstract 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.

Abstract Maximal lifespan of mammalian species, even if closely related, may differ more than 10-fold, however the nature of the mechanisms that determine this variability is unresolved. Here, we assess the relationship between maximal lifespan duration and concentrations of more than 20,000 lipid compounds, measured in 669 tissue samples from 6 tissues of 35 species representing three mammalian clades: primates, rodents and bats. We identify lipids associated with species’ longevity across the three clades, uncoupled from other parameters, such as basal metabolic rate, body size, or body temperature. These lipids clustered in specific lipid classes and pathways, and enzymes linked to them display signatures of greater stabilizing selection in long-living species, and cluster in functional groups related to signaling and protein-modification processes. These findings point towards the existence of defined molecular mechanisms underlying variation in maximal lifespan among mammals.

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.

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.

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.

Metabolite concentrations reflect the physiological states of tissues and cells. However, the role of metabolic changes in species evolution is currently unknown. Here, we present a study of metabolome evolution conducted in three brain regions and two non-neural tissues from humans, chimpanzees, macaque monkeys, and mice based on over 10,000 hydrophilic compounds. While chimpanzee, macaque, and mouse metabolomes diverge following the genetic distances among species, we detect remarkable acceleration of metabolome evolution in human prefrontal cortex and skeletal muscle affecting neural and energy metabolism pathways. These metabolic changes could not be attributed to environmental conditions and were confirmed against the expression of their corresponding enzymes. We further conducted muscle strength tests in humans, chimpanzees, and macaques. The results suggest that, while humans are characterized by superior cognition, their muscular performance might be markedly inferior to that of chimpanzees and macaque monkeys.

Caloric restriction (CR) and chemical agents, such as resveratrol and rapamycin that partially mimic the CR effect, can delay morbidity and mortality across a broad range of species. In humans, however, the effects of CR or other life-extending agents have not yet been investigated systematically. Human maximal lifespan is already substantially greater compared to that of closely related primate species. It is therefore possible that humans have acquired genetic mutations that mimic the CR effect. Here, we tested this notion by comparing transcriptome differences between humans and other primates, with the transcriptome changes observed in mice subjected to CR. We show that the human transcriptome state, relative to other primate transcriptomes, does not match that of the CR mice or mice treated with resveratrol, but resembles the transcriptome state of ad libitum fed mice. At the same time, the transcriptome changes induced by CR in mice are enriched among genes showing age-related changes in primates, concentrated in specific expression patterns, and can be linked with specific functional pathways, including insulin signalling, cancer, and the immune response. These findings indicate that the evolution of human longevity was likely independent of CR-induced lifespan extension mechanisms. Consequently, application of CR or CR-mimicking agents may yet offer a promising direction for the extension of healthy human lifespan.

While splicing differences between tissues, sexes and species are well documented, little is known about the extent and the nature of splicing changes that take place during human or mammalian development and aging. Here, using high‐throughput transcriptome sequencing, we have characterized splicing changes that take place during whole human lifespan in two brain regions: prefrontal cortex and cerebellum. Identified changes were confirmed using independent human and rhesus macaque RNA‐seq data sets, exon arrays and PCR, and were detected at the protein level using mass spectrometry. Splicing changes across lifespan were abundant in both of the brain regions studied, affecting more than a third of the genes expressed in the human brain. Approximately 15% of these changes differed between the two brain regions. Across lifespan, splicing changes followed discrete patterns that could be linked to neural functions, and associated with the expression profiles of the corresponding splicing factors. More than 60% of all splicing changes represented a single splicing pattern reflecting preferential inclusion of gene segments potentially targeting transcripts for nonsense‐mediated decay in infants and elderly.

BackgroundHumans have a widely different diet from other primate species, and are dependent on its high nutritional content. The molecular mechanisms responsible for adaptation to the human diet are currently unknown. Here, we addressed this question by investigating whether the gene expression response observed in mice fed human and chimpanzee diets involves the same regulatory mechanisms as expression differences between humans and chimpanzees.ResultsUsing mouse and primate transcriptomic data, we identified the transcription factor EGR1 (early growth response 1) as a putative regulator of diet-related differential gene expression between human and chimpanzee livers. Specifically, we predict that EGR1 regulates the response to the high caloric content of human diets. However, we also show that close to 90% of the dietary response to the primate diet found in mice, is not observed in primates. This might be explained by changes in tissue-specific gene expression between taxa.ConclusionOur results suggest that the gene expression response to the nutritionally rich human diet is partially mediated by the transcription factor EGR1. While this EGR1-driven response is conserved between mice and primates, the bulk of the mouse response to human and chimpanzee dietary differences is not observed in primates. This result highlights the rapid evolution of diet-related expression regulation and underscores potential limitations of mouse models in dietary studies.

Recent studies demonstrate that the organization of the chromatin within the nuclear space might play a crucial role in the regulation of gene expression. The ongoing progress in determination of the 3D structure of the nuclear chromatin allows one to study correlations between spatial proximity of genome domains and their epigenetic state. We combined the data on three-dimensional architecture of the whole human genome with results of high-throughput studies of the chromatin functional state and observed that fragments of different chromosomes that are spatially close tend to have similar patterns of histone modifications, methylation state, DNAse sensitivity, expression level, and chromatin states in general. Moreover, clustering of genome regions by spatial proximity produced compact clusters characterized by the high level of histone modifications and DNAse sensitivity and low methylation level, and loose clusters with the opposite characteristics. We also associated the spatial proximity data with previously detected chimeric transcripts and the results of RNA-seq experiments and observed that the frequency of formation of chimeric transcripts from fragments of two different chromosomes is higher among spatially proximal genome domains. A fair fraction of these chimeric transcripts seems to arise post-transcriptionally via trans-splicing.

While multiple studies have reported the accelerated evolution of brain gene expression in the human lineage, the mechanisms underlying such changes are unknown. Here, we address this issue from a developmental perspective, by analyzing mRNA and microRNA (miRNA) expression in two brain regions within macaques, chimpanzees, and humans throughout their lifespan. We find that constitutive gene expression divergence (species differences independent of age) is comparable between humans and chimpanzees. However, humans display a 3-5 times faster evolutionary rate in divergence of developmental patterns, compared to chimpanzees. Such accelerated evolution of human brain developmental patterns (i) cannot be explained by life-history changes among species, (ii) is twice as pronounced in the prefrontal cortex than the cerebellum, (iii) preferentially affects neuron-related genes, and (iv) unlike constitutive divergence does not depend on cis-regulatory changes, but might be driven by human-specific changes in expression of trans-acting regulators. We show that developmental profiles of miRNAs, as well as their target genes, show the fastest rates of human-specific evolutionary change, and using a combination of computational and experimental methods, we identify miR-92a, miR-454, and miR-320b as possible regulators of human-specific neural development. Our results suggest that different mechanisms underlie adaptive and neutral transcriptome divergence, and that changes in the expression of a few key regulators may have been a major driving force behind rapid evolution of the human brain.

Among other factors, changes in gene expression on the human evolutionary lineage have been suggested to play an important role in the establishment of human-specific phenotypes. However, the molecular mechanisms underlying these expression changes are largely unknown. Here, we have explored the role of microRNA (miRNA) in the regulation of gene expression divergence among adult humans, chimpanzees, and rhesus macaques, in two brain regions: prefrontal cortex and cerebellum. Using a combination of high-throughput sequencing, miRNA microarrays, and Q-PCR, we have shown that up to 11% of the 325 expressed miRNA diverged significantly between humans and chimpanzees and up to 31% between humans and macaques. Measuring mRNA and protein expression in human and chimpanzee brains, we found a significant inverse relationship between the miRNA and the target genes expression divergence, explaining 2%-4% of mRNA and 4%-6% of protein expression differences. Notably, miRNA showing human-specific expression localize in neurons and target genes that are involved in neural functions. Enrichment in neural functions, as well as miRNA-driven regulation on the human evolutionary lineage, was further confirmed by experimental validation of predicted miRNA targets in two neuroblastoma cell lines. Finally, we identified a signature of positive selection in the upstream region of one of the five miRNA with human-specific expression, miR-34c-5p. This suggests that miR-34c-5p expression change took place after the split of the human and the Neanderthal lineages and had adaptive significance. Taken together these results indicate that changes in miRNA expression might have contributed to evolution of human cognitive functions.

Transcription is the first step connecting genetic information with an organism's phenotype. While expression of annotated genes in the human brain has been characterized extensively, our knowledge about the scope and the conservation of transcripts located outside of the known genes' boundaries is limited. Here, we use high-throughput transcriptome sequencing (RNA-Seq) to characterize the total non-ribosomal transcriptome of human, chimpanzee, and rhesus macaque brain. In all species, only 20-28% of non-ribosomal transcripts correspond to annotated exons and 20-23% to introns. By contrast, transcripts originating within intronic and intergenic repetitive sequences constitute 40-48% of the total brain transcriptome. Notably, some repeat families show elevated transcription. In non-repetitive intergenic regions, we identify and characterize 1,093 distinct regions highly expressed in the human brain. These regions are conserved at the RNA expression level across primates studied and at the DNA sequence level across mammals. A large proportion of these transcripts (20%) represents 3'UTR extensions of known genes and may play roles in alternative microRNA-directed regulation. Finally, we show that while transcriptome divergence between species increases with evolutionary time, intergenic transcripts show more expression differences among species and exons show less. Our results show that many yet uncharacterized evolutionary conserved transcripts exist in the human brain. Some of these transcripts may play roles in transcriptional regulation and contribute to evolution of human-specific phenotypic traits.

With few exceptions, current methods for short read mapping make use of simple seed heuristics to speed up the search. Most of the underlying matching models neglect the necessity to allow not only mismatches, but also insertions and deletions. Current evaluations indicate, however, that very different error models apply to the novel high-throughput sequencing methods. While the most frequent error-type in Illumina reads are mismatches, reads produced by 454's GS FLX predominantly contain insertions and deletions (indels). Even though 454 sequencers are able to produce longer reads, the method is frequently applied to small RNA (miRNA and siRNA) sequencing. Fast and accurate matching in particular of short reads with diverse errors is therefore a pressing practical problem. We introduce a matching model for short reads that can, besides mismatches, also cope with indels. It addresses different error models. For example, it can handle the problem of leading and trailing contaminations caused by primers and poly-A tails in transcriptomics or the length-dependent increase of error rates. In these contexts, it thus simplifies the tedious and error-prone trimming step. For efficient searches, our method utilizes index structures in the form of enhanced suffix arrays. In a comparison with current methods for short read mapping, the presented approach shows significantly increased performance not only for 454 reads, but also for Illumina reads. Our approach is implemented in the software segemehl available at http://www.bioinf.uni-leipzig.de/Software/segemehl/.

Although the human diet is markedly different from the diets of closely related primate species, the influence of diet on phenotypic and genetic differences between humans and other primates is unknown. In this study, we analyzed gene expression in laboratory mice fed diets typical of humans and of chimpanzees. The effects of human diets were found to be significantly different from that of a chimpanzee diet in the mouse liver, but not in the brain. Importantly, 10% of the genes that differ in their expression between humans and chimpanzee livers differed also between the livers of mice fed the human and chimpanzee diets. Furthermore, both the promoter sequences and the amino acid sequences of these diet-related genes carry more differences between humans and chimpanzees than random genes. Our results suggest that the mouse can be used to study at least some aspects of human-specific traits.

Even though mRNA expression levels are commonly used as a proxy for estimating functional differences that occur at the protein level, the relation between mRNA and protein expression is not well established. Further, no study to date has tested whether the evolutionary differences in mRNA expression observed between species reflect those observed in protein expression. Since a large proportion of mRNA expression differences observed between mammalian species appears to have no functional consequences for the phenotype, it is conceivable that many or most mRNA expression differences are not reflected at the protein level. If this is true, then differences in protein expression may largely reflect functional adaptations observed in species phenotypes. In this paper, we present the first direct comparison of mRNA and protein expression differences seen between humans and chimpanzees. We reproducibly find a significant positive correlation between mRNA expression and protein expression differences. This correlation is comparable in magnitude to that found between mRNA and protein expression changes at different developmental stages or in different physiological conditions within one species. Noticeably, this correlation is mainly due to genes with large expression differences between species. Our study opens the door to a new level of understanding of regulatory evolution and poses many new questions that remain to be answered.

Although a large proportion of human transcription occurs outside the boundaries of known genes, the functional significance of this transcription remains unknown. We have compared the expression patterns of known genes as well as intergenic transcripts within the ENCODE regions between humans and chimpanzees in brain, heart, testis, and lymphoblastoid cell lines. We find that intergenic transcripts show patterns of tissue-specific conservation of their expression, which are comparable to exonic transcripts of known genes. This suggests that intergenic transcripts are subject to functional constraints that restrict their rate of evolutionary change as well as putative positive selection to an extent comparable to that of classical protein-coding genes. In brain and testis, we find that part of this intergenic transcription is caused by widespread use of alternative promoters. Further, we find that about half of the expression differences between humans and chimpanzees are due to intergenic transcripts.

It is well established that gene expression levels in many organisms change during the aging process, and the advent of DNA microarrays has allowed genome-wide patterns of transcriptional changes associated with aging to be studied in both model organisms and various human tissues. Understanding the effects of aging on gene expression in the human brain is of particular interest, because of its relation to both normal and pathological neurodegeneration. Here we show that human cerebral cortex, human cerebellum, and chimpanzee cortex each undergo different patterns of age-related gene expression alterations. In humans, many more genes undergo consistent expression changes in the cortex than in the cerebellum; in chimpanzees, many genes change expression with age in cortex, but the pattern of changes in expression bears almost no resemblance to that of human cortex. These results demonstrate the diversity of aging patterns present within the human brain, as well as how rapidly genome-wide patterns of aging can evolve between species; they may also have implications for the oxidative free radical theory of aging, and help to improve our understanding of human neurodegenerative diseases.

Retroviral infections of the germline have the potential to episodically alter gene function and genome structure during the course of evolution. Horizontal transmissions between species have been proposed, but little evidence exists for such events in the human/great ape lineage of evolution. Based on analysis of finished BAC chimpanzee genome sequence, we characterize a retroviral element (Pan troglodytes endogenous retrovirus 1 [PTERV1]) that has become integrated in the germline of African great ape and Old World monkey species but is absent from humans and Asian ape genomes. We unambiguously map 287 retroviral integration sites and determine that approximately 95.8% of the insertions occur at non-orthologous regions between closely related species. Phylogenetic analysis of the endogenous retrovirus reveals that the gorilla and chimpanzee elements share a monophyletic origin with a subset of the Old World monkey retroviral elements, but that the average sequence divergence exceeds neutral expectation for a strictly nuclear inherited DNA molecule. Within the chimpanzee, there is a significant integration bias against genes, with only 14 of these insertions mapping within intronic regions. Six out of ten of these genes, for which there are expression data, show significant differences in transcript expression between human and chimpanzee. Our data are consistent with a retroviral infection that bombarded the genomes of chimpanzees and gorillas independently and concurrently, 3-4 million years ago. We speculate on the potential impact of such recent events on the evolution of humans and great apes.

Microarray technologies allow the identification of large numbers of expression differences within and between species. Although environmental and physiological stimuli are clearly responsible for changes in the expression levels of many genes, it is not known whether the majority of changes of gene expression fixed during evolution between species and between various tissues within a species are caused by Darwinian selection or by stochastic processes. We find the following: (1) expression differences between species accumulate approximately linearly with time; (2) gene expression variation among individuals within a species correlates positively with expression divergence between species; (3) rates of expression divergence between species do not differ significantly between intact genes and expressed pseudogenes; (4) expression differences between brain regions within a species have accumulated approximately linearly with time since these regions emerged during evolution. These results suggest that the majority of expression differences observed between species are selectively neutral or nearly neutral and likely to be of little or no functional significance. Therefore, the identification of gene expression differences between species fixed by selection should be based on null hypotheses assuming functional neutrality. Furthermore, it may be possible to apply a molecular clock based on expression differences to infer the evolutionary history of tissues.