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 approaches to the study of HIV‐1 acquisition.

Host genome studies are increasingly available for the study of infectious disease susceptibility. Current technologies include large-scale genotyping, genome-wide screens such as transcriptome and silencing (silencing RNA) studies, and increasingly, the possibility to sequence complete genomes. These approaches are of interest for the study of individuals who remain uninfected despite documented exposure to human immunodeficiency virus type 1. The main limitation remains the ascertainment of exposure and establishing large cohorts of informative individuals. The pattern of enrichment for CCR5 Δ32 homozygosis should serve as the standard for assessing the extent to which a given cohort (of white subjects) includes a large proportion of exposed uninfected individuals.

CNVs and genetic medicine (excitement and consequences of a rediscovery).

The extensive variability of individual human genomes contributes to phenotypic variability. Structural genomic variants, and copy number variants (CNVs) in particular, have recently been rediscovered as contributors to the genomic plasticity and evolution and as pathoetiologic elements for both monogenic and complex traits. Herein we review some of the consequences of CNVs in the context of human inherited diseases.

The genomic impact of 100 million years of social evolution in seven ant species.

Ants (Hymenoptera, Formicidae) represent one of the most successful eusocial taxa in terms of both their geographic distribution and species number. The publication of seven ant genomes within the past year was a quantum leap for socio- and ant genomics. The diversity of social organization in ants makes them excellent model organisms to study the evolution of social systems. Comparing the ant genomes with those of the honeybee, a lineage that evolved eusociality independently from ants, and solitary insects suggests that there are significant differences in key aspects of genome organization between social and solitary insects, as well as among ant species. Altogether, these seven ant genomes open exciting new research avenues and opportunities for understanding the genetic basis and regulation of social species, and adaptive complex systems in general.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (.100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

Genome Biology and Evolution of the Animal Kingdom

Estudi realitzat a partir d’una estada a la Institut J.W. Jenkinson Laboratory for Evolution and Development of the University of Oxford, Regne Unit, entre 2010 i 2012. He estat membre del laboratori del Professor Peter W.H. Holland com a becari post-doctoral Beatriu de Pinós des de setembre de 2010 al setembre de 2012. El nostre projecte de recerca se centra en l'anàlisi genòmic comparatiu del Regne Animal, tot explorant el contingut dels genomes a través de totes les branques de l'arbre dels animals. Totes les referències a les meves publicacions durant aquest post-doc es poden trobar a http://about.me/jordi_paps. Crec que el nombre i la qualitat dels resultats del meu post-doc, un total de 8 publicacions incloent dos articles a la prestigiosa revista Nature, són prova de l'èxit d'aquest post-doc. Prof Peter W. H. Holland (Departament de Zoologia de la Universitat d'Oxford) i jo som coautors de tres articles de genòmica comparativa, resultats directes d'aquest projecte: 1) comparació de families gèniques entre vertebrats invertebrats (Briefings in Functional Genomics), 2) el genoma de l'ostra (publicat a la revista Nature), i 3) els genomes de 6 platihelmints paràsits (acceptat també a Nature). A més, tenim altres 2 treballs en preparació. Un d'ells analitza l'evolució, expressió i funció dels gens Hox al a la tènia Hymenolepis. El perfil fi d'aquests gens clau del desenvolupament esclareix els canvis d'estil de vida dels organismes. A més, durant aquest últim post-doc he participat en diverses col•laboracions, incloent anàlisi de gens d'envelliment a cucs plans, un estudi sobre la filogènia del grup Gastrotricha, una revisió de l'evolució phylum Platyhelminthes, així com un capítol d'un llibre sobre l'evolució dels animals bilaterals. Finalment, gràcies a la beca Beatriu de Pinós, el Prof. Peter W.H. Holland m'ha convidat a formar part del seu equip com un investigador post-doctoral en el seu projecte ERC Advance actual sobre duplicacions genòmiques.

Targeted Gene Deletion and In Vivo Analysis of Putative Virulence Gene Function in the Pathogenic Dermatophyte Arthroderma benhamiae.

Dermatophytes cause the majority of superficial mycoses in humans and animals. However, little is known about the pathogenicity of this specialized group of filamentous fungi, for which molecular research has been limited thus far. During experimental infection of guinea pigs by the human pathogenic dermatophyte Arthroderma benhamiae, we recently detected the activation of the fungal gene encoding malate synthase AcuE, a key enzyme of the glyoxylate cycle. By the establishment of the first genetic system for A. benhamiae, specific ΔacuE mutants were constructed in a wild-type strain and, in addition, in a derivative in which we inactivated the nonhomologous end-joining pathway by deletion of the A. benhamiae KU70 gene. The absence of AbenKU70 resulted in an increased frequency of the targeted insertion of linear DNA by homologous recombination, without notably altering the monitored in vitro growth abilities of the fungus or its virulence in a guinea pig infection model. Phenotypic analyses of ΔacuE mutants and complemented strains depicted that malate synthase is required for the growth of A. benhamiae on lipids, major constituents of the skin. However, mutant analysis did not reveal a pathogenic role of the A. benhamiae enzyme in guinea pig dermatophytosis or during epidermal invasion of the fungus in an in vitro model of reconstituted human epidermis. The presented efficient system for targeted genetic manipulation in A. benhamiae, paired with the analyzed infection models, will advance the functional characterization of putative virulence determinants in medically important dermatophytes.

Presence of Chlamydiales DNA in ticks and fleas suggests that ticks are carriers of Chlamydiae.

The Chlamydiales order includes the Chlamydiaceae, Parachlamydiaceae, Waddliaceae, Simkaniaceae, Criblamydiaceae, Rhabdochlamydiaceae, Clavichlamydiaceae, and Piscichlamydiaceae families. Members of the Chlamydiales order are obligate intracellular bacteria that replicate within eukaryotic cells of different origins including humans, animals, and amoebae. Many of these bacteria are pathogens or emerging pathogens of both humans and animals, but their true diversity is largely underestimated, and their ecology remains to be investigated. Considering their potential threat on human health, it is important to expand our knowledge on the diversity of Chlamydiae, but also to define the host range colonized by these bacteria. Thus, using a new pan-Chlamydiales PCR, we analyzed the prevalence of Chlamydiales DNA in ticks and fleas, which are important vectors of several viral and bacterial infectious diseases. To conduct this study, 1340 Ixodes ricinus ticks prepared in 192 pools were collected in Switzerland and 55 other ticks belonging to different tick species and 97 fleas belonging to different flea species were harvested in Algeria. In Switzerland, the prevalence of Chlamydiales DNA in the 192 pools was equal to 28.1% (54/192) which represents an estimated prevalence in the 1340 individual ticks of between 4.0% and 28.4%. The pan-Chlamydiales qPCR was positive for 45% (25/55) of tick samples collected in Algeria. The sequencing of the positive qPCR amplicons revealed a high diversity of Chlamydiales species. Most of them belonged to the Rhabdochlamydiaceae and Parachlamydiaceae families. Thus, ticks may carry Chlamydiales and should thus be considered as possible vectors for Chlamydiales propagation to both humans and animals.

Network-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanoma.

Cancer genomes frequently contain somatic copy number alterations (SCNA) that can significantly perturb the expression level of affected genes and thus disrupt pathways controlling normal growth. In melanoma, many studies have focussed on the copy number and gene expression levels of the BRAF, PTEN and MITF genes, but little has been done to identify new genes using these parameters at the genome-wide scale. Using karyotyping, SNP and CGH arrays, and RNA-seq, we have identified SCNA affecting gene expression ('SCNA-genes') in seven human metastatic melanoma cell lines. We showed that the combination of these techniques is useful to identify candidate genes potentially involved in tumorigenesis. Since few of these alterations were recurrent across our samples, we used a protein network-guided approach to determine whether any pathways were enriched in SCNA-genes in one or more samples. From this unbiased genome-wide analysis, we identified 28 significantly enriched pathway modules. Comparison with two large, independent melanoma SCNA datasets showed less than 10% overlap at the individual gene level, but network-guided analysis revealed 66% shared pathways, including all but three of the pathways identified in our data. Frequently altered pathways included WNT, cadherin signalling, angiogenesis and melanogenesis. Additionally, our results emphasize the potential of the EPHA3 and FRS2 gene products, involved in angiogenesis and migration, as possible therapeutic targets in melanoma. Our study demonstrates the utility of network-guided approaches, for both large and small datasets, to identify pathways recurrently perturbed in cancer.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are “genomic fossils” valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome’s structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction (∼80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Comparative gene prediction in human and mouse

The completion of the sequencing of the mouse genome promises to help predict human genes with greater accuracy. While current ab initio gene prediction programs are remarkably sensitive (i.e., they predict at least a fragment of most genes), their specificity is often low, predicting a large number of false-positive genes in the human genome. Sequence conservation at the protein level with the mouse genome can help eliminate some of those false positives. Here we describe SGP2, a gene prediction program that combines ab initio gene prediction with TBLASTX searches between two genome sequences to provide both sensitive and specific gene predictions. The accuracy of SGP2 when used to predict genes by comparing the human and mouse genomes is assessed on a number of data sets, including single-gene data sets, the highly curated human chromosome 22 predictions, and entire genome predictions from ENSEMBL. Results indicate that SGP2 outperforms purely ab initio gene prediction methods. Results also indicate that SGP2 works about as well with 3x shotgun data as it does with fully assembled genomes. SGP2 provides a high enough specificity that its predictions can be experimentally verified at a reasonable cost. SGP2 was used to generate a complete set of gene predictions on both the human and mouse by comparing the genomes of these two species. Our results suggest that another few thousand human and mouse genes currently not in ENSEMBL are worth verifying experimentally.

Tandem chimeric transcription as a means to increase protein diversity in the human genome

The “one-gene, one-protein” rule, coined by Beadle and Tatum, has been fundamental to molecular biology. The rule implies that the genetic complexity of an organism depends essentially on its gene number. The discovery, however, that alternative gene splicing and transcription are widespread phenomena dramatically altered our understanding of the genetic complexity of higher eukaryotic organisms; in these, a limited number of genes may potentially encode a much larger number of proteins. Here we investigate yet another phenomenon that may contribute to generate additional protein diversity. Indeed, by relying on both computational and experimental analysis, we estimate that at least 4%–5% of the tandem gene pairs in the human genome can be eventually transcribed into a single RNA sequence encoding a putative chimeric protein. While the functional significance of most of these chimeric transcripts remains to be determined, we provide strong evidence that this phenomenon does not correspond to mere technical artifacts and that it is a common mechanism with the potential of generating hundreds of additional proteins in the human genome.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.