Is ectopic expression caused by deregulatory mutations or due to gene-regulation leaks with evolutionary potential?

Pax3 Hypomorphs Reveal Hidden Pax7 Functional Genetic Compensation in Utero

Pax3 and Pax7 transcription factors are paralogs within the Pax gene family that that are expressed in early embryos in partially overlapping expression domains and have distinct functions. Significantly, mammalian development is largely unaffected by Pax7 systemic deletion but systemic Pax3 deletion results in defects in neural tube closure, neural crest emigration, cardiac outflow tract septation, muscle hypoplasia and in utero lethality by E14. However, we previously demonstrated that Pax3 hypomorphs expressing only 20% functional Pax3 protein levels exhibit normal neural tube and heart development, but myogenesis is selectively impaired. To determine why only some Pax3-expressing cell lineages are affected and to further titrate Pax3 threshold levels required for neural tube and heart development, we generated hypomorphs containing both a hypomorphic and a null Pax3 allele. This resulted in mutants only expressing 10% functional Pax3 protein with exacerbated neural tube, neural crest and muscle defects, but still a normal heart. To examine why the cardiac neural crest appears resistant to very low Pax3 levels, we examined its paralog Pax7. Significantly, Pax7 expression is both ectopically expressed in Pax3-expressing dorsal neural tube cells and is also upregulated in the Pax3-expressing lineages. To test whether this compensatory Pax7 expression is functional, we deleted Pax7 both systemically and lineage-specifically in hypomorphs expressing only 10% Pax3. Removal of one Pax7 allele resulted in partial outflow tract defects, and complete loss of Pax7 resulted in full penetrance outflow tract defects and in utero lethality. Moreover, combinatorial loss of Pax3 and Pax7 resulted in severe craniofacial defects and a total block of neural crest cell emigration from the neural tube. Pax7^Cre lineage mapping revealed ectopic labeling of Pax3-derived neural crest tissues and within the outflow tract of the heart, experimentally confirming the observation of ectopic activation of Pax7 in 10% Pax3 hypomorphs. Finally, genetic cell ablation of Pax7^Cre-marked cells is sufficient to cause outflow tract defects in hypomorphs expressing only 10% Pax3, confirming that ectopic and induced Pax7 can play an overlapping functional genetic compensational role in both cardiac neural crest lineage and during craniofacial development, which is normally masked by the dominant role of Pax3.

Small Molecule-Induced Pancreatic β-Like Cell Development: Mechanistic Approaches and Available Strategies

Diabetes is a metabolic disease which affects not only glucose metabolism but also lipid and protein metabolism. It encompasses two major types: type 1 and 2 diabetes. Despite the different etiologies of type 1 and 2 diabetes mellitus (T1DM and T2DM, respectively), the defining features of the two forms are insulin deficiency and resistance, respectively. Stem cell therapy is an efficient method for the treatment of diabetes, which can be achieved by differentiating pancreatic β-like cells. The consistent generation of glucose-responsive insulin releasing cells remains challenging. In this review article, we present basic concepts of pancreatic organogenesis, which intermittently provides a basis for engineering differentiation procedures, mainly based on the use of small molecules. Small molecules are more auspicious than any other growth factors, as they have unique, valuable properties like cell-permeability, as well as a nonimmunogenic nature; furthermore, they offer immense benefits in terms of generating efficient functional beta-like cells. We also summarize advances in the generation of stem cell-derived pancreatic cell lineages, especially endocrine β-like cells or islet organoids. The successful induction of stem cells depends on the quantity and quality of available stem cells and the efficient use of small molecules.

Human Olfactory Receptors: Novel Cellular Functions Outside of the Nose

Olfactory receptors (ORs) are not exclusively expressed in the olfactory sensory neurons; they are also observed outside of the olfactory system in all other human tissues tested to date, including the testis, lung, intestine, skin, heart, and blood. Within these tissues, certain ORs have been determined to be exclusively expressed in only one tissue, whereas other ORs are more widely distributed in many different tissues throughout the human body. For most of the ectopically expressed ORs, limited data are available for their functional roles. They have been shown to be involved in the modulation of cell-cell recognition, migration, proliferation, the apoptotic cycle, exocytosis, and pathfinding processes. Additionally, there is a growing body of evidence that they have the potential to serve as diagnostic and therapeutic tools, as ORs are highly expressed in different cancer tissues. Interestingly, in addition to the canonical signaling pathways activated by ORs in olfactory sensory neurons, alternative pathways have been demonstrated in nonolfactory tissues. In this review, the existing data concerning the expression, as well as the physiological and pathophysiological functions, of ORs outside of the nose are highlighted to provide insights into future lines of research.

The developmental and genetic trajectory of coloration in the guppy (Poecilia reticulata)

Examining the association between trait variation and development is crucial for understanding the evolution of phenotypic differences. Male guppy ornamental caudal fin coloration is one trait that shows a striking degree of variation within and between guppy populations. Males initially have no caudal fin coloration, then gradually develop it as they reach sexual maturity. For males, there is a trade-off between female preference for caudal fin coloration and increased visibility to predators. This trade-off may reach unique endpoints in males from different predation regimes. Caudal fin coloration includes black melanin, orange/yellow pteridines or carotenoids, and shimmering iridescence. This study examined the phenotypic trajectory and genetics associated with color development. We found that black coloration always developed first, followed by orange/yellow, then iridescence. The ordering and timing of color appearance was the same regardless of predation regime. The increased expression of melanin synthesis genes correlated well with the visual appearance of black coloration, but there was no correlation between carotenoids or pteridine synthesis gene expression and the appearance of orange/yellow. The lack of orange/yellow coloration in earlier male caudal fin developmental stages may be due to reduced expression of genes underlying the development of orange/yellow xanthophores.

The White-Knight Hypothesis, or Does the Environment Limit Innovations?

Organisms often harbor latent traits that are byproducts of other adaptations. Such latent traits are not themselves adaptive but can become adaptive in the right environment. Here I discuss several examples of such traits. Their abundance suggests that environmental change rather than new mutations might often limit the origin of evolutionary adaptations and innovations. This is important, because environments can change much faster than new mutations arise. I introduce a conceptual model that distinguishes between mutation-limited and environment-limited trait origins and suggest how experiments could help discriminate between them. Wherever latent traits are plentiful, ecology rather than genetics might determine how fast new adaptations originate and thus how fast adaptive Darwinian evolution proceeds.

Proteome Profiling Outperforms Transcriptome Profiling for Coexpression Based Gene Function Prediction

Coexpression of mRNAs under multiple conditions is commonly used to infer cofunctionality of their gene products despite well-known limitations of this "guilt-by-association" (GBA) approach. Recent advancements in mass spectrometry-based proteomic technologies have enabled global expression profiling at the protein level; however, whether proteome profiling data can outperform transcriptome profiling data for coexpression based gene function prediction has not been systematically investigated. Here, we address this question by constructing and analyzing mRNA and protein coexpression networks for three cancer types with matched mRNA and protein profiling data from The Cancer Genome Atlas (TCGA) and the Clinical Proteomic Tumor Analysis Consortium (CPTAC). Our analyses revealed a marked difference in wiring between the mRNA and protein coexpression networks. Whereas protein coexpression was driven primarily by functional similarity between coexpressed genes, mRNA coexpression was driven by both cofunction and chromosomal colocalization of the genes. Functionally coherent mRNA modules were more likely to have their edges preserved in corresponding protein networks than functionally incoherent mRNA modules. Proteomic data strengthened the link between gene expression and function for at least 75% of Gene Ontology (GO) biological processes and 90% of KEGG pathways. A web application Gene2Net (http://cptac.gene2net.org) developed based on the three protein coexpression networks revealed novel gene-function relationships, such as linking ERBB2 (HER2) to lipid biosynthetic process in breast cancer, identifying PLG as a new gene involved in complement activation, and identifying AEBP1 as a new epithelial-mesenchymal transition (EMT) marker. Our results demonstrate that proteome profiling outperforms transcriptome profiling for coexpression based gene function prediction. Proteomics should be integrated if not preferred in gene function and human disease studies.

Functionalities of expressed messenger RNAs revealed from mutant phenotypes

Total messenger RNAs mRNAs that are produced from a given gene under a certain set of conditions include both functional and nonfunctional transcripts. The high prevalence of nonfunctional mRNAs that have been detected in cells has raised questions regarding the functional implications of mRNA expression patterns and divergences. Phenotypes that result from the mutagenesis of protein-coding genes have provided the most straightforward descriptions of gene functions, and such data obtained from model organisms have facilitated investigations of the functionalities of expressed mRNAs. Mutant phenotype data from mouse tissues have revealed various attributes of functional mRNAs, including tissue-specificity, strength of expression, and evolutionary conservation. In addition, the role that mRNA expression evolution plays in driving morphological evolution has been revealed from studies designed to exploit morphological and physiological phenotypes of mouse mutants. Investigations into yeast essential genes (defined by an absence of colony growth after gene deletion) have further described gene regulatory strategies that reduce protein expression noise by mediating the rates of transcription and translation. In addition to the functional significance of expressed mRNAs as described in the abovementioned findings, the functionalities of other type of RNAs (i.e., noncoding RNAs) remain to be characterized with systematic mutations and phenotyping of the DNA regions that encode these RNA molecules. WIREs RNA 2016, 7:416-427. doi: 10.1002/wrna.1329 For further resources related to this article, please visit the WIREs website.

Contrasting gene expression patterns induced by levodopa and pramipexole treatments in the rat model of Parkinson's disease

Whether the treatment of Parkinson's disease has to be initiated with levodopa or a D2 agonist like pramipexole remains debatable. Levodopa is more potent against symptoms than D2 agonists, but D2 agonists are less prone to induce motor complications and may have neuroprotective effects. Although regulation of plastic changes in striatal circuits may be the key to their different therapeutic potential, the gene expression patterns induced by de novo treatments with levodopa or D2 agonists are currently unknown. By studying the whole striatal transcriptome in a rodent model of early stage Parkinson's disease, we have identified the gene expression patterns underlying therapeutically comparable chronic treatments with levodopa or pramipexole. Despite the overall relatively small size of mRNA expression changes at the level of individual transcripts, our data show a robust and complete segregation of the transcript expression patterns induced by both treatments. Moreover, transcripts related to oxidative metabolism and mitochondrial function were enriched in levodopa-treated compared to vehicle-treated and pramipexole-treated animals, whereas transcripts related to olfactory transduction pathways were enriched in both treatment groups compared to vehicle-treated animals. Thus, our data reveal the plasticity of genetic striatal networks possibly contributing to the therapeutic effects of the most common initial treatments for Parkinson's disease, suggesting a role for oxidative stress in the long term complications induced by levodopa and identifying previously overlooked signaling cascades as potentially new therapeutic targets.

Unraveling the association between mRNA expressions and mutant phenotypes in a genome-wide assessment of mice

High-throughput gene expression profiling has revealed substantial leaky and extraneous transcription of eukaryotic genes, challenging the perceptions that transcription is strictly regulated and that changes in transcription have phenotypic consequences. To assess the functional implications of mRNA transcription directly, we analyzed mRNA expression data derived from microarrays, RNA-sequencing, and in situ hybridization, together with phenotype data of mouse mutants as a proxy of gene function at the tissue level. The results indicated that despite the presence of widespread ectopic transcription, mRNA expression and mutant phenotypes of mammalian genes or tissues remain associated. The expression-phenotype association at the gene level was particularly strong for tissue-specific genes, and the association could be underestimated due to data insufficiency and incomprehensive phenotyping of mouse mutants; the strength of expression-phenotype association at the tissue level depended on tissue functions. Mutations on genes expressed at higher levels or expressed at earlier embryonic stages more often result in abnormal phenotypes in the tissues where they are expressed. The mRNA expression profiles that have stronger associations with their phenotype profiles tend to be more evolutionarily conserved, indicating that the evolution of transcriptome and the evolution of phenome are coupled. Therefore, mutations resulting in phenotypic aberrations in expressed tissues are more likely to occur in highly transcribed genes, tissue-specific genes, genes expressed during early embryonic stages, or genes with evolutionarily conserved mRNA expression profiles.

Olfactory receptors in non-chemosensory tissues

Olfactory receptors (ORs) detect volatile chemicals that lead to the initial perception of smell in the brain. The olfactory receptor (OR) is the first protein that recognizes odorants in the olfactory signal pathway and it is present in over 1,000 genes in mice. It is also the largest member of the G protein-coupled receptors (GPCRs). Most ORs are extensively expressed in the nasal olfactory epithelium where they perform the appropriate physiological functions that fit their location. However, recent whole-genome sequencing shows that ORs have been found outside of the olfactory system, suggesting that ORs may play an important role in the ectopic expression of non-chemosensory tissues. The ectopic expressions of ORs and their physiological functions have attracted more attention recently since MOR23 and testicular hOR17-4 have been found to be involved in skeletal muscle development, regeneration, and human sperm chemotaxis, respectively. When identifying additional expression profiles and functions of ORs in non-olfactory tissues, there are limitations posed by the small number of antibodies available for similar OR genes. This review presents the results of a research series that identifies ectopic expressions and functions of ORs in non-chemosensory tissues to provide insight into future research directions. [BMB Reports 2012; 45(11): 612-622]

Extensive divergence of yeast stress responses through transitions between induced and constitutive activation

Closely related species show a high degree of differences in gene expression, but the functional significance of these differences remains unclear. Similarly, stress responses in yeast typically involve differential expression of numerous genes, and it is unclear how many of these are functionally significant. To address these issues, we compared the expression programs of four yeast species under different growth conditions, and found that the response of these species to stress has diverged extensively. On an individual gene basis, most transcriptional responses are not conserved in any pair of species, and there are very limited common responses among all four species. We present evidence that many evolutionary changes in stress responses are compensated either (i) by the response of related genes or (ii) by changes in the basal expression levels of the genes whose responses have diverged. Thus, stress-related genes are often induced upon stress in some species but maintain high levels even in the absence of stress at other species, indicating a transition between induced and constitutive activation. In addition, ~15% of the stress responses are specific to only one of the four species, with no evidence for compensating effects or stress-related annotations, and these may reflect fortuitous regulation that is unimportant for the stress response (i.e., biological noise). Frequent compensatory changes and biological noise may explain how diverged expression responses support similar physiological responses.

The Vein Patterning 1 (VEP1) gene family laterally spread through an ecological network

Lateral gene transfer (LGT) is a major evolutionary mechanism in prokaryotes. Knowledge about LGT— particularly, multicellular— eukaryotes has only recently started to accumulate. A widespread assumption sees the gene as the unit of LGT, largely because little is yet known about how LGT chances are affected by structural/functional features at the subgenic level. Here we trace the evolutionary trajectory of VEin Patterning 1, a novel gene family known to be essential for plant development and defense. At the subgenic level VEP1 encodes a dinucleotide-binding Rossmann-fold domain, in common with members of the short-chain dehydrogenase/reductase (SDR) protein family. We found: i) VEP1 likely originated in an aerobic, mesophilic and chemoorganotrophic α-proteobacterium, and was laterally propagated through nets of ecological interactions, including multiple LGTs between phylogenetically distant green plant/fungi-associated bacteria, and five independent LGTs to eukaryotes. Of these latest five transfers, three are ancient LGTs, implicating an ancestral fungus, the last common ancestor of land plants and an ancestral trebouxiophyte green alga, and two are recent LGTs to modern embryophytes. ii) VEP1's rampant LGT behavior was enabled by the robustness and broad utility of the dinucleotide-binding Rossmann-fold, which provided a platform for the evolution of two unprecedented departures from the canonical SDR catalytic triad. iii) The fate of VEP1 in eukaryotes has been different in different lineages, being ubiquitous and highly conserved in land plants, whereas fungi underwent multiple losses. And iv) VEP1-harboring bacteria include non-phytopathogenic and phytopathogenic symbionts which are non-randomly distributed with respect to the type of harbored VEP1 gene. Our findings suggest that VEP1 may have been instrumental for the evolutionary transition of green plants to land, and point to a LGT-mediated ‘Trojan Horse’ mechanism for the evolution of bacterial pathogenesis against plants. VEP1 may serve as tool for revealing microbial interactions in plant/fungi-associated environments.

Systematic analysis of posttranscriptional gene expression

Recent systems studies of gene expression have begun to dissect the layers of regulation that underlie the eukaryotic transcriptome, the combined consequence of transcriptional and posttranscriptional events. Among the regulatory layers of the transcriptome are those of the ribonome, a highly dynamic environment of ribonucleoproteins in which RNA-binding proteins (RBPs), noncoding regulatory RNAs (ncRNAs) and messenger RNAs (mRNAs) interact. While multiple mRNAs are coordinated together in groups within the ribonome of a eukaryotic cell, each individual type of mRNA consists of multiple copies, each of which has an opportunity to be a member of more than one modular group termed a posttranscriptional RNA operon or regulon (PTRO). The mRNAs associated with each PTRO encode functionally related proteins and are coordinated at the levels of RNA stability and translation by the actions of the specific RBPs and noncoding regulatory RNAs. This article examines the methods that led to the elucidation of PTROs and the coordinating mechanisms that appear to regulate the RNA components of PTROs. Moreover, the article considers the characteristics of the dynamic systems that drive PTROs and how mRNA components are bound collectively in physical 'states' to respond to cellular perturbations and diseases. In conclusion, these studies have challenged the extent to which cellular mRNA abundance can inform investigators of the functional status of a biological system. We argue that understanding the ribonome has greater potential for illuminating the underlying coordination principles of growth, differentiation, and disease.

The transcriptional landscape of cross-specific hybrids and its possible link with growth in brook charr (Salvelinus fontinalis Mitchill)

The genetic mechanisms underlying hybridization are poorly understood despite their potentially important roles in speciation processes, adaptative evolution, and agronomical innovation. In this study, transcription profiles were compared among three populations of brook charr and their hybrids using microarrays to assess the influence of hybrid origin on modes of transcription regulation inheritance and on the mechanisms underlying growth. We found that twice as many transcripts were differently expressed between the domestic population and the two wild populations (Rupert and Laval) than between wild ones, despite their deeper genetic distance. This could reflect the consequence of artificial selection during domestication. We detected that hybrids exhibited strikingly different patterns of mode of transcription regulation, being mostly additive (94%) for domestic × Rupert, and nonadditive for Laval × domestic (45.7%) and Rupert × Laval hybrids (37.5%). Both heterosis and outbreeding depression for growth were observed among the crosses. Our results indicated that prevalence of dominance in transcription regulation seems related to growth heterosis, while prevalence of transgressive transcription regulation may be more related to outbreeding depression. Our study clearly shows, for the first time in vertebrates, that the consequences of hybridization on both the transcriptome level and the phenotype are highly dependent on the specific genetic architectures of crossed populations and therefore hardly predictable.

Thoc1 deficiency compromises gene expression necessary for normal testis development in the mouse

Accumulating evidence suggests that regulation of RNA processing through an RNP-driven mechanism is important for coordinated gene expression. This hypothesis predicts that defects in RNP biogenesis will adversely affect the elaboration of specific gene expression programs. To explore the role of RNP biogenesis on mammalian development, we have characterized the phenotype of mice hypomorphic for Thoc1. Thoc1 encodes an essential component of the evolutionarily conserved TREX complex. TREX accompanies the elongating RNA polymerase II and facilitates RNP assembly and recruitment of RNA processing factors. Hypomorphic Thoc1 mice are viable despite significantly reduced Thoc1 expression in the tissues examined. While most tissues of Thoc1-deficient mice appear to develop and function normally, gametogenesis is severely compromised. Male infertility is associated with a loss in spermatocyte viability and abnormal endocrine signaling. We suggest that loss of spermatocyte viability is a consequence of defects in the expression of genes required for normal differentiation of cell types within the testes. A number of the genes affected appear to be direct targets for regulation by Thoc1. These findings support the notion that Thoc1-mediated RNP assembly contributes to the coordinated expression of genes necessary for normal differentiation and development in vivo.

Determination and validation of principal gene products

MOTIVATION

Alternative splicing has the potential to generate a wide range of protein isoforms. For many computational applications and for experimental research, it is important to be able to concentrate on the isoform that retains the core biological function. For many genes this is far from clear.

RESULTS

We have combined five methods into a pipeline that allows us to detect the principal variant for a gene. Most of the methods were based on conservation between species, at the level of both gene and protein. The five methods used were the conservation of exonic structure, the detection of non-neutral evolution, the conservation of functional residues, the existence of a known protein structure and the abundance of vertebrate orthologues. The pipeline was able to determine a principal isoform for 83% of a set of well-annotated genes with multiple variants.

Hierarchical tree snipping: clustering guided by prior knowledge

MOTIVATION

Hierarchical clustering is widely used to cluster genes into groups based on their expression similarity. This method first constructs a tree. Next this tree is partitioned into subtrees by cutting all edges at some level, thereby inducing a clustering. Unfortunately, the resulting clusters often do not exhibit significant functional coherence.

RESULTS

To improve the biological significance of the clustering, we develop a new framework of partitioning by snipping--cutting selected edges at variable levels. The snipped edges are selected to induce clusters that are maximally consistent with partially available background knowledge such as functional classifications. Algorithms for two key applications are presented: functional prediction of genes, and discovery of functionally enriched clusters of co-expressed genes. Simulation results and cross-validation tests indicate that the algorithms perform well even when the actual number of clusters differs considerably from the requested number. Performance is improved compared with a previously proposed algorithm.

AVAILABILITY

A java package is available at http://www.cs.bgu.ac.il/~dotna/ TreeSnipping

RNA regulons: coordination of post-transcriptional events

Recent findings demonstrate that multiple mRNAs are co-regulated by one or more sequence-specific RNA-binding proteins that orchestrate their splicing, export, stability, localization and translation. These and other observations have given rise to a model in which mRNAs that encode functionally related proteins are coordinately regulated during cell growth and differentiation as post-transcriptional RNA operons or regulons, through a ribonucleoprotein-driven mechanism. Here I describe several recently discovered examples of RNA operons in budding yeast, fruitfly and mammalian cells, and their potential importance in processes such as immune response, oxidative metabolism, stress response, circadian rhythms and disease. I close by considering the evolutionary wiring and rewiring of these combinatorial post-transcriptional gene-expression networks.

Biological clocks and the coordination theory of RNA operons and regulons

One of the regulatory models of circadian rhythms involves the oscillation of transcription and translation. Although transcription factors have been widely examined during circadian processes, posttranscriptional mechanisms are less well-studied. Several laboratories have used microarrays to detect changes in mRNA expression throughout the circadian cycle and have found that mRNAs encoding the RNA-binding proteins (RBPs) nocturnin and butyrate response factor (BRF1) undergo rhythmic changes. Nocturnin is a deadenylation enzyme that removes poly(A) from the 3' ends of mRNAs, whereas BRF1 destabilizes mRNAs encoding early response gene (ERG) transcripts that contain AU-rich sequences in their 3'-untranslated regions (UTRs). Moroni and coworkers proposed that BRF1 functions as an oscillating posttranscriptional RNA operon (PTRO) that diurnally degrades ERG transcripts in peripheral organs (Keene and Tenenbaum 2002; Benjamin et al. 2006). The PTRO model posits that mRNAs can be members of one or more discrete functionally related subsets of mRNAs as determined by cis elements in mRNA and trans-acting RBPs or microRNAs that collectively recognize these cis elements (Keene 2007). This chapter describes the basis of posttranscriptional coordination by RNA operons and their potential for horizontal transfer among cells and discusses the potential for RBPs and microRNAs to participate in coordinating circadian rhythms and other biological clocks.

Origins and Evolution of Spliceosomal Introns

Research into the origins of introns is at a critical juncture in the resolution of theories on the evolution of early life (which came first, RNA or DNA?), the identity of LUCA (the last universal common ancestor, was it prokaryotic- or eukaryotic-like?), and the significance of noncoding nucleotide variation. One early notion was that introns would have evolved as a component of an efficient mechanism for the origin of genes. But alternative theories emerged as well. From the debate between the "introns-early" and "introns-late" theories came the proposal that introns arose before the origin of genetically encoded proteins and DNA, and the more recent "introns-first" theory, which postulates the presence of introns at that early evolutionary stage from a reconstruction of the "RNA world." Here we review seminal and recent ideas about intron origins. Recent discoveries about the patterns and causes of intron evolution make this one of the most hotly debated and exciting topics in molecular evolutionary biology today.

Analysing the developing brain transcriptome with the GenePaint platform

We discuss technical means by which the complexity of gene and protein signalling cascades can be projected onto the complex structure of the mammalian brain. We argue that this requires both robotics and novel computational tools to register images of gene expression, annotate expression patterns and quantify gene expression. When sufficiently enriched and detailed, such gene expression/neuroanatomical atlases are hypothesis-generating tools and contain in themselves much of the information needed to investigate function in normal and genetically or otherwise modified brains. To be successful and useful, data-rich and comprehensive gene expression/neuroanatomical atlases have to be web accessible and structured in a way that allows the application of data exploration and mining tools.

Widespread ectopic expression of olfactory receptor genes

BACKGROUND

Olfactory receptors (ORs) are the largest gene family in the human genome. Although they are expected to be expressed specifically in olfactory tissues, some ectopic expression has been reported, with special emphasis on sperm and testis. The present study systematically explores the expression patterns of OR genes in a large number of tissues and assesses the potential functional implication of such ectopic expression.

RESULTS

We analyzed the expression of hundreds of human and mouse OR transcripts, via EST and microarray data, in several dozens of human and mouse tissues. Different tissues had specific, relatively small OR gene subsets which had particularly high expression levels. In testis, average expression was not particularly high, and very few highly expressed genes were found, none corresponding to ORs previously implicated in sperm chemotaxis. Higher expression levels were more common for genes with a non-OR genomic neighbor. Importantly, no correlation in expression levels was detected for human-mouse orthologous pairs. Also, no significant difference in expression levels was seen between intact and pseudogenized ORs, except for the pseudogenes of subfamily 7E which has undergone a human-specific expansion.

CONCLUSION

The OR superfamily as a whole, show widespread, locus-dependent and heterogeneous expression, in agreement with a neutral or near neutral evolutionary model for transcription control. These results cannot reject the possibility that small OR subsets might play functional roles in different tissues, however considerable care should be exerted when offering a functional interpretation for ectopic OR expression based only on transcription information.

Similar gene expression profiles do not imply similar tissue functions

Although similarities in gene expression among tissues are commonly inferred to reflect functional constraints, this has never been formally tested. Furthermore, it is unclear which evolutionary processes are responsible for the observed similarities. When examining genome-wide expression data in mouse, we found that patterns of expression similarity between tissues extend to genes that are unlikely to function in the tissues. Thus, ectopic expression can seem coordinated across tissues. This indicates that knowledge of gene expression patterns per se is insufficient to infer gene function. Ectopic expression is possibly explained as expression leakage, caused by spreading of chromatin modifications or the transcription apparatus into neighboring genes.

Models of spliceosomal intron proliferation in the face of widespread ectopic expression

It is now certain that today living organisms can acquire new spliceosomal introns in their genes. The proposed sources of spliceosomal introns are exons, transposons, and other introns, including spliceosomal and group II self-splicing introns. Spliceosomal introns are thought to be the most likely source, because the inserted sequence would immediately be endowed with the essential set of intron recognition sequences, thereby preventing the deleterious effects associated with incorrect splicing. The most obvious spliceosomal intron duplication pathways involve an RNA transcript intermediate step. Therefore, for a spliceosomal intron to be originated by duplication, either the source gene from which the novel intron is derived, or that gene and the recipient gene, which contains the novel intron, would need to be expressed in the germ line. Intron proliferation surveys indicate that putative intron duplicate-containing genes do not always match detectable expression in the germ line, which casts doubt on the generality of the duplication model. However, judging mechanisms of intron gain (or loss) from present-day gene expression profiles could be erroneous, if expression patterns were different at the time the introns arose. In fact, this may likely be so in most cases. Ectopic expression, i.e., the expression of genes at times and locations where the target gene is not known to have a function, is a much more common phenomenon than previously realized. We conclude with a speculation on a possible interplay between spliceosomal introns and ectopic expression at the origin of multicellularity.

Genomic fossils as a snapshot of the human transcriptome

Processed pseudogenes (PPGs) are cDNA sequences that were generated through reverse transcription of mature, spliced mRNAs and have subsequently been reinserted at a new genomic location. These cDNA sequences are usually no longer transcribed and are considered "dead on arrival." Here we show that PPGs can be used to generate a map of the transcriptome. By analyzing thousands of human PPGs, we were able to discover hundreds of transcript variants so far unidentified. An experimental verification of a subset of these variants by RT-PCR indicates that most of them are still active in the human transcriptome. Furthermore, we demonstrate that PPGs can enable the identification of ancient splice variants that were expressed ancestrally but are now extinct. Our results show that the genome itself carries a "virtual cDNA library" that can readily be used to analyze both present and ancestral transcripts. Our approach can be applied to sequenced metazoan genomes to computationally annotate splicing variation even when expressed sequences are unavailable.

Evolutionary Conservation of Expression Profiles Between Human and Mouse Orthologous Genes

Mouse models are often used to study human genes because it is believed that the expression and function are similar for the majority of orthologous genes between the two species. However, recent comparisons of microarray data from thousands of orthologous human and mouse genes suggested rapid evolution of gene expression profiles under minimal or no selective constraint. These findings appear to contradict non-array-based observations from many individual genes and imply the uselessness of mouse models for studying human genes. Because absolute levels of gene expression are not comparable between species when the data are generated by species-specific microarrays, use of relative mRNA abundance among tissues (RA) is preferred to that of absolute expression signals. We thus reanalyze human and mouse genome-wide gene expression data generated by oligonucleotide microarrays. We show that the mean correlation coefficient among expression profiles detected by different probe sets of the same gene is only 0.38 for humans and 0.28 for mice, indicating that current measures of expression divergence are flawed because the large estimation error (discrepancy in expression signal detected by different probe sets of the same gene) is mistakenly included in the between-species divergence. When this error is subtracted, 84% of human-mouse orthologous gene pairs show significantly lower expression divergence than that of random gene pairs. In contrast to a previous finding, but consistent with the common sense, expression profiles of orthologous tissues between species are more similar to each other than to those of nonorthologous tissues. Furthermore, the evolutionary rate of expression divergence and that of coding sequence divergence are found to be weakly, but significantly positively correlated, when RA and the Euclidean distance are used to measure expression-profile divergence. These results highlight the importance of proper consideration of various estimation errors in comparing the microarray data between species.