Conservation of core gene expression in vertebrate tissues

Gene expression supports a single origin of horns and antlers in hoofed mammals

Horns, antlers, and other bony cranial appendages of even-toed hoofed mammals (ruminant artiodactyls) challenge traditional morphological homology assessments. Cranial appendages all share a permanent bone portion with family-specific integument coverings, but homology determination depends on whether the integument covering is an essential component or a secondary elaboration of each structure. To enhance morphological homology assessments, we tested whether juvenile cattle horn bud transcriptomes share homologous gene expression patterns with deer antlers relative to pig outgroup tissues, treating the integument covering as a secondary elaboration. We uncovered differentially expressed genes that support horn and antler homology, potentially distinguish them from non-cranial-appendage bone and other tissues, and highlight the importance of phylogenetic outgroups in homology assessments. Furthermore, we found differentially expressed genes that could support a shared cranial neural crest origin for horns and antlers and expression patterns that refine our understanding of the timing of horn and antler differentiation.

Meta-analysis of RNA-Seq datasets highlights novel genes/pathways involved in fat deposition in fat-tail of sheep

Introduction

Fat-tail in sheep is considered as an important energy reservoir to provide energy as a survival buffer during harsh challenges. However, fat-tail is losing its importance in modern sheep industry systems and thin-tailed breeds are more desirable. Using comparative transcriptome analysis to compare fat-tail tissue between fat- and thin-tailed sheep breeds provides a valuable approach to study the complex genetic factors associated with fat-tail development. However, transcriptomic studies often suffer from issues with reproducibility, which can be improved by integrating multiple studies based on a meta-analysis.

Methods

Hence, for the first time, an RNA-Seq meta-analysis on sheep fat-tail transcriptomes was performed using six publicly available datasets.

Results and discussion

A total of 500 genes (221 up-regulated, 279 down-regulated) were identified as differentially expressed genes (DEGs). A jackknife sensitivity analysis confirmed the robustness of the DEGs. Moreover, QTL and functional enrichment analysis reinforced the importance of the DEGs in the underlying molecular mechanisms of fat deposition. Protein-protein interactions (PPIs) network analysis revealed the functional interactions among the DEGs and the subsequent sub-network analysis led to identify six functional sub-networks. According to the results of the network analysis, down-regulated DEGs in green and pink sub-networks (like collagen subunits IV, V, and VI, integrins 1 and 2, SCD, SCD5, ELOVL6, ACLY, SLC27A2, and LPIN1) may impair lipolysis or fatty acid oxidation and cause fat accumulation in tail. On the other hand, up-regulated DEGs, especially those are presented in green and pink sub-networks (like IL6, RBP4, LEPR, PAI-1, EPHX1, HSD11B1, and FMO2), might contribute to a network controlling fat accumulation in the tail of sheep breed through mediating adipogenesis and fatty acid biosynthesis. Our results highlighted a set of known and novel genes/pathways associated with fat-tail development, which could improve the understanding of molecular mechanisms behind fat deposition in sheep fat-tail.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation

Longevity is a defining, heritable trait that varies dramatically between species. To resolve the genetic regulation of this trait, we have mined genomic variation in rockfishes, which range in longevity from 11 to over 205 years. Multiple shifts in rockfish longevity have occurred independently and in a short evolutionary time frame, thus empowering convergence analyses. Our analyses reveal a common network of genes under convergent evolution, encompassing established aging regulators such as insulin signaling, yet also identify flavonoid (aryl-hydrocarbon) metabolism as a pathway modulating longevity. The selective pressures on these pathways indicate the ancestral state of rockfishes was long lived and that the changes in short-lived lineages are adaptive. These pathways were also used to explore genome-wide association studies of human longevity, identifying the aryl-hydrocarbon metabolism pathway to be significantly associated with human survival to the 99th percentile. This evolutionary intersection defines and cross-validates a previously unappreciated genetic architecture that associates with the evolution of longevity across vertebrates.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Non-targeted proteomics reveals altered immune response in geographically distinct populations of green sea turtles (Chelonia mydas)

All seven species of sea turtle are facing increasing pressures from human activities that are impacting their health. Changes in circulating blood proteins of an individual, or all members of a population, can provide an early indicator of adverse health outcomes. Non-targeted measurement of all detectable proteins in a blood sample can indicate physiological changes. In the context of wildlife toxicology, this technique can provide a powerful tool for discovering biomarkers of chemical exposure and effect. This study presents a non-targeted examination of the protein abundance in sea turtle plasma obtained from three geographically distinct foraging populations of green turtles (Chelonia mydas) on the Queensland coast. Relative changes in protein expression between sites were compared, and potential markers of contaminant exposure were investigated. Blood plasma protein profiles were distinct between populations, with 85 out of the 116 identified proteins differentially expressed (p < 0.001). The most strongly dysregulated proteins were predominantly acute phase proteins, suggestive of differing immune status between the populations. The highest upregulation of known markers of immunotoxicity, such as pentraxin fusion and complement factor h, was observed in the Moreton Bay turtles. Forty-five different organohalogens were also measured in green turtle plasma samples as exposure to some organohalogens (e.g., polychlorinated biphenyls) has previously been identified as a cause for immune dysregulation in marine animals. The few detected organohalogens were at very low (pg/mL) concentrations in turtles from all sites, and are unlikely to be the cause of the proteome differences observed. However, the changes in protein expression may be indicative of exposure to other chemicals or environmental stressors. The results of this study provide important information about differences in protein expression between different populations of turtles, and guide future toxicological and health studies on east-Australian green sea turtles.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Extensive Adaptive Variation in Gene Expression within and between Closely Related Horseshoe Bats (Chiroptera, Rhinolophus) Revealed by Three Organs

In the process of species differentiation and adaption, the relative influence of natural selection on gene expression variation often remains unclear (especially its impact on phenotypic divergence). In this study, we used differentially expressed genes from brain, cochlea, and liver samples collected from two species of bats to determine the gene expression variation forced by natural selection when comparing at the interspecific (Rhinolophus siamensis and R. episcopus episcopus) and the intraspecific (R. e. episcopus and R. episcopus spp.) levels. In both cases, gene expression variation was extensively adaptive (>66.0%) and mainly governed by directional selection, followed by stabilizing selection, and finally balancing selection. The expression variation related to acoustic signals (resting frequency, RF) and body size (forearm length, FA) was also widely governed by natural selection (>69.1%). Different functional patterns of RF- or FA-related adaptive expression variation were found between the two comparisons, which manifested as abundant immune-related regulations between subspecies (indicating a relationship between immune response and phenotypic adaption). Our study verifies the extensive adaptive expression variation between both species and subspecies and provides insight into the effects of natural selection on species differentiation and adaptation as well as phenotypic divergence at the expression level.

Cardiac ion channel expression in the equine model - In-silico prediction utilising RNA sequencing data from mixed tissue samples

Understanding cardiomyocyte ion channel expression is crucial to understanding normal cardiac electrophysiology and underlying mechanisms of cardiac pathologies particularly arrhythmias. Hitherto, equine cardiac ion channel expression has rarely been investigated. Therefore, we aim to predict equine cardiac ion channel gene expression. Raw RNAseq data from normal horses from 9 datasets was retrieved from ArrayExpress and European Nucleotide Archive and reanalysed. The normalised (FPKM) read counts for a gene in a mix of tissue were hypothesised to be the average of the expected expression in each tissue weighted by the proportion of the tissue in the mix. The cardiac-specific expression was predicted by estimating the mean expression in each other tissues. To evaluate the performance of the model, predicted gene expression values were compared to the human cardiac gene expression. Cardiac-specific expression could be predicted for 91 ion channels including most expressed Na+ channels, K+ channels and Ca2+ -handling proteins. These revealed interesting differences from what would be expected based on human studies. These differences included predominance of NaV 1.4 rather than NaV 1.5 channel, and RYR1, SERCA1 and CASQ1 rather than RYR2, SERCA2, CASQ2 Ca2+ -handling proteins. Differences in channel expression not only implicate potentially different regulatory mechanisms but also pathological mechanisms of arrhythmogenesis.

regCNN: identifying Drosophila genome-wide cis-regulatory modules via integrating the local patterns in epigenetic marks and transcription factor binding motifs

Transcription regulation in metazoa is controlled by the binding events of transcription factors (TFs) or regulatory proteins on specific modular DNA regulatory sequences called cis-regulatory modules (CRMs). Understanding the distributions of CRMs on a genomic scale is essential for constructing the metazoan transcriptional regulatory networks that help diagnose genetic disorders. While traditional reporter-assay CRM identification approaches can provide an in-depth understanding of functions of some CRM, these methods are usually cost-inefficient and low-throughput. It is generally believed that by integrating diverse genomic data, reliable CRM predictions can be made. Hence, researchers often first resort to computational algorithms for genome-wide CRM screening before specific experiments. However, current existing in silico methods for searching potential CRMs were restricted by low sensitivity, poor prediction accuracy, or high computation time from TFBS composition combinatorial complexity. To overcome these obstacles, we designed a novel CRM identification pipeline called regCNN by considering the base-by-base local patterns in TF binding motifs and epigenetic profiles. On the test set, regCNN shows an accuracy/auROC of 84.5%/92.5% in CRM identification. And by further considering local patterns in epigenetic profiles and TF binding motifs, it can accomplish 4.7% (92.5%–87.8%) improvement in the auROC value over the average value-based pure multi-layer perceptron model. We also demonstrated that regCNN outperforms all currently available tools by at least 11.3% in auROC values. Finally, regCNN is verified to be robust against its resizing window hyperparameter in dealing with the variable lengths of CRMs. The model of regCNN can be downloaded athttp://cobisHSS0.im.nuk.edu.tw/regCNN/.

Comparative Analyses of Gene Co-expression Networks: Implementations and Applications in the Study of Evolution

Similarities and differences in the associations of biological entities among species can provide us with a better understanding of evolutionary relationships. Often the evolution of new phenotypes results from changes to interactions in pre-existing biological networks and comparing networks across species can identify evidence of conservation or adaptation. Gene co-expression networks (GCNs), constructed from high-throughput gene expression data, can be used to understand evolution and the rise of new phenotypes. The increasing abundance of gene expression data makes GCNs a valuable tool for the study of evolution in non-model organisms. In this paper, we cover motivations for why comparing these networks across species can be valuable for the study of evolution. We also review techniques for comparing GCNs in the context of evolution, including local and global methods of graph alignment. While some protein-protein interaction (PPI) bioinformatic methods can be used to compare co-expression networks, they often disregard highly relevant properties, including the existence of continuous and negative values for edge weights. Also, the lack of comparative datasets in non-model organisms has hindered the study of evolution using PPI networks. We also discuss limitations and challenges associated with cross-species comparison using GCNs, and provide suggestions for utilizing co-expression network alignments as an indispensable tool for evolutionary studies going forward.

Transcription Factor Reprogramming in the Inner Ear: Turning on Cell Fate Switches to Regenerate Sensory Hair Cells

Non-mammalian vertebrates can restore their auditory and vestibular hair cells naturally by triggering the regeneration of adjacent supporting cells. The transcription factor ATOH1 is a key regulator of hair cell development and regeneration in the inner ear. Following the death of hair cells, supporting cells upregulate ATOH1 and give rise to new hair cells. However, in the mature mammalian cochlea, such natural regeneration of hair cells is largely absent. Transcription factor reprogramming has been used in many tissues to convert one cell type into another, with the long-term hope of achieving tissue regeneration. Reprogramming transcription factors work by altering the transcriptomic and epigenetic landscapes in a target cell, resulting in a fate change to the desired cell type. Several studies have shown that ATOH1 is capable of reprogramming cochlear non-sensory tissue into cells resembling hair cells in young animals. However, the reprogramming ability of ATOH1 is lost with age, implying that the potency of individual hair cell-specific transcription factors may be reduced or lost over time by mechanisms that are still not clear. To circumvent this, combinations of key hair cell transcription factors have been used to promote hair cell regeneration in older animals. In this review, we summarize recent findings that have identified and studied these reprogramming factor combinations for hair cell regeneration. Finally, we discuss the important questions that emerge from these findings, particularly the feasibility of therapeutic strategies using reprogramming factors to restore human hearing in the future.

Juxtapose: a gene-embedding approach for comparing co-expression networks

BACKGROUND

Gene co-expression networks (GCNs) are not easily comparable due to their complex structure. In this paper, we propose a tool, Juxtapose, together with similarity measures that can be utilized for comparative transcriptomics between a set of organisms. While we focus on its application to comparing co-expression networks across species in evolutionary studies, Juxtapose is also generalizable to co-expression network comparisons across tissues or conditions within the same species.

METHODS

A word embedding strategy commonly used in natural language processing was utilized in order to generate gene embeddings based on walks made throughout the GCNs. Juxtapose was evaluated based on its ability to embed the nodes of synthetic structures in the networks consistently while also generating biologically informative results. Evaluation of the techniques proposed in this research utilized RNA-seq datasets from GTEx, a multi-species experiment of prefrontal cortex samples from the Gene Expression Omnibus, as well as synthesized datasets. Biological evaluation was performed using gene set enrichment analysis and known gene relationships in literature.

RESULTS

We show that Juxtapose is capable of globally aligning synthesized networks as well as identifying areas that are conserved in real gene co-expression networks without reliance on external biological information. Furthermore, output from a matching algorithm that uses cosine distance between GCN embeddings is shown to be an informative measure of similarity that reflects the amount of topological similarity between networks.

CONCLUSIONS

Juxtapose can be used to align GCNs without relying on known biological similarities and enables post-hoc analyses using biological parameters, such as orthology of genes, or conserved or variable pathways.

AVAILABILITY

A development version of the software used in this paper is available at https://github.com/klovens/juxtapose.

Evolution of dosage compensation does not depend on genomic background

Organisms have evolved various mechanisms to cope with the differences in the gene copy numbers between sexes caused by degeneration of Y and W sex chromosomes. Complete dosage compensation or at least expression balance between sexes has been reported predominantly in XX/XY systems, but rarely in ZZ/ZW systems. However, this often-reported pattern is based on comparisons of lineages where sex chromosomes evolved from nonhomologous genomic regions, potentially differing in sensitivity to differences in gene copy numbers. Here we document that two reptilian lineages (XX/XY iguanas and ZZ/ZW softshell turtles), which independently co-opted the same ancestral genomic region for the function of sex chromosomes, evolved different gene dose regulatory mechanisms. The independent co-option of the same genomic region for the role of sex chromosomes as in the iguanas and the softshell turtles offers great opportunity for testing evolutionary scenarios on sex chromosome evolution under the explicit control of the genomic background and gene identity. We show that the parallel loss of functional genes from the Y chromosome of the green anole and the W chromosome of the Florida softshell turtle led to different dosage compensation mechanisms. Our approach controlling for genetic background thus does not support that the variability in the regulation of gene dose differences is a consequence of ancestral autosomal gene content.

Gene expression dynamics during rapid organismal diversification in African cichlid fishes

Changes in gene expression play a fundamental role in phenotypic evolution. Transcriptome evolutionary dynamics have so far mainly been compared among distantly related species and remain largely unexplored during rapid organismal diversification, in which gene regulatory changes have been suggested as particularly effective drivers of phenotypic divergence. Here we studied gene expression evolution in a model system of adaptive radiation, the cichlid fishes of African Lake Tanganyika. By comparing gene expression profiles of 6 different organs in 74 cichlid species representing all subclades of this radiation, we demonstrate that the rate of gene expression evolution varies among organs, transcriptome parts and the subclades of the radiation, indicating different strengths of selection. We found that the noncoding part of the transcriptome evolved more rapidly than the coding part, and that the gonadal transcriptomes evolved more rapidly than the somatic ones, with the exception of liver. We further show that the rate of gene expression change was not constant over the course of the radiation but accelerated at its later phase. Finally, we show that-at the per-gene level-the evolution of expression patterns is dominated by stabilizing selection.

A Synthesis Concerning Conservation and Divergence of Cell Types across Epithelia

Advances in single-cell RNA-seq (scRNA-seq) and computational analysis have enabled the systematic interrogation of the cellular composition of tissues. Combined with tools from developmental biology, cell biology, and genetics, these approaches are revealing fundamental aspects of tissue geometry and physiology, including the distribution, origins, and inferred functions of specialized cell types, and the dynamics of cellular turnover and differentiation. By comparing different tissues, such studies can delineate shared and specialized features of cell types and their lineage. Here, we compare two developmentally related murine epithelia, the airway and the small intestinal epithelia, which are both derived from the embryonic endodermal gut tube. We examine how airway and intestine generate and functionalize common archetypal cell types to fulfill similar shared physiologic functionalities. We point to cases in which similar cell types are repurposed to accommodate each tissue's unique physiologic role, and highlight tissue-specific cells whose specializations contribute to the distinct functional roles of each organ. We discuss how archetypal and unique cell types are incorporated within a cellular lineage, and how the regulation of the proportions of these cell types enables tissue-level organization to meet functional demands and maintain homeostasis.

Gene expression variability in human and chimpanzee populations share common determinants

Inter-individual variation in gene expression has been shown to be heritable and is often associated with differences in disease susceptibility between individuals. Many studies focused on mapping associations between genetic and gene regulatory variation, yet much less attention has been paid to the evolutionary processes that shape the observed differences in gene regulation between individuals in humans or any other primate. To begin addressing this gap, we performed a comparative analysis of gene expression variability and expression quantitative trait loci (eQTLs) in humans and chimpanzees, using gene expression data from primary heart samples. We found that expression variability in both species is often determined by non-genetic sources, such as cell-type heterogeneity. However, we also provide evidence that inter-individual variation in gene regulation can be genetically controlled, and that the degree of such variability is generally conserved in humans and chimpanzees. In particular, we found a significant overlap of orthologous genes associated with eQTLs in both species. We conclude that gene expression variability in humans and chimpanzees often evolves under similar evolutionary pressures.

Comparative Transcriptome Analysis of Pine Trees Treated with Resistance-Inducing Substances against the Nematode Bursaphelenchus xylophilus

The pinewood nematode (PWN) Bursaphelenchus xylophilus causes pine wilt disease, which results in substantial economic and environmental losses across pine forests worldwide. Although systemic acquired resistance (SAR) is effective in controlling PWN, the detailed mechanisms underlying the resistance to PWN are unclear. Here, we treated pine samples with two SAR elicitors, acibenzolar-S-methyl (ASM) and methyl salicylic acid (MeSA) and constructed an in vivo transcriptome of PWN-infected pines under SAR conditions. A total of 252 million clean reads were obtained and mapped onto the reference genome. Compared with untreated pines, 1091 and 1139 genes were differentially upregulated following the ASM and MeSA treatments, respectively. Among these, 650 genes showed co-expression patterns in response to both SAR elicitors. Analysis of these patterns indicated a functional linkage among photorespiration, peroxisome, and glycine metabolism, which may play a protective role against PWN infection-induced oxidative stress. Further, the biosynthesis of flavonoids, known to directly control parasitic nematodes, was commonly upregulated under SAR conditions. The ASM- and MeSA-specific expression patterns revealed functional branches for myricetin and quercetin production in flavonol biosynthesis. This study will enhance the understanding of the dynamic interactions between pine hosts and PWN under SAR conditions.

Comparative transcriptomics of primary cells in vertebrates

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Cross-Species Single-Cell Analysis Reveals Divergence of the Primate Microglia Program

Microglia, the brain-resident immune cells, are critically involved in many physiological and pathological brain processes, including neurodegeneration. Here we characterize microglia morphology and transcriptional programs across ten species spanning more than 450 million years of evolution. We find that microglia express a conserved core gene program of orthologous genes from rodents to humans, including ligands and receptors associated with interactions between glia and neurons. In most species, microglia show a single dominant transcriptional state, whereas human microglia display significant heterogeneity. In addition, we observed notable differences in several gene modules of rodents compared with primate microglia, including complement, phagocytic, and susceptibility genes to neurodegeneration, such as Alzheimer's and Parkinson's disease. Our study provides an essential resource of conserved and divergent microglia pathways across evolution, with important implications for future development of microglia-based therapies in humans.

Animal Use in Neurobiological Research

The fact that neurobiological research is reliant upon laboratory-reared rodents is well known. The following paper discusses this topic broadly, but also aims to highlight other species used in the study of the nervous system and the evolution of animal species usage from the end of World War II through recent investigations. Attention is drawn to the dramatic reduction in the diversity of species used in neuroscience, with a significant shift toward two species, the mouse (Mus musculus) and rat (Rattus norvegicus). Such a limitation in animal species causes many difficulties in the development of new therapies for various neuropsychiatric diseases. Based on numerous scientific publications, the advantages of using a greater diversity of species in neuroscience and the disadvantages of focusing on mice and rats are presented.

Comparative Transcriptome Analyses Reveal the Role of Conserved Function in Electric Organ Convergence Across Electric Fishes

The independent origins of multiple electric organs (EOs) of fish are fascinating examples of convergent evolution. However, comparative transcriptomics of different electric fish lineages are scarce. In this study, we found that the gene expression of EOs and skeletal muscles from three lineages (Mormyroidea, Siluriformes, and Gymnotiformes) tended to cluster together based on the species of origin, irrespective of the organ from which they are derived. A pairwise comparison of differentially expressed genes (DEGs) revealed that no less than half of shared DEGs exhibited parallel expression differentiation, indicating conserved directionality of differential expression either in or between lineages, but only a few shared DEGs were identified across all focal species. Nevertheless, the functional enrichment analysis of DEGs indicated that there were more parallel gene expression changes at the level of pathways and biological functions. Therefore, we may conclude that there is no parallel evolution of the entire transcriptomes of EOs among different lineages. Further, our results support the hypothesis that it is not different genes but conserved biological functions that play a crucial role in the convergence of complex phenotypes. This study provides insight into the genetic basis underlying the EO convergent evolution; however, more studies in different cases will be needed to demonstrate whether this pattern can be extended to other cases to derive a general rule for convergent evolution.

Evidence for Weak Selective Constraint on Human Gene Expression

Variation in human complex traits is connected to variation in gene expression, and selection on complex traits can be reflected in selection on gene expression. Here, Glassberg and Gao et al. analyze polymorphic....

Gene expression variation is a major contributor to phenotypic variation in human complex traits. Selection on complex traits may therefore be reflected in constraint on gene expression. Here, we explore the effects of stabilizing selection on cis-regulatory genetic variation in humans. We analyze patterns of expression variation at copy number variants and find evidence for selection against large increases in gene expression. Using allele-specific expression (ASE) data, we further show evidence of selection against smaller-effect variants. We estimate that, across all genes, singletons in a sample of 122 individuals have ∼2.2× greater effects on expression variation than the average variant across allele frequencies. Despite their increased effect size relative to common variants, we estimate that singletons in the sample studied explain, on average, only 5% of the heritability of gene expression from cis-regulatory variants. Finally, we show that genes depleted for loss-of-function variants are also depleted for cis-eQTLs and have low levels of allelic imbalance, confirming tighter constraint on the expression levels of these genes. We conclude that constraint on gene expression is present, but has relatively weak effects on most cis-regulatory variants, thus permitting high levels of gene-regulatory genetic variation.

A quantitative framework for characterizing the evolutionary history of mammalian gene expression

The evolutionary history of a gene helps predict its function and relationship to phenotypic traits. Although sequence conservation is commonly used to decipher gene function and assess medical relevance, methods for functional inference from comparative expression data are lacking. Here, we use RNA-seq across seven tissues from 17 mammalian species to show that expression evolution across mammals is accurately modeled by the Ornstein–Uhlenbeck process, a commonly proposed model of continuous trait evolution. We apply this model to identify expression pathways under neutral, stabilizing, and directional selection. We further demonstrate novel applications of this model to quantify the extent of stabilizing selection on a gene's expression, parameterize the distribution of each gene's optimal expression level, and detect deleterious expression levels in expression data from individual patients. Our work provides a statistical framework for interpreting expression data across species and in disease.

Prediction of gene regulatory enhancers across species reveals evolutionarily conserved sequence properties

Genomic regions with gene regulatory enhancer activity turnover rapidly across mammals. In contrast, gene expression patterns and transcription factor binding preferences are largely conserved between mammalian species. Based on this conservation, we hypothesized that enhancers active in different mammals would exhibit conserved sequence patterns in spite of their different genomic locations. To investigate this hypothesis, we evaluated the extent to which sequence patterns that are predictive of enhancers in one species are predictive of enhancers in other mammalian species by training and testing two types of machine learning models. We trained support vector machine (SVM) and convolutional neural network (CNN) classifiers to distinguish enhancers defined by histone marks from the genomic background based on DNA sequence patterns in human, macaque, mouse, dog, cow, and opossum. The classifiers accurately identified many adult liver, developing limb, and developing brain enhancers, and the CNNs outperformed the SVMs. Furthermore, classifiers trained in one species and tested in another performed nearly as well as classifiers trained and tested on the same species. We observed similar cross-species conservation when applying the models to human and mouse enhancers validated in transgenic assays. This indicates that many short sequence patterns predictive of enhancers are largely conserved. The sequence patterns most predictive of enhancers in each species matched the binding motifs for a common set of TFs enriched for expression in relevant tissues, supporting the biological relevance of the learned features. Thus, despite the rapid change of active enhancer locations between mammals, cross-species enhancer prediction is often possible. Our results suggest that short sequence patterns encoding enhancer activity have been maintained across more than 180 million years of mammalian evolution.

The Grayling Genome Reveals Selection on Gene Expression Regulation after Whole-Genome Duplication

Whole-genome duplication (WGD) has been a major evolutionary driver of increased genomic complexity in vertebrates. One such event occurred in the salmonid family ∼80 Ma (Ss4R) giving rise to a plethora of structural and regulatory duplicate-driven divergence, making salmonids an exemplary system to investigate the evolutionary consequences of WGD. Here, we present a draft genome assembly of European grayling (Thymallus thymallus) and use this in a comparative framework to study evolution of gene regulation following WGD. Among the Ss4R duplicates identified in European grayling and Atlantic salmon (Salmo salar), one-third reflect nonneutral tissue expression evolution, with strong purifying selection, maintained over ∼50 Myr. Of these, the majority reflect conserved tissue regulation under strong selective constraints related to brain and neural-related functions, as well as higher-order protein–protein interactions. A small subset of the duplicates have evolved tissue regulatory expression divergence in a common ancestor, which have been subsequently conserved in both lineages, suggestive of adaptive divergence following WGD. These candidates for adaptive tissue expression divergence have elevated rates of protein coding- and promoter-sequence evolution and are enriched for immune- and lipid metabolism ontology terms. Lastly, lineage-specific duplicate divergence points toward underlying differences in adaptive pressures on expression regulation in the nonanadromous grayling versus the anadromous Atlantic salmon. Our findings enhance our understanding of the role of WGD in genome evolution and highlight cases of regulatory divergence of Ss4R duplicates, possibly related to a niche shift in early salmonid evolution.

Systems level expression correlation of Ras GTPase regulators

BACKGROUND

Proteins of the ubiquitously expressed core proteome are quantitatively correlated across multiple eukaryotic species. In addition, it was found that many protein paralogues exhibit expression anticorrelation, suggesting that the total level of protein with a given functionality must be kept constant.

METHODS

We performed Spearman's rank correlation analyses of gene expression levels for the RAS GTPase subfamily and their regulatory GEF and GAP proteins across tissues and across individuals for each tissue. A large set of published data for normal tissues from a wide range of species, human cancer tissues and human cell lines was analysed.

RESULTS

We show that although the multidomain regulatory proteins of Ras GTPases exhibit considerable tissue and individual gene expression variability, their total amounts are balanced in normal tissues. In a given tissue, the sum of activating (GEFs) and deactivating (GAPs) domains of Ras GTPases can vary considerably, but each person has balanced GEF and GAP levels. This balance is impaired in cell lines and in cancer tissues for some individuals.

CONCLUSIONS

Our results are relevant for critical considerations of knock out experiments, where functionally related homologs may compensate for the down regulation of a protein.

Similar Evolutionary Trajectories for Retrotransposon Accumulation in Mammals

The factors guiding retrotransposon insertion site preference are not well understood. Different types of retrotransposons share common replication machinery and yet occupy distinct genomic domains. Autonomous long interspersed elements accumulate in gene-poor domains and their nonautonomous short interspersed elements accumulate in gene-rich domains. To determine genomic factors that contribute to this discrepancy we analyzed the distribution of retrotransposons within the framework of chromosomal domains and regulatory elements. Using comparative genomics, we identified large-scale conserved patterns of retrotransposon accumulation across several mammalian genomes. Importantly, retrotransposons that were active after our sample-species diverged accumulated in orthologous regions. This suggested a similar evolutionary interaction between retrotransposon activity and conserved genome architecture across our species. In addition, we found that retrotransposons accumulated at regulatory element boundaries in open chromatin, where accumulation of particular retrotransposon types depended on insertion size and local regulatory element density. From our results, we propose a model where density and distribution of genes and regulatory elements canalize retrotransposon accumulation. Through conservation of synteny, gene regulation and nuclear organization, mammalian genomes with dissimilar retrotransposons follow similar evolutionary trajectories.

Complexity and conservation of regulatory landscapes underlie evolutionary resilience of mammalian gene expression

To gain insight into how mammalian gene expression is controlled by rapidly evolving regulatory elements, we jointly analysed promoter and enhancer activity with downstream transcription levels in liver samples from 15 species. Genes associated with complex regulatory landscapes generally exhibit high expression levels that remain evolutionarily stable. While the number of regulatory elements is the key driver of transcriptional output and resilience, regulatory conservation matters: elements active across mammals most effectively stabilize gene expression. In contrast, recently evolved enhancers typically contribute weakly, consistent with their high evolutionary plasticity. These effects are observed across the entire mammalian clade and are robust to potential confounders, such as the gene expression level. Using liver as a representative somatic tissue, our results illuminate how the evolutionary stability of gene expression is profoundly entwined with both the number and conservation of surrounding promoters and enhancers.

A Method for Cross-Species Visualization and Analysis of RNA-Sequence Data

In this methods article, I describe a computational workflow for cross-species visualization and comparison of mRNA-seq transcriptome profiling data. The workflow is based on gene set variation analysis (GSVA) and is illustrated using commands in the R programming language. I provide a complete step-by-step procedure for the workflow using mRNA-seq data sets from dog and human bladder cancer as an example.

Large-scale mapping of mammalian transcriptomes identifies conserved genes associated with different cell states

Distinguishing cell states based only on gene expression data remains a challenging task. This is true even for analyses within a species. In cross-species comparisons, the results obtained by different groups have varied widely. Here, we integrate RNA-seq data from more than 40 cell and tissue types of four mammalian species to identify sets of associated genes as indicators for specific cell states in each species. We employ a statistical method, TROM, to identify both protein-coding and non-coding indicators. Next, we map the cell states within each species and also between species using these indicator genes. We recapitulate known phenotypic similarity between related cell and tissue types and reveal molecular basis for their similarity. We also report novel associations between several tissues and cell types with functional support. Moreover, our identified conserved associated genes are found to be a good resource for studying cell differentiation and reprogramming. Lastly, long non-coding RNAs can serve well as associated genes to indicate cell states. We further infer the biological functions of those non-coding associated genes based on their co-expressed protein-coding genes. This study demonstrates that combining statistical modeling with public RNA-seq data can be powerful for improving our understanding of cell identity control.

Genomic dissection of conserved transcriptional regulation in intestinal epithelial cells

The intestinal epithelium serves critical physiologic functions that are shared among all vertebrates. However, it is unknown how the transcriptional regulatory mechanisms underlying these functions have changed over the course of vertebrate evolution. We generated genome-wide mRNA and accessible chromatin data from adult intestinal epithelial cells (IECs) in zebrafish, stickleback, mouse, and human species to determine if conserved IEC functions are achieved through common transcriptional regulation. We found evidence for substantial common regulation and conservation of gene expression regionally along the length of the intestine from fish to mammals and identified a core set of genes comprising a vertebrate IEC signature. We also identified transcriptional start sites and other putative regulatory regions that are differentially accessible in IECs in all 4 species. Although these sites rarely showed sequence conservation from fish to mammals, surprisingly, they drove highly conserved IEC expression in a zebrafish reporter assay. Common putative transcription factor binding sites (TFBS) found at these sites in multiple species indicate that sequence conservation alone is insufficient to identify much of the functionally conserved IEC regulatory information. Among the rare, highly sequence-conserved, IEC-specific regulatory regions, we discovered an ancient enhancer upstream from her6/HES1 that is active in a distinct population of Notch-positive cells in the intestinal epithelium. Together, these results show how combining accessible chromatin and mRNA datasets with TFBS prediction and in vivo reporter assays can reveal tissue-specific regulatory information conserved across 420 million years of vertebrate evolution. We define an IEC transcriptional regulatory network that is shared between fish and mammals and establish an experimental platform for studying how evolutionarily distilled regulatory information commonly controls IEC development and physiology.

The epithelium lining the intestine is an ancient animal tissue that serves as a primary site of nutrient absorption and interaction with microbiota. Its formation and function require complex patterns of gene transcription that vary along the intestine and in specialized intestinal epithelial cell (IEC) subtypes. However, it is unknown how the underlying transcriptional regulatory mechanisms have changed over the course of vertebrate evolution. Here, we used genome-wide profiling of mRNA levels and chromatin accessibility to identify conserved IEC genes and regulatory regions in 4 vertebrate species (zebrafish, stickleback, mouse, and human) separated from a common ancestor by 420 million years. We identified substantial similarities in genes expressed along the vertebrate intestine. These data disclosed putative conserved transcription factor binding sites (TFBS) enriched in accessible chromatin near IEC genes and in regulatory sites with accessibility restricted to IECs. Fluorescent reporter assays in transparent zebrafish showed that these regions, which frequently lacked sequence conservation, were still capable of driving conserved expression patterns. We also found a highly conserved region near mammalian and fish hes1 sufficient to drive expression in a specific population of IECs with active Notch signaling. These results establish a platform to define the conserved transcriptional networks underlying vertebrate IEC physiology.

Variable cardiac α-actin (Actc1) expression in early adult skeletal muscle correlates with promoter methylation

Different genes encode the α-actin isoforms that are predominantly expressed in heart and skeletal muscle. Mutations in the skeletal muscle α-actin gene (ACTA1) cause muscle diseases that are mostly lethal in the early postnatal period. We previously demonstrated that the disease phenotype of ACTA1 mouse models could be rescued by transgenic over-expression of cardiac α-actin (ACTC1). ACTC1 is the predominant striated α-actin isoform in the heart but is also expressed in developing skeletal muscle. To develop a translatable therapy, we investigated the genetic regulation of Actc1 expression. Using strains from The Collaborative Cross (CC) genetic resource, we found that Actc1 varies in expression by up to 24-fold in skeletal muscle. We defined significant expression quantitative trait loci (eQTL) associated with early adult Actc1 expression in soleus and heart. eQTL in both heart and soleus mapped to the Actc1 locus and replicate an eQTL mapped for Actc1 in BXD heart and quadriceps. We built on this previous work by analysing genes within the eQTL peak regions to prioritise likely candidates for modifying Actc1 expression. Additionally we interrogated the CC founder haplotype contributions to enable prioritisation of genetic variants for functional analyses. Methylation around the Actc1 transcriptional start site in early adult skeletal muscle negatively correlated with Actc1 expression in a strain-dependent manner, while other marks of regulatory potential (histone modification and chromatin accessibility) were unaltered. This study provides novel insights into the complex genetic regulation of Actc1 expression in early adult skeletal muscles.

Genomic and phenotypic signatures of climate adaptation in an Anolis lizard

Integrated knowledge on phenotype, physiology, and genomic adaptations is required to understand the effects of climate on evolution. The functional genomic basis of organismal adaptation to changes in the abiotic environment, its phenotypic consequences, and its possible convergence across vertebrates are still understudied. In this study, we use a comparative approach to verify predicted gene functions for vertebrate thermal adaptation with observed functions underlying repeated genomic adaptations in response to elevation in the lizard Anolis cybotes. We establish a direct link between recurrently evolved phenotypes and functional genomics of altitude-related climate adaptation in three highland and lowland populations in the Dominican Republic. We show that across vertebrates, genes contained in this interactome are expressed within the brain, the endocrine system, and during development. These results are relevant to elucidate the effect of global climate change across vertebrates and might aid in furthering insight into gene-environment relationships under disturbances to homeostasis.

Comparative transcriptomics in human and mouse

Cross-species comparisons of genomes, transcriptomes and gene regulation are now feasible at unprecedented resolution and throughput, enabling the comparison of human and mouse biology at the molecular level. Insights have been gained into the degree of conservation between human and mouse at the level of not only gene expression but also epigenetics and inter-individual variation. However, a number of limitations exist, including incomplete transcriptome characterization and difficulties in identifying orthologous phenotypes and cell types, which are beginning to be addressed by emerging technologies. Ultimately, these comparisons will help to identify the conditions under which the mouse is a suitable model of human physiology and disease, and optimize the use of animal models.

Putting the "Biology" Back into "Neurobiology": The Strength of Diversity in Animal Model Systems for Neuroscience Research

Current trends in neuroscience research have moved toward a reliance on rodent animal models to study most aspects of brain function. Such laboratory-reared animals are highly inbred, have been disengaged from their natural environments for generations and appear to be of limited predictive value for successful clinical outcomes. In this Perspective article, we argue that research on a rich diversity of animal model systems is fundamental to new discoveries in evolutionarily conserved core physiological and molecular mechanisms that are the foundation of human brain function. Analysis of neural circuits across phyla will reveal general computational solutions that form the basis for adaptive behavioral responses. Further, we stress that development of ethoexperimental approaches to improve our understanding of behavioral nuance will help to realign our research strategies with therapeutic goals and improve the translational validity of specific animal models. Finally, we suggest that neuroscience has a role in environmental conservation of habitat and fauna that will preserve and protect the ecological settings that drive species-specific behavioral adaptations. A rich biodiversity will enhance our understanding of human brain function and lead in unpredicted directions for development of therapeutic treatments for neurological disorders.

A guide to translation of research results from model organisms to human

A new study helps resolve a controversy about determinants of gene expression variability and might facilitate the effective translation of research results across species.

Gene-specific patterns of expression variation across organs and species

BACKGROUND

A comparison of transcriptional profiles derived from different tissues in a given species or among different species assumes that commonalities reflect evolutionarily conserved programs and that differences reflect species or tissue responses to environmental conditions or developmental program staging. Apparently conflicting results have been published regarding whether organ-specific transcriptional patterns dominate over species-specific patterns, or vice versa, making it unclear to what extent the biology of a given organism can be extrapolated to another. These studies have in common that they treat the transcriptomes monolithically, implicitly ignoring that each gene is likely to have a specific pattern of transcriptional variation across organs and species.

RESULTS

We use linear models to quantify this pattern. We find a continuum in the spectrum of expression variation: the expression of some genes varies considerably across species and little across organs, and simply reflects evolutionary distance. At the other extreme are genes whose expression varies considerably across organs and little across species; these genes are much more likely to be associated with diseases than are genes whose expression varies predominantly across species.

CONCLUSIONS

Whether transcriptomes, when considered globally, cluster preferentially according to one component or the other may not be a property of the transcriptomes, but rather a consequence of the dominant behavior of a subset of genes. Therefore, the values of the components of the variance of expression for each gene could become a useful resource when planning, interpreting, and extrapolating experimental data from mouse to humans.

Conservation of uORF repressiveness and sequence features in mouse, human and zebrafish

Upstream open reading frames (uORFs) are ubiquitous repressive genetic elements in vertebrate mRNAs. While much is known about the regulation of individual genes by their uORFs, the range of uORF-mediated translational repression in vertebrate genomes is largely unexplored. Moreover, it is unclear whether the repressive effects of uORFs are conserved across species. To address these questions, we analyse transcript sequences and ribosome profiling data from human, mouse and zebrafish. We find that uORFs are depleted near coding sequences (CDSes) and have initiation contexts that diminish their translation. Linear modelling reveals that sequence features at both uORFs and CDSes modulate the translation of CDSes. Moreover, the ratio of translation over 5′ leaders and CDSes is conserved between human and mouse, and correlates with the number of uORFs. These observations suggest that the prevalence of vertebrate uORFs may be explained by their conserved role in repressing CDS translation.

Upstream open reading frames (uORFs) can repress gene expression. Here, Guo-Liang Chew and colleagues use bioinformatics approaches to show that conservation of uORF-mediated translational repression is mediated by sequence features in human, mouse and zebrafish genomes.

Dissecting the expression relationships between RNA-binding proteins and their cognate targets in eukaryotic post-transcriptional regulatory networks

RNA-binding proteins (RBPs) are pivotal in orchestrating several steps in the metabolism of RNA in eukaryotes thereby controlling an extensive network of RBP-RNA interactions. Here, we employed CLIP (cross-linking immunoprecipitation)-seq datasets for 60 human RBPs and RIP-ChIP (RNP immunoprecipitation-microarray) data for 69 yeast RBPs to construct a network of genome-wide RBP- target RNA interactions for each RBP. We show in humans that majority (~78%) of the RBPs are strongly associated with their target transcripts at transcript level while ~95% of the studied RBPs were also found to be strongly associated with expression levels of target transcripts when protein expression levels of RBPs were employed. At transcript level, RBP - RNA interaction data for the yeast genome, exhibited a strong association for 63% of the RBPs, confirming the association to be conserved across large phylogenetic distances. Analysis to uncover the features contributing to these associations revealed the number of target transcripts and length of the selected protein-coding transcript of an RBP at the transcript level while intensity of the CLIP signal, number of RNA-Binding domains, location of the binding site on the transcript, to be significant at the protein level. Our analysis will contribute to improved modelling and prediction of post-transcriptional networks.

Divergence in gene expression within and between two closely related flycatcher species

Relatively little is known about the character of gene expression evolution as species diverge. It is for instance unclear if gene expression generally evolves in a clock‐like manner (by stabilizing selection or neutral evolution) or if there are frequent episodes of directional selection. To gain insights into the evolutionary divergence of gene expression, we sequenced and compared the transcriptomes of multiple organs from population samples of collared (Ficedula albicollis) and pied flycatchers (F. hypoleuca), two species which diverged less than one million years ago. Ordination analysis separated samples by organ rather than by species. Organs differed in their degrees of expression variance within species and expression divergence between species. Variance was negatively correlated with expression breadth and protein interactivity, suggesting that pleiotropic constraints reduce gene expression variance within species. Variance was correlated with between‐species divergence, consistent with a pattern expected from stabilizing selection and neutral evolution. Using an expression P_ST approach, we identified genes differentially expressed between species and found 16 genes uniquely expressed in one of the species. For one of these, DPP7, uniquely expressed in collared flycatcher, the absence of expression in pied flycatcher could be associated with a ≈20‐kb deletion including 11 of 13 exons. This study of a young vertebrate speciation model system expands our knowledge of how gene expression evolves as natural populations become reproductively isolated.

On the genome base composition of teleosts: the effect of environment and lifestyle

BACKGROUND

The DNA base composition is well known to be highly variable among organisms. Bio-physic studies on the effect of the GC increments on the DNA structure have shown that GC-richer DNA sequences are more bendable. The result was the keystone of the hypothesis proposing the metabolic rate as the major force driving the GC content variability, since an increased resistance to the torsion stress is mainly required during the transcription process to avoid DNA breakage. Hence, the aim of the present work is to test if both salinity and migration, suggested to affect the metabolic rate of teleostean fishes, affect the average genomic GC content as well. Moreover, since the gill surface has been reported to be a major morphological expression of metabolic rate, this parameter was also analyzed in the light of the above hypothesis.

RESULTS

Teleosts living in different environments (freshwater and seawater) and with different lifestyles (migratory and non-migratory) were analyzed studying three variables: routine metabolic rate, gill area and genomic GC-content, none of them showing a phylogenetic signal among fish species. Routine metabolic rate, specific gill area and average genomic GC were higher in seawater than freshwater species. The same trend was observed comparing migratory versus non-migratory species. Crossing salinity and lifestyle, the active migratory species living in seawater show coincidentally the highest routine metabolic rate, the highest specific gill area and the highest average genomic GC content.

CONCLUSIONS

The results clearly highlight that environmental factors (salinity) and lifestyle (migration) affect not only the physiology (i.e. the routine metabolic rate), and the morphology (i.e. gill area) of teleosts, but also basic genome feature (i.e. the GC content), thus opening to an interesting liaison among the three variables in the light of the metabolic rate hypothesis.

Divergent MLS1 Promoters Lie on a Fitness Plateau for Gene Expression

Qualitative patterns of gene activation and repression are often conserved despite an abundance of quantitative variation in expression levels within and between species. A major challenge to interpreting patterns of expression divergence is knowing which changes in gene expression affect fitness. To characterize the fitness effects of gene expression divergence, we placed orthologous promoters from eight yeast species upstream of malate synthase (MLS1) in Saccharomyces cerevisiae. As expected, we found these promoters varied in their expression level under activated and repressed conditions as well as in their dynamic response following loss of glucose repression. Despite these differences, only a single promoter driving near basal levels of expression caused a detectable loss of fitness. We conclude that the MLS1 promoter lies on a fitness plateau whereby even large changes in gene expression can be tolerated without a substantial loss of fitness.

Meta-analysis of RNA-seq expression data across species, tissues and studies

BACKGROUND

Differences in gene expression drive phenotypic differences between species, yet major organs and tissues generally have conserved gene expression programs. Several comparative transcriptomic studies have observed greater similarity in gene expression between homologous tissues from different vertebrate species than between diverse tissues of the same species. However, a recent study by Lin and colleagues reached the opposite conclusion. These studies differed in the species and tissues analyzed, and in technical details of library preparation, sequencing, read mapping, normalization, gene sets, and clustering methods.

RESULTS

To better understand gene expression evolution we reanalyzed data from four studies, including that of Lin, encompassing 6-13 tissues each from 11 vertebrate species using standardized mapping, normalization, and clustering methods. An analysis of independent data showed that the set of tissues chosen by Lin et al. were more similar to each other than those analyzed by previous studies. Comparing expression in five common tissues from the four studies, we observed that samples clustered exclusively by tissue rather than by species or study, supporting conservation of organ physiology in mammals. Furthermore, inter-study distances between homologous tissues were generally less than intra-study distances among different tissues, enabling informative meta-analyses. Notably, when comparing expression divergence of tissues over time to expression variation across 51 human GTEx tissues, we could accurately predict the clustering of expression for arbitrary pairs of tissues and species.

CONCLUSIONS

These results provide a framework for the design of future evolutionary studies of gene expression and demonstrate the utility of comparing RNA-seq data across studies.

Evidence for a common evolutionary rate in metazoan transcriptional networks

Genome sequences diverge more rapidly in mammals than in other animal lineages, such as birds or insects. However, the effect of this rapid divergence on transcriptional evolution remains unclear. Recent reports have indicated a faster divergence of transcription factor binding in mammals than in insects, but others found the reverse for mRNA expression. Here, we show that these conflicting interpretations resulted from differing methodologies. We performed an integrated analysis of transcriptional network evolution by examining mRNA expression, transcription factor binding and cis-regulatory motifs across >25 animal species, including mammals, birds and insects. Strikingly, we found that transcriptional networks evolve at a common rate across the three animal lineages. Furthermore, differences in rates of genome divergence were greatly reduced when restricting comparisons to chromatin-accessible sequences. The evolution of transcription is thus decoupled from the global rate of genome sequence evolution, suggesting that a small fraction of the genome regulates transcription.

DOI:http://dx.doi.org/10.7554/eLife.11615.001

The genetic information that makes each individual unique is encoded in DNA molecules. Cells read this molecular instruction manual by a process called transcription, in which proteins called transcription factors bind to DNA in specific places and regulate which sections of the DNA will be expressed. These 'transcripts' are active molecules that determine the cell’s – and ultimately the individual’s – characteristics. However, it is not well understood how alterations in the DNA of different individuals or species can lead to changes in where the transcription factors bind, and in which transcripts are expressed.

Carvunis, Wang, Skola et al. set out to determine if there is a relationship between how often DNA changes and how often transcription changes during the evolution of animals. The experiments examined the abundance of transcripts in the cells of a variety of animal species with close or distant evolutionary relationships. For example, the house mouse was compared to a close relative called the Algerian mouse, to another species of rodent (rat) and to humans.

The experiments show that the changes in transcript abundances are happening at similar rates in mammals, birds and insects, even though DNA changes at very different rates in these groups of animals. This similarity was also observed for other aspects of transcription, such as in changes to where transcription factors bind to DNA. The next challenges are to find out what makes transcription evolve at such similar rates in these groups of animals, and whether these findings extend to other species and to other processes in cells.

DOI:http://dx.doi.org/10.7554/eLife.11615.002

A Computational Pipeline for Cross-Species Analysis of RNA-seq Data Using R and Bioconductor

RNA sequencing (RNA-seq) has revolutionized transcriptome analysis through profiling the expression of thousands of genes at the same time. Systematic analysis of orthologous transcripts across species is critical for understanding the evolution of gene expression and uncovering important information in animal models of human diseases. Several computational methods have been published for analyzing gene expression between species, but they often lack crucial details and therefore cannot serve as a practical guide. Here, we present the first step-by-step protocol for cross-species RNA-seq analysis with a concise workflow that is largely based on the free open-source R language and Bioconductor packages. This protocol covers the entire process from short-read mapping, gene expression quantification, differential expression analysis to pathway enrichment. Many useful utilities for data visualization are included. This complete and easy-to-follow protocol provides hands-on guidance for users who are new to cross-species gene expression analysis.

A comparison of human and mouse gene co-expression networks reveals conservation and divergence at the tissue, pathway and disease levels

Background

A deeper understanding of differences and similarities in transcriptional regulation between species can uncover important information about gene functions and the role of genes in disease. Deciphering such patterns between mice and humans is especially important since mice play an essential role in biomedical research.

Results

Here, in order to characterize evolutionary changes between humans and mice, we compared gene co-expression maps to evaluate the conservation of co-expression. We show that the conservation of co-expression connectivity of homologous genes is negatively correlated with molecular evolution rates, as expected. Then we investigated evolutionary aspects of gene sets related to functions, tissues, pathways and diseases. Genes expressed in the testis, eye and skin, and those associated with regulation of transcription, olfaction, PI3K signalling, response to virus and bacteria were more divergent between mice and humans in terms of co-expression connectivity. Surprisingly, a deeper investigation of the PI3K signalling cascade revealed that its divergence is caused by the most crucial genes of this pathway, such as mTOR and AKT2. On the other hand, our analysis revealed that genes expressed in the brain and in the bone, and those associated with cell adhesion, cell cycle, DNA replication and DNA repair are most strongly conserved in terms of co-expression network connectivity as well as having a lower rate of duplication events. Genes involved in lipid metabolism and genes specific to blood showed a signature of increased co-expression connectivity in the mouse. In terms of diseases, co-expression connectivity of genes related to metabolic disorders is the most strongly conserved between mice and humans and tumor-related genes the most divergent.

Conclusions

This work contributes to discerning evolutionary patterns between mice and humans in terms of gene interactions. Conservation of co-expression is a powerful approach to identify gene targets and processes with potential similarity and divergence between mice and humans, which has implications for drug testing and other studies employing the mouse as a model organism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0534-7) contains supplementary material, which is available to authorized users.

Complex transcriptional regulation and independent evolution of fungal-like traits in a relative of animals

Cell-type specification through differential genome regulation is a hallmark of complex multicellularity. However, it remains unclear how this process evolved during the transition from unicellular to multicellular organisms. To address this question, we investigated transcriptional dynamics in the ichthyosporean Creolimax fragrantissima, a relative of animals that undergoes coenocytic development. We find that Creolimax utilizes dynamic regulation of alternative splicing, long inter-genic non-coding RNAs and co-regulated gene modules associated with animal multicellularity in a cell-type specific manner. Moreover, our study suggests that the different cell types of the three closest animal relatives (ichthyosporeans, filastereans and choanoflagellates) are the product of lineage-specific innovations. Additionally, a proteomic survey of the secretome reveals adaptations to a fungal-like lifestyle. In summary, the diversity of cell types among protistan relatives of animals and their complex genome regulation demonstrates that the last unicellular ancestor of animals was already capable of elaborate specification of cell types.

DOI:http://dx.doi.org/10.7554/eLife.08904.001

All living animals are descended from a single-celled ancestor, and understanding how these ancestors became the first multicellular animals remains a major challenge in the field of evolutionary biology. An early breakthrough towards this goal was the realization that, even though they’re mostly single-celled organisms, the closest living relatives of animals share most of the basic gene toolkit that animals use to support their multicellular lifestyles.

This shared toolkit also includes the genes that allow each specialized cell type in an animal (for example, a skin cell or liver cell) to express the subset of genes that it needs to fulfil its specific role. Discovering how the single-celled relatives of animals regulate these and other “multicellularity-related” genes during their life cycles is the next crucial step towards understanding how animals became multicellular.

Creolimax fragrantissima is a single-celled relative of animals. One stage in this organism’s life cycle involves its nucleus (which contains its genetic material) replicating multiple times without the cell itself dividing. After this stage of development, new cells are formed, each receiving with a single nucleus, and released to live freely in the environment. Characterizing how C. fragrantissima regulates which genes are expressed during these two very different stages of development could shed new light on how multicellular animals evolved to regulate their genes in specific cell types. However, little is known about these processes in C. fragrantissima.

Now, de Mendoza et al. have both sequenced C. fragrantissima’s genome and analysed which genes are expressed during the stages of its life cycle. This analysis reveals that this organism regulates its gene expression in several ways that are more commonly associated with gene regulation in multicellular animals. Furthermore, when compared to two other living relatives of animals that have brief multicellular stages in their life cycles, de Mendoza et al. found that the three organisms expressed similar genes during these similar life cycle stages.

Furthermore, like fungi, C. fragrantissima digests its food externally and then absorbs the nutrients. Using a range of techniques, de Mendoza et al. identified the proteins involved in these processes and discovered that many had evolved independently from their counterparts in fungi. Furthermore, in some cases, the genes for these proteins had actually been acquired from bacteria via a process called lateral gene transfer.

Together these findings suggest that it was likely that the last single-celled ancestor of multicellular animals already had the biological ability to create different cell types. Understanding if the cell types found in single-celled species resemble cell types from simple animals, such as sponges and comb jellies, at a molecular level is the next step towards determining what the ancestor of living animals looked like.

DOI:http://dx.doi.org/10.7554/eLife.08904.002

On the evolutionary relationship between chondrocytes and osteoblasts

Vertebrates are the only animals that produce bone, but the molecular genetic basis for this evolutionary novelty remains obscure. Here, we synthesize information from traditional evolutionary and modern molecular genetic studies in order to generate a working hypothesis on the evolution of the gene regulatory network (GRN) underlying bone formation. Since transcription factors are often core components of GRNs (i.e., kernels), we focus our analyses on Sox9 and Runx2. Our argument centers on three skeletal tissues that comprise the majority of the vertebrate skeleton: immature cartilage, mature cartilage, and bone. Immature cartilage is produced during early stages of cartilage differentiation and can persist into adulthood, whereas mature cartilage undergoes additional stages of differentiation, including hypertrophy and mineralization. Functionally, histologically, and embryologically, these three skeletal tissues are very similar, yet unique, suggesting that one might have evolved from another. Traditional studies of the fossil record, comparative anatomy and embryology demonstrate clearly that immature cartilage evolved before mature cartilage or bone. Modern molecular approaches show that the GRNs regulating differentiation of these three skeletal cell fates are similar, yet unique, just like the functional and histological features of the tissues themselves. Intriguingly, the Sox9 GRN driving cartilage formation appears to be dominant to the Runx2 GRN of bone. Emphasizing an embryological and evolutionary transcriptomic view, we hypothesize that the Runx2 GRN underlying bone formation was co-opted from mature cartilage. We discuss how modern molecular genetic experiments, such as comparative transcriptomics, can test this hypothesis directly, meanwhile permitting levels of constraint and adaptation to be evaluated quantitatively. Therefore, comparative transcriptomics may revolutionize understanding of not only the clade-specific evolution of skeletal cells, but also the generation of evolutionary novelties, providing a modern paradigm for the evolutionary process.