Search

Abstract Background The arrangement of regulatory motifs in gene promoters, or promoter architecture, is the result of mutation and selection processes that have operated over many millions of years. In mammals, tissue-specific transcriptional regulation is related to the presence of specific protein-interacting DNA motifs in gene promoters. However, little is known about the relative location and spacing of these motifs. To fill this gap, we have performed a systematic search for motifs that show significant bias at specific promoter locations in a large collection of housekeeping and tissue-specific genes. Results We observe that promoters driving housekeeping gene expression are enriched in particular motifs with strong positional bias, such as YY1, which are of little relevance in promoters driving tissue-specific expression. We also identify a large number of motifs that show positional bias in genes expressed in a highly tissue-specific manner. They include well-known tissue-specific motifs, such as HNF1 and HNF4 motifs in liver, kidney and small intestine, or RFX motifs in testis, as well as many potentially novel regulatory motifs. Based on this analysis, we provide predictions for 559 tissue-specific motifs in mouse gene promoters. Conclusion The study shows that motif positional bias is an important feature of mammalian proximal promoters and that it affects both general and tissue-specific motifs. Motif positional constraints define very distinct promoter architectures depending on breadth of expression and type of tissue.

Large double-stranded DNA viruses deploy multiple strategies to subvert host immune defenses. Some of these tactics are mediated by viral gene products acquired by horizontal gene transfer from the corresponding hosts and shaped throughout evolution. The programmed death-1 (PD-1) receptor and its ligands, PD-L1 and PD-L2, play a pivotal role attenuating T-cell responses and regulating immune tolerance. In this study, we report the first functional PD-L1 homolog gene (De2) found in a pathogen. De2, captured by a γ-herpesvirus from its host during co-evolution around 50 million years ago, encodes a cell-surface glycoprotein that interacts with high affinity and stability with host PD-1. We also find that mutations evolved by the viral protein result in a significant loss of its ability to interact in cis with CD80, an interaction that for PD-L1:CD80 has been reported to block PD-1 inhibitory pathways. Furthermore, we demonstrate that the viral protein strongly inhibits T-cell signaling. Our observations suggest that PD-L1 homologs may enable viruses to evade T cell responses, favor their replication, and prevent excessive tissue damage. Altogether, our findings reveal a novel viral immunosuppressive strategy and highlight the importance of the modulation of the PD-1/PD-L1 axis during viral infections., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Martínez-Vicente, Poblador, Leitner, Farré, Steinberger, Engel and Angulo.)

The genesis of gene families by the capture of host genes and their subsequent duplication is a crucial process in the evolution of large DNA viruses. CD48 is a cell surface molecule that interacts via its N-terminal immunoglobulin (Ig) domain with the cell surface receptor 2B4 (CD244), regulating leukocyte cytotoxicity. We previously reported the presence of five CD48 homologs (vCD48s) in two related cytomegaloviruses, and demonstrated that one of them, A43, binds 2B4 and acts as a soluble CD48 decoy receptor impairing NK cell function. Here, we have characterized the rest of these vCD48s. We show that they are highly glycosylated proteins that display remarkably distinct features: divergent biochemical properties, cellular locations, and temporal expression kinetics. In contrast to A43, none of them interacts with 2B4. Consistent with this, molecular modeling of the N-terminal Ig domains of these vCD48s evidences notable changes as compared to CD48, suggesting that they interact with alternative targets. Accordingly, we demonstrate that one of them, S30, tightly binds CD2, a crucial T- and NK-cell adhesion and costimulatory molecule. Thus, our findings show how a key host immune receptor gene captured by a virus can be subsequently remodeled to evolve new immunoevasins with altered binding properties.

Throughout evolution, cytomegaloviruses (CMVs) have been capturing genes from their hosts, employing the derived proteins to evade host immune defenses. We have recently reported the presence of a number of CD48 homologs (vCD48s) encoded by different pathogenic viruses, including several CMVs. However, their properties and biological relevance remain as yet unexplored. CD48, a cosignaling molecule expressed on the surface of most hematopoietic cells, modulates the function of natural killer (NK) and other cytotoxic cells by binding to its natural ligand 2B4 (CD244). Here, we have characterized A43, the vCD48 exhibiting the highest amino acid sequence identity with host CD48. A43, which is encoded by owl monkey CMV, is a soluble molecule released from the cell after being proteolytically processed through its membrane proximal region. A43 is expressed with immediate-early kinetics, yielding a protein that is rapidly detected in the supernatant of infected cells. Remarkably, surface plasmon resonance assays revealed that this viral protein binds to host 2B4 with high affinity and slow dissociation rates. We demonstrate that soluble A43 is capable to abrogate host CD48:2B4 interactions. Moreover, A43 strongly binds to human 2B4 and prevents 2B4-mediated NK-cell adhesion to target cells, therefore reducing the formation of conjugates and the establishment of immunological synapses between human NK cells and CD48-expressing target cells. Furthermore, in the presence of this viral protein, 2B4-mediated cytotoxicity and IFN-γ production by NK cells are severely impaired. In summary, we propose that A43 may serve as a functional soluble CD48 decoy receptor by binding and masking 2B4, thereby impeding effective NK cell immune control during viral infections. Thus, our findings provide a novel example of the immune evasion strategies developed by viruses., Competing Interests: The authors have declared that no competing interests exist.

Human cytomegalovirus (HCMV) persistence in infected individuals relies on a plethora of mechanisms to efficiently reduce host immune responses. To that end, HCMV uses a variety of gene products, some of which have not been identified yet. Here we characterized the UL8 gene, which consists of two exons, sharing the first with the HCMV RL11 family member UL7 UL8 is a transmembrane protein with an N-terminal immunoglobulin (Ig)-like domain in common with UL7 but with an extended stalk and a distinctive cytoplasmic tail. The UL8 open reading frame gives rise to a heavily glycosylated protein predominantly expressed on the cell surface, from where it can be partially endocytosed and subsequently degraded. Infections with UL8-tagged viruses indicated that UL8 was synthesized with late-phase kinetics. By virtue of its highly conserved Ig-like domain, this viral protein interacted with a surface molecule present on activated neutrophils. Notably, when ectopically expressed in THP-1 myeloid cells, UL8 was able to significantly reduce the production of a variety of proinflammatory cytokines. Mutations in UL8 indicated that this functional effect was mediated by the cell surface expression of its Ig-like domain. To investigate the impact of the viral protein in the infection context, we engineered HCMVs lacking the UL8 gene and demonstrated that UL8 decreases the release of a large number of proinflammatory factors at late times after infection of THP-1 cells. Our data indicate that UL8 may exert an immunosuppressive role key for HCMV survival in the host. IMPORTANCE HCMV is a major pathogen that causes life-threatening diseases and disabilities in infected newborns and immunocompromised individuals. Containing one of the largest genomes among all reported human viruses, HCMV encodes an impressive repertoire of gene products. However, the functions of a large proportion of them still remain unknown, a fact that complicates the design of new therapeutic approaches to prevent or treat HCMV-associated diseases. In this report, we have conducted an extensive study of UL8 , one of the previously uncharacterized HCMV open reading frames. We found that the UL8 protein is expressed at late times postinfection and utilized by HCMV to reduce the production of proinflammatory factors by infected myeloid cells. Thus, the work presented here points to a key role of UL8 as a novel HCMV immune modulator capable of restraining host antiviral defenses., (Copyright © 2018 American Society for Microbiology.)

Human cytomegalovirus (HCMV) infection promotes the differentiation and persistent expansion of a mature NK cell subset, which displays high surface levels of the activating CD94/NKG2C NK cell receptor, together with additional distinctive phenotypic and functional features. The mechanisms underlying the development of adaptive NK cells remain uncertain but some observations support the involvement of a cognate interaction of CD94/NKG2C with ligand(s) displayed by HCMV-infected cells. To approach this issue, the heterodimer and its adaptor (DAP12) were expressed in the human Jurkat leukemia T cell line; signaling was detected by transfection of a reporter plasmid encoding for Luciferase (Luc) under NFAT/AP1-dependent control. Engagement of the receptor by solid-phase bound CD94- or NKG2C-specific monoclonal antibodies (mAbs) triggered Luc expression. Moreover, reporter activation was detectable upon interaction with HLA-E+ 721.221 (.221-AEH) cells, as well as with 721.221 cells incubated with synthetic peptides, which stabilized surface expression of endogenous HLA-E; the response was specifically antagonized by soluble NKG2C- and HLA-E-specific mAbs. By contrast, activation of Jurkat-NKG2C+ was undetectable upon interaction with Human Fetal Foreskin Fibroblasts (HFFF) infected with HCMV laboratory strains (i.e., AD169, Towne), regardless of their differential ability to preserve surface HLA-E expression. On the other hand, infection with two clinical isolates or with the endotheliotropic TB40/E strain triggered Jurkat-NKG2C+ activation; yet, this response was not inhibited by blocking mAbs and was independent of CD94/NKG2C expression. The results are discussed in the framework of previous observations supporting the hypothetical existence of specific ligand(s) for CD94/NKG2C in HCMV-infected cells.

Pathogens have developed a plethora of strategies to undermine host immune defenses in order to guarantee their survival. For large DNA viruses, these immune evasion mechanisms frequently rely on the expression of genes acquired from host genomes. Horizontally transferred genes include members of the immunoglobulin superfamily, whose products constitute the most diverse group of proteins of vertebrate genomes. Their promiscuous immunoglobulin domains, which comprise the building blocks of these molecules, are involved in a large variety of functions mediated by ligand-binding interactions. The flexible structural nature of the immunoglobulin domains makes them appealing targets for viral capture due to their capacity to generate high functional diversity. Here, we present an up-to-date review of immunoglobulin superfamily gene homologs encoded by herpesviruses, poxviruses, and adenoviruses, that include CD200, CD47, Fc receptors, interleukin-1 receptor 2, interleukin-18 binding protein, CD80, carcinoembryonic antigen-related cell adhesion molecules, and signaling lymphocyte activation molecules. We discuss their distinct structural attributes, binding properties, and functions, shaped by evolutionary pressures to disarm specific immune pathways. We include several novel genes identified from extensive genome database surveys. An understanding of the properties and modes of action of these viral proteins may guide the development of novel immune-modulatory therapeutic tools., (© 2017 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Since the discovery of the high abundance of Alu elements in the human genome, the interest for the functional significance of these retrotransposons has been increasing. Primate Alu and rodent Alu-like elements are retrotransposed by a mechanism driven by the LINE1 (L1) encoded proteins, the same machinery that generates the L1 repeats, the processed pseudogenes (PPs), and other retroelements. Apart from free Alu RNAs, Alus are also transcribed and retrotranscribed as part of cellular gene transcripts, generally embedded inside 3' untranslated regions (UTRs). Despite different proposed hypotheses, the functional implication of the presence of Alus inside 3'UTRs remains elusive. In this study we hypothesized that Alu elements in 3'UTRs could be involved in the genesis of PPs. By analyzing human genome data we discovered that the existence of 3'UTR-embedded Alu elements is overrepresented in genes source of PPs. In contrast, the presence of other retrotransposable elements in 3'UTRs does not show this PP linked overrepresentation. This research was extended to mouse and rat genomes and the results accordingly reveal overrepresentation of 3'UTR-embedded B1 (Alu-like) elements in PP parent genes. Interestingly, we also demonstrated that the overrepresentation of 3'UTR-embedded Alus is particularly significant in PP parent genes with low germline gene expression level. Finally, we provide data that support the hypothesis that the L1 machinery is also the system that herpesviruses, and possibly other large DNA viruses, use to capture host genes expressed in germline or somatic cells. Altogether our results suggest a novel role for Alu or Alu-like elements inside 3'UTRs as facilitators of the genesis of PPs, particularly in lowly expressed genes. Moreover, we propose that this L1-driven mechanism, aided by the presence of 3'UTR-embedded Alus, may also be exploited by DNA viruses to incorporate host genes to their viral genomes., Competing Interests: The authors have declared that no competing interests exist.

Radioligand binding assays to rat striatal dopamine D1 receptors showed that brain lateralization of the dopaminergic system were not due to changes in expression but in agonist affinity. D1 receptor-mediated striatal imbalance resulted from a significantly higher agonist affinity in the left striatum. D1 receptors heteromerize with dopamine D3 receptors, which are considered therapeutic targets for dyskinesia in parkinsonian patients. Expression of both D3 and D1-D3 receptor heteromers were increased in samples from 6-hydroxy-dopamine-hemilesioned rats rendered dyskinetic by treatment with 3, 4-dihydroxyphenyl-L-alanine (L-DOPA). Similar findings were obtained using striatal samples from primates. Radioligand binding studies in the presence of a D3 agonist led in dyskinetic, but not in lesioned or L-DOPA-treated rats, to a higher dopamine sensitivity. Upon D3-receptor activation, the affinity of agonists for binding to the right striatal D1 receptor increased. Excess dopamine coming from L-DOPA medication likely activates D3 receptors thus making right and left striatal D1 receptors equally responsive to dopamine. These results show that dyskinesia occurs concurrently with a right/left striatal balance in D1 receptor-mediated neurotransmission.

Unlabelled: Throughout evolution, large DNA viruses have been usurping genes from their hosts to equip themselves with proteins that restrain host immune defenses. Signaling lymphocytic activation molecule (SLAM) family (SLAMF) receptors are involved in the regulation of both innate and adaptive immunity, which occurs upon engagement with their ligands via homotypic or heterotypic interactions. Here we report a total of seven SLAMF genes encoded by the genomes of two cytomegalovirus (CMV) species, squirrel monkey CMV (SMCMV) and owl monkey CMV (OMCMV), that infect New World monkeys. Our results indicate that host genes were captured by retrotranscription at different stages of the CMV-host coevolution. The most recent acquisition led to S1 in SMCMV. S1 is a SLAMF6 homolog with an amino acid sequence identity of 97% to SLAMF6 in its ligand-binding N-terminal Ig domain. We demonstrate that S1 is a cell surface glycoprotein capable of binding to host SLAMF6. Furthermore, the OMCMV genome encodes A33, an LY9 (SLAMF3) homolog, and A43, a CD48 (SLAMF2) homolog, two soluble glycoproteins which recognize their respective cellular counterreceptors and thus are likely to be viral SLAMF decoy receptors. In addition, distinct copies of further divergent CD48 homologs were found to be encoded by both CMV genomes. Remarkably, all these molecules display a number of unique features, including cytoplasmic tails lacking characteristic SLAMF signaling motifs. Taken together, our findings indicate a novel immune evasion mechanism in which incorporation of host SLAMF receptors that retain their ligand-binding properties enables viruses to interfere with SLAMF functions and to supply themselves with convenient structural molds for expanding their immunomodulatory repertoires., Importance: The way in which viruses shape their genomes under the continual selective pressure exerted by the host immune system is central for their survival. Here, we report that New World monkey cytomegaloviruses have broadly captured and duplicated immune cell receptors of the signaling lymphocyte activation molecule (SLAM) family during host-virus coevolution. Notably, we demonstrate that several of these viral SLAMs exhibit exceptional preservation of their N-terminal immunoglobulin domains, which results in maintenance of their ligand-binding capacities. At the same time, these molecules present distinctive structural properties which include soluble forms and the absence of typical SLAM signaling motifs in their cytoplasmic domains, likely reflecting the evolutionary adaptation undergone to efficiently interfere with host SLAM family activities. The observation that the genomes of other large DNA viruses might bear SLAM family homologs further underscores the importance of these molecules as a novel class of immune regulators and as convenient scaffolds for viral evolution., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

The molecular basis of priming for L-DOPA-induced dyskinesias in Parkinson's disease (PD), which depends on the indirect pathway of motor control, is not known. In rodents, the indirect pathway contains striatopallidal GABAergic neurons that express heterotrimers composed of A(2A) adenosine, CB(1) cannabinoid and D(2) dopamine receptors that regulate dopaminergic neurotransmission. The present study was designed to investigate the expression of these heteromers in the striatum of a primate model of Parkinson's disease and to determine whether their expression and pharmacological properties are altered upon L-DOPA treatment. By using the recently developed in situ proximity ligation assay and by identification of a biochemical fingerprint, we discovered a regional distribution of A(2A)/CB(1) /D(2) receptor heteromers that predicts differential D(2)-mediated neurotransmission in the caudate-putamen of Macaca fascicularis. Whereas heteromers were abundant in the caudate nucleus of both naïve and MPTP-treated monkeys, L-DOPA treatment blunted the biochemical fingerprint and led to weak heteromer expression. These findings constitute the first evidence of altered receptor heteromer expression in pathological conditions and suggest that drugs targeting A(2A)-CB(1) -D(2) receptor heteromers may be successful to either normalize basal ganglia output or prevent L-DOPA-induced side effects., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Receptors of the signalling lymphocyte-activation molecules (SLAM) family are involved in the functional regulation of a variety of immune cells upon engagement through homotypic or heterotypic interactions amongst them. Here we show that murine cytomegalovirus (MCMV) dampens the surface expression of several SLAM receptors during the course of the infection of macrophages. By screening a panel of MCMV deletion mutants, we identified m154 as an immunoevasin that effectively reduces the cell-surface expression of the SLAM family member CD48, a high-affinity ligand for natural killer (NK) and cytotoxic T cell receptor CD244. m154 is a mucin-like protein, expressed with early kinetics, which can be found at the cell surface of the infected cell. During infection, m154 leads to proteolytic degradation of CD48. This viral protein interferes with the NK cell cytotoxicity triggered by MCMV-infected macrophages. In addition, we demonstrate that an MCMV mutant virus lacking m154 expression results in an attenuated phenotype in vivo, which can be substantially restored after NK cell depletion in mice. This is the first description of a viral gene capable of downregulating CD48. Our novel findings define m154 as an important player in MCMV innate immune regulation.

Long-term therapy with L-3,4-dihydroxyphenylalanine (L-DOPA), still the most effective treatment in Parkinson's disease (PD), is associated with severe motor complications such as dyskinesia. Experimental and clinical data have indicated that adenosine A2A receptor antagonists can provide symptomatic improvement by potentiating L-DOPA efficacy and minimizing its side effects. It is known that the G-protein-coupled adenosine A2A, cannabinoid CB1 and dopamine D2 receptors may interact and form functional A2A-CB1-D2 receptor heteromers in co-transfected cells as well as in rat striatum. These data suggest that treatment with a combination of drugs or a single compound selectively acting on A2A-CB1-D2 heteromers may represent an alternative therapeutic treatment of PD. We investigated the expression of A2A-CB1-D2 receptor heteromers in the striatum of both naïve and hemiparkinsonian rats (HPD-rats) bearing a unilateral 6-hydroxydopamine (6-OHDA) lesion, and assessed how receptor heteromer expression and biochemical properties were affected by L-DOPA treatment. Radioligand binding data showed that A2A-CB1-D2 receptor heteromers are present in the striatum of both naïve and HPD-rats. However, behavioral results indicated that the combined administration of A2A (MSX-3 or SCH58261) and CB1 (rimonabant) receptor antagonists, in the presence of L-DOPA does not produce a response different from administration of the A2A receptor antagonist alone. These behavioral results prompted identification of heteromers in L-DOPA-treated animals. Interestingly, the radioligand binding results in samples from lesioned animals suggest that the heteromer is lost following acute or chronic treatment with L-DOPA., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Recently there has been much interest in the Regulators of Calcineurin (RCAN) proteins which are important endogenous modulators of the calcineurin-NFATc signalling pathway. They have been shown to have a crucial role in cellular programmes such as the immune response, muscle fibre remodelling and memory, but also in pathological processes such as cardiac hypertrophy and neurodegenerative diseases. In vertebrates, the RCAN family form a functional subfamily of three members RCAN1, RCAN2 and RCAN3 whereas only one RCAN is present in the rest of Eukarya. In addition, RCAN genes have been shown to collocate with RUNX and CLIC genes in ACD clusters (ACD21, ACD6 and ACD1). How the RCAN genes and their clustering in ACDs evolved is still unknown. After analysing RCAN gene family evolution using bioinformatic tools, we propose that the three RCAN vertebrate genes within the ACD clusters, which evolved from single copy genes present in invertebrates and lower eukaryotes, are the result of two rounds of whole genome duplication, followed by a segmental duplication. This evolutionary scenario involves the loss or gain of some RCAN genes during evolution. In addition, we have analysed RCAN gene structure and identified the existence of several characteristic features that can be involved in RCAN evolution and gene expression regulation. These included: several transposable elements, CpG islands in the 5' region of the genes, the existence of antisense transcripts (NAT) associated with the three human genes, and considerable evidence for bidirectional promoters that regulate RCAN gene expression. Furthermore, we show that the CpG island associated with the RCAN3 gene promoter is unmethylated and transcriptionally active. All these results provide timely new insights into the molecular mechanisms underlying RCAN function and a more in depth knowledge of this gene family whose members are obvious candidates for the development of future therapies.

Under normal conditions the brain maintains a delicate balance between inputs of reward seeking controlled by neurons containing the D1-like family of dopamine receptors and inputs of aversion coming from neurons containing the D2-like family of dopamine receptors. Cocaine is able to subvert these balanced inputs by altering the cell signaling of these two pathways such that D1 reward seeking pathway dominates. Here, we provide an explanation at the cellular and biochemical level how cocaine may achieve this. Exploring the effect of cocaine on dopamine D2 receptors function, we present evidence of σ1 receptor molecular and functional interaction with dopamine D2 receptors. Using biophysical, biochemical, and cell biology approaches, we discovered that D2 receptors (the long isoform of the D2 receptor) can complex with σ1 receptors, a result that is specific to D2 receptors, as D3 and D4 receptors did not form heteromers. We demonstrate that the σ1-D2 receptor heteromers consist of higher order oligomers, are found in mouse striatum and that cocaine, by binding to σ1 -D2 receptor heteromers, inhibits downstream signaling in both cultured cells and in mouse striatum. In contrast, in striatum from σ1 knockout animals these complexes are not found and this inhibition is not seen. Taken together, these data illuminate the mechanism by which the initial exposure to cocaine can inhibit signaling via D2 receptor containing neurons, destabilizing the delicate signaling balance influencing drug seeking that emanates from the D1 and D2 receptor containing neurons in the brain.

The enzyme adenosine deaminase (ADA) is a multifunctional protein that can both degrade adenosine and bind extracellularly to adenosine receptors, acting as an allosteric modulator regulating the hormonal effects of adenosine. The molecular regions of ADA responsible for the latter are unknown. In this work, alanine scanning mutagenesis of various ADA amino acid stretches, selected through in silico docking experiments, allowed us to identify regions of the enzyme responsible for modulating both its catalytic activity and its ability to modulate agonist binding to A and A adenosine receptors (AR and AR). The combination of computational and in vitro experiments show that the structural gate to the catalytic site; i.e., the α-1 helix containing residues L58-I72 and the loop containing residues A184-I188 of ADA, were important to maintain both the catalytic efficiency of the enzyme and its action as an allosteric modulator of the adenosine receptors. These data are consistent with a predicted supramolecular assembly, in which ADA bridges AR and CD26 and are in line with the notion that the interaction of ADA with adenosine receptors has an important role in the immunosynapse. We propose that it is the ADA open form, but not the closed one, that is responsible for the functional interaction with A₁R and A₂AR.

'Orang-utan' is derived from a Malay term meaning 'man of the forest' and aptly describes the southeast Asian great apes native to Sumatra and Borneo. The orang-utan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orang-utan draft genome assembly and short read sequence data from five Sumatran and five Bornean orang-utan genomes. Our analyses reveal that, compared to other primates, the orang-utan genome has many unique features. Structural evolution of the orang-utan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe a primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orang-utan genome structure. Orang-utans have extremely low energy usage for a eutherian mammal, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400,000 years ago, is more recent than most previous studies and underscores the complexity of the orang-utan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (N(e)) expanded exponentially relative to the ancestral N(e) after the split, while Bornean N(e) declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts.

Gene duplication is a major mechanism for molecular evolutionary innovation. Young gene duplicates typically exhibit elevated rates of protein evolution and, according to a number of recent studies, increased expression divergence. However, the nature of these changes is still poorly understood. To gain novel insights into the functional consequences of gene duplication, we have undertaken an in-depth analysis of a large data set of gene families containing primate- and/or rodent-specific gene duplicates. We have found a clear tendency toward an increase in protein, promoter, and expression divergence with increasing number of duplication events undergone by each gene since the human-mouse split. In addition, gene duplication is significantly associated with a reduction in expression breadth and intensity. Interestingly, it is possible to identify three main groups regarding the evolution of gene expression following gene duplication. The first group, which comprises around 25% of the families, shows patterns compatible with tissue-expression partitioning. The second and largest group, comprising 33-53% of the families, shows broad expression of one of the gene copies and reduced, overlapping, expression of the other copy or copies. This can be attributed, in most cases, to loss of expression in several tissues of one or more gene copies. Finally, a substantial number of families, 19-35%, maintain a very high level of tissue-expression overlap (>0.8) after tens of millions of years of evolution. These families may have been subject to selection for increased gene dosage.

Unlabelled: Many DNA functional motifs tend to accumulate or cluster at specific gene locations. These locations can be detected, in a group of gene sequences, as high frequency 'peaks' with respect to a reference position, such as the transcription start site (TSS). We have developed a web tool for the identification of regions containing significant motif peaks. We show, by using different yeast gene datasets, that peak regions are strongly enriched in experimentally-validated motifs and contain potentially important novel motifs., Availability: http://genomics.imim.es/peaks

Background: Understanding the constraints that operate in mammalian gene promoter sequences is of key importance to understand the evolution of gene regulatory networks. The level of promoter conservation varies greatly across orthologous genes, denoting differences in the strength of the evolutionary constraints. Here we test the hypothesis that the number of tissues in which a gene is expressed is related in a significant manner to the extent of promoter sequence conservation., Results: We show that mammalian housekeeping genes, expressed in all or nearly all tissues, show significantly lower promoter sequence conservation, especially upstream of position -500 with respect to the transcription start site, than genes expressed in a subset of tissues. In addition, we evaluate the effect of gene function, CpG island content and protein evolutionary rate on promoter sequence conservation. Finally, we identify a subset of transcription factors that bind to motifs that are specifically over-represented in housekeeping gene promoters., Conclusion: This is the first report that shows that the promoters of housekeeping genes show reduced sequence conservation with respect to genes expressed in a more tissue-restricted manner. This is likely to be related to simpler gene expression, requiring a smaller number of functional cis-regulatory motifs.

Information about the genomic coordinates and the sequence of experimentally identified transcription factor binding sites is found scattered under a variety of diverse formats. The availability of standard collections of such high-quality data is important to design, evaluate and improve novel computational approaches to identify binding motifs on promoter sequences from related genes. ABS (http://genome.imim.es/datasets/abs2005/index.html) is a public database of known binding sites identified in promoters of orthologous vertebrate genes that have been manually curated from bibliography. We have annotated 650 experimental binding sites from 68 transcription factors and 100 orthologous target genes in human, mouse, rat or chicken genome sequences. Computational predictions and promoter alignment information are also provided for each entry. A simple and easy-to-use web interface facilitates data retrieval allowing different views of the information. In addition, the release 1.0 of ABS includes a customizable generator of artificial datasets based on the known sites contained in the collection and an evaluation tool to aid during the training and the assessment of motif-finding programs.

In this paper we present several web-based tools to identify conserved patterns in sequences. In particular we present details on the functionality of PROMO version 2.0, a program for the prediction of transcription factor binding site in a single sequence or in a group of related sequences and, of MALGEN, a tool to visualize sequence correspondences among long DNA sequences. The web tools and associated documentation can be accessed at http://www.lsi.upc.es/~alggen (RESEARCH link).

We have developed a set of tools to construct positional weight matrices from known transcription factor binding sites in a species or taxon-specific manner, and to search for matches in DNA sequences.