Bioinformatic study of selection in animal genomes

Les gènes orthologues divergent sur plusieurs aspects durant l'évolution. Après une revue de la littérature cherchant à montrer de la divergence entre les orthologues de l'humain et de la souris, j'ai souligné les différentes causes de cette divergence. En comparant les gènes qui divergent en fonction, je n'ai pas trouvé de lien avec la divergence des séquences, pour cette raison je me suis penché sur l'étude de l'expression. Notamment, j'ai étudié le niveau, la spécificité ainsi que la présence/absence d'expression des orthologues humain-souris liés aux maladies Mendéliennes. Malgré les similarités trouvées entre l'humain et la souris, j'ai détecté une différence d'expression spécifique à une des deux espèces liée a un phénotype précis (gène essentiel/non-essentiel). Cela m'a permis de conclure que la différence sur le plan phénotypique entre l'humain et la souris est mieux expliquée par les patrons d'expression plutôt que le niveau d'expression ou la sélection. J'ai été également intéressé par l'évolution des séquences d'ADN codantes pour des protéines, en particulier sur le rôle de la sélection. J'ai commencé par une étude sur la fiabilité de détection de la sélection positive en comparant des séquences divergentes. J'ai trouvé, en utilisant le model de branche-site que la sélection peut être détectée sur des séquences qui ont divergé il y a plus de 500 millions d'années. J'ai analysé le biais de GC entres les séquences sans trouver une influence sur l'estimation de la sélection positive. Finalement, Je crois que ce travail est une première étape dans l'établissement d'un lien entre la sélection et les patrons d'expression des gènes chez les vertébrés.

Low-frequency drug-resistant HIV-1 and risk of virological failure to first-line NNRTI-based ART: a multicohort European case-control study using centralized ultrasensitive 454 pyrosequencing.

OBJECTIVES: It is still debated if pre-existing minority drug-resistant HIV-1 variants (MVs) affect the virological outcomes of first-line NNRTI-containing ART. METHODS: This Europe-wide case-control study included ART-naive subjects infected with drug-susceptible HIV-1 as revealed by population sequencing, who achieved virological suppression on first-line ART including one NNRTI. Cases experienced virological failure and controls were subjects from the same cohort whose viraemia remained suppressed at a matched time since initiation of ART. Blinded, centralized 454 pyrosequencing with parallel bioinformatic analysis in two laboratories was used to identify MVs in the 1%-25% frequency range. ORs of virological failure according to MV detection were estimated by logistic regression. RESULTS: Two hundred and sixty samples (76 cases and 184 controls), mostly subtype B (73.5%), were used for the analysis. Identical MVs were detected in the two laboratories. 31.6% of cases and 16.8% of controls harboured pre-existing MVs. Detection of at least one MV versus no MVs was associated with an increased risk of virological failure (OR = 2.75, 95% CI = 1.35-5.60, P = 0.005); similar associations were observed for at least one MV versus no NRTI MVs (OR = 2.27, 95% CI = 0.76-6.77, P = 0.140) and at least one MV versus no NNRTI MVs (OR = 2.41, 95% CI = 1.12-5.18, P = 0.024). A dose-effect relationship between virological failure and mutational load was found. CONCLUSIONS: Pre-existing MVs more than double the risk of virological failure to first-line NNRTI-based ART.

Genes associated with retinitis pigmentosa and allied diseases are frequently mutated in the general population.

BACKGROUND: Retinitis pigmentosa and other hereditary retinal degenerations (HRD) are rare genetic diseases leading to progressive blindness. Recessive HRD are caused by mutations in more than 100 different genes. Laws of population genetics predict that, on a purely theoretical ground, such a high number of genes should translate into an extremely elevated frequency of unaffected carriers of mutations. In this study we estimate the proportion of these individuals within the general population, via the analyses of data from whole-genome sequencing. METHODOLOGY/PRINCIPAL FINDINGS: We screened complete and high-quality genome sequences from 46 control individuals from various world populations for HRD mutations, using bioinformatic tools developed in-house. All mutations detected in silico were validated by Sanger sequencing. We identified clear-cut, null recessive HRD mutations in 10 out of the 46 unaffected individuals analyzed (∼22%). CONCLUSIONS/SIGNIFICANCE: Based on our data, approximately one in 4-5 individuals from the general population may be a carrier of null mutations that are responsible for HRD. This would be the highest mutation carrier frequency so far measured for a class of Mendelian disorders, especially considering that missenses and other forms of pathogenic changes were not included in our assessment. Among other things, our results indicate that the risk for a consanguineous couple of generating a child with a blinding disease is particularly high, compared to other genetic conditions.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

Segmental copy number variation shapes tissue transcriptomes

Abstract : Copy number variation (CNV) of DNA segments has recently gained considerable interest as a source of genetic variation likely to play a role in phenotypic diversity and evolution. Much effort has been put into the identification and mapping of regions that vary in copy number among seemingly normal individuals, both in humans and in a number of model organisms, using both bioinformatic and hybridization-based methods. Synteny studies suggest the existence of CNV hotspots in mammalian genomes, often in connection with regions of segmental duplication. CNV alleles can be in equilibrium within a population, but can also arise de novo between generations, illustrating the highly dynamic nature of these regions. A small number of studies have assessed the effect of CNV on single loci, however, at the genome-wide scale, the functional impact of CNV remains poorly studied. We have explored the influence of CNV on gene expression, first using the Williams-Beuren syndrome (WBS) associated deletion as a model, and second at the genome-wide scale in inbred mouse strains. We found that the WBS deletion influences the expression levels not only of the hemizygous genes, but also affects the euploid genes mapping nearby. Consistently, on a genome wide scale we observe that CNV genes are expressed at more variable levels than genes that do not vary in copy number. Likewise, CNVs influence the relative expression levels of genes that map to the flank of the genome rearrangements, thus globally influencing tissue transcriptomes. Further studies are warranted to complete cataloguing and fine mapping of CNV regions, as well as to elucidate the different mechanisms by which CNVs influence gene expression. Résumé : La variation en nombre de copies (copy number variation ou CNV) de segments d'ADN suscite un intérêt en tant que variation génétique susceptible de jouer un r81e dans la diversité phénotypique et l'évolution. Les régions variables en nombre de copies parmi des individus apparemment normaux ont été cartographiées et cataloguées au moyen de puces à ADN et d'analyse bioinformatique. L'étude de la synténie entre plusieurs espèces de mammifères laisse supposer l'existence de régions à haut taux de variation, souvent liées à des duplications segmentaires. Les allèles CNV peuvent être en équilibre au sein d'une population ou peuvent apparaître de novo. Ces faits illustrent la nature hautement dynamique de ces régions. Quelques études se sont penchées sur l'effet de la variation en nombre de copies de loci isolés, cependant l'impact de ce phénomène n'a pas été étudié à l'échelle génomique. Nous avons examiné l'influence des CNV sur l'expression des gènes. Dans un premier temps nous avons utilisé la délétion associée au syndrome de Williams-Beuren (WBS), puis, dans un second temps, nous avons poursuivi notre étude à l'échelle du génome, dans des lignées consanguines de souris. Nous avons établi que la délétion WBS influence l'expression non seulement des gènes hémizygotes, mais également celle des gènes euploïdes voisins. A l'échelle génomique, nous observons des phénomènes concordants. En effet, l'expression des gènes variant en nombre de copies est plus variable que celles des gènes ne variant pas. De plus, à l'instar de la délétion WBS, les CNV influencent l'expression des gènes adjacents, exerçant ainsi un impact global sur les profils d'expression dans les tissus. Résumé pour un large public : De nombreuses maladies ont pour cause un défaut génétique. Parmi les types de mutations, on compte la disparition (délétion) d'une partie de notre génome ou sa duplication. Bien que l'on connaisse les anomalies associées à certaines maladies, les mécanismes moléculaires par lesquels ces réarrangements de notre matériel génétique induisent les maladies sont encore méconnus. C'est pourquoi nous nous sommes intéressés à la régulation des gènes dans les régions susceptibles à délétion ou duplication. Dans ce travail, nous avons démontré que les délétions et les duplications influencent la régulation des gènes situés à proximité, et que ces changements interviennent dans plusieurs organes.

'Next generation' biogeography: towards understanding the drivers of species diversification and persistence

The drivers of species diversification and persistence are of great interest to current biogeography, especially in those global biodiversity hotspots' harbouring most of Earth's animal and plant life. Classical multispecies biogeographical work has yielded fascinating insights into broad-scale patterns of diversification, and DNA-based intraspecific phylogeographical studies have started to complement this picture at much finer temporal and spatial scales. The advent of novel next-generation sequencing (NGS) technologies provides the opportunity to greatly scale up the numbers of individuals, populations and species sampled, potentially merging intraspecific and interspecific approaches to biogeographical inference. Here, we outline these prospects and issues by using the example of an undisputed hotspot, the Cape of southern Africa. We outline the current state of knowledge on the biogeography of species diversification within the Cape, review the literature for phylogeographical evidence of its likely drivers and mechanisms, and suggest possible ways forward based on NGS approaches. We demonstrate the potential of these methods and current bioinformatic issues with the help of restriction-site-associated DNA (RAD) sequencing data for three highly divergent species of the Restionaceae, an important plant radiation in the Cape. A thorough understanding of the mechanisms that facilitate species diversification and persistence in spatially structured, species-rich environments will require the adoption of novel genomic and bioinformatic tools in biogeographical studies.

Pathogenesis-related proteins and the jasmonate signaling pathway

Résumé Etant une importante source d'énergie, les plantes sont constamment attaquées par des pathogènes. Ne pouvant se mouvoir, elles ont développé des systèmes de défense sophistiqués afin de lutter contre ces prédateurs. Parmi ces systèmes, les voies de signalisation mettant en jeu des éliciteurs endog8nes tels que les jasmonates permettent d'induire la production de protéines de défense telles que les protéines dites "liées à la pathogénèse". Les gènes codant pour ces protéines appartiennent à des familles multigéniques. Le premier but de cette thèse est d'évaluer le nombre de ces gènes dans le génome d'Arabidopsis thaliana et d'estimer la part de ce système de défense, dépendant de la voie de signalisation des jasmonates. Nous avons défini un cluster de seulement 1S gènes sur 266, "liés à la pathogénèse", exclusivement régulés par les jasmonates. De multiples membres des familles des lectines de type jacaline et des inhibiteurs de trypsines semblent dépendre du jasmonate. Présente dans tous les systèmes immunitaires des eucaryotes, la famille des défensines est une famille très intéressante. Chez Arabidopsis thaliana, 317 protéines similaires aux défensines ont été définies, cependant seulement 15 défensines (PDF) sont bien annotées. Ces 15 défensines sont séparées en deux groupes dont un semble avoir évolué plus récemment. Le second but de cette thèse est d'étudier ce groupe de défensines à l'aide de la bioinformatique et des techniques de biologie moléculaire (gêne rapporteur, PCR en temps réel). Nous avons montré que ce groupe contenait une défensine acide intéressante, PDF1.5, qui semblait avoir subi une sélection positive. Cette protéine n'avait encore jamais été étudiée. Contrairement à ce que nous pensions, nous avons établi que cette protéine pouvait avoir une activité biologique liée à la défense. Ce travail de thèse a permis de préciser le nombre de gènes "liées à la pathogénèse" induits par la voie des jasmonates et d'apporter des éléments de réponse sur la question de la redondance des gènes de défense. En conclusion, même si de nombreuses familles de gènes intervenant dans la défense sont bien définies chez Arabidopsis, il reste encore de nombreuses études à faire sur chacun de ces membres. Abstract Being an important source of energy, plants are constantly attacked by herbivores and pathogens. As sessile organisms, they have developed sophisticated defense responses to cope with attack. Among these responses, signalling pathways, using endogenous elicitors including jasmonates (JA), allow the plant to induce the production of defense proteins such as pathogenesis-related (PR) proteins. The genes encoding these proteins belong to multigenic families. The first goal of this thesis was to evaluate the number of PR genes in the genome of Arabidopsis thaliana and estimate how much of this plant defense system was dependent on the jasmonate signaling pathway in leaves. Surprisingly a cluster of only 1S genes out of 2ó6 PR genes was exclusively regulated by JA. Multiple members of the jacalin lectin and trypsin inhibitor gene families were shown to be regulated by JA. Present in all eukaryotic immune systems, defensins are an attractive PR family to study. In Arabidopsis thaliana, 317 defensin-related proteins have been found but just 1S defensins (i.e. PDF family) are well annotated. These defensins are split into 2 groups. One of these groups may have appeared and diversified recently. The second goal of this thesis was to study this defensin gene group combining bioinformatic, reporter gene and quantitative PCR techniques. We have shown that this group contains an interesting acidic defensin, PDF1.S, which seems to have undergone positive selection. No information was known on this protein. We have established that this protein may have a biological activity in plant defense. This thesis allowed us to define the number of PR genes induced by the jasmonate pathway and gave initial leads to explain the redundancy of the PR genes in the genome of Arabidopsis. In conclusion, even if many defense gene families are already defined in the Arabidopsis genome, much work remains to be done on individual members.

Identification of biologic markers of the premature rupture of fetal membranes: proteomic approach.

In obstetrics, premature rupture of the membranes (PROM) is a frequent observation which is responsible for many premature deliveries. PROM is also associated with an increased risk of fetal and maternal infections. Early diagnosis is mandatory in order to decrease such complications. Despite that current biological tests allowing the diagnosis of PROM are both sensitive and specific, contamination of the samples by maternal blood can induce false positive results. Therefore, in order to identify new potential markers of PROM (present only in amniotic blood, and absent in maternal blood), proteomic studies were undertaken on samples collected from six women at terms (pairs of maternal plasma and amniotic fluid) as well as on four samples of amniotic fluid collected from other women at the 17(th) week of gestation. All samples (N = 16) were analyzed by two-dimensional (2-D) high-resolution electrophoresis, followed by sensitive silver staining. The gel images were studied using bioinformatic tools. Analyses were focused on regions corresponding to pI between 4.5 and 7 and to molecular masses between 20 and 50 kDa. In this area, 646 +/- 113 spots were detected, and 27 spots appeared to be present on the gels of amniotic fluid, but were absent on those of maternal plasma. Nine out of these 27 spots were also observed on the gels of the four samples of amniotic fluids collected at the 17(th) week of pregnancy. Five of these 9 spots were unambiguously detected on preparative 2-D gels stained by Coomassie blue, and were identified by mass spectrometry analyses. Three spots corresponded to fragments of plasma proteins, and 2 appeared to be fragments of proteins not known to be present in plasma. These 2 proteins were agrin (SWISS-PROT: O00468) and perlecan (SWISS-PROT: P98160). Our results show that proteomics is a valuable approach to identify new potential biological markers for future PROM diagnosis.

Altered expression of beta-galactosidase-1-like protein 3 (Glb1l3) in the retinal pigment epithelium (RPE)-specific 65-kDa protein knockout mouse model of Leber's congenital amaurosis

Purpose: In this study, we investigated the expression of the gene encoding beta-galactosidase (Glb)-1-like protein 3 (Glb1l3), a member of the glycosyl hydrolase 35 family, during retinal degeneration in the retinal pigment epithelium (RPE)-specific 65-kDa protein knockout (Rpe65(-/-)) mouse model of Leber congenital amaurosis (LCA). Additionally, we assessed the expression of the other members of this protein family, including beta-galactosidase-1 (Glb1), beta-galactosidase-1-like (Glb1l), and beta-galactosidase-1-like protein 2 (Glb1l2).Methods: The structural features of Glb1l3 were assessed using bioinformatic tools. mRNA expression of Glb-related genes was investigated by oligonucleotide microarray, real-time PCR, and reverse transcription (RT) -PCR. The localized expression of Glb1l3 was assessed by combined in situ hybridization and immunohistochemistry.Results: Glb1l3 was the only Glb-related member strongly downregulated in Rpe65(-/-) retinas before the onset and during progression of the disease. Glb1l3 mRNA was only expressed in the retinal layers and the RPE/choroid. The other Glb-related genes were ubiquitously expressed in different ocular tissues, including the cornea and lens. In the healthy retina, expression of Glb1l3 was strongly induced during postnatal retinal development; age-related increased expression persisted during adulthood and aging.Conclusions: These data highlight early-onset downregulation of Glb1l3 in Rpe65-related disease. They further indicate that impaired expression of Glb1l3 is mostly due to the absence of the chromophore 11-cis retinal, suggesting that Rpe65 deficiency may have many metabolic consequences in the underlying neuroretina.

Identification of new MHC class I-restricted human tumor antigens for immunotherapy

Summary Cancer is a leading cause of morbidity and mortality in Western countries (as an example, colorectal cancer accounts for about 300'000 new cases and 200'000 deaths each year in Europe and in the USA). Despite that many patients with cancer have complete macroscopic clearance of their disease after resection, radiotherapy and/or chemotherapy, many of these patients develop fatal recurrence. Vaccination with immunogenic peptide tumor antigens has shown encouraging progresses in the last decade; immunotherapy might therefore constitute a fourth therapeutic option in the future. We dissect here and critically evaluate the numerous steps of reverse immunology, a forecast procedure to identify antigenic peptides from the sequence of a gene of interest. Bioinformatic algorithms were applied to mine sequence databases for tumor-specific transcripts. A quality assessment of publicly available sequence databanks allowed defining strengths and weaknesses of bioinformatics-based prediction of colon cancer-specific alternative splicing: new splice variants could be identified, however cancer-restricted expression could not be significantly predicted. Other sources of target transcripts were quantitatively investigated by polymerase chain reactions, as cancer-testis genes or reported overexpressed transcripts. Based on the relative expression of a defined set of housekeeping genes in colon cancer tissues, we characterized a precise procedure for accurate normalization and determined a threshold for the definition of significant overexpression of genes in cancers versus normal tissues. Further steps of reverse immunology were applied on a splice variant of the Melan¬A gene. Since it is known that the C-termini of antigenic peptides are directly produced by the proteasome, longer precursor and overlapping peptides encoded by the target sequence were synthesized chemically and digested in vitro with purified proteasome. The resulting fragments were identified by mass spectroscopy to detect cleavage sites. Using this information and based on the available anchor motifs for defined HLA class I molecules, putative antigenic peptides could be predicted. Their relative affinity for HLA molecules was confirmed experimentally with functional competitive binding assays and they were used to search patients' peripheral blood lymphocytes for the presence of specific cytolytic T lymphocytes (CTL). CTL clones specific for a splice variant of Melan-A could be isolated; although they recognized peptide-pulsed cells, they failed to lyse melanoma cells in functional assays of antigen recognition. In the conclusion, we discuss advantages and bottlenecks of reverse immunology and compare the technical aspects of this approach with the more classical procedure of direct immunology, a technique introduced by Boon and colleagues more than 10 years ago to successfully clone tumor antigens.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

Reverse immunology approach for the identification of CD8 T-cell-defined antigens: advantages and hurdles.

One of the challenges of tumour immunology remains the identification of strongly immunogenic tumour antigens for vaccination. Reverse immunology, that is, the procedure to predict and identify immunogenic peptides from the sequence of a gene product of interest, has been postulated to be a particularly efficient, high-throughput approach for tumour antigen discovery. Over one decade after this concept was born, we discuss the reverse immunology approach in terms of costs and efficacy: data mining with bioinformatic algorithms, molecular methods to identify tumour-specific transcripts, prediction and determination of proteasomal cleavage sites, peptide-binding prediction to HLA molecules and experimental validation, assessment of the in vitro and in vivo immunogenic potential of selected peptide antigens, isolation of specific cytolytic T lymphocyte clones and final validation in functional assays of tumour cell recognition. We conclude that the overall low sensitivity and yield of every prediction step often requires a compensatory up-scaling of the initial number of candidate sequences to be screened, rendering reverse immunology an unexpectedly complex approach.

Functional analysis of a conserved DNA methyltransferase in caulobacter crescentus

CcrM is a DNA methyltransferase that methylates the adenine in GANTC motifs in the chromo-some of the bacterial model Caulobacter crescentus. The loss of the CcrM homolog is lethal in C. crescentus and in several other species of Alphaproteobacteria. In this research, we used different experimental and bioinformatic approaches to determine why CcrM is so critical to the physiology of C. crescentus. We first showed that CcrM is a resident orphan DNA methyltransferase in non-Rickettsiales Alphaproteobacteria and that its gene is strictly conserved in this clade (with only one ex¬ception among the genomes sequenced so far). In C. crescentus, cells depleted in CcrM in rich medium quickly lose viability and present an elongated phenotype characteristic of an im¬pairment in cell division. Using minimal medium instead of rich medium as selective and main¬tenance substrate, we could generate a AccrM mutant that presents a viability comparable to the wild type strain and only mild morphological defects. On the basis of a transcriptomic ap¬proach, we determined that several genes essential for cell division were downregulated in the AccrM strain in minimal medium. We offered decisive arguments to support that the efficient transcription of two of these genes, ftsZ and mipZ, coding respectively for the Z-ring forming GTPase FtsZ and an inhibitor of FtsZ polymerization needed for the correct positioning of the Z- ring at mid-cell, requires the methylation of an adenine in a conserved GANTC motif located in their core promoter region. We propose a model, according to which the genome of C. crescentus encodes a transcriptional activator that requires a methylated adenine in a GANTC context to bind to DNA and suggest that this transcriptional regulator might be the global cell-cycle regulator GcrA. In addition, combining a classic genetic approach and in vitro evolution experiments, we showed that the mortality and cell division defects of the AccrM strain in rich medium are mainly due to limiting intracellular levels of the FtsZ protein. We also studied the dynamics of GANTC methylation in C. crescentus using the SMRT technol¬ogy developed by Pacific Biosciences. Our findings support the commonly accepted model, accord¬ing to which the methylation state of GANTC motifs varies during the cell cycle of C. crescentus: before the initiation of DNA replication, the GANTC motifs are fully-methylated (methylated on both strands); when the DNA gets replicated, the GANTC motifs become hemi-methylated (methyl¬ated on one strand only) and this occurs at different times during replication for different loci along the chromosome depending on their position relative to the origin of replication; the GANTC mo¬tifs are only remethylated after DNA replication has finished as a consequence of the massive and short-lived expression of CcrM in predivisional cells. About 30 GANTC motifs in the C. crescentus chromosome were found to be undermethylated in most of the bacterial population; these might be protected from CcrM activity by DNA binding proteins and some of them could be involved in methylation-based bistable transcriptional switches. - CcrM est une ADN méthyltransférase qui méthyle les adénines dans le contexte GANTC dans le génome de la bactérie modèle Caulobacter crescentus. La perte de l'homologue de CcrM chez C. crescentus et chez plusieurs autres espèces d'Alphaproteobactéries est létale. Dans le courant de cette recherche, nous tentons de déterminer pourquoi la protéine CcrM est cruciale pour la survie de C. crescentus. Nous démontrons d'abord que CcrM est une adénine méthyltransférase orpheline résidente, dont le gène fait partie du génome minimal partagé par les Alphaprotéobactéries non-Rickettsiales (à une exception près). Lorsqu'une souche de C. crescentus est privée de CcrM, sa viabilité décroît rapi¬dement et ses cellules présentent une morphologie allongée qui suggère que la division cellulaire est inhibée. Nous sommes parvenus à créer une souche AccrM en utilisant un milieu minimum, au lieu du milieu riche classiquement employé, comme milieu de sélection et de maintenance pour la souche. Lorsque nous avons étudié le transcriptome de cette souche de C. crescentus privée de CcrM, nous avons pu constater que plusieurs gènes essentiels pour le bon déroulement de la division cellulaire bactérienne étaient réprimés. En particulier, l'expression adéquate des gènes ftsZ et mipZ - qui codent, respectivement, pour FtsZ, la protéine qui constitue, au milieu de la cellule, un anneau protéique qui initie le processus de division et pour MipZ, un inhibiteur de la polymérisation de FtsZ qui est indispensable pour le bon positionnement de l'anneau FtsZ - est dépendante de la présence d'une adénine méthylée dans un motif GANTC conservé situé dans leur région promotrice. Nous présentons un modèle selon lequel le génome de C. crescentus code pour un facteur de transcription qui exige la présence d'une adénine méthylée dans un contexte GANTC pour s'attacher à l'ADN et nous suggérons qu'il pourrait s'agir du régulateur global du cycle cellulaire GcrA. En outre, nous montrons, en combinant la génétique classique et une approche basée sur l'évolution expérimentale, que la mortalité et l'inhibition de la division cellulaire caractéristiques de la souche àccrMeη milieu riche sont dues à des niveaux excessivement bas de protéine FtsZ. Nous avons aussi étudié la dynamique de la méthylation du chromosome de C. crescentus sur la base de la technologie SMRT développée par Pacific Biosciences. Nous confirmons le modèle communément accepté, qui affirme que l'état de méthylation des motifs GANTC change durant le cycle cellulaire de C. crescentus: les motifs GANTC sont complètement méthylés (méthylés sur les deux brins) avant de début de la réplication de l'ADN; ils deviennent hémi-méthylés (méthylés sur un brin seulement) une fois répliqués, ce qui arrive à différents moments durant la réplication pour différents sites le long du chromosome en fonction de leur position par rapport à l'origine de répli-cation; finalement, les motifs GANTC sont reméthylés après la fin de la réplication du chromosome lorsque la protéine CcrM est massivement, mais très transitoirement, produite. Par ailleurs, nous identifions dans le chromosome de C. crescentus environ 30 motifs GANTC qui restent en perma-nence non-méthylés dans une grande partie de la population bactérienne; ces motifs sont probable-ment protégés de l'action de CcrM par des protéines qui s'attachent à l'ADN et certains d'entre eux pourraient être impliqués dans des mécanismes de régulation générant une transcription bistable.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Domain shuffling in a sensor protein contributed to the evolution of insect pathogenicity in plant-beneficial Pseudomonas protegens.

Pseudomonas protegens is a biocontrol rhizobacterium with a plant-beneficial and an insect pathogenic lifestyle, but it is not understood how the organism switches between the two states. Here, we focus on understanding the function and possible evolution of a molecular sensor that enables P. protegens to detect the insect environment and produce a potent insecticidal toxin specifically during insect infection but not on roots. By using quantitative single cell microscopy and mutant analysis, we provide evidence that the sensor histidine kinase FitF is a key regulator of insecticidal toxin production. Our experimental data and bioinformatic analyses indicate that FitF shares a sensing domain with DctB, a histidine kinase regulating carbon uptake in Proteobacteria. This suggested that FitF has acquired its specificity through domain shuffling from a common ancestor. We constructed a chimeric DctB-FitF protein and showed that it is indeed functional in regulating toxin expression in P. protegens. The shuffling event and subsequent adaptive modifications of the recruited sensor domain were critical for the microorganism to express its potent insect toxin in the observed host-specific manner. Inhibition of the FitF sensor during root colonization could explain the mechanism by which P. protegens differentiates between the plant and insect host. Our study establishes FitF of P. protegens as a prime model for molecular evolution of sensor proteins and bacterial pathogenicity.