SGP-1: prediction and validation of homologous genes based on sequence alignments

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of SGP-1 depends little on species-specific properties such as codon usage or the nucleotide distribution. SGP-1 may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

Pattern formation during mesophase growth in a homologous series

Remarkable differences in the shape of the nematic-smectic-B interface in a quasi-two-dimensional geometry have been experimentally observed in three liquid crystals of very similar molecular structure, i.e., neighboring members of a homologous series. In the thermal equilibrium of the two mesophases a faceted rectanglelike shape was observed with considerably different shape anisotropies for the three homologs. Various morphologies such as dendritic, dendriticlike, and faceted shapes of the rapidly growing smectic-B germ were also observed for the three substances. Experimental results were compared with computer simulations based on the phase field model. The pattern forming behavior of a binary mixture of two homologs was also studied.

Interaccions arbre-arbust al límit superior del bosc i efectes de gradients bioclimàtics en la vegetació

Durant el periode d’elaboració d’aquesta tesi hem aprofundit en el coneixement dels factors que controlen les dinàmiques espacio-temporals del límit superior del bosc. Aquest ecotò se situa entre el límit superior del bosc i els prats alpins i és susceptible a canvis ambientals, fet que provoca que fluctuï altitudinalment i latitudinalment en funció d’aquests canvis. Els motius d’aquesta dinàmica s’ha estudiat sovint des d’un punt de vista climàtic, però mai fins ara s’havia estudiat des d’un punt de vista de les interaccions entre organismes. Per aquest fet hem estat evaluant l’efecte de les interaccions planta-planta en la regulació de la dinàmica supraforestal. L’estudi l’hem emmarcat en un context alpí (als Pirineus Catalans) i en un context subàrtic (Lapònia, Suècia), fet que ens ha permès fer un estudi comparatiu en dos ecotons contrastats però homòlegs ecològicament. Hem desenvolupat una sèrie d’experiments considerant diversos factors (augment de temperatura, quantitat de nutrients, presència d’arbust, posició en l’ecotò); en les dues zones d’estudi hem fet una plantació de plançons dels arbres formadors del límit del bosc en les diverses situacions derivades de la combinació d’aquests factors, i hem fet el seguiment fenològic dels plançons durant tres periodes de creixement. Els resultats dels experiments ens han permès veure que les interaccions entre organismes tenen una gran importància en la regulació de la dinàmica supraforestal, tant als Pirineus com a Lapònia. Les interaccions planta-planta i planta-herbívors determinen el reclutament de plançons i per tant l’estructuració de les comunitats supraforestals. Per altra banda, la posició en l’ecotò evidencia la presència d’un gradient bioclimàtic; les manipulacions ambientals de temperatura i nutrients originen una resposta generalment positiva en el desenvolupament dels plançons, indicant que canvis en aquestes variables pot suposar alteracions notables de l’estructura forestal del límit del bosc. Per altra banda en aquest projecte també hem aprofundit en temes relacionats amb l'efecte dels gradients altitudinals en la distribució de plantes vasculars als Pirineus Catalans.

Interaccions arbre-arbust al límit superior del bosc i efectes de gradients bioclimàtics en la vegetació

Durant el periode d’elaboració d’aquesta tesi hem aprofundit en el coneixement dels factors que controlen les dinàmiques espacio-temporals del límit superior del bosc. Aquest ecotò se situa entre el límit superior del bosc i els prats alpins i és susceptible a canvis ambientals, fet que provoca que fluctuï altitudinalment i latitudinalment en funció d’aquests canvis. Els motius d’aquesta dinàmica s’ha estudiat sovint des d’un punt de vista climàtic, però mai fins ara s’havia estudiat des d’un punt de vista de les interaccions entre organismes. Per aquest fet hem estat evaluant l’efecte de les interaccions planta-planta en la regulació de la dinàmica supraforestal. L’estudi l’hem emmarcat en un context alpí (als Pirineus Catalans) i en un context subàrtic (Lapònia, Suècia), fet que ens ha permès fer un estudi comparatiu en dos ecotons contrastats però homòlegs ecològicament. Hem desenvolupat una sèrie d’experiments considerant diversos factors (augment de temperatura, quantitat de nutrients, presència d’arbust, posició en l’ecotò); en les dues zones d’estudi hem fet una plantació de plançons dels arbres formadors del límit del bosc en les diverses situacions derivades de la combinació d’aquests factors, i hem fet el seguiment fenològic dels plançons durant tres periodes de creixement. Els resultats dels experiments ens han permès veure que les interaccions entre organismes tenen una gran importància en la regulació de la dinàmica supraforestal, tant als Pirineus com a Lapònia. Les interaccions planta-planta i planta-herbívors determinen el reclutament de plançons i per tant l’estructuració de les comunitats supraforestals. Per altra banda, la posició en l’ecotò evidencia la presència d’un gradient bioclimàtic; les manipulacions ambientals de temperatura i nutrients originen una resposta generalment positiva en el desenvolupament dels plançons, indicant que canvis en aquestes variables pot suposar alteracions notables de l’estructura forestal del límit del bosc. Per altra banda en aquest projecte també hem aprofundit en temes relacionats amb l'efecte dels gradients altitudinals en la distribució de plantes vasculars als Pirineus Catalans.

Identificació i caracterització funcional dels complexes proteics dels receptors BRL1 i BRL3 en el teixit vascular d’Arabidopsis Thaliana

Els esteroids juguen papers clau en el creixement I el desenvolupament d’eucariotes multicel•lulars. En plantes, aquestes hormones, anomenades Brassinosteroides (BRs), estan involucrades en una gran varietat de processos biològics essencials per a les plantes. S’han descrit anteriorment dos receptors de BRs del tipus Leucine Rich Repeat Receptor Like Kinase LRR-RLK, BRASSINOSTEROID RECEPTOR LIKE 1 i 3 (BRL1 i BRL3 respectivalemt) que són homòlegs al receptor principal BRI1 i són necessaris pel desenvolupament vascular. Tot i que els principals components de la senyal ja han estat identificats pel seu homòleg més pròxim, el receptor BRI1, els complexes de BRL1 i BRL3 juntament amb els candidats co-receptors així com els components de la ruta de senyalització encara no han sigut identificats. Per tal d’entendre millor la funció molecular d’aquests receptors de BRs en la planta aquesta tesis doctoral planteja dues aproximacions: com a primera aproximació, vaig realitzar un estudi fenotípic del desenvolupament del teixit vascular a la planta model Arabidopsis thaliana (Arabidopsis). Disposant d'una amplia bateria de mutants de síntesis de la hormona i senyalització del receptor BRI1, vam analitzar quantitativament el seu patró vascular a la tija d'Arabidopsis. Vam establir els paràmetres en les plantes silvestres [Col-0 wild type, (WT)] i els vam analitzar a tots i cadascun dels mutants. Això conjuntament amb una col•laboració amb la Dr. Marta Ibañes, física de la Universitat de Barcelona que va construir un model matemàtic per simular la formació del patró vascular ens va permetre el•laborar una hipòtesis que vam demostrar experimentalment i va ser publicada a la revista PNAS. Posteriorment vam observar que les plantes knock-out d'aquests dos receptors BRL1 y BRL3 a diferència de BRI1, no tenien cap fenotip obvi en el teixit vascular de la planta adulta. Així, a continuació, per entendre quina necessitat té la planta de disposar de tres receptors tant altament homòlegs que poden percebre la mateixa hormona, vam utilitzar una aproximació bioquímica en col•laboració amb el Prof. de Vries de la Universitat de Wageningen (Holanda) per tal de purificar els complexes dels receptors in vivo i els seus interactors. Això ens ha permès entendre millor el paper funcional d'aquests receptors en la planta. Els resultats d’aquests experiments estan resumits en un article en preparació que aviat estarà en revisió.

Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

Nematode selenoproteome: the use of the selenocysteine insertion system to decode one codon in an animal genome?

Selenocysteine (Sec) is co-translationally inserted into selenoproteins in response to codon UGA with the help of the selenocysteine insertion sequence (SECIS) element. The number of selenoproteins in animals varies, with humans having 25 and mice having 24 selenoproteins. To date, however, only one selenoprotein, thioredoxin reductase, has been detected in Caenorhabditis elegans, and this enzyme contains only one Sec. Here, we characterize the selenoproteomes of C.elegans and Caenorhabditis briggsae with three independent algorithms, one searching for pairs of homologous nematode SECIS elements, another searching for Cys- or Sec-containing homologs of potential nematode selenoprotein genes and the third identifying Sec-containing homologs of annotated nematode proteins. These methods suggest that thioredoxin reductase is the only Sec-containing protein in the C.elegans and C.briggsae genomes. In contrast, we identified additional selenoproteins in other nematodes. Assuming that Sec insertion mechanisms are conserved between nematodes and other eukaryotes, the data suggest that nematode selenoproteomes were reduced during evolution, and that in an extreme reduction case Sec insertion systems probably decode only a single UGA codon in C.elegans and C.briggsae genomes. In addition, all detected genes had a rare form of SECIS element containing a guanosine in place of a conserved adenosine present in most other SECIS structures, suggesting that in organisms with small selenoproteomes SECIS elements may change rapidly.

Prediction of enzyme function by combining sequence similarity and protein interactions

Background: A number of studies have used protein interaction data alone for protein function prediction. Here, we introduce a computational approach for annotation of enzymes, based on the observation that similar protein sequences are more likely to perform the same function if they share similar interacting partners. Results: The method has been tested against the PSI-BLAST program using a set of 3,890 protein sequences from which interaction data was available. For protein sequences that align with at least 40% sequence identity to a known enzyme, the specificity of our method in predicting the first three EC digits increased from 80% to 90% at 80% coverage when compared to PSI-BLAST. Conclusion: Our method can also be used in proteins for which homologous sequences with known interacting partners can be detected. Thus, our method could increase 10% the specificity of genome-wide enzyme predictions based on sequence matching by PSI-BLAST alone.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

Background: Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms. Results: While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs. Conclusions: From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.

Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system

Chemoreception is a biological process essential for the survival of animals, as it allows the recognition of important volatile cues for the detection of food, egg-laying substrates, mates or predators, among other purposes. Furthermore, its role in pheromone detection may contribute to evolutionary processes such as reproductive isolation and speciation. This key role in several vital biological processes makes chemoreception a particularly interesting system for studying the role of natural selection in molecular adaptation. Two major gene families are involved in the perireceptor events of the chemosensory system: the odorant-binding protein (OBP) and chemosensory protein (CSP) families. Here, we have conducted an exhaustive comparative genomic analysis of these gene families in twenty Arthropoda species. We show that the evolution of the OBP and CSP gene families is highly dynamic, with a high number of gains and losses of genes, pseudogenes and independent origins of subfamilies. Taken together, our data clearly support the birth-and-death model for the evolution of these gene families with an overall high gene-turnover rate. Moreover, we show that the genome organization of the two families is significantly more clustered than expected by chance and, more important, that this pattern appears to be actively maintained across the Drosophila phylogeny. Finally, we suggest the homologous nature of the OBP and CSP gene families, dating back their MRCA (most recent common ancestor) to 380¿420 Mya, and we propose a scenario for the origin and diversification of these families.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Contribution of S6K1/MAPK signaling pathways in the response to oxidative stress: activation of RSK and MSK by hydrogen peroxide

Cells respond to different kind of stress through the coordinated activation of signaling pathways such as MAPK or p53. To find which molecular mechanisms are involved, we need to understand their cell adaptation. The ribosomal protein, S6 kinase 1 (S6K1), is a common downstream target of signaling by hormonal or nutritional stress. Here, we investigated the initial contribution of S6K1/MAPK signaling pathways in the cell response to oxidative stress produced by hydrogen peroxide (H2O2). To analyze S6K1 activation, we used the commercial anti-phospho-Thr389-S6K1 antibody most frequently mentioned in the bibliography. We found that this antibody detected an 80-90 kDa protein that was rapidly phosphorylated in response to H2O2 in several human cells. Unexpectedly, this phosphorylation was insensitive to both mTOR and PI3K inhibitors, and knock-down experiments showed that this protein was not S6K1. RSK and MSK proteins were candidate targets of this phosphorylation. We demonstrated that H2O2 stimulated phosphorylation of RSK and MSK kinases at residues that are homologous to Thr389 in S6K1. This phosphorylation required the activity of either p38 or ERK MAP kinases. Kinase assays showed activation of RSK and MSK by H2O2. Experiments with mouse embryonic fibroblasts from p38 animals" knockout confirmed these observations. Altogether, these findings show that the S6K1 signaling pathway is not activated under these conditions, clarify previous observations probably misinterpreted by non-specific detection of proteins RSK and MSK by the anti-phospho-Thr389-S6K1 antibody, and demonstrate the specific activation of MAPK signaling pathways through ERK/p38/RSK/MSK by H2O2.

PhylomeDB v4: zooming into the plurality of evolutionary histories of a genome

Phylogenetic trees representing the evolutionary relationships of homologous genes are the entry point for many evolutionary analyses. For instance, the use of a phylogenetic tree can aid in the inference of orthology and paralogy relationships, and in the detection of relevant evolutionary events such as gene family expansions and contractions, horizontal gene transfer, recombination or incomplete lineage sorting. Similarly, given the plurality of evolutionary histories among genes encoded in a given genome, there is a need for the combined analysis of genome-wide collections of phylogenetic trees (phylomes). Here, we introduce a new release of PhylomeDB (http://phylomedb.org), a public repository of phylomes. Currently, PhylomeDB hosts 120 public phylomes, comprising >1.5 million maximum likelihood trees and multiple sequence alignments. In the current release, phylogenetic trees are annotated with taxonomic, protein-domain arrangement, functional and evolutionary information. PhylomeDB is also a major source for phylogeny-based predictions of orthology and paralogy, covering >10 million proteins across 1059 sequenced species. Here we describe newly implemented PhylomeDB features, and discuss a benchmark of the orthology predictions provided by the database, the impact of proteome updates and the use of the phylome approach in the analysis of newly sequenced genomes and transcriptomes.

Analysis of meiotic recombination in 22q11.2, a region that frequently undergoes deletions and duplications

Background: The 22q11.2 deletion syndrome is the most frequent genomic disorder with an estimated frequency of 1/4000 live births. The majority of patients (90%) have the same deletion of 3 Mb (Typically Deleted Region, TDR) that results from aberrant recombination at meiosis between region specific low-copy repeats (LCRs). Methods: As a first step towards the characterization of recombination rates and breakpoints within the 22q11.2 region we have constructed a high resolution recombination breakpoint map based on pedigree analysis and a population-based historical recombination map based on LD analysis. Results: Our pedigree map allows the location of recombination breakpoints with a high resolution (potential recombination hotspots), and this approach has led to the identification of 5 breakpoint segments of 50 kb or less (8.6 kb the smallest), that coincide with historical hotspots. It has been suggested that aberrant recombination leading to deletion (and duplication) is caused by low rates of Allelic Homologous Recombination (AHR) within the affected region. However, recombination rate estimates for 22q11.2 region show that neither average recombination rates in the 22q11.2 region or within LCR22-2 (the LCR implicated in most deletions and duplications), are significantly below chromosome 22 averages. Furthermore, LCR22-2, the repeat most frequently implicated in rearrangements, is also the LCR22 with the highest levels of AHR. In addition, we find recombination events in the 22q11.2 region to cluster within families. Within this context, the same chromosome recombines twice in one family; first by AHR and in the next generation by NAHR resulting in an individual affected with the del22q11.2 syndrome. Conclusion: We show in the context of a first high resolution pedigree map of the 22q11.2 region that NAHR within LCR22 leading to duplications and deletions cannot be explained exclusively under a hypothesis of low AHR rates. In addition, we find that AHR recombination events cluster within families. If normal and aberrant recombination are mechanistically related, the fact that LCR22s undergo frequent AHR and that we find familial differences in recombination rates within the 22q11.2 region would have obvious health-related implications.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

Background: Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms. Results: While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs. Conclusions: From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.