vHOG, a multispecies vertebrate ontology of homologous organs groups.

MOTIVATION: Most anatomical ontologies are species-specific, whereas a framework for comparative studies is needed. We describe the vertebrate Homologous Organs Groups ontology, vHOG, used to compare expression patterns between species.¦RESULTS: vHOG is a multispecies anatomical ontology for the vertebrate lineage. It is based on the HOGs used in the Bgee database of gene expression evolution. vHOG version 1.4 includes 1184 terms, follows OBO principles and is based on the Common Anatomy Reference Ontology (CARO). vHOG only describes structures with historical homology relations between model vertebrate species. The mapping to species-specific anatomical ontologies is provided as a separate file, so that no homology hypothesis is stated within the ontology itself. Each mapping has been manually reviewed, and we provide support codes and references when available. Availability and implementation: vHOG is available from the Bgee download site (http://bgee.unil.ch/), as well as from the OBO Foundry and the NCBO Bioportal websites.¦CONTACT: bgee@isb-sib.ch; frederic.bastian@unil.ch.

Lymphodepletion by short-term chemotherapy does not alter the highly stable and persistent EBV specific CD8 T cell repertoire

Protective T cell responses againstpersistent viruses like Epstein-Barrvirus in healthy individuals are characterizedby a remarkable stability ofthe T cell receptor (TCR) clonotypicrepertoire, with highly preservedclonotype selection and persistenceover time. Here, we extended recentwork to the study of EBV-specificCD8 T cell responses in melanomapatients treated by short-term chemotherapyfor transient lymphodepletion,followed by adoptive cell transfer(ACT) and immune reconstitutionfor cancer therapy. After this treatment,we observed increased proportionsof virus-specific T cells in 3/5patients, accompanied by a more differentiatedphenotype (EMRA/CD28neg), compared to specific cells ofhealthy individuals. Yet, similarly tohealthy donors, clonotype selectionand composition of virus-specific Tcells varied along the pathway of celldifferentiation, with some clonotypesthat were selected with late differentiation,while others were not. Aftertreatment, we did not observe noveldominant clonotypes, likely related toabsence of EBV reactivation measuredas viral load levels by quantitativePCR in PBMCs and antibody levelsin plasma samples. Furthermore,public TCR BV signatures were frequentlyfound within T cell clonotypesthat dominated the repertoiresof patients, in line with those observedin healthy individuals. Ourfindings indicate that even in situationswhere the immune system isstrongly challenged such as followinglymphodepletion and homeostatic repopulation,cytotoxic T cells specificfor EBV remain strikingly stable interms of clonotype selection and com-position along T cell differentiation.We are currently characterizing theclonotype selection and gene expressionprofiles of single EBV-specificCD8 T lymphocytes sorted ex-vivo inone patient who underwent two cyclesof lymphodepletion with escaladingdoses of chemotherapy overone-year interval. Observations madefrom this setting will provide additionalinsight into the degree of stabilityof virus specific T cells, and changesin the expression levels of genesimportant for cytolytic function andlong-term survival of T cells. Thiswork is supported by the Swiss NationalCenter of Competence in Research(NCCR) Molecular Oncology,and the Swiss National Science Foundation.

Substantial deletion overlap among divergent Arabidopsis genomes revealed by intersection of short reads and tiling arrays.

ABSTRACT: Identification of small polymorphisms from next generation sequencing short read data is relatively easy, but detection of larger deletions is less straightforward. Here, we analyzed four divergent Arabidopsis accessions and found that intersection of absent short read coverage with weak tiling array hybridization signal reliably flags deletions. Interestingly, individual deletions were frequently observed in two or more of the accessions examined, suggesting that variation in gene content partly reflects a common history of deletion events.

HAMAP in 2015: updates to the protein family classification and annotation system.

HAMAP (High-quality Automated and Manual Annotation of Proteins-available at http://hamap.expasy.org/) is a system for the automatic classification and annotation of protein sequences. HAMAP provides annotation of the same quality and detail as UniProtKB/Swiss-Prot, using manually curated profiles for protein sequence family classification and expert curated rules for functional annotation of family members. HAMAP data and tools are made available through our website and as part of the UniRule pipeline of UniProt, providing annotation for millions of unreviewed sequences of UniProtKB/TrEMBL. Here we report on the growth of HAMAP and updates to the HAMAP system since our last report in the NAR Database Issue of 2013. We continue to augment HAMAP with new family profiles and annotation rules as new protein families are characterized and annotated in UniProtKB/Swiss-Prot; the latest version of HAMAP (as of 3 September 2014) contains 1983 family classification profiles and 1998 annotation rules (up from 1780 and 1720). We demonstrate how the complex logic of HAMAP rules allows for precise annotation of individual functional variants within large homologous protein families. We also describe improvements to our web-based tool HAMAP-Scan which simplify the classification and annotation of sequences, and the incorporation of an improved sequence-profile search algorithm.

Repercussion of a deficiency in mitochondrial ß-oxidation on the carbon flux of short-chain fatty acids to the peroxisomal ß-oxidation cycle in Aspergillus nidulans.

The fungus Aspergillus nidulans contains both a mitochondrial and peroxisomal ß-oxidation pathway. This work was aimed at studying the influence of mutations in the foxA gene, encoding a peroxisomal multifunctional protein, or in the scdA/echA genes, encoding a mitochondrial short-chain dehydrogenase and an enoyl-CoA hydratase, respectively, on the carbon flux to the peroxisomal ß-oxidation pathway. A. nidulans transformed with a peroxisomal polyhydroxyalkanoate (PHA) synthase produced PHA from the polymerization of 3-hydroxyacyl-CoA intermediates derived from the peroxisomal ß-oxidation of external fatty acids. PHA produced from erucic acid or heptadecanoic acid contained a broad spectrum of monomers, ranging from 5 to 14 carbons, revealing that the peroxisomal ß-oxidation cycle can handle both long and short-chain intermediates. While the ∆foxA mutant grown on erucic acid or oleic acid synthesized 10-fold less PHA compared to wild type, the same mutant grown on octanoic acid or heptanoic acid produced 3- to 6-fold more PHA. Thus, while FoxA has an important contribution to the degradation of long-chain fatty acids, the flux of short-chain fatty acids to peroxisomal ß-oxidation is actually enhanced in its absence. While no change in PHA was observed in the ∆scdA∆echA mutant grown on erucic acid or oleic acid compared to wild type, there was a 2- to 4-fold increased synthesis of PHA in ∆scdA∆echA cells grown in octanoic acid or heptanoic acid. These results reveal that a compensatory mechanism exists in A. nidulans that increases the flux of short-chain fatty acids towards the peroxisomal ß-oxidation cycle when the mitochondrial ß-oxidation pathway is defective.

A comparison of random sequence reads versus 16S rDNA sequences for estimating the biodiversity of a metagenomic library

The construction of metagenomic libraries has permitted the study of microorganisms resistant to isolation and the analysis of 16S rDNA sequences has been used for over two decades to examine bacterial biodiversity. Here, we show that the analysis of random sequence reads (RSRs) instead of 16S is a suitable shortcut to estimate the biodiversity of a bacterial community from metagenomic libraries. We generated 10,010 RSRs from a metagenomic library of microorganisms found in human faecal samples. Then searched them using the program BLASTN against a prokaryotic sequence database to assign a taxon to each RSR. The results were compared with those obtained by screening and analysing the clones containing 16S rDNA sequences in the whole library. We found that the biodiversity observed by RSR analysis is consistent with that obtained by 16S rDNA. We also show that RSRs are suitable to compare the biodiversity between different metagenomic libraries. RSRs can thus provide a good estimate of the biodiversity of a metagenomic library and, as an alternative to 16S, this approach is both faster and cheaper.

Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

A short motif in Drosophila SECIS Binding Protein 2 provides differential binding affinity to SECIS RNA hairpins

Selenoproteins contain the amino acid selenocysteine which is encoded by a UGA Sec codon. Recoding UGA Sec requires a complex mechanism, comprising the cis-acting SECIS RNA hairpin in the 3′UTR of selenoprotein mRNAs, and trans-acting factors. Among these, the SECIS Binding Protein 2 (SBP2) is central to the mechanism. SBP2 has been so far functionally characterized only in rats and humans. In this work, we report the characterization of the Drosophila melanogaster SBP2 (dSBP2). Despite its shorter length, it retained the same selenoprotein synthesis-promoting capabilities as the mammalian counterpart. However, a major difference resides in the SECIS recognition pattern: while human SBP2 (hSBP2) binds the distinct form 1 and 2 SECIS RNAs with similar affinities, dSBP2 exhibits high affinity toward form 2 only. In addition, we report the identification of a K (lysine)-rich domain in all SBP2s, essential for SECIS and 60S ribosomal subunit binding, differing from the well-characterized L7Ae RNA-binding domain. Swapping only five amino acids between dSBP2 and hSBP2 in the K-rich domain conferred reversed SECIS-binding properties to the proteins, thus unveiling an important sequence for form 1 binding.

Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Increasing the carbon flux toward synthesis of short-chain-length--medium-chain-length polyhydroxyalkanoate in the peroxisome of Saccharomyces cerevisiae through modification of the beta-oxidation cycle.

Short-chain-length-medium-chain-length polyhydroxyalkanoates were synthesized in Saccharomyces cerevisiae from intermediates of the beta-oxidation cycle by expressing the polyhydroxyalkanoate synthases from Aeromonas caviae and Ralstonia eutropha in the peroxisomes. The quantity of polymer produced was increased by using a mutant of the beta-oxidation-associated multifunctional enzyme with low dehydrogenase activity toward R-3-hydroxybutyryl coenzyme A.

Use of a short hPDE6b promoter for rod gene transfer in a model of severe retinal cystrophy, the Rd10 mouse

Purpose: Gene therapy of severe retinal dystrophies directly affecting photoreceptor is still a challenge in terms of clinical application. One of the main hurdles is to generate high transgene expression specifically in rods or cones. In the present study, we are investigating the possibility to drive hPDE6b expression in the Rd10 mouse retina using a specific sequence of the human PDE6b promoter. Methods: Two 5' flanking fragments of the human PDE6b gene: (-93 to +53 (146 bp) and -297 to +53 (350 bp, see Di Polo and Farber, 1995) were cloned in different plasmids in order to check their expression in vitro and in vivo. These elements drove the activity of either luciferase (pGL3 plasmids) or EGFP (AAV2/8 backbone). Then, an AAV2/8 vector carrying the PDE6b cDNA was tested with subretinal injections at P9 in the Rd10 eyes. Eye fundus, OCT, ERG recordings and histological investigations were performed to assess the efficacy of the gene transfer. Results: The short PDE6b promoter containing 146bp (-93 to +53) showed the highest activity in the Y-79 cells, as described previously (Di Polo and Farber, 1995). Subretinal administrations of AAV2/8-PDE6bpromoter-EGFP allowed a rapid expression specifically in rods and not in cones. The expression is faster than a vector containing the CMV promoter. The AAV2/8-PDE6bpromoter-PDE6b and the control vector were injected at P9 in the Rd10 mouse retina and investigated 5 weeks post-injection. Out of 14 eyes, 6 presented an increased rod sensitivity of about 300 fold, and increased a- and b-wave responses in ERG recordings. Flicker stimulations revealed that cones are also functional. OCT images and histological analyses revealed an increased ONL size in the injected area. The retina treated with the therapeutic vector presented 4-6 rows of photoreceptors with outersegments containing PDE6b. In the control eyes, only 2-4 rows of photoreceptors with almost no OS were observed . Conclusions: The 146 bp promoter sequence (-93 to + 53) is the shortest regulatory element described to date which allows to obtain efficient rod-specific expression in the context of somatic gene transfer. This first result is of great interest for AAV vector design in general allowing more space for the accommodation of transgenes of interest and good expression in rods. Moreover we showed the proof of principle of the efficacy of AAV2/8-PDE6bp-PDE6b vector in the Rd10 mouse model of severe photoreceptor degeneration without using neither AAV mutated capsids, nor self-complementary vectors.

The mitochondrial targeting function of randomly generated peptide sequences correlates with predicted helical amphiphilicity.

A pool of oligonucleotides encoding a start methionine and nine random amino acids was inserted at the 5'-end of the gene for the yeast cytochrome oxidase subunit IV lacking its own mitochondrial targeting sequence. Approximately one-quarter of the randomly generated sequences targeted subunit IV to its correct intramitochondrial location in vivo. Sequence analysis of 89 randomly generated sequences showed that their efficiencies as mitochondrial targeting signals correlated with the potential to fold into an amphiphilic alpha-helix. Functional targeting sequences were enriched in arginine and isoleucine residues but contained few aspartate, glutamate, and proline residues. Nonfunctional sequences predicted to have significant helical amphiphilicity often had at least one acidic or multiple helix-breaking residues that would be expected to interfere with targeting functioning. These results support the hypothesis that the signal for targeting a protein into the mitochondrial matrix is usually a positively charged amphiphilic helix.

Contribution pollinique à l'étude du pléistocène de la région lémanique (Suisse) : paléoenvironnement des derniers 800'000 ans

Abstact : Pollinic contribution to a Pleistocene study of the Lake Geneva region (Switzerland) Palaeoenvironment of the last 800,000 years. Although the "Plateau Romand" (tableland of French-speaking Switzerland) has been marked by the repetitive action of different Quaternary glacial thrusts, numerous preserved outcrops here and there, thanks to favourable topographical conditions, provide sufficient pollinic content to reconstitute the significant botanical fluctuations and allow bio- and chrono-stratigraphical correlations with the long European sequences. The present study examines around ten Pleistocene deposits rich in pollens, whose presumed ages cover a range of 800,000 years. - The Early and Early-Middle Pleistocene are examined at Ecoteaux (VD), undeniably the oldest site in this study. The "Lower Formation of Ecoteaux", whose low pollinic content reveals an environment of the cold desert type, shows an inverse palaeomagnetic remanence. An age earlier than 780,000 years BP (Matuyama period or earlier) is proposed. The "Upper Formation of Ecoteaux" contains over eight major fluctuations in the vegetal covering, of which at least four temporal phases and one climatic optimum. The presence of botanical genera, relics of the Tertiary in particular Pterocarya and Carya, associated with the vegetal dynamic made up of short oscillations, contrives to link this formation to the Cromerian Complex. The thermomeres recorded could correspond to one or several interglacial periods from Middle Pleistocene. Possible correlations with marine isotopic stages 15, 17 and 19 are discussed. - The Middle and Late Pleistocene is approached through the sediments of Onnens (VD), Creux d'Enfer (FR), Port-Valais (VS), La Dénériaz (VD) and Cortaillod (NE). Each site provides a temperate interglacial flora, whose pollinic content is evaluated according to the different conceivable bio-stratigraphic correlations. - Finally the Würm is studied thanks to the outcrop of Marly (FR), together with the excavations at Villars-sous-Yens (VD), Versoix (GE) and Ollon (VD). Together, these observations tend to show that the west of the Swiss Plateau is a zone of transition between the meridian vegetation of the lower Rhone valley and that of beyond the Jura. At the outlet of the Rhone glacier, and on the borders of the prealpine and Jurassian glaciers, the region is situated at the crossroads of post-glacial botanical migrations, between the major axis of the Rhone valley and the Eastern refuge zones; this situation confers particular importance on it. Although believed to be deeply scraped by ice, one notices today that the "Plateau Romand" harbours overtwenty interglacial and interstadial sites. The abrasive action of the Rhone glacier has thus been more modest than previously estimated. The effective extension of glacial tongues, together with atmospheric circulation, are doubtless varying from one interglacial to another. Finally, the speeds of glacial movements, and the unbelievable vegetal vitality, are important parameters, capable of leading to deposits, very early enriched in pollens. RÉSUMÉS : Contribution pollinique à l'étude du Pléistocène de la région lémanique (Suisse) - Paléoenvironnement des derniers 800'000 ans Bien que le Plateau romand soit marqué par l'action répétitive des différentes poussées glaciaires quaternaires, de nombreux affleurements préservés ça et là, à la faveur de conditions topographiques favorables, livrent un contenu pollinique suffisant pour reconstituer des fluctuations botaniques significatives et autoriser des corrélations bioet chronostratigraphiques avec les longues séquences européennes. Le présent travail étudie une dizaine de dépôts pléistocènes riches en pollens, dont les âges présumés balaient un éventail de plus de 800'000 ans. - Le Pléistocène Ancien et Moyen inférieur, sont examinés à Ecoteaux (VD); indéniablement le site le plus ancien de cette étude. La "Formation inférieure d'Ecoteaux", dont le maigre contenu pollinique traduit un environnement de type désertique froid, présente une rémanence magnétique inverse. Un âge antérieur à 780'000 ans BP (période de Matuyama ou antérieure) est proposé. La "Formation supérieure d'Ecoteaux" comprend plus de huit fluctuations majeures du couvert végétal, parmi lesquelles on compte au moins quatre épisodes tempérés et un optimum climatique. La présence de genres botaniques reliques du Tertiaire, en particulier Pterocarya et Carya, associée à la dynamique végétale faite de courtes oscillations, concourent à lier cette formation au Complexe Cromérien. Les thermomères enregistrées pourraient correspondre à un- ou plusieurs interglaciaires du Pléistocène moyen. Les corrélations possibles avec les stades isotopiques marins, MIS 15, 17 et 19 sont discutées. - Le Pléistocène Moyen et Récent est abordé à travers les sédiments d'Onnens (VD), du Creux d'Enfer (FR), de Port-Valais (VS), de La Dénériaz (VD) et de Cortaillod (NE). Chaque site livre une flore tempérée interglaciaire, dont le contenu pollinique est évalué en fonction des différentes corrélations bio-stratigraphiques envisageables. - Enfin le Würm est étudié grâce à l'affleurement de Marly (FR), accompagné des forages de Villars-sous-Yens (VD), Versoix (GE) et Ollon (VD). L'ensemble de ces observations tend à montrer que l'ouest du Plateau suisse est une zone charnière entre la végétation méridionale de la basse vallée du Rhône et celle d'au-delà du Jura. Au débouché du glacier du Rhône et aux confins des glaciers préalpins et jurassiens, cette région est située au carrefour des migrations botaniques post-glaciaires, entre l'axe majeur de la vallée du Rhône et les zones refuges de l'Est. Cette situation lui confère une importance toute particulière. Alors qu'on pensait le Plateau romand profondément raboté par les glaces, on constate aujourd'hui qu'il recèle plus d'une vingtaine de sites interglaciaires et interstadiaires. L'action abrasive du glacier rhodanien a donc été plus modeste qu'envisagée jusqu'ici. L'extension effective des langues glaciaires, tout comme la circulation atmosphérique, diffèrent sans doute d'un interglaciaire à l'autre. Enfin la vitesse des mouvements glaciaires et l'incroyable vitalité végétale sont des paramètres importants, susceptibles d'avoir entraîné des dépôts, très précocement enrichis en pollens.