Convergent evolution of [D-Leucine1] microcystin-LR in taxonomically disparate cyanobacteria

Abstract Background Many important toxins and antibiotics are produced by non-ribosomal biosynthetic pathways. Microcystins are a chemically diverse family of potent peptide toxins and the end-products of a hybrid NRPS and PKS secondary metabolic pathway. They are produced by a variety of cyanobacteria and are responsible for the poisoning of humans as well as the deaths of wild and domestic animals around the world. The chemical diversity of the microcystin family is attributed to a number of genetic events that have resulted in the diversification of the pathway for microcystin assembly. Results Here, we show that independent evolutionary events affecting the substrate specificity of the microcystin biosynthetic pathway have resulted in convergence on a rare [D-Leu1] microcystin-LR chemical variant. We detected this rare microcystin variant from strains of the distantly related genera Microcystis, Nostoc, and Phormidium. Phylogenetic analysis performed using sequences of the catalytic domains within the mcy gene cluster demonstrated a clear recombination pattern in the adenylation domain phylogenetic tree. We found evidence for conversion of the gene encoding the McyA2 adenylation domain in strains of the genera Nostoc and Phormidium. However, point mutations affecting the substrate-binding sequence motifs of the McyA2 adenylation domain were associated with the change in substrate specificity in two strains of Microcystis. In addition to the main [D-Leu1] microcystin-LR variant, these two strains produced a new microcystin that was identified as [Met1] microcystin-LR. Conclusions Phylogenetic analysis demonstrated that both point mutations and gene conversion result in functional mcy gene clusters that produce the same rare [D-Leu1] variant of microcystin in strains of the genera Microcystis, Nostoc, and Phormidium. Engineering pathways to produce recombinant non-ribosomal peptides could provide new natural products or increase the activity of known compounds. Our results suggest that the replacement of entire adenylation domains could be a more successful strategy to obtain higher specificity in the modification of the non-ribosomal peptides than point mutations.

Design and implementation of bioinformatics tools for large scale genome annotation

The continuous increase of genome sequencing projects produced a huge amount of data in the last 10 years: currently more than 600 prokaryotic and 80 eukaryotic genomes are fully sequenced and publically available. However the sole sequencing process of a genome is able to determine just raw nucleotide sequences. This is only the first step of the genome annotation process that will deal with the issue of assigning biological information to each sequence. The annotation process is done at each different level of the biological information processing mechanism, from DNA to protein, and cannot be accomplished only by in vitro analysis procedures resulting extremely expensive and time consuming when applied at a this large scale level. Thus, in silico methods need to be used to accomplish the task. The aim of this work was the implementation of predictive computational methods to allow a fast, reliable, and automated annotation of genomes and proteins starting from aminoacidic sequences. The first part of the work was focused on the implementation of a new machine learning based method for the prediction of the subcellular localization of soluble eukaryotic proteins. The method is called BaCelLo, and was developed in 2006. The main peculiarity of the method is to be independent from biases present in the training dataset, which causes the over‐prediction of the most represented examples in all the other available predictors developed so far. This important result was achieved by a modification, made by myself, to the standard Support Vector Machine (SVM) algorithm with the creation of the so called Balanced SVM. BaCelLo is able to predict the most important subcellular localizations in eukaryotic cells and three, kingdom‐specific, predictors were implemented. In two extensive comparisons, carried out in 2006 and 2008, BaCelLo reported to outperform all the currently available state‐of‐the‐art methods for this prediction task. BaCelLo was subsequently used to completely annotate 5 eukaryotic genomes, by integrating it in a pipeline of predictors developed at the Bologna Biocomputing group by Dr. Pier Luigi Martelli and Dr. Piero Fariselli. An online database, called eSLDB, was developed by integrating, for each aminoacidic sequence extracted from the genome, the predicted subcellular localization merged with experimental and similarity‐based annotations. In the second part of the work a new, machine learning based, method was implemented for the prediction of GPI‐anchored proteins. Basically the method is able to efficiently predict from the raw aminoacidic sequence both the presence of the GPI‐anchor (by means of an SVM), and the position in the sequence of the post‐translational modification event, the so called ω‐site (by means of an Hidden Markov Model (HMM)). The method is called GPIPE and reported to greatly enhance the prediction performances of GPI‐anchored proteins over all the previously developed methods. GPIPE was able to predict up to 88% of the experimentally annotated GPI‐anchored proteins by maintaining a rate of false positive prediction as low as 0.1%. GPIPE was used to completely annotate 81 eukaryotic genomes, and more than 15000 putative GPI‐anchored proteins were predicted, 561 of which are found in H. sapiens. In average 1% of a proteome is predicted as GPI‐anchored. A statistical analysis was performed onto the composition of the regions surrounding the ω‐site that allowed the definition of specific aminoacidic abundances in the different considered regions. Furthermore the hypothesis that compositional biases are present among the four major eukaryotic kingdoms, proposed in literature, was tested and rejected. All the developed predictors and databases are freely available at: BaCelLo http://gpcr.biocomp.unibo.it/bacello eSLDB http://gpcr.biocomp.unibo.it/esldb GPIPE http://gpcr.biocomp.unibo.it/gpipe

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

Large scale variation in Enterococcus faecalis illustrated by the genome analysis of strain OG1RF.

BACKGROUND: Enterococcus faecalis has emerged as a major hospital pathogen. To explore its diversity, we sequenced E. faecalis strain OG1RF, which is commonly used for molecular manipulation and virulence studies. RESULTS: The 2,739,625 base pair chromosome of OG1RF was found to contain approximately 232 kilobases unique to this strain compared to V583, the only publicly available sequenced strain. Almost no mobile genetic elements were found in OG1RF. The 64 areas of divergence were classified into three categories. First, OG1RF carries 39 unique regions, including 2 CRISPR loci and a new WxL locus. Second, we found nine replacements where a sequence specific to V583 was substituted by a sequence specific to OG1RF. For example, the iol operon of OG1RF replaces a possible prophage and the vanB transposon in V583. Finally, we found 16 regions that were present in V583 but missing from OG1RF, including the proposed pathogenicity island, several probable prophages, and the cpsCDEFGHIJK capsular polysaccharide operon. OG1RF was more rapidly but less frequently lethal than V583 in the mouse peritonitis model and considerably outcompeted V583 in a murine model of urinary tract infections. CONCLUSION: E. faecalis OG1RF carries a number of unique loci compared to V583, but the almost complete lack of mobile genetic elements demonstrates that this is not a defining feature of the species. Additionally, OG1RF's effects in experimental models suggest that mediators of virulence may be diverse between different E. faecalis strains and that virulence is not dependent on the presence of mobile genetic elements.

Pdsg1 and Pdsg2, novel proteins involved in developmental genome remodelling in Paramecium

The epigenetic influence of maternal cells on the development of their progeny has long been studied in various eukaryotes. Multicellular organisms usually provide their zygotes not only with nutrients but also with functional elements required for proper development, such as coding and non-coding RNAs. These maternally deposited RNAs exhibit a variety of functions, from regulating gene expression to assuring genome integrity. In ciliates, such as Paramecium these RNAs participate in the programming of large-scale genome reorganization during development, distinguishing germline-limited DNA, which is excised, from somatic-destined DNA. Only a handful of proteins playing roles in this process have been identified so far, including typical RNAi-derived factors such as Dicer-like and Piwi proteins. Here we report and characterize two novel proteins, Pdsg1 and Pdsg2 (Paramecium protein involved in Development of the Somatic Genome 1 and 2), involved in Paramecium genome reorganization. We show that these proteins are necessary for the excision of germline-limited DNA during development and the survival of sexual progeny. Knockdown of PDSG1 and PDSG2 genes affects the populations of small RNAs known to be involved in the programming of DNA elimination (scanRNAs and iesRNAs) and chromatin modification patterns during development. Our results suggest an association between RNA-mediated trans-generational epigenetic signal and chromatin modifications in the process of Paramecium genome reorganization.

GCN5 maintains genome stability during development

The histone acetyltransferase, GCN5, is essential for survival of mice during embryogenesis. GCN5 null embryos die early during development due to increased apoptosis. We have demonstrated that the increased apoptosis in associated with increased p53 protein levels. Loss of p53 rescues the embryonic apoptosis in the GCN5 null embryos. These results raised the question of what molecular trigger leads to p53 stabilization and cell death in the absence of GCN5. p53 is generally referred to as the gatekeeper of the cell, monitoring cellular responses to DNA damage, genotoxic stress, and other unfavorable conditions in the cell. Therefore, we examined individual cells in wild type and mutant embryos for gross chromosomal aberrations that might trigger a genome integrity checkpoint. Karyotype analysis indicates that approximately 30% of the cells in an E8.5 GCN5 null embryo display chromosomal aberrations, predominantly chromosomal end adhesions and associations. In wild type E8.5 embryos, only 6% of the cells have chromosomal aberrations. Recent data using telomeric FISH demonstrates that cells from GCN5 null embryos have a decreased telomeric signal. Telomere maintenance is essential for maintaining genome integrity. Telomeric defects are associated with loss of chromosomes and chromosomal rearrangements that can lead to detrimental gene fusions involved in many types of cancers. Little is known about the chromatin structures present near the telomeric ends, or whether any of the telomere-associated proteins are subject to post-translational modification such as acetylation. Our results are the first data to demonstrate the involvement of a histone acetyltransferase, GCN5, in maintaining genome integrity through telomere maintenance and/or capping. ^

Pathway analysis of genome-wide association study data in primary biliary cirrhosisa

Genome-wide association studies (GWAS) have successfully identified several genetic loci associated with inherited predisposition to primary biliary cirrhosis (PBC), the most common autoimmune disease of the liver. Pathway-based tests constitute a novel paradigm for GWAS analysis. By evaluating genetic variation across a biological pathway (gene set), these tests have the potential to determine the collective impact of variants with subtle effects that are individually too weak to be detected in traditional single variant GWAS analysis. To identify biological pathways associated with the risk of development of PBC, GWAS of PBC from Italy (449 cases and 940 controls) and Canada (530 cases and 398 controls) were independently analyzed. The linear combination test (LCT), a recently developed pathway-level statistical method was used for this analysis. For additional validation, pathways that were replicated at the P <0.05 level of significance in both GWAS on LCT analysis were also tested for association with PBC in each dataset using two complementary GWAS pathway approaches. The complementary approaches included a modification of the gene set enrichment analysis algorithm (i-GSEA4GWAS) and Fisher's exact test for pathway enrichment ratios. Twenty-five pathways were associated with PBC risk on LCT analysis in the Italian dataset at P<0.05, of which eight had an FDR<0.25. The top pathway in the Italian dataset was the TNF/stress related signaling pathway (p=7.38×10 -4, FDR=0.18). Twenty-six pathways were associated with PBC at the P<0.05 level using the LCT in the Canadian dataset with the regulation and function of ChREBP in liver pathway (p=5.68×10-4, FDR=0.285) emerging as the most significant pathway. Two pathways, phosphatidylinositol signaling system (Italian: p=0.016, FDR=0.436; Canadian: p=0.034, FDR=0.693) and hedgehog signaling (Italian: p=0.044, FDR=0.636; Canadian: p=0.041, FDR=0.693), were replicated at LCT P<0.05 in both datasets. Statistically significant association of both pathways with PBC genetic susceptibility was confirmed in the Italian dataset on i-GSEA4GWAS. Results for the phosphatidylinositol signaling system were also significant in both datasets on applying Fisher's exact test for pathway enrichment ratios. This study identified a combination of known and novel pathway-level associations with PBC risk. If functionally validated, the findings may yield fresh insights into the etiology of this complex autoimmune disease with possible preventive and therapeutic application.^

Genome-Wide Analysis of the Response of Dickeya dadantii 3937 to Plant Antimicrobial Peptides

Antimicrobial peptides constitute an important factor in the defense of plants against pathogens, and bacterial resistance to these peptides have previously been shown to be an important virulence factor in Dickeya dadantii, the causal agent of soft-rot disease of vegetables. In order to understand the bacterial response to antimicrobial pep- tides, a transcriptional microarray analysis was performed upon treatment with sub-lethal concentration of thionins, a widespread plant peptide. In all, 36 genes were found to be overexpressed, and were classified according to their deduced function as i) transcriptional regulators, ii) transport, and iii) modification of the bacterial membrane. One gene encoding a uricase was found to be repressed. The majority of these genes are known to be under the control of the PhoP/PhoQ system. Five genes representing the different functions induced were selected for further analysis. The results obtained indicate that the presence of antimicrobial peptides induces a complex response which includes peptide-specific elements and general stress-response elements contributing differentially to the virulence in different hosts.

BAC-VAC, a novel generation of (DNA) vaccines: A bacterial artificial chromosome (BAC) containing a replication-competent, packaging-defective virus genome induces protective immunity against herpes simplex virus 1

This study aimed to exploit bacterial artificial chromosomes (BAC) as large antigen-capacity DNA vaccines (BAC-VAC) against complex pathogens, such as herpes simplex virus 1 (HSV-1). The 152-kbp HSV-1 genome recently has been cloned as an F-plasmid-based BAC in Escherichia coli (fHSV), which can efficiently produce infectious virus progeny upon transfection into mammalian cells. A safe modification of fHSV, fHSVΔpac, does not give rise to progeny virus because the signals necessary to package DNA into virions have been excluded. However, in mammalian cells fHSVΔpac DNA can still replicate, express the HSV-1 genes, cause cytotoxic effects, and produce virus-like particles. Because these functions mimic the lytic cycle of the HSV-1 infection, fHSVΔpac was expected to stimulate the immune system as efficiently as a modified live virus vaccine. To test this hypothesis, mice were immunized with fHSVΔpac DNA applied intradermally by gold-particle bombardment, and the immune responses were compared with those induced by infection with disabled infectious single cycle HSV-1. Immunization with either fHSVΔpac or disabled infectious single cycle HSV-1 induced the priming of HSV-1-specific cytotoxic T cells and the production of virus-specific antibodies and conferred protection against intracerebral injection of wild-type HSV-1 at a dose of 200 LD50. Protection probably was cell-mediated, as transfer of serum from immunized mice did not protect naive animals. We conclude that BAC-VACs per se, or in combination with genetic elements that support replicative amplification of the DNA in the cell nucleus, represent a useful new generation of DNA-based vaccination strategies for many viral and nonviral antigens.

Identification of an ATP-binding cassette transporter involved in bicarbonate uptake in the cyanobacterium Synechococcus sp. strain PCC 7942

Exposure of cells of cyanobacteria (blue–green algae) grown under high-CO2 conditions to inorganic C-limitation induces transcription of particular genes and expression of high-affinity CO2 and HCO3− transport systems. Among the low-CO2-inducible transcription units of Synechococcus sp. strain PCC 7942 is the cmpABCD operon, encoding an ATP-binding cassette transporter similar to the nitrate/nitrite transporter of the same cyanobacterium. A nitrogen-regulated promoter was used to selectively induce expression of the cmpABCD genes by growth of transgenic cells on nitrate under high CO2 conditions. Measurements of the initial rate of HCO3− uptake after onset of light, and of the steady-state rate of HCO3− uptake in the light, showed that the controlled induction of the cmp genes resulted in selective expression of high-affinity HCO3− transport activity. The forced expression of cmpABCD did not significantly increase the CO2 uptake capabilities of the cells. These findings demonstrated that the cmpABCD genes encode a high-affinity HCO3− transporter. A deletion mutant of cmpAB (M42) retained low CO2-inducible activity of HCO3− transport, indicating the occurrence of HCO3− transporter(s) distinct from the one encoded by cmpABCD. HCO3− uptake by low-CO2-induced M42 cells showed lower affinity for external HCO3− than for wild-type cells under the same conditions, showing that the HCO3− transporter encoded by cmpABCD has the highest affinity for HCO3− among the HCO3− transporters present in the cyanobacterium. This appears to be the first unambiguous identification and description of a primary active HCO3− transporter.

Nucleotide binding and autophosphorylation of the clock protein KaiC as a circadian timing process of cyanobacteria

A negative feedback control of kaiC expression by KaiC protein has been proposed to generate a basic oscillation of the circadian clock in the cyanobacterium Synechococcus sp. PCC 7942. KaiC has two P loops or Walker's motif As, that are potential ATP-/GTP-binding motifs and DXXG motifs conserved in various GTP-binding proteins. Herein, we demonstrate that in vitro KaiC binds ATP and, with lower affinity, GTP. Point mutation by site-directed mutagenesis of P loop 1 completely nullified the circadian rhythm of kaiBC expression and markedly reduced ATP-binding activity. Moreover, KaiC can be autophosphorylated in vitro. These results suggest that the nucleotide-binding activity of KaiC plays important roles in the generation of circadian oscillation in cyanobacteria.

Comparative genomics of the restriction-modification systems in Helicobacter pylori

Helicobacter pylori is a Gram-negative bacterial pathogen with a small genome of 1.64–1.67 Mb. More than 20 putative DNA restriction-modification (R-M) systems, comprising more than 4% of the total genome, have been identified in the two completely sequenced H. pylori strains, 26695 and J99, based on sequence similarities. In this study, we have investigated the biochemical activities of 14 Type II R-M systems in H. pylori 26695. Less than 30% of the Type II R-M systems in 26695 are fully functional, similar to the results obtained from strain J99. Although nearly 90% of the R-M genes are shared by the two H. pylori strains, different sets of these R-M genes are functionally active in each strain. Interestingly, all strain-specific R-M genes are active, whereas most shared genes are inactive. This agrees with the notion that strain-specific genes have been acquired more recently through horizontal transfer from other bacteria and selected for function. Thus, they are less likely to be impaired by random mutations. Our results also show that H. pylori has extremely diversified R-M systems in different strains, and that the diversity may be maintained by constantly acquiring new R-M systems and by inactivating and deleting the old ones.

Circadian gating of cell division in cyanobacteria growing with average doubling times of less than 24 hours.

To ascertain whether the circadian oscillator in the prokaryotic cyanobacterium Synechococcus PCC 7942 regulates the timing of cell division in rapidly growing cultures, we measured the rate of cell division, DNA content, cell size, and gene expression (monitored by luminescence of the PpsbAI::luxAB reporter) in cultures that were continuously diluted to maintain an approximately equal cell density. We found that populations dividing at rates as rapid as once per 10 h manifest circadian gating of cell division, since phases in which cell division slows or stops recur with a circadian periodicity. The data clearly show that Synechococcus cells growing with doubling times that are considerably faster than once per 24 h nonetheless express robust circadian rhythms of cell division and gene expression. Apparently Synechococcus cells are able to simultaneously sustain two timing circuits that express significantly different periods.

The notion of the Cambrian pananimalia genome.

The toil by photosynthesizing cyanobacteria and blue-green algae of nearly three billion years appeared to have finally resulted in the sufficient accumulation of molecular oxygen. So, the stage was set for the emergence, at the ocean bottom, of diverse animals that were consumers of molecular oxygen. It now appears that this Cambrian explosion, during which nearly all the extant animal phyla have emerged, was of an astonishingly short duration, lasting only 6-10 million years. Inasmuch as only a 1% DNA base sequence change is expected in 10 million years under the standard spontaneous mutation rate, I propose that all those diverse animals of the early Cambrian period, some 550 million years ago, were endowed with nearly identical genomes, with differential usage of the same set of genes accounting for the extreme diversities of body forms. Some of the more pertinent genes that are thought to be included in the Cambrian pananimalia genome are as follows. (i) A gene for lysyloxidase that, in the presence of molecular oxygen, crosslinked collagen triple helices to produce ligaments and tendons, thus contributing to the stout bodies of the Cambrian animals. (ii) Genes for hemoglobin; these internal transporters of molecular oxygen are today seen sporadically in members of diverse animal phyla. (iii) The Pax-6 gene for eye formation; the eyes of a ribbon worm to a human are organized by this gene. In animals without eyes, the same gene organizes other sensory systems and organs. (iv) A series of Hox genes for the anterior-posterior (cranio-caudal) body plans: these genes are also present in all phyla of the kingdom Animalia.