EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

Comparative gene finding in chicken indicates that we are closing in on the set of multi-exonic widely-expressed human genes

The recent availability of the chicken genome sequence poses the question of whether there are human protein-coding genes conserved in chicken that are currently not included in the human gene catalog. Here, we show, using comparative gene finding followed by experimental verification of exon pairs by RT–PCR, that the addition to the multi-exonic subset of this catalog could be as little as 0.2%, suggesting that we may be closing in on the human gene set. Our protocol, however, has two shortcomings: (i) the bioinformatic screening of the predicted genes, applied to filter out false positives, cannot handle intronless genes; and (ii) the experimental verification could fail to identify expression at a specific developmental time. This highlights the importance of developing methods that could provide a reliable estimate of the number of these two types of genes.

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.

Gene finding in the chicken genome

Background: Despite the continuous production of genome sequence for a number of organisms,reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularlytrue for genomes for which there is not a large collection of known gene sequences, such as therecently published chicken genome. We used the chicken sequence to test comparative andhomology-based gene-finding methods followed by experimental validation as an effective genomeannotation method.Results: We performed experimental evaluation by RT-PCR of three different computational genefinders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram wascomputed and each component of it was evaluated. The results showed that de novo comparativemethods can identify up to about 700 chicken genes with no previous evidence of expression, andcan correctly extend about 40% of homology-based predictions at the 5' end.Conclusions: De novo comparative gene prediction followed by experimental verification iseffective at enhancing the annotation of the newly sequenced genomes provided by standardhomology-based methods.

Genome-wide screen of genes required for caffeine tolerance in fission yeast

Background: An excess of caffeine is cytotoxic to all eukaryotic cell types. We aim to study how cells become tolerant to atoxic dose of this drug, and the relationship between caffeine and oxidative stress pathways.Methodology/Principal Findings: We searched for Schizosaccharomyces pombe mutants with inhibited growth on caffeinecontainingplates. We screened a collection of 2,700 haploid mutant cells, of which 98 were sensitive to caffeine. The genes mutated in these sensitive clones were involved in a number of cellular roles including the H2O2-induced Pap1 and Sty1 stress pathways, the integrity and calcineurin pathways, cell morphology and chromatin remodeling. We have investigated the role of the oxidative stress pathways in sensing and promoting survival to caffeine. The Pap1 and the Sty1 pathways are both required for normal tolerance to caffeine, but only the Sty1 pathway is activated by the drug. Cells lacking Pap1 aresensitive to caffeine due to the decreased expression of the efflux pump Hba2. Indeed, ?hba2 cells are sensitive to caffeine, and constitutive activation of the Pap1 pathway enhances resistance to caffeine in an Hba2-dependent manner. Conclusions/Significance: With our caffeine-sensitive, genome-wide screen of an S. pombe deletion collection, we havedemonstrated the importance of some oxidative stress pathway components on wild-type tolerance to the drug.

Topological comparison of methods for predicting transcriptional cooperativity in yeast

Background: The cooperative interaction between transcription factors has a decisive role in the control of the fate of the eukaryotic cell. Computational approaches for characterizing cooperative transcription factors in yeast, however, are based on different rationales and provide a low overlap between their results. Because the wealth of information contained in protein interaction networks and regulatory networks has proven highly effective in elucidating functional relationships between proteins, we compared different sets of cooperative transcription factor pairs (predicted by four different computational methods) within the frame of those networks. Results: Our results show that the overlap between the sets of cooperative transcription factors predicted by the different methods is low yet significant. Cooperative transcription factors predicted by all methods are closer and more clustered in the protein interaction network than expected by chance. On the other hand, members of a cooperative transcription factor pair neither seemed to regulate each other nor shared similar regulatory inputs, although they do regulate similar groups of target genes. Conclusion: Despite the different definitions of transcriptional cooperativity and the different computational approaches used to characterize cooperativity between transcription factors, the analysis of their roles in the framework of the protein interaction network and the regulatory network indicates a common denominator for the predictions under study. The knowledge of the shared topological properties of cooperative transcription factor pairs in both networks can be useful not only for designing better prediction methods but also for better understanding the complexities of transcriptional control in eukaryotes.

Draft genome of the red harvester ant Pogonomyrmex barbatus.

We report the draft genome sequence of the red harvester ant, Pogonomyrmex barbatus. The genome was sequenced using 454 pyrosequencing, and the current assembly and annotation were completed in less than 1 y. Analyses of conserved gene groups (more than 1,200 manually annotated genes to date) suggest a high-quality assembly and annotation comparable to recently sequenced insect genomes using Sanger sequencing. The red harvester ant is a model for studying reproductive division of labor, phenotypic plasticity, and sociogenomics. Although the genome of P. barbatus is similar to other sequenced hymenopterans (Apis mellifera and Nasonia vitripennis) in GC content and compositional organization, and possesses a complete CpG methylation toolkit, its predicted genomic CpG content differs markedly from the other hymenopterans. Gene networks involved in generating key differences between the queen and worker castes (e.g., wings and ovaries) show signatures of increased methylation and suggest that ants and bees may have independently co-opted the same gene regulatory mechanisms for reproductive division of labor. Gene family expansions (e.g., 344 functional odorant receptors) and pseudogene accumulation in chemoreception and P450 genes compared with A. mellifera and N. vitripennis are consistent with major life-history changes during the adaptive radiation of Pogonomyrmex spp., perhaps in parallel with the development of the North American deserts.

Effects of epigenetic regulatory elements on telomeric gene expression

Transgene expression in eukaryotic cells strongly depends on the locus of integration in the host genome and often results in limited transcription level because of unfavorable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which are believed to act on the chromatin structure and may protect transgenes from this so-called position effect. Despite being extensively used to increase transgene expression, the mechanism of action of many of these elements remains largely unknown. Here we evaluated different epigenetic regulatory DNA elements for their ability to protect transgene transcription at telomeres, a defined chromatin environment associated to low or inconsistent expression caused by the Telomere Position Effect (TPE). For the assessment of the effects of epigenetic regulators at telomeres, a novel dual reporter system had to be designed. Telomeric integration of the newly-developed dual reporter system carrying different epigenetic regulators showed that MARs (Matrix Attachment Regions), a UCOE (Ubiquitous Chromatin-Opening Element) or the chicken cHS4 insulator have strong barrier activity which prevented TPE from spreading toward the centromere, resulting in stable and in some cases increased expression of a telomeric-distal reporter gene. In addition, MARs and STAR element 40 resulted in an increase of cells expressing the telomeric-proximal reporter gene, suggesting also an anti-silencing effect. Chromatin immunoprecipitation assays revealed that at telomeres MARs actively promote the deposition of euchromatic histone marks, especially acetylation of both histone H3 and H4, which might be involved in MARs' barrier and transcriptional activator activities. Differently, the chromatin in proximity of the UCOE element was depleted of several repressive chromatin marks, such as trimethylated lysine 9 and lysine 27 on histone H3 and trimethylated lysine 20 of histone H4, possibly favoring the preservation of an open chromatin structure at the integration site. We conclude that epigenetic regulatory elements that may be used to enhance and sustain transgene expression have all a specific epigenetic signature which might be at the basis of their mechanism of action, and that a combination of different classes of epigenetic regulators might be advantageous when high levels of protein expression are required. - L'expression des transgènes dans les cellules eucaryotes est fortement influencée par leur site d'intégration dans le génome. En effet, une structure chromatinienne défavorable au niveau du site d'intégration peut fortement limiter le degré d'expression d'un transgène. Il existe toutefois des séquences d'ADN qui, en agissant sur la structure de la chromatine, permettent de limiter cet effet de position et, par conséquent, de promouvoir l'expression soutenue d'un transgène. Ces éléments génomiques, connus comme régulateurs épigénétiques, sont largement utilisés dans plusieurs domaines où une expression élevée et soutenue est requise, malgré un mode de fonctionnement parfois méconnu. Dans cette étude, j'ai évalué la capacité de différents régulateurs épigénétiques à protéger la transcription de transgènes au niveau des télomères, régions chromatiniennes bien définies qui ont été associées à un fort effet de silençage, connu comme «effet de position télomérique». Pour cela, un nouveau système à deux gènes rapporteurs a été développé. Lorsque des MARs (Matrix Attachment Regions, séquences d'ADN pouvant s'associer à la matrice nucléaire), un UCOE (Ubiquitous Chromatin-Opening Element, élément pouvant ouvrir la chromatine) ou l'isolateur génétique cHS4 (dérivé du locus de la β-globine de poulet) sont placés entre les deux gènes rapporteurs, une forte activité barrière bloquant la propagation de la chromatine répressive télomérique est observée, résultant en un plus grand nombre de cellules exprimant le gène télomérique-distal. D'autre part, une augmentation du nombre de cellules exprimant le gène télomérique-proximal, observée en présence des éléments MAR et STAR 40 (Stabilizing Anti-Repressor element 40, un élément pouvant prévenir la répression génique), suggère aussi un faible effet anti-silençage pour ces éléments. Des expériences d'immunoprécipitation de la chromatine démontrent qu'au télomère, les MARs favorisent l'assemblage de marqueurs de la chromatine active, surtout l'acétylation des histones H3 et H4, qui pourraient être à la base de l'activité barrière et de celle d'activateur transcriptionel. Différemment, la chromatine à proximité de l'élément UCOE est particulièrement pauvre en marqueurs de la chromatine silencieuse, comme la trimethylation des lysines 9 et 27 de l'histone H3, ainsi que la trimethylation de la lysine 20 de l'histone H4. Cela suggère que UCOE pourrait préserver une structure chromatinienne ouverte au site d'intégration, favorisant l'expression des gènes à sa proximité. En conclusion, les régulateurs épigénétiques analysés lors de cette étude ont tous montré une signature épigénétique spécifique qui pourrait être à la base de leurs mécanismes de fonctionnement, suggérant aussi qu'une utilisation d'éléments épigénétiques de classe différente dans un même vecteur d'expression pourrait être avantageuse lorsque de hauts et soutenus niveaux d'expression sont nécessaires.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin.

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

The SM protein Sly1 accelerates assembly of the ER-Golgi SNARE complex.

Soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) and Sec1/Munc18 (SM) proteins constitute the core of an ancient vesicle fusion machine that diversified into distinct sets that now function in different trafficking steps in eukaryotic cells. Deciphering their precise mode of action has proved challenging. SM proteins are thought to act primarily through one type of SNARE protein, the syntaxins. Despite high structural similarity, however, contrasting binding modes have been found for different SM proteins and syntaxins. Whereas the secretory SM protein Munc18 binds to the ‟closed conformation" of syntaxin 1, the ER-Golgi SM protein Sly1 interacts only with the N-peptide of Sed5. Recent findings, however, indicate that SM proteins might interact simultaneously with both syntaxin regions. In search for a common mechanism, we now reinvestigated the Sly1/Sed5 interaction. We found that individual Sed5 adopts a tight closed conformation. Sly1 binds to both the closed conformation and the N-peptide of Sed5, suggesting that this is the original binding mode of SM proteins and syntaxins. In contrast to Munc18, however, Sly1 facilitates SNARE complex formation by loosening the closed conformation of Sed5.

Analysis of LEDGF/p75 expression regulation

Background To replicate, retroviruses must insert DNA copies of their RNA genomes into the host genome. This integration process is catalyzed by the viral integrase protein. The site of viral integration has been shown to be non-random and retrovirus-specific. LEDGF/p75, a splice variant encoded by PSIP1 gene and described as a general transcription coactivator, was identified as a tethering factor binding both to chromatin and to lentiviral integrases, thereby affecting integration efficiency as well as integration site selection. LEDGF/p75 is still a poorly characterized protein, and its cellular endogenous function has yet to be fully determined. In order to start unveiling the roles of LEDGF/p75 in the cell, we started to investigate the mechanisms involved in the regulation of LEDGF/p75. Materials and methods To identify PSIP1 minimal promoter and associated regulatory elements, we cloned a region starting 5 kb upstream the transcription start site (TSS, +1 reference position) to the ATG start codon (+816), as well as systematic truncations, in a plasmid containing the firefly luciferase reporter gene. These constructs were co-transfected into HEK293 cells with a plasmid encoding the Renilla luciferase under the pTK promoter as an internal control for transfection efficiency. Both luciferase activities were assessed by luminescence as an indicator of promoter activity. Results Luciferase assays identified regions -76 to +1 and +1 to +94 as two independent minimal promoters showing respectively a 3.7x and 2.3x increase in luciferase activity. These two independent minimal promoters worked synergistically increasing luciferase activity up to 16.3x as compared to background. Moreover, we identified five regulatory blocks which modulated luciferase activity depending on the DNA region tested, three enhancers (- 2007 to -1159, -284 to -171 and +94 to +644) and two silencers (-171 to -76 and +796 to +816). However, the silencing effect of the region -171 to -76 is dependent on the presence of the +94 to +644 region, ruling out the enhancer activity of the latter. Computational analysis of PSIP1 promoter revealed the absence of TATA box and initiator (INR) sequences, classifying this promoter as nonconventional. TATA-less and INR-less promoters are characterized by multiple Sp1 binding sites, involved in the recruitment of the RNA pol II complex. Consistent with this, PSIP1 promoter contains multiple putative Sp1 binding sequences in regions -76 to +1 and +1 to +94.

Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

Analysis of susceptibility to human immunodeficiency virus 1 (HIV-1) by comparative genetics

Summary [résumé français voir ci-dessous] From the beginning of the 20th century the world population has been confronted with the human immune deficiency virus 1 (HIV-1). This virus has the particularity to mutate fast, and could thus evade and adapt to the human host. Our closest evolutionary related organisms, the non-human primates, are less susceptible to HIV-1. In a broader sense, primates are differentially susceptible to various retrovirus. Species specificity may be due to genetic differences among primates. In the present study we applied evolutionary and comparative genetic techniques to characterize the evolutionary pattern of host cellular determinants of HIV-1 pathogenesis. The study of the evolution of genes coding for proteins participating to the restriction or pathogenesis of HIV-1 may help understanding the genetic basis of modern human susceptibility to infection. To perform comparative genetics analysis, we constituted a collection of primate DNA and RNA to allow generation of de novo sequence of gene orthologs. More recently, release to the public domain of two new primate complete genomes (bornean orang-utan and common marmoset) in addition of the three previously available genomes (human, chimpanzee and Rhesus monkey) help scaling up the evolutionary and comparative genome analysis. Sequence analysis used phylogenetic and statistical methods for detecting molecular adaptation. We identified different selective pressures acting on host proteins involved in HIV-1 pathogenesis. Proteins with HIV-1 restriction properties in non-human primates were under strong positive selection, in particular in regions of interaction with viral proteins. These regions carried key residues for the antiviral activity. Proteins of the innate immunity presented an evolutionary pattern of conservation (purifying selection) but with signals of relaxed constrain if we compared them to the average profile of purifying selection of the primate genomes. Large scale analysis resulted in patterns of evolutionary pressures according to molecular function, biological process and cellular distribution. The data generated by various analyses served to guide the ancestral reconstruction of TRIM5a a potent antiviral host factor. The resurrected TRIM5a from the common ancestor of Old world monkeys was effective against HIV-1 and the recent resurrected hominoid variants were more effective against other retrovirus. Thus, as the result of trade-offs in the ability to restrict different retrovirus, human might have been exposed to HIV-1 at a time when TRIM5a lacked the appropriate specific restriction activity. The application of evolutionary and comparative genetic tools should be considered for the systematical assessment of host proteins relevant in viral pathogenesis, and to guide biological and functional studies. Résumé La population mondiale est confrontée depuis le début du vingtième siècle au virus de l'immunodéficience humaine 1 (VIH-1). Ce virus a un taux de mutation particulièrement élevé, il peut donc s'évader et s'adapter très efficacement à son hôte. Les organismes évolutivement le plus proches de l'homme les primates nonhumains sont moins susceptibles au VIH-1. De façon générale, les primates répondent différemment aux rétrovirus. Cette spécificité entre espèces doit résider dans les différences génétiques entre primates. Dans cette étude nous avons appliqué des techniques d'évolution et de génétique comparative pour caractériser le modèle évolutif des déterminants cellulaires impliqués dans la pathogenèse du VIH- 1. L'étude de l'évolution des gènes, codant pour des protéines impliquées dans la restriction ou la pathogenèse du VIH-1, aidera à la compréhension des bases génétiques ayant récemment rendu l'homme susceptible. Pour les analyses de génétique comparative, nous avons constitué une collection d'ADN et d'ARN de primates dans le but d'obtenir des nouvelles séquences de gènes orthologues. Récemment deux nouveaux génomes complets ont été publiés (l'orang-outan du Bornéo et Marmoset commun) en plus des trois génomes déjà disponibles (humain, chimpanzé, macaque rhésus). Ceci a permis d'améliorer considérablement l'étendue de l'analyse. Pour détecter l'adaptation moléculaire nous avons analysé les séquences à l'aide de méthodes phylogénétiques et statistiques. Nous avons identifié différentes pressions de sélection agissant sur les protéines impliquées dans la pathogenèse du VIH-1. Des protéines avec des propriétés de restriction du VIH-1 dans les primates non-humains présentent un taux particulièrement haut de remplacement d'acides aminés (sélection positive). En particulier dans les régions d'interaction avec les protéines virales. Ces régions incluent des acides aminés clé pour l'activité de restriction. Les protéines appartenant à l'immunité inné présentent un modèle d'évolution de conservation (sélection purifiante) mais avec des traces de "relaxation" comparé au profil général de sélection purifiante du génome des primates. Une analyse à grande échelle a permis de classifier les modèles de pression évolutive selon leur fonction moléculaire, processus biologique et distribution cellulaire. Les données générées par les différentes analyses ont permis la reconstruction ancestrale de TRIM5a, un puissant facteur antiretroviral. Le TRIM5a ressuscité, correspondant à l'ancêtre commun entre les grands singes et les groupe des catarrhiniens, est efficace contre le VIH-1 moderne. Les TRIM5a ressuscités plus récents, correspondant aux ancêtres des grands singes, sont plus efficaces contre d'autres rétrovirus. Ainsi, trouver un compromis dans la capacité de restreindre différents rétrovirus, l'homme aurait été exposé au VIH-1 à une période où TRIM5a manquait d'activité de restriction spécifique contre celui-ci. L'application de techniques d'évolution et de génétique comparative devraient être considérées pour l'évaluation systématique de protéines impliquées dans la pathogenèse virale, ainsi que pour guider des études biologiques et fonctionnelles

Role of myosin V in actin cable organization in fission yeast

Functional specialization is tightly linked to the ability of eukaryotic cells to acquire a particular shape. Cell morphogenesis, in turn, relies on the capacity to establish and maintain cell "polarity", which is achieved by orienting the trafficking of signaling molecules and organelles towards specific cellular locations and/or membrane domains. The "oriented" transport is based upon cytoskeletal polymers, microtubules and actin filaments, which serve as tracks for molecular motors. These latter generate motion that is translated either into pulling forces or directed transport. Fission yeast, a rod-like unicellular eukaryote, shapes itself by restricting growth at cell tips through the concerted activity of microtubules and actin cables. Microtubules, which assemble into 2-6 bundles and run parallel to the long axis of the cell, serve to orient growth to the tips. Growth is supported by the actin cytoskeleton, which provides tracks, the cables, for motor-based transport of secretory vesicles. The molecular motors, which bind cargos and deliver them to the tips along cables, are also known as type V myosins (hereafter indicated as myosin V). How the bundles of parallel actin filaments, i.e. the cables, extend from the tips through the cell and whether they serve any other purpose, besides providing tracks, is poorly understood. It is also unclear how the crosstalk between the two cytoskeletal systems is achieved. These are the basic questions I addressed during my PhD. The first part of the thesis work (Chapter two) suggests that the sole function of actin cables in polarized growth is to serve as tracks for motors. The data indicate that cells may have evolved two cytoskeletal systems to provide robustness to the polarization process but in principle a unique cytoskeleton might have been able to direct and support polarized growth. How actin cables are organized within the cell to optimize cargo transport is addressed later on (Chapter three). The major finding, based on the actin cable defect of cells lacking myosin Vs, is that actin filaments self-organize through the activity of the transport motors. In fact, by delivering cargos to cell tips and exerting physical pulling forces on actin filaments, Myosin Vs contribute not only to polarize cargo transport but also actin tracks. Among the cargos transported by Myosin V, which may be relevant to its function in organizing cables, there is likely the endoplasmic reticulum (ER). Actin cables, which run parallel to cortical ER, may serve as tracks for Myosin V. Myosin V-driven displacement, in turn, may account for the dynamic expansion and organization of ER during polarized growth as suggested in Chapter four. The last part of the work (Chapter five) highlights the existence of a crosstalk between actin and microtubules. In absence of myosin V, indeed, microtubules contribute to actin cable organization, likely playing a scaffolding/tethering function. Whether or not the kinesin 1, Klp3, plays any role in such process has to be demonstrated. In conclusion the work proposes a novel role for myosin Vs in actin organization, besides its transport function, and provides molecular tools to further dissect the role of this type of myosin in fission yeast. - La spécialisation fonctionnelle est étroitement connectée à la capacité des cellules eucaryotes d'acquérir une forme particulière. La morphogenèse cellulaire à son tour, est basée sur la capacité d'établir et de maintenir la polarité cellulaire, polarité réalisée en orientant le trafic des molécules signales et des organelles vers des zones cellulaires spécifiques. Ce transport directionnel dépend des polymères du cytosquelette, microtubules et microfilaments, qui servent comme des voies pour les moteurs moléculaires. Ces derniers engendrent du mouvement, traduit soit en force de traction soit en transport directionnel. La levure fissipare, un eucaryote unicellulaire en forme de bâtonnet, acquière sa forme en limitant sa croissance aux extrémités par l'action concertée des microtubules et de l'actine. Les microtubules, qui s'assemblent de façon antiparallèle et parcourent la cellule parallèlement à l'axe longitudinal, servent à orienter la croissance aux extrémités. Cette croissance est permise par le cytosquelette d'actine, fournissant des voies, les câbles, pour le transport actif des vésicules de sécrétion. Les moteurs moléculaires, responsables de ce transport actif sont aussi appelés myosines de type V (par la suite appelés myosines V). La manière dont ces câbles s'étendent depuis l'extrémité jusqu'à l'intérieur de la cellule est peu connue. De plus, on ignore également si ces câbles présentent une fonction autre que le transport. L'interaction entre les deux cytosquelettes est également obscure. Ce sont ces questions de base auxquelles j'ai tenté de répondre lors de ma thèse. La première partie de cette thèse (chapitre II) suggère que les câbles d'actine, pendant la croissance polarisée, fonctionnent uniquement comme des voies pour les moteurs moléculaires. Les données indiqueraient que les cellules ont fait évoluer deux systèmes de cytosquelette pour assurer plus de robustesse au processus de polarisation, bien que, comme nous le verrons, un système unique est suffisant. Au chapitre III, nous verrons comment les câbles d'actine sont organisés à l'intérieur de la cellule afin d'optimiser le transport des cargo. La découverte majeure, réalisée en observant des cellules dont la myosine V fait défaut, est que ces filaments d'actine s'auto organisent grâce au passage des moteurs moléculaires le long de ces voies. En réalité, en délivrant les cargos aux extrémités de la cellule et en exerçant des forces de traction sur les câbles, les myosines V contribuent non seulement à polariser le transport mais également à polariser les voies elles mêmes. Nous verrons également au chapitre IV, que parmi les cargos importants pour l'organisation des câbles, il y aurait le réticulum endoplasmique (RE). En effet, les câbles d'actine, qui s'étalent parallèlement au RE cortical, pourraient servir comme voie pour la myosine V. Cette dernière en retour pourrait être responsable de l'expansion dynamique et de l'organisation du RE pendant la croissance polarisée.