Mixed-Model Reanalysis of Primate Data Suggests Tissue and Species Biases in Oligonucleotide-Based Gene Expression Profiles

Transcriptome analysis revealed potential mechanisms of channel catfish growth advantage over blue catfish in a tank culture environment

Channel catfish (Ictalurus punctatus) and blue catfish (Ictalurus furcatus) are two economically important freshwater aquaculture species in the United States, with channel catfish contributing to nearly half of the country's aquaculture production. While differences in economic traits such as growth rate and disease resistance have been noted, the extent of transcriptomic variance across various tissues between these species remains largely unexplored. The hybridization of female channel catfish with male blue catfish has led to the development of superior hybrid catfish breeds that exhibit enhanced growth rates and improved disease resistance, which dominate more than half of the total US catfish production. While hybrid catfish have significant growth advantages in earthen ponds, channel catfish were reported to grow faster in tank culture environments. In this study, we confirmed channel fish's superiority in growth over blue catfish in 60-L tanks at 10.8 months of age (30.3 g and 11.6 g in this study, respectively; p < 0.001). In addition, we conducted RNA sequencing experiments and established transcriptomic resources for the heart, liver, intestine, mucus, and muscle of both species. The number of expressed genes varied across tissues, ranging from 5,036 in the muscle to over 20,000 in the mucus. Gene Ontology analysis has revealed the functional specificity of differentially expressed genes within their respective tissues, with significant pathway enrichment in metabolic pathways, immune activity, and stress responses. Noteworthy tissue-specific marker genes, including lrrc10, fabp2, myog, pth1a, hspa9, cyp21a2, agt, and ngtb, have been identified. This transcriptome resource is poised to support future investigations into the molecular mechanisms underlying environment-dependent heterosis and advance genetic breeding efforts of hybrid catfish.

High rates of evolution preceded shifts to sex-biased gene expression in Leucadendron, the most sexually dimorphic angiosperms

Differences between males and females are usually more subtle in dioecious plants than animals, but strong sexual dimorphism has evolved convergently in the South African Cape plant genus Leucadendron. Such sexual dimorphism in leaf size is expected largely to be due to differential gene expression between the sexes. We compared patterns of gene expression in leaves among 10 Leucadendron species across the genus. Surprisingly, we found no positive association between sexual dimorphism in morphology and the number or the percentage of sex-biased genes (SBGs). Sex bias in most SBGs evolved recently and was species specific. We compared rates of evolutionary change in expression for genes that were sex biased in one species but unbiased in others and found that SBGs evolved faster in expression than unbiased genes. This greater rate of expression evolution of SBGs, also documented in animals, might suggest the possible role of sexual selection in the evolution of gene expression. However, our comparative analysis clearly indicates that the more rapid rate of expression evolution of SBGs predated the origin of bias, and shifts towards bias were depleted in signatures of adaptation. Our results are thus more consistent with the view that sex bias is simply freer to evolve in genes less subject to constraints in expression level.

Experimental ischaemic stroke induces transient cardiac atrophy and dysfunction

BACKGROUND

Stroke can lead to cardiac dysfunction in patients, but the mechanisms underlying the interaction between the injured brain and the heart are poorly understood. The objective of the study is to investigate the effects of experimental murine stroke on cardiac function and molecular signalling in the heart.

METHODS AND RESULTS

Mice were subjected to filament-induced left middle cerebral artery occlusion for 30 or 60 min or sham surgery and underwent repetitive micro-echocardiography. Left ventricular contractility was reduced early (24-72 h) but not late (2 months) after brain ischaemia. Cardiac dysfunction was accompanied by a release of high-sensitive cardiac troponin (hsTNT (ng/ml): d1: 7.0 ± 1.0 vs. 25.0 ± 3.2*; d3: 7.3 ± 1.1 vs. 52.2 ± 16.7*; d14: 5.7 ± 0.8 vs. 5.2 ± 0.3; sham vs. 60 min. MCAO; mean ± SEM; *p < 0.05); reduced heart weight (heart weight/tibia length ratio: d1: 6.9 ± 0.2 vs. 6.4 ± 0.1*; d3: 6.7 ± 0.2 vs. 5.8 ± 0.1*; d14: 6.7 ± 0.2 vs. 6.4 ± 03; sham vs. 60 min. MCAO; mean ± SEM; *p < 0.05); resulting from cardiomyocyte atrophy (cardiomyocyte size: d1: 12.8% ± 0.002**; d3: 13.5% ± 0.002**; 14d: 6.3% ± 0.003*; 60 min. MCAO vs. sham; mean ± SEM; **p < 0.01; *p < 0.05), accompanied by increased atrogin-1 and the E3 ubiquitin ligase murf-1. Net norepinephrine but not synthesis was increased, suggesting a reduced norepinephrine release or an increase of norepinephrine re-uptake, resulting in a functional denervation. Transcriptome analysis in cardiac tissue identified the transcription factor peroxisome proliferator-activated receptor gamma as a potential mediator of stroke-induced transcriptional dysregulation involved in cardiac atrophy.

CONCLUSIONS

Stroke induces a complex molecular response in the heart muscle with immediate but transient cardiac atrophy and dysfunction.

Comparative Hippocampal Synaptic Proteomes of Rodents and Primates: Differences in Neuroplasticity-Related Proteins

Key to the human brain’s unique capacities are a myriad of neural cell types, specialized molecular expression signatures, and complex patterns of neuronal connectivity. Neurons in the human brain communicate via well over a quadrillion synapses. Their specific contribution might be key to the dynamic activity patterns that underlie primate-specific cognitive function. Recently, functional differences were described in transmission capabilities of human and rat synapses. To test whether unique expression signatures of synaptic proteins are at the basis of this, we performed a quantitative analysis of the hippocampal synaptic proteome of four mammalian species, two primates, human and marmoset, and two rodents, rat and mouse. Abundance differences down to 1.15-fold at an FDR-corrected p-value of 0.005 were reliably detected using SWATH mass spectrometry. The high measurement accuracy of SWATH allowed the detection of a large group of differentially expressed proteins between individual species and rodent vs. primate. Differentially expressed proteins between rodent and primate were found highly enriched for plasticity-related proteins.

RNA sequencing (RNA-Seq) of lymph node, spleen, and thymus transcriptome from wild Peninsular Malaysian cynomolgus macaque (Macaca fascicularis)

The cynomolgus macaque (Macaca fascicularis) is an extensively utilised nonhuman primate model for biomedical research due to its biological, behavioural, and genetic similarities to humans. Genomic information of cynomolgus macaque is vital for research in various fields; however, there is presently a shortage of genomic information on the Malaysian cynomolgus macaque. This study aimed to sequence, assemble, annotate, and profile the Peninsular Malaysian cynomolgus macaque transcriptome derived from three tissues (lymph node, spleen, and thymus) using RNA sequencing (RNA-Seq) technology. A total of 174,208,078 paired end 70 base pair sequencing reads were obtained from the Illumina Hi-Seq 2500 sequencer. The overall mapping percentage of the sequencing reads to the M. fascicularis reference genome ranged from 53–63%. Categorisation of expressed genes to Gene Ontology (GO) and KEGG pathway categories revealed that GO terms with the highest number of associated expressed genes include Cellular process, Catalytic activity, and Cell part, while for pathway categorisation, the majority of expressed genes in lymph node, spleen, and thymus fall under the Global overview and maps pathway category, while 266, 221, and 138 genes from lymph node, spleen, and thymus were respectively enriched in the Immune system category. Enriched Immune system pathways include Platelet activation pathway, Antigen processing and presentation, B cell receptor signalling pathway, and Intestinal immune network for IgA production. Differential gene expression analysis among the three tissues revealed 574 differentially expressed genes (DEG) between lymph and spleen, 5402 DEGs between lymph and thymus, and 7008 DEGs between spleen and thymus. Venn diagram analysis of expressed genes revealed a total of 2,630, 253, and 279 tissue-specific genes respectively for lymph node, spleen, and thymus tissues. This is the first time the lymph node, spleen, and thymus transcriptome of the Peninsular Malaysian cynomolgus macaque have been sequenced via RNA-Seq. Novel transcriptomic data will further enrich the present M. fascicularis genomic database and provide future research potentials, including novel transcript discovery, comparative studies, and molecular markers development.

Tailoring high-density oligonucleotide arrays for transcript profiling of different Arabidopsis thaliana accessions using a sequence-based approach

Excluding polymorphic probes from GeneChip ^® transcript profiling experiments via a sequence-based approach results in improved detection of differentially expressed genes in developing seeds of Arabidopsis thaliana accessions Col-0 and C24. GeneChip^® arrays represent a powerful tool for transcript profiling experiments. The ATH1 GeneChip^® has been designed based on the sequence of the Arabidopsis thaliana reference genome Col-0, hence the features on the array exactly match the sequences of Col-0 transcripts. In contrast, transcripts of other A. thaliana accessions or related species may show nucleotide differences and/or insertions/deletions when compared to the corresponding Col-0 transcripts, therefore, comparisons of transcript abundance involving different A. thaliana accessions or related species may be compromised for a certain number of transcripts. To tackle this limitation, a sequence-based strategy was developed. Only features on the array that were identical in sequence for the specimen to be compared were considered for transcript profiling. The impact of the proposed strategy was evaluated for transcript profiles that were established for developing seeds of A. thaliana accessions Col-0 and C24.

Evolution of the Human Nervous System Function, Structure, and Development

The nervous system-in particular, the brain and its cognitive abilities-is among humans' most distinctive and impressive attributes. How the nervous system has changed in the human lineage and how it differs from that of closely related primates is not well understood. Here, we consider recent comparative analyses of extant species that are uncovering new evidence for evolutionary changes in the size and the number of neurons in the human nervous system, as well as the cellular and molecular reorganization of its neural circuits. We also discuss the developmental mechanisms and underlying genetic and molecular changes that generate these structural and functional differences. As relevant new information and tools materialize at an unprecedented pace, the field is now ripe for systematic and functionally relevant studies of the development and evolution of human nervous system specializations.

Global analysis of regulatory divergence in the evolution of mouse alternative polyadenylation

Alternative polyadenylation (APA), which is regulated by both cis‐elements and trans‐factors, plays an important role in post‐transcriptional regulation of eukaryotic gene expression. However, comparing to the extensively studied transcription and alternative splicing, the extent of APA divergence during evolution and the relative cis‐ and trans‐contribution remain largely unexplored. To directly address these questions for the first time in mammals, by using deep sequencing‐based methods, we measured APA divergence between C57BL/6J and SPRET/EiJ mouse strains as well as allele‐specific APA pattern in their F1 hybrids. Among the 24,721 polyadenylation sites (pAs) from 7,271 genes expressing multiple pAs, we identified 3,747 pAs showing significant divergence between the two strains. After integrating the allele‐specific data from F1 hybrids, we demonstrated that these events could be predominately attributed to cis‐regulatory effects. Further systematic sequence analysis of the regions in proximity to cis‐divergent pAs revealed that the local RNA secondary structure and a poly(U) tract in the upstream region could negatively modulate the pAs usage.

Transcriptional co-factor Transducin beta-like (TBL) 1 acts as a checkpoint in pancreatic cancer malignancy

Pancreatic ductal adenocarcinoma (PDAC) is the fourth leading cause of cancer fatalities in Western societies, characterized by high metastatic potential and resistance to chemotherapy. Critical molecular mechanisms of these phenotypical features still remain unknown, thus hampering the development of effective prognostic and therapeutic measures in PDAC. Here, we show that transcriptional co-factor Transducin beta-like (TBL) 1 was over-expressed in both human and murine PDAC. Inactivation of TBL1 in human and mouse pancreatic cancer cells reduced cellular proliferation and invasiveness, correlating with diminished glucose uptake, glycolytic flux, and oncogenic PI3 kinase signaling which in turn could rescue TBL1 deficiency-dependent phenotypes. TBL1 deficiency both prevented and reversed pancreatic tumor growth, mediated transcriptional PI3 kinase inhibition, and increased chemosensitivity of PDAC cells in vivo. As TBL1 mRNA levels were also found to correlate with PI3 kinase levels and overall survival in a cohort of human PDAC patients, TBL1 was identified as a checkpoint in the malignant behavior of pancreatic cancer and its expression may serve as a novel molecular target in the treatment of human PDAC.

Tandem repeat variation in human and great ape populations and its impact on gene expression divergence

Tandem repeats (TRs) are stretches of DNA that are highly variable in length and mutate rapidly. They are thus an important source of genetic variation. This variation is highly informative for population and conservation genetics. It has also been associated with several pathological conditions and with gene expression regulation. However, genome-wide surveys of TR variation in humans and closely related species have been scarce due to technical difficulties derived from short-read technology. Here we explored the genome-wide diversity of TRs in a panel of 83 human and nonhuman great ape genomes, in a total of six different species, and studied their impact on gene expression evolution. We found that population diversity patterns can be efficiently captured with short TRs (repeat unit length, 1-5 bp). We examined the potential evolutionary role of TRs in gene expression differences between humans and primates by using 30,275 larger TRs (repeat unit length, 2-50 bp). Genes that contained TRs in the promoters, in their 3' untranslated region, in introns, and in exons had higher expression divergence than genes without repeats in the regions. Polymorphic small repeats (1-5 bp) had also higher expression divergence compared with genes with fixed or no TRs in the gene promoters. Our findings highlight the potential contribution of TRs to human evolution through gene regulation.

Tandem repeat variation in human and great ape populations and its impact on gene expression divergence

Tandem repeats (TRs) are stretches of DNA that are highly variable in length and mutate rapidly. They are thus an important source of genetic variation. This variation is highly informative for population and conservation genetics. It has also been associated with several pathological conditions and with gene expression regulation. However, genome-wide surveys of TR variation in humans and closely related species have been scarce due to technical difficulties derived from short-read technology. Here we explored the genome-wide diversity of TRs in a panel of 83 human and nonhuman great ape genomes, in a total of six different species, and studied their impact on gene expression evolution. We found that population diversity patterns can be efficiently captured with short TRs (repeat unit length, 1–5 bp). We examined the potential evolutionary role of TRs in gene expression differences between humans and primates by using 30,275 larger TRs (repeat unit length, 2–50 bp). Genes that contained TRs in the promoters, in their 3′ untranslated region, in introns, and in exons had higher expression divergence than genes without repeats in the regions. Polymorphic small repeats (1–5 bp) had also higher expression divergence compared with genes with fixed or no TRs in the gene promoters. Our findings highlight the potential contribution of TRs to human evolution through gene regulation.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Browning of white adipose tissue uncouples glucose uptake from insulin signaling

Presence of thermogenically active adipose tissue in adult humans has been inversely associated with obesity and type 2 diabetes. While it had been shown that insulin is crucial for the development of classical brown fat, its role in development and function of inducible brown-in-white (brite) adipose tissue is less clear. Here we show that insulin deficiency impaired differentiation of brite adipocytes. However, adrenergic stimulation almost fully induced the thermogenic program under these settings. Although brite differentiation of adipocytes as well as browning of white adipose tissue entailed substantially elevated glucose uptake by adipose tissue, the capacity of insulin to stimulate glucose uptake surprisingly was not higher in the brite state. Notably, in line with the insulin-independent stimulation of glucose uptake, our data revealed that brite recruitment results in induction of solute carrier family 2 (GLUT-1) expression in adipocytes and inguinal WAT. These results for the first time demonstrate that insulin signaling is neither essential for brite recruitment, nor is it improved in cells or tissues upon browning.

Transcriptomic insights into human brain evolution: acceleration, neutrality, heterochrony

Primate brain transcriptome comparisons within the last 12 years have yielded interesting but contradictory observations on how the transcriptome evolves, and its adaptive role in human cognitive evolution. Since the human-chimpanzee common ancestor, the human prefrontal cortex transcriptome seems to have evolved more than that of the chimpanzee. But at the same time, most expression differences among species, especially those observed in adults, appear as consequences of neutral evolution at cis-regulatory sites. Adaptive expression changes in the human brain may be rare events involving timing shifts, or heterochrony, in specific neurodevelopmental processes. Disentangling adaptive and neutral expression changes, and associating these with human-specific features of the brain require improved methods, comparisons across more species, and further work on comparative development.

Dynamics of transcriptome evolution in the model eukaryote Neurospora

Mounting evidence indicates that changes in the transcriptome contribute significantly to the phenotypic differentiation of closely related species. Nonetheless, further genome-wide studies, spanning a broad range of organisms, are needed to decipher the factors driving transcriptome evolution. The model Neurospora (Ascomycota) comprises a simple system for empirically studying the evolutionary dynamics of the transcriptome. Here, we studied the evolution of gene expression in Neurospora crassa and Neurospora tetrasperma and show that patterns of transcriptome evolution are connected to genome evolution, tissue type and sexual identity (mating types, mat A and mat a) in these eukaryotes. Based on the comparisons of inter- and intraspecies expression divergence, our data reveal that rapid expression divergence is more apt to occur in sexual/female (SF) than vegetative/male (VM) tissues. In addition, interspecies gene expression and protein sequence divergence were strongly correlated for SF, but not VM, tissue. A correlation between transcriptome and protein evolution parallels findings from certain animals, but not yeast, and add support for the theory that expression evolution differs fundamentally among multicellular and unicellular eukaryotes. Finally, we found that sexual identity in these hermaphroditic Neurospora species is connected to interspecies expression divergence in a tissue-dependent manner: rapid divergence occurred for mat A- and mat a-biased genes from SF and VM tissues, respectively. Based on these findings, it is hypothesized that rapid interspecies transcriptome evolution is shifting the mating types of Neurospora towards distinct female and male phenotypes, that is, sexual dimorphism.

Tempo and mode of regulatory evolution in Drosophila

Genetic changes affecting gene expression contribute to phenotypic divergence; thus, understanding how regulatory networks controlling gene expression change over time is critical for understanding evolution. Prior studies of expression differences within and between species have identified properties of regulatory divergence, but technical and biological differences among these studies make it difficult to assess the generality of these properties or to understand how regulatory changes accumulate with divergence time. Here, we address these issues by comparing gene expression among strains and species of Drosophila with a range of divergence times and use F₁ hybrids to examine inheritance patterns and disentangle cis- and trans-regulatory changes. We find that the fixation of compensatory changes has caused the regulation of gene expression to diverge more rapidly than gene expression itself. Specifically, we observed that the proportion of genes with evidence of cis-regulatory divergence has increased more rapidly with divergence time than the proportion of genes with evidence of expression differences. Surprisingly, the amount of expression divergence explained by cis-regulatory changes did not increase steadily with divergence time, as was previously proposed. Rather, one species (Drosophila sechellia) showed an excess of cis-regulatory divergence that we argue most likely resulted from positive selection in this lineage. Taken together, this work reveals not only the rate at which gene expression evolves, but also the molecular and evolutionary mechanisms responsible for this evolution.

Microarray RNA transcriptional profiling: Part I. Platforms, experimental design and standardization

This review summarizes, in a balanced and comprehensive manner, the various components of microarrays and their types, substrate architecture, platforms for microarray probe implementation, standardizations and confounders. The review is intended to familiarize the beginner with the principles of experimental design and the selection of an appropriate microarray platform. This parallel technology has revolutionized transcriptomic approaches to data profiling and has a major role in the identification of expressed genes, classification and diagnosis studies. The technology is still evolving and guidelines for standardization and reporting have been developed and are being improved.

Modeling gene expression evolution with an extended Ornstein-Uhlenbeck process accounting for within-species variation

Much of the phenotypic variation observed between even closely related species may be driven by differences in gene expression levels. The current availability of reliable techniques like RNA-Seq, which can quantify expression levels across species, has enabled comparative studies. Ornstein-Uhlenbeck (OU) processes have been proposed to model gene expression evolution as they model both random drift and stabilizing selection and can be extended to model changes in selection regimes. The OU models provide a statistical framework that allows comparisons of specific hypotheses of selective regimes, including random drift, constrained drift, and expression level shifts. In this way, inferences may be made about the mode of selection acting on the expression level of a gene. We augment this model to include within-species expression variance, allowing for modeling of nonevolutionary expression variance that could be caused by individual genetic, environmental, or technical variation. Through simulations, we explore the reliability of parameter estimates and the extent to which different selective regimes can be distinguished using phylogenies of varying size using both the typical OU model and our extended model. We find that if individual variation is not accounted for, nonevolutionary expression variation is often mistaken for strong stabilizing selection. The methods presented in this article are increasingly relevant as comparative expression data becomes more available and researchers turn to expression as a primary evolving phenotype.

TECHNICAL ADVANCES: A microarray for large-scale genomic and transcriptional analyses of the zebra finch (Taeniopygia guttata) and other passerines

The microarray technology has revolutionized biological research in the last decade. By monitoring the expression of many genes simultaneously, microarrays can elucidate gene function, as well as scan entire genomes for candidate genes encoding complex traits. However, because of high costs of sequencing and design, microarrays have largely been restricted to a few model species. Cross-species microarray (CSM) analyses, where microarrays are used for other species than the one they were designed for, have had varied success. We have conducted a CSM analysis by hybridizing genomic DNA from the common whitethroat (Sylvia communis) on a newly developed Affymetrix array designed for the zebra finch (Taeniopygia guttata), the Lund-zf array. The results indicate a very high potential for the zebra finch array to act as a CSM utility in other passerine birds. When hybridizing zebra finch genomic DNA, 98% of the gene representatives had higher signal intensities than the background cut-off, and for the common whitethroat, we found the equivalent proportion to be as high as 96%. This was surprising given the fact that finches and warblers diverged 25-50 million years ago, but may be explained by a relatively low sequence divergence between passerines (89-93%). Passerine birds are widely used in studies of ecology and evolution, and a zebra finch array that can be used for many species may have a large impact on future research directions.

'maskBAD'--a package to detect and remove Affymetrix probes with binding affinity differences

BACKGROUND

Hybridization differences caused by target sequence differences can be a confounding factor in analyzing gene expression on microarrays, lead to false positives and reduce power to detect real expression differences. We prepared an R Bioconductor compatible package to detect, characterize and remove such probes in Affymetrix 3'IVT and exon-based arrays on the basis of correlation of signal intensities from probes within probe sets.

RESULTS

Using completely mouse genomes we determined type 1 (false negatives) and type 2 (false positives) errors with high accuracy and we show that our method routinely outperforms previous methods. When detecting 76.2% of known SNP/indels in mouse expression data, we obtain at most 5.5% false positives. At the same level of false positives, best previous method detected 72.6%. We also show that probes with differing binding affinity both hinder differential expression detection and introduce artifacts in cancer-healthy tissue comparison.

CONCLUSIONS

Detection and removal of such probes should be a routine step in Affymetrix data preprocessing. We prepared a user friendly R package, compatible with Bioconductor, that allows the filtering and improving of data from Affymetrix microarrays experiments.

Inhibition of VCAM-1 expression in endothelial cells by CORM-3: the role of the ubiquitin-proteasome system, p38, and mitochondrial respiration

Carbon monoxide (CO) abrogates TNF-α-mediated inflammatory responses in endothelial cells, yet the underlying mechanism thereof is still elusive. We have previously shown that the anti-inflammatory effect of CO-releasing molecule-3 (CORM-3) is not completely mediated via deactivation of the NF-κB pathway. In this study, we sought to explore other potential mechanisms by which CORM-3 downregulates VCAM-1 expression on TNF-α-stimulated HUVECs. By genome-wide gene expression profiling and pathway analysis we studied the relevance of particular pathways for the anti-inflammatory effect of CORM-3. In CORM-3-stimulated HUVECs significant changes in expression were found for genes implicated in the proteasome and porphyrin pathways. Although proteasome activities were increased by CORM-3, proteasome inhibitors did not abolish the effect of CORM-3. Likewise, heme oxygenase-1 inhibitors did not abrogate the ability of CORM-3 to downregulate VCAM-1 expression. Interestingly, CORM-3 inhibited MAPK p38, and the p38 inhibitor SB203580 downregulated VCAM-1 expression. However, downregulation of VCAM-1 by CORM-3 occurred only at concentrations that partly inhibit ATP production and sodium azide and oligomycin paralleled the effect of CORM-3 in this regard. Our results indicate that CORM-3-induced downregulation of VCAM-1 is mediated via p38 inhibition and mitochondrial respiration, whereas the ubiquitin-proteasome system seems not to be involved.

Synergistic activation by p38MAPK and glucocorticoid signaling mediates induction of M2-like tumor-associated macrophages expressing the novel CD20 homolog MS4A8A

Tumor-associated macrophages (TAMs) represent alternatively activated (M2) macrophages that support tumor growth. Previously, we have described a special LYVE-1(+) M2 TAM subset in vitro and in vivo; gene profiling of this TAM subset identified MS4A8A as a novel TAM molecule expressed in vivo by TAM in mammary carcinoma and malignant melanoma. In vitro, Ms4a8a mRNA and MS4A8A protein expression was strongly induced in bone marrow-derived macrophages (BMDMs) by combining M2 mediators (IL-4, glucocorticoids) and tumor-conditioned media (TCM). Admixture of MS4A8A(+) TCM/IL-4/GC-treated BMDM significantly enhanced the tumor growth rate of subcutaneously transplanted TS/A mammary carcinomas. Upon forced overexpression of MS4A8A, Raw 264.7 macrophage-like cells displayed a special gene signature. Admixture of these MS4A8A(+) Raw 264.7 cells also significantly enhanced the tumor growth rate of subcutaneously transplanted mammary carcinomas. To identify the signaling pathways involved in synergistic induction of MS4A8A, the major signaling cascades with known functions in TAM were analyzed. Although inhibitors of NF-κB activation and of the MAPK JNK and ERK did not show relevant effects, the p38α/β MAPK inhibitor SB203580 strongly and highly significantly (p > 0.001) inhibited MS4A8A expression on mRNA and protein level. In addition, MS4A8A expression was restricted in M2 BMDM from mice with defective GC receptor (GR) dimerization indicating that classical GR gene regulation is mandatory for MS4A8A induction. In conclusion, expression of MS4A8A within the complex signal integration during macrophage immune responses may act to fine tune gene regulation. Furthermore, MS4A8A(+) TAM may serve as a novel cellular target for selective cancer therapy.

Cyst formation in the PKD2 (1-703) transgenic rat precedes deregulation of proliferation-related pathways

Background

Polycystic Kidney Disease is characterized by the formation of large fluid-filled cysts that eventually destroy the renal parenchyma leading to end-stage renal failure. Although remarkable progress has been made in understanding the pathologic mechanism of the disease, the precise orchestration of the early events leading to cyst formation is still unclear. Abnormal cellular proliferation was traditionally considered to be one of the primary irregularities leading to cyst initiation and growth. Consequently, many therapeutic interventions have focused on targeting this abnormal proliferation, and some have even progressed to clinical trials. However, the role of proliferation in cyst development was primarily examined at stages where cysts are already visible in the kidneys and therefore at later stages of disease development.

Methods

In this study we focused on the cystic phenotype since birth in an attempt to clarify the temporal contribution of cellular proliferation in cyst development. Using a PKD2 transgenic rat model (PKD2 (1-703)) of different ages (0-60 days after birth) we performed gene expression profiling and phenotype analysis by measuring various kidney parameters.

Results

Phenotype analysis demonstrated that renal cysts appear immediately after birth in the PKD2 transgenic rat model (PKD2 (1-703)). On the other hand, abnormal proliferation occurs at later stages of the disease as identified by gene expression profiling. Interestingly, other pathways appear to be deregulated at early stages of the disease in this PKD model. Specifically, gene expression analysis demonstrated that at day 0 the RAS system is involved. This is altered at day 6, when Wnt signaling and focal adhesion pathways are affected. However, at and after 24 days, proliferation, apoptosis, altered ECM signaling and many other factors become involved.

Conclusions

Our data suggest that cystogenesis precedes deregulation of proliferation-related pathways, suggesting that proliferation abnormalities may contribute in cyst growth rather than cyst formation.

Liver sinusoidal endothelium: a microenvironment-dependent differentiation program in rat including the novel junctional protein liver endothelial differentiation-associated protein-1

Liver sinusoidal endothelium (LSEC) is a prime example of organ-specific microvascular differentiation and functions. Disease-associated capillarization of LSEC in vivo and dedifferentiation of LSEC in vitro indicate the importance of the hepatic microenvironment. To identify the LSEC-specific molecular differentiation program in the rat we used a two-sided gene expression profiling approach comparing LSEC freshly isolated ex vivo with both lung microvascular endothelial cells (LMEC) and with LSEC cultured for 42 hours. The LSEC signature consisted of 48 genes both down-regulated in LMEC and in LSEC upon culture (fold change >7 in at least one comparison); quantitative reverse-transcription polymerase chain reaction confirmation of these genes included numerous family members and signaling pathway-associated molecules. The LSEC differentiation program comprised distinct sets of growth (Wnt2, Fzd4, 5, 9, Wls, vascular endothelial growth factors [VEGFR] 1, 2, 3, Nrp2) and transcription factors (Gata4, Lmo3, Tcfec, Maf) as well as endocytosis-related (Stabilin-1/2, Lyve1, and Ehd3) and cytoskeleton-associated molecules (Rnd3/RhoE). Specific gene induction in cultured LSEC versus freshly isolated LSEC as well as LMEC (Esm-1, Aatf) and up-regulation of gene expression to LMEC levels (CXCR4, Apelin) confirmed true transdifferentiation of LSEC in vitro. In addition, our analysis identified a novel 26-kDa single-pass transmembrane protein, liver endothelial differentiation-associated protein (Leda)-1, that was selectively expressed in all liver endothelial cells and preferentially localized to the abluminal cell surface. Upon forced overexpression in MDCK cells, Leda-1 was sorted basolaterally to E-cadherin-positive adherens junctions, suggesting functional involvement in cell adhesion and polarity.

CONCLUSION

Comparative microvascular analysis in rat identified a hepatic microenvironment-dependent LSEC-specific differentiation program including the novel junctional molecule Leda-1.

Regulatory divergence in Drosophila revealed by mRNA-seq

The regulation of gene expression is critical for organismal function and is an important source of phenotypic diversity between species. Understanding the genetic and molecular mechanisms responsible for regulatory divergence is therefore expected to provide insight into evolutionary change. Using deep sequencing, we quantified total and allele-specific mRNA expression levels genome-wide in two closely related Drosophila species (D. melanogaster and D. sechellia) and their F(1) hybrids. We show that 78% of expressed genes have divergent expression between species, and that cis- and trans-regulatory divergence affects 51% and 66% of expressed genes, respectively, with 35% of genes showing evidence of both. This is a relatively larger contribution of trans-regulatory divergence than was expected based on prior studies, and may result from the unique demographic history of D. sechellia. Genes with antagonistic cis- and trans-regulatory changes were more likely to be misexpressed in hybrids, consistent with the idea that such regulatory changes contribute to hybrid incompatibilities. In addition, cis-regulatory differences contributed more to divergent expression of genes that showed additive rather than nonadditive inheritance. A correlation between sequence similarity and the conservation of cis-regulatory activity was also observed that appears to be a general feature of regulatory evolution. Finally, we examined regulatory divergence that may have contributed to the evolution of a specific trait--divergent feeding behavior in D. sechellia. Overall, this study illustrates the power of mRNA sequencing for investigating regulatory evolution, provides novel insight into the evolution of gene expression in Drosophila, and reveals general trends that are likely to extend to other species.

Modulation of brain dead induced inflammation by vagus nerve stimulation

Because the vagus nerve is implicated in control of inflammation, we investigated if brain death (BD) causes impairment of the parasympathetic nervous system, thereby contributing to inflammation. BD was induced in rats. Anaesthetised ventilated rats (NBD) served as control. Heart rate variability (HRV) was assessed by ECG. The vagus nerve was electrically stimulated (BD + STIM) during BD. Intestine, kidney, heart and liver were recovered after 6 hours. Affymetrix chip-analysis was performed on intestinal RNA. Quantitative PCR was performed on all organs. Serum was collected to assess TNFalpha concentrations. Renal transplantations were performed to address the influence of vagus nerve stimulation on graft outcome. HRV was significantly lower in BD animals. Vagus nerve stimulation inhibited the increase in serum TNFalpha concentrations and resulted in down-regulation of a multiplicity of pro-inflammatory genes in intestinal tissue. In renal tissue vagal stimulation significantly decreased the expression of E-selectin, IL1beta and ITGA6. Renal function was significantly better in recipients that received a graft from a BD + STIM donor. Our study demonstrates impairment of the parasympathetic nervous system during BD and inhibition of serum TNFalpha through vagal stimulation. Vagus nerve stimulation variably affected gene expression in donor organs and improved renal function in recipients.

FOXO3 modulates endothelial gene expression and function by classical and alternative mechanisms

FOXO transcription factors represent targets of the phosphatidylinositol 3-kinase/protein kinase B survival pathway controlling important biological processes, such as cell cycle progression, apoptosis, vascular remodeling, stress responses, and metabolism. Recent studies suggested the existence of alternative mechanisms of FOXO-dependent gene expression beyond classical binding to a FOXO-responsive DNA-binding element (FRE). Here we analyzed the relative contribution of those mechanisms to vascular function by comparing the transcriptional and cellular responses to conditional activation of FOXO3 and a corresponding FRE-binding mutant in human primary endothelial cells. We demonstrate that FOXO3 controls expression of vascular remodeling genes in an FRE-dependent manner. In contrast, FOXO3-induced cell cycle arrest and apoptosis occurs independently of FRE binding, albeit FRE-dependent gene expression augments the proapoptotic response. These findings are supported by bioinformatical analysis, which revealed a statistical overrepresentation of cell cycle regulators and apoptosis-related genes in the group of co-regulated genes. Molecular analysis of FOXO3-induced endothelial apoptosis excluded modulators of the extrinsic death receptor pathway and demonstrated important roles for the BCL-2 family members BIM and NOXA in this process. Although NOXA essentially contributed to FRE-dependent apoptosis, BIM was effectively induced in the absence of FRE-binding, and small interfering RNA-mediated BIM depletion could rescue apoptosis induced by both FOXO3 mutants. These data suggest BIM as a critical cell type-specific mediator of FOXO3-induced endothelial apoptosis, whereas NOXA functions as an amplifying factor. Our study provides the first comprehensive analysis of alternatively regulated FOXO3 targets in relevant primary cells and underscores the importance of such genes for endothelial function and integrity.

On theoretical models of gene expression evolution with random genetic drift and natural selection

BACKGROUND

The relative contributions of natural selection and random genetic drift are a major source of debate in the study of gene expression evolution, which is hypothesized to serve as a bridge from molecular to phenotypic evolution. It has been suggested that the conflict between views is caused by the lack of a definite model of the neutral hypothesis, which can describe the long-run behavior of evolutionary change in mRNA abundance. Therefore previous studies have used inadequate analogies with the neutral prediction of other phenomena, such as amino acid or nucleotide sequence evolution, as the null hypothesis of their statistical inference.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we introduced two novel theoretical models, one based on neutral drift and the other assuming natural selection, by focusing on a common property of the distribution of mRNA abundance among a variety of eukaryotic cells, which reflects the result of long-term evolution. Our results demonstrated that (1) our models can reproduce two independently found phenomena simultaneously: the time development of gene expression divergence and Zipf's law of the transcriptome; (2) cytological constraints can be explicitly formulated to describe long-term evolution; (3) the model assuming that natural selection optimized relative mRNA abundance was more consistent with previously published observations than the model of optimized absolute mRNA abundances.

CONCLUSIONS/SIGNIFICANCE

The models introduced in this study give a formulation of evolutionary change in the mRNA abundance of each gene as a stochastic process, on the basis of previously published observations. This model provides a foundation for interpreting observed data in studies of gene expression evolution, including identifying an adequate time scale for discriminating the effect of natural selection from that of random genetic drift of selectively neutral variations.

Inherited variation in gene expression

Variation in gene expression constitutes an important source of biological variability within and between populations that is likely to contribute significantly to phenotypic diversity. Recent conceptual, technical, and methodological advances have enabled the genome-scale dissection of transcriptional variation. Here, we outline common approaches for detecting gene expression quantitative trait loci, and summarize the insights gleaned from these studies regarding the genetic architecture of transcriptional variation and the nature of regulatory alleles. Particular emphasis is placed on human studies, and we discuss experimental designs that ensure that increasingly large and complex studies continue to advance our understanding of gene expression variation. We conclude by discussing the evolution of gene expression levels, and we explore prospects for leveraging new technological developments to investigate inherited variation in gene expression in even greater depth.

On the reproducibility of results of pathway analysis in genome-wide expression studies of colorectal cancers

One of the major problems in genomics and medicine is the identification of gene networks and pathways deregulated in complex and polygenic diseases, like cancer. In this paper, we address the problem of assessing the variability of results of pathways analysis identified in different and independent genome wide expression studies, in which the same phenotypic conditions are assayed. To this end, we assessed the deregulation of 1891 curated gene sets in four independent gene expression data sets of subjects affected by colorectal cancer (CRC). In this comparison we used two well-founded statistical models for evaluating deregulation of gene networks. We found that the results of pathway analysis in expression studies are highly reproducible. Our study revealed 53 pathways identified by the two methods in all the four data sets analyzed with high statistical significance and strong biological relevance with the pathology examined. This set of pathways associated to single markers as well as to whole biological processes altered constitutes a signature of the disease which sheds light on the genetics bases of CRC.

Genetical toxicogenomics in Drosophila identifies master-modulatory loci that are regulated by developmental exposure to lead

The genetics of gene expression in recombinant inbred lines (RILs) can be mapped as expression quantitative trait loci (eQTLs). So-called "genetical genomics" studies have identified locally acting eQTLs (cis-eQTLs) for genes that show differences in steady-state RNA levels. These studies have also identified distantly acting master-modulatory trans-eQTLs that regulate tens or hundreds of transcripts (hotspots or transbands). We expand on these studies by performing genetical genomics experiments in two environments in order to identify trans-eQTL that might be regulated by developmental exposure to the neurotoxin lead. Flies from each of 75 RIL were raised from eggs to adults on either control food (made with 250 microM sodium acetate), or lead-treated food (made with 250 microM lead acetate, PbAc). RNA expression analyses of whole adult male flies (5-10 days old) were performed with Affymetrix DrosII whole genome arrays (18,952 probesets). Among the 1389 genes with cis-eQTL, there were 405 genes unique to control flies and 544 genes unique to lead-treated ones (440 genes had the same cis-eQTLs in both samples). There are 2396 genes with trans-eQTL which mapped to 12 major transbands with greater than 95 genes. Permutation analyses of the strain labels but not the expression data suggests that the total number of eQTL and the number of transbands are more important criteria for validation than the size of the transband. Two transbands, one located on the 2nd chromosome and one on the 3rd chromosome, co-regulate 33 lead-induced genes, many of which are involved in neurodevelopmental processes. For these 33 genes, rather than allelic variation at one locus exerting differential effects in two environments, we found that variation at two different loci are required for optimal effects on lead-induced expression.

Integrating multiple microarray data for cancer pathway analysis using bootstrapping K-S test

Previous applications of microarray technology for cancer research have mostly focused on identifying genes that are differentially expressed between a particular cancer and normal cells. In a biological system, genes perform different molecular functions and regulate various biological processes via interactions with other genes thus forming a variety of complex networks. Therefore, it is critical to understand the relationship (e.g., interactions) between genes across different types of cancer in order to gain insights into the molecular mechanisms of cancer. Here we propose an integrative method based on the bootstrapping Kolmogorov-Smirnov test and a large set of microarray data produced with various types of cancer to discover common molecular changes in cells from normal state to cancerous state. We evaluate our method using three key pathways related to cancer and demonstrate that it is capable of finding meaningful alterations in gene relations.

The effects of probe binding affinity differences on gene expression measurements and how to deal with them

MOTIVATION

When comparing gene expression levels between species or strains using microarrays, sequence differences between the groups can cause false identification of expression differences. Our simulated dataset shows that a sequence divergence of only 1% between species can lead to falsely reported expression differences for >50% of the transcripts-similar levels of effect have been reported previously in comparisons of human and chimpanzee expression. We propose a method for identifying probes that cause such false readings, using only the microarray data, so that problematic probes can be excluded from analysis. We then test the power of the method to detect sequence differences and to correct for falsely reported expression differences. Our method can detect 70% of the probes with sequence differences using human and chimpanzee data, while removing only 18% of probes with no sequence differences. Although only 70% of the probes with sequence differences are detected, the effect of removing probes on falsely reported expression differences is more dramatic: the method can remove 98% of the falsely reported expression differences from a simulated dataset. We argue that the method should be used even when sequence data are available.

CONTACT

lachmann@eva.mpg.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolutionary processes acting on candidate cis-regulatory regions in humans inferred from patterns of polymorphism and divergence

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.

It has been suggested that changes in gene expression may have played a more important role in the evolution of modern humans than changes in protein-coding sequences. In order to identify signatures of natural selection on candidate cis-regulatory regions, we examined single nucleotide polymorphisms obtained from the complete re-sequencing of conserved non-coding sites (CNCs) in the flanking regions of over 15,000 genes in 35 humans. Patterns of allele frequencies in CNCs indicate the presence of both positive and negative selection acting on standing variation within these candidate cis-regulatory regions, particularly for the 5′ and 3′ UTRs of genes. Gene-specific tests comparing levels of polymorphism and divergence identify several genes with strong signatures of selection on candidate cis-regulatory regions and suggest that the biological characteristics of genes subject to selection are different between coding and candidate cis-regulatory regions with respect to gene expression and function. For example, we find stronger signatures of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, which we do not observe in a concurrent analysis on protein-coding regions. Our results suggest that both positive and negative selection have acted on candidate cis-regulatory regions and that the evolution of non-coding DNA has played an important role throughout hominid evolution.

A systems biology-based gene expression classifier of glioblastoma predicts survival with solid tumors

Accurate prediction of survival of cancer patients is still a key open problem in clinical research. Recently, many large-scale gene expression clusterings have identified sets of genes reportedly predictive of prognosis; however, those gene sets shared few genes in common and were poorly validated using independent data. We have developed a systems biology-based approach by using either combined gene sets and the protein interaction network (Method A) or the protein network alone (Method B) to identify common prognostic genes based on microarray gene expression data of glioblastoma multiforme and compared with differential gene expression clustering (Method C). Validations of prediction performance show that the 23-prognostic gene classifier identified by Method A outperforms other gene classifiers identified by Methods B and C or previously reported for gliomas on 17 of 20 independent sample cohorts across five tumor types. We also find that among the 23 genes are 21 related to cellular proliferation and two related to response to stress/immune response. We further find that the increased expression of the 21 genes and the decreased expression of the other two genes are associated with poorer survival, which is supportive with the notion that cellular proliferation and immune response contribute to a significant portion of predictive power of prognostic classifiers. Our results demonstrate that the systems biology-based approach enables to identify common survival-associated genes.

SNEP: Simultaneous detection of nucleotide and expression polymorphisms using Affymetrix GeneChip

BACKGROUND

High-density short oligonucleotide microarrays are useful tools for studying biodiversity, because they can be used to investigate both nucleotide and expression polymorphisms. However, when different strains (or species) produce different signal intensities after mRNA hybridization, it is not easy to determine whether the signal intensities were affected by nucleotide or expression polymorphisms. To overcome this difficulty, nucleotide and expression polymorphisms are currently examined separately.

RESULTS

We have developed SNEP, a new method that allows simultaneous detection of both nucleotide and expression polymorphisms. SNEP involves a robust statistical procedure based on the idea that a nucleotide polymorphism observed at the probe level can be regarded as an outlier, because the nucleotide polymorphism can reduce the hybridization signal intensity. To investigate the performance of SNEP, we used three species: barley, rice and mice. In addition to the publicly available barley data, we obtained new rice and mouse data from the strains with available genome sequences. The sensitivity and false positive rate of nucleotide polymorphism detection were estimated based on the sequence information. The robustness of expression polymorphism detection against nucleotide polymorphisms was also investigated.

CONCLUSION

SNEP performed well regardless of the genome size and showed a better performance for nucleotide polymorphism detection, when compared with other previously proposed methods. The R-software 'SNEP' is available at http://www.ism.ac.jp/~fujisawa/SNEP/.

Effects of a spaceflight environment on heritable changes in wheat gene expression

Once it was established that the spaceflight environment was not a drastic impediment to plant growth, a remaining space biology question was whether long-term spaceflight exposure could cause changes in subsequent generations, even if they were returned to a normal Earth environment. In this study, we used a genomic approach to address this question. We tested whether changes in gene expression patterns occur in wheat plants that are several generations removed from growth in space, compared to wheat plants with no spaceflight exposure in their lineage. Wheat flown on Mir for 167 days in 1991 formed viable seeds back on Earth. These seeds were grown on the ground for three additional generations. Gene expression of fourth-generation Mir flight leaves was compared to that of the control leaves by using custom-made wheat microarrays. The data were evaluated using analysis of variance, and transcript abundance of each gene was contrasted among samples with t-tests. After corrections were made for multiple tests, none of the wheat genes represented on the microarrays showed a statistically significant difference in expression between wheat that has spaceflight exposure in their lineage and plants with no spaceflight exposure. This suggests that exposure to the spaceflight environment in low Earth orbit space stations does not cause significant, heritable changes in gene expression patterns in plants.

Strain-dependent differences in several reproductive traits are not accompanied by early postmating transcriptome changes in female Drosophila melanogaster

Upon mating, Drosophila melanogaster females undergo numerous alterations in their behavior and reproductive physiology that are accompanied by small-magnitude transcript-level changes in up to 1700 genes. Many of these postmating transcriptome changes are the direct result of the sperm and seminal fluid proteins (Acps) that females receive from their mates. To begin to determine if the genetic background of the female's mate contributes to the previously described gene expression changes, we assessed whether interactions between the genotypes of two commonly used laboratory strains of D. melanogaster (Canton-S and Oregon R) influence the female's postmating transcriptome as well as several pre- and postcopulatory phenotypes. We find negligible differences in the female's transcriptome at 1-3 hr postmating regardless of the strain of the male with whom she mated. However, a male x female genotype interaction significantly influenced mate selection, and, in some cases, fecundity, fertility, and hatchability. Our data support previous work suggesting that many of the early postmating changes observed in D. melanogaster females are not caused by large modifications of transcript levels. Instead, early postmating phenotypes result from preexisting receptors or pathways that are already in place upon sexual maturity.

Differences in genetic and environmental influences on the human cerebral cortex associated with development during childhood and adolescence

In this report, we present the first regional quantitative analysis of age-related differences in the heritability of cortical thickness using anatomic MRI with a large pediatric sample of twins, twin siblings, and singletons (n = 600, mean age 11.1 years, range 5-19). Regions of primary sensory and motor cortex, which develop earlier, both phylogenetically and ontologically, show relatively greater genetic effects earlier in childhood. Later developing regions within the dorsal prefrontal cortex and temporal lobes conversely show increasingly prominent genetic effects with maturation. The observation that regions associated with complex cognitive processes such as language, tool use, and executive function are more heritable in adolescents than children is consistent with previous studies showing that IQ becomes increasingly heritable with maturity(Plomin et al. 1997: Psychol Sci 8:442-447). These results suggest that both the specific cortical region and the age of the population should be taken into account when using cortical thickness as an intermediate phenotype to link genes, environment, and behavior.

Gene regulation in primates evolves under tissue-specific selection pressures

Regulatory changes have long been hypothesized to play an important role in primate evolution. To identify adaptive regulatory changes in humans, we performed a genome-wide survey for genes in which regulation has likely evolved under natural selection. To do so, we used a multi-species microarray to measure gene expression levels in livers, kidneys, and hearts from six humans, chimpanzees, and rhesus macaques. This comparative gene expression data allowed us to identify a large number of genes, as well as specific pathways, whose inter-species expression profiles are consistent with the action of stabilizing or directional selection on gene regulation. Among the latter set, we found an enrichment of genes involved in metabolic pathways, consistent with the hypothesis that shifts in diet underlie many regulatory adaptations in humans. In addition, we found evidence for tissue-specific selection pressures, as well as lower rates of protein evolution for genes in which regulation evolves under natural selection. These observations are consistent with the notion that adaptive circumscribed changes in gene regulation have fewer deleterious pleiotropic effects compared with changes at the protein sequence level.

SLEPR: a sample-level enrichment-based pathway ranking method -- seeking biological themes through pathway-level consistency

Analysis of microarray and other high throughput data often involves identification of genes consistently up or down-regulated across samples as the first step in extraction of biological meaning. This gene-level paradigm can be limited as a result of valid sample fluctuations and biological complexities. In this report, we describe a novel method, SLEPR, which eliminates this limitation by relying on pathway-level consistencies. Our method first selects the sample-level differentiated genes from each individual sample, capturing genes missed by other analysis methods, ascertains the enrichment levels of associated pathways from each of those lists, and then ranks annotated pathways based on the consistency of enrichment levels of individual samples from both sample classes. As a proof of concept, we have used this method to analyze three public microarray datasets with a direct comparison with the GSEA method, one of the most popular pathway-level analysis methods in the field. We found that our method was able to reproduce the earlier observations with significant improvements in depth of coverage for validated or expected biological themes, but also produced additional insights that make biological sense. This new method extends existing analyses approaches and facilitates integration of different types of HTP data.

Single-species microarrays and comparative transcriptomics

BACKGROUND

Prefabricated expression microarrays are currently available for only a few species but methods have been proposed to extend their application to comparisons between divergent genomes.

METHODOLOGY/PRINCIPAL FINDINGS

Here we demonstrate that the hybridization intensity of genomic DNA is a poor basis on which to select unbiased probes on Affymetrix expression arrays for studies of comparative transcriptomics, and that doing so produces spurious results. We used the Affymetrix Xenopus laevis microarray to evaluate expression divergence between X. laevis, X. borealis, and their F1 hybrids. When data are analyzed with probes that interrogate only sequences with confirmed identity in both species, we recover results that differ substantially analyses that use genomic DNA hybridizations to select probes.

CONCLUSIONS/SIGNIFICANCE

Our findings have implications for the experimental design of comparative expression studies that use single-species microarrays, and for our understanding of divergent expression in hybrid clawed frogs. These findings also highlight important limitations of single-species microarrays for studies of comparative transcriptomics of polyploid species.

Evolution of primate gene expression: drift and corrective sweeps?

Changes in gene expression play an important role in species' evolution. Earlier studies uncovered evidence that the effect of mutations on expression levels within the primate order is skewed, with many small downregulations balanced by fewer but larger upregulations. In addition, brain-expressed genes appeared to show an increased rate of evolution on the branch leading to human. However, the lack of a mathematical model adequately describing the evolution of gene expression precluded the rigorous establishment of these observations. Here, we develop mathematical tools that allow us to revisit these earlier observations in a model-testing and inference framework. We introduce a model for skewed gene-expression evolution within a phylogenetic tree and use a separate model to account for biological or experimental outliers. A Bayesian Markov chain Monte Carlo inference procedure allows us to infer the phylogeny and other evolutionary parameters, while quantifying the confidence in these inferences. Our results support previous observations; in particular, we find strong evidence for a sustained positive skew in the distribution of gene-expression changes in primate evolution. We propose a "corrective sweep" scenario to explain this phenomenon.

The environmental contribution to gene expression profiles

Microarray analysis provides a bridge between the molecular genetic analysis of model organisms in laboratory settings and studies of physiology, development, and adaptation in the wild. By sampling species across a range of environments, it is possible to gain a broad picture of the genomic response to environmental perturbation. Incorporating estimates of genetic relationships into study designs will facilitate genomic analysis of environmental plasticity by aiding the identification of major regulatory loci in natural populations.

Post-mating gene expression profiles of female Drosophila melanogaster in response to time and to four male accessory gland proteins

In Drosophila melanogaster, the genetic and molecular bases of post-mating changes in the female's behavior and physiology are poorly understood. However, DNA microarray studies have demonstrated that, shortly after mating, transcript abundance of >1700 genes is altered in the female's reproductive tract as well as in other tissues. Many of these changes are elicited by sperm and seminal fluid proteins (Acps) that males transfer to females. To further dissect the transcript-level changes that occur following mating, we examined gene expression profiles of whole female flies at four time points following copulation. We found that, soon after copulation ends, a large number of small-magnitude transcriptional changes occurred in the mated female. At later time points, larger magnitude changes were seen, although these occurred in a smaller number of genes. We then explored how four individual Acps (ovulin, Acp36DE, Acp29AB, and Acp62F) with unique functions independently affected gene expression in females shortly after mating. Consistent with their early and possibly local action within the female, ovulin and Acp36DE caused relatively few gene expression changes in whole bodies of mated females. In contrast, Acp29AB and Acp62F modulated a large number of transcriptional changes shortly after mating.

Cascading transcriptional effects of a naturally occurring frameshift mutation in Saccharomyces cerevisiae

Gene-expression variation in natural populations is widespread, and its phenotypic effects can be acted upon by natural selection. Only a few naturally segregating genetic differences associated with expression variation have been identified at the molecular level. We have identified a single nucleotide insertion in a vineyard isolate of Saccharomyces cerevisiae that has cascading effects through the gene-expression network. This allele is responsible for about 45% (103/230) of the genes that show differential gene expression among the homozygous diploid progeny produced by a vineyard isolate. Using isogenic laboratory strains, we confirm that this allele causes dramatic differences in gene-expression levels of key genes involved in amino acid biosynthesis. The mutation is a frameshift mutation in a mononucleotide run of eight consecutive T's in the coding region of the gene SSY1, which encodes a key component of a plasma-membrane sensor of extracellular amino acids. The potentially high rate of replication slippage of this mononucleotide repeat, combined with its relatively mild effects on growth rate in heterozygous genotypes, is sufficient to account for the persistence of this phenotype at low frequencies in natural populations.

Reconstructing phylogenies and phenotypes: a molecular view of human evolution

This review broadly summarizes how molecular biology has contributed to our understanding of human evolution. Molecular anthropology began in the 1960s with immunological comparisons indicating that African apes and humans were closely related and, indeed, shared a common ancestor as recently as 5 million years ago. Although initially dismissed, this finding has proven robust and numerous lines of molecular evidence now firmly place the human-ape divergence at 4-8 Ma. Resolving the trichotomy among humans, chimpanzees and gorillas took a few more decades. Despite the readily apparent physical similarities shared by African apes to the exclusion of modern humans (body hair, knuckle-walking, thin tooth enamel), the molecular support for a human-chimpanzee clade is now overwhelming. More recently, whole genome sequencing and gene mapping have shifted the focus of molecular anthropology from phylogenetic analyses to phenotypic reconstruction and functional genomics. We are starting to identify the genetic basis of the morphological, physiological and behavioural traits that distinguish modern humans from apes and apes from other primates. Most notably, recent comparative genomic analyses strongly indicate that the marked differences between modern humans and chimpanzees are likely due more to changes in gene regulation than to modifications of the genes themselves, an idea first proposed over 30 years ago. Almost weekly, press releases describe newly identified genes and regulatory elements that seem to have undergone strong positive selection along the human lineage. Loci involved in speech (e.g. FOXP2), brain development (e.g. ASPM), and skull musculature (e.g. MYH16) have been of particular interest, but some surprising candidate loci (e.g. those involved in auditory capabilities) have emerged as well. Exciting new research avenues, such as the Neanderthal Genome Project, promise that molecular analyses will continue to provide novel insights about our evolution. Ultimately, however, these molecular findings can only be understood in light of data from field sites, morphology labs, and museum collections. Indeed, molecular anthropology depends on these sources for calibrating molecular clocks and placing genetic data within the context of key morphological and ecological transitions in human evolution.

Sharing and reusing gene expression profiling data in neuroscience

As public availability of gene expression profiling data increases, it is natural to ask how these data can be used by neuroscientists. Here we review the public availability of high-throughput expression data in neuroscience and how it has been reused, and tools that have been developed to facilitate reuse. There is increasing interest in making expression data reuse a routine part of the neuroscience tool-kit, but there are a number of challenges. Data must become more readily available in public databases; efforts to encourage investigators to make data available are important, as is education on the benefits of public data release. Once released, data must be better-annotated. Techniques and tools for data reuse are also in need of improvement. Integration of expression profiling data with neuroscience-specific resources such as anatomical atlases will further increase the value of expression data.

Constraint and turnover in sex-biased gene expression in the genus Drosophila

Both genome content and deployment contribute to phenotypic differences between species. Sex is the most important difference between individuals in a species and has long been posited to be rapidly evolving. Indeed, in the Drosophila genus, traits such as sperm length, genitalia, and gonad size are the most obvious differences between species. Comparative analysis of sex-biased expression should deepen our understanding of the relationship between genome content and deployment during evolution. Using existing and newly assembled genomes, we designed species-specific microarrays to examine sex-biased expression of orthologues and species-restricted genes in D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. pseudoobscura, D. virilis and D. mojavensis. We show that averaged sex-biased expression changes accumulate monotonically over time within the genus. However, different genes contribute to expression variance within species groups compared to between groups. We observed greater turnover of species-restricted genes with male-biased expression, indicating that gene formation and extinction may play a significant part in species differences. Genes with male-biased expression also show the greatest expression and DNA sequence divergence. This higher divergence and turnover of genes with male-biased expression may be due to high transcription rates in the male germline, greater functional pleiotropy of genes expressed in females, and/or sexual competition.