Characterization of mitochondrial genomes in bivalve species with doubly uniparental inheritance of mitochondria

Many bivalve species possess two distinct mtDNA lineages, called F and M, respectively inherited maternally and paternally: this system is called doubly uniparental inheritance (DUI). The main experimental project of my PhD was the quantification of the two mtDNAs during the development of the DUI species Ruditapes philippinarum, from early embryos to sub-adults, using Real-Time qPCR. I identified the time interval in which M mtDNA is lost from female individuals, while it is retained in males (which are heteroplasmic through all of their life cycle). The results also suggested absence of mtDNA replication during early embryogenesis, a process constituting a bottleneck that highly reduces the copy number of mtDNA molecules in cells of developing larvae. In males this bottleneck may produce cells homoplasmic for M mtDNA, and could be considered as a first step of the segregation of M in the male germ line. Another finding was the characterization, in young clams approaching the first reproductive season, of a significant boost in copy number of F mtDNA in females and of M in males. Given the age of animals in which this mtDNA-specific growth was observed, the finding could probably be the outcome of the first round of gonads and gametes production. Other lines of research included the characterization of the unassigned regions in mt genomes of DUI bivalves. These regions can harbor signals involved in the control of replication and/or transcription of the mtDNA molecule, as well as additional open reading frames (ORFs) not related to oxidative phosphorylation. These features in DUI species could be associated to the maintenance of separate inheritance routes for the two mtDNAs. Additional ORFs are also found in other animal mt genomes: I summarized the presence of gene duplications as a co-author in a review focusing on animal mt genomes with unusual gene content.

Identification of Drosophila heart-specific Cis-Regulatory Modules under Hox control

Cardiac morphogenesis is a complex process governed by evolutionarily conserved transcription factors and signaling molecules. The Drosophila cardiac tube is linear, made of 52 pairs of cardiomyocytes (CMs), which express specific transcription factor genes that have human homologues implicated in Congenital Heart Diseases (CHDs) (NKX2-5, GATA4 and TBX5). The Drosophila cardiac tube is linear and composed of a rostral portion named aorta and a caudal one called heart, distinguished by morphological and functional differences controlled by Hox genes, key regulators of axial patterning. Overexpression and inactivation of the Hox gene abdominal-A (abd-A), which is expressed exclusively in the heart, revealed that abd-A controls heart identity. The aim of our work is to isolate the heart-specific cisregulatory sequences of abd-A direct target genes, the realizator genes granting heart identity. In each segment of the heart, four pairs of cardiomyocytes (CMs) express tinman (tin), homologous to NKX2-5, and acquire strong contractile and automatic rhythmic activities. By tyramide amplified FISH, we found that seven genes, encoding ion channels, pumps or transporters, are specifically expressed in the Tin-CMs of the heart. We initially used online available tools to identify their heart-specific cisregutatory modules by looking for Conserved Non-coding Sequences containing clusters of binding sites for various cardiac transcription factors, including Hox proteins. Based on these data we generated several reporter gene constructs and transgenic embryos, but none of them showed reporter gene expression in the heart. In order to identify additional abd-A target genes, we performed microarray experiments comparing the transcriptomes of aorta versus heart and identified 144 genes overexpressed in the heart. In order to find the heart-specific cis-regulatory regions of these target genes we developed a new bioinformatic approach where prediction is based on pattern matching and ordered statistics. We first retrieved Conserved Noncoding Sequences from the alignment between the D.melanogaster and D.pseudobscura genomes. We scored for combinations of conserved occurrences of ABD-A, ABD-B, TIN, PNR, dMEF2, MADS box, T-box and E-box sites and we ranked these results based on two independent strategies. On one hand we ranked the putative cis-regulatory sequences according to best scored ABD-A biding sites, on the other hand we scored according to conservation of binding sites. We integrated and ranked again the two lists obtained independently to produce a final rank. We generated nGFP reporter construct flies for in vivo validation. We identified three 1kblong heart-specific enhancers. By in vivo and in vitro experiments we are determining whether they are direct abd-A targets, demonstrating the role of a Hox gene in the realization of heart identity. The identified abd-A direct target genes may be targets also of the NKX2-5, GATA4 and/or TBX5 homologues tin, pannier and Doc genes, respectively. The identification of sequences coregulated by a Hox protein and the homologues of transcription factors causing CHDs, will provide a mean to test whether these factors function as Hox cofactors granting cardiac specificity to Hox proteins, increasing our knowledge on the molecular mechanisms underlying CHDs. Finally, it may be investigated whether these Hox targets are involved in CHDs.

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

Computational methods for genome screening

Motivation An actual issue of great interest, both under a theoretical and an applicative perspective, is the analysis of biological sequences for disclosing the information that they encode. The development of new technologies for genome sequencing in the last years, opened new fundamental problems since huge amounts of biological data still deserve an interpretation. Indeed, the sequencing is only the first step of the genome annotation process that consists in the assignment of biological information to each sequence. Hence given the large amount of available data, in silico methods became useful and necessary in order to extract relevant information from sequences. The availability of data from Genome Projects gave rise to new strategies for tackling the basic problems of computational biology such as the determination of the tridimensional structures of proteins, their biological function and their reciprocal interactions. Results The aim of this work has been the implementation of predictive methods that allow the extraction of information on the properties of genomes and proteins starting from the nucleotide and aminoacidic sequences, by taking advantage of the information provided by the comparison of the genome sequences from different species. In the first part of the work a comprehensive large scale genome comparison of 599 organisms is described. 2,6 million of sequences coming from 551 prokaryotic and 48 eukaryotic genomes were aligned and clustered on the basis of their sequence identity. This procedure led to the identification of classes of proteins that are peculiar to the different groups of organisms. Moreover the adopted similarity threshold produced clusters that are homogeneous on the structural point of view and that can be used for structural annotation of uncharacterized sequences. The second part of the work focuses on the characterization of thermostable proteins and on the development of tools able to predict the thermostability of a protein starting from its sequence. By means of Principal Component Analysis the codon composition of a non redundant database comprising 116 prokaryotic genomes has been analyzed and it has been showed that a cross genomic approach can allow the extraction of common determinants of thermostability at the genome level, leading to an overall accuracy in discriminating thermophilic coding sequences equal to 95%. This result outperform those obtained in previous studies. Moreover, we investigated the effect of multiple mutations on protein thermostability. This issue is of great importance in the field of protein engineering, since thermostable proteins are generally more suitable than their mesostable counterparts in technological applications. A Support Vector Machine based method has been trained to predict if a set of mutations can enhance the thermostability of a given protein sequence. The developed predictor achieves 88% accuracy.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

Mitochondrial dysfunction in hereditary optic neuropathies

MITOCHONDRIAL DYSFUNCTION IN HEREDITARY OPTIC NEUROPATHIES Mitochondrial pathologies are a heterogeneous group of clinical manifestations characterized by oxidative phosphorylation impairment. At the beginning of their recognition mitochondrial pathologies were regarded as rare disorders but indeed they are more frequent than originally thought. Due to the unique mitochondria peculiarities mitochondrial pathologies can be caused by mutations in both mitochondrial and nuclear genomes. The poor knowledge of pathologic mechanism of these disorders has not allowed a real development of the “mitochondrial medicine”, that is currently limited to symptoms mitigation. Leber hereditary optic neuropathy (LHON) was the first pathology to be linked to a point mutation in the mtDNA. The mechanism by which point mutations in mitochondrial gene encoding Complex I subunits leads to optic nerve degeneration is still unknown, although is well accepted that other genetic or environmental factors are involved in the modulation of pathology, where a pivotal role is certainly played by oxidative stress. We studied the relationship between the Ala16Val dimorphism in the mitochondrial targeting sequence of nuclear gene SOD2 and the 3460/ND1 LHON mutation. Our results show that, in control population, the heterozygous SOD2 genotype is associated to a higher activity and quantity of MnSOD, particularly with respect to Val homozygotes. Furthermore, we demonstrated that LHON patients harboring at least one Ala allele are characterized by an increased MnSOD activity with respect to relative control population. Since the ATP synthesis rate – severely reduced in LHON patients lymphocytes - is not affected by the SOD2 genotype, we concluded that SOD2 gene could modulate the pathogenicity of LHON mutations through a mechanism associated to an increase of reactive oxygen species production. Autosomal dominant optic atrophy (ADOA) is a pathology linked to mutations in nuclear gene encoding Opa1, a dynamin-related protein localized in the mitochondrial matrix. Although the clinical course is slightly different, the endpoint of ADOA is exactly the same of LHON: optic nerve degeneration with specific involvement of retinal ganglion cells. Opa1 is a relatively new protein, whose major role is the regulation of mitochondrial fusion. Mitochondrial morphology is the results of the equilibrium between two opposite force: fusion and fission, two processes that have to be finely regulated in order to preserve mitochondrial and cellular physiology. We studied fibroblasts deriving from ADOA patients characterized by a new deletion in the GTPase domain of the OPA1 gene. The biochemical characterization of ADOA and control fibroblasts has concerned the evaluation of ATP synthesis rate, mitochondrial membrane potential in different metabolic conditions and the morphological status of mitochondria. Regarding ATP synthesis rate we did not find significant differences between ADOA and control fibroblasts even though a trend toward increased reduction in ADOA samples is observed when fibroblasts are grown in absence of glucose or in the medium containing gramicidin. Furthermore, we found that also in ADOA fibroblasts membrane potential is actively maintained by proton pumping of fully functional respiratory chain complexes. Our results indicate that the mutation found in the pedigree analyzed acts primary impairing the mitochondrial fusion without affecting the energy production, supporting the notion that cell function is tightly linked to mitochondrial morphology. Mitochondrial dysfunctions are acquiring great attention because of their recognized relevance not only in aging but also in age-related pathologies including cancer, cardiovascular disease, type II diabetes, and neurodegenerative disorders. The involvement of mitochondria in such detrimental pathologies that, currently, have become so common enhances the necessity of standardization of therapeutic strategies capable of rescuing the normal mitochondrial function. In order to propose an alternative treatment for energy deficiency-disorders we tested the effect of substrates capable to stimulate the substrate-level phosphorylation on viability and energy availability in different experimental models grown under different metabolic conditions. In fibroblasts, the energy defect was achieved by culturing cells in presence of oligomycin, an inhibitor of ATP synthase complex. NARP cybrids have been used as model of mitochondrial pathology. Cell viability and ATP content have been considered as parameters to assay the capability of exogenous substrate to rescue energy failure. Our results suggest that patients suffering for some forms of ATP synthase deficiency, or characterized by a deficiency in energy production, might benefit from dietary or pharmacological treatment based on supplementation of α-ketoglutarate and aspartate.

Meccanismi evolutivi dei parapoxvirus: caratterizzazione genomica di pseudocowpoxvirus e messa a punto di sistemi per lo studio delle ricombinazioni

The Poxviruses are a family of double stranded DNA (dsDNA) viruses that cause disease in many species, both vertebrate and invertebrate. Their genomes range in size from 135 to 365 kbp and show conservation in both organization and content. In particular, the central genomic regions of the chordopoxvirus subfamily (those capable of infecting vertebrates) contain 88 genes which are present in all the virus species characterised to date and which mostly occur in the same order and orientation. In contrast, however, the terminal regions of the genomes frequently contain genes that are species or genera-specific and that are not essential for the growth of the virus in vitro but instead often encode factors with important roles in vivo including modulation of the host immune response to infection and determination of the host range of the virus. The Parapoxviruses (PPV), of which Orf virus is the prototypic species, represent a genus within the chordopoxvirus subfamily of Poxviridae and are characterised by their ability to infect ruminants and humans. The genus currently contains four recognised species of virus, bovine papular stomatitis virus (BPSV) and pseudocowpox virus (PCPV) both of which infect cattle, orf virus (OV) that infects sheep and goats, and parapoxvirus of red deer in New Zealand (PVNZ). The ORFV genome has been fully sequenced, as has that of BPSV, and is ~138 kb in length encoding ~132 genes. The vast majority of these genes allow the virus to replicate in the cytoplasm of the infected host cell and therefore encode proteins involved in replication, transcription and metabolism of nucleic acids. These genes are well conserved between all known genera of poxviruses. There is however another class of genes, located at either end of the linear dsDNA genome, that encode proteins which are non-essential for replication and generally dictate host range and virulence of the virus. The non-essential genes are often the most variable within and between species of virus and therefore are potentially useful for diagnostic purposes. Given their role in subverting the host-immune response to infection they are also targets for novel therapeutics. The function of only a relatively small number of these proteins has been elucidated and there are several genes whose function still remains obscure principally because there is little similarity between them and proteins of known function in current sequence databases. It is thought that by selectively removing some of the virulence genes, or at least neutralising the proteins in some way, current vaccines could be improved. The evolution of poxviruses has been proposed to be an adaptive process involving frequent events of gene gain and loss, such that the virus co-evolves with its specific host. Gene capture or horizontal gene transfer from the host to the virus is considered an important source of new viral genes including those likely to be involved in host range and those enabling the virus to interfere with the host immune response to infection. Given the low rate of nucleotide substitution, recombination can be seen as an essential evolutionary driving force although it is likely underestimated. Recombination in poxviruses is intimately linked to DNA replication with both viral and cellular proteins participate in this recombination-dependent replication. It has been shown, in other poxvirus genera, that recombination between isolates and perhaps even between species does occur, thereby providing another mechanism for the acquisition of new genes and for the rapid evolution of viruses. Such events may result in viruses that have a selective advantage over others, for example in re-infections (a characteristic of the PPV), or in viruses that are able to jump the species barrier and infect new hosts. Sequence data related to viral strains isolated from goats suggest that possible recombination events may have occurred between OV and PCPV (Ueda et al. 2003). The recombination events are frequent during poxvirus replication and comparative genomic analysis of several poxvirus species has revealed that recombinations occur frequently on the right terminal region. Intraspecific recombination can occur between strains of the same PPV species, but also interspecific recombination can happen depending on enough sequence similarity to enable recombination between distinct PPV species. The most important pre-requisite for a successful recombination is the coinfection of the individual host by different virus strains or species. Consequently, the following factors affecting the distribution of different viruses to shared target cells need to be considered: dose of inoculated virus, time interval between inoculation of the first and the second virus, distance between the marker mutations, genetic homology. At present there are no available data on the replication dynamics of PPV in permissive and non permissive hosts and reguarding co-infetions there are no information on the interference mechanisms occurring during the simultaneous replication of viruses of different species. This work has been carried out to set up permissive substrates allowing the replication of different PPV species, in particular keratinocytes monolayers and organotypic skin cultures. Furthermore a method to isolate and expand ovine skin stem cells was has been set up to indeep further aspects of viral cellular tropism during natural infection. The study produced important data to elucidate the replication dynamics of OV and PCPV virus in vitro as well as the mechanisms of interference that can arise during co-infection with different viral species. Moreover, the analysis carried on the genomic right terminal region of PCPV 1303/05 contributed to a better knowledge of the viral genes involved in host interaction and pathogenesis as well as to locate recombination breakpoints and genetic homologies between PPV species. Taken together these data filled several crucial gaps for the study of interspecific recombinations of PPVs which are thought to be important for a better understanding of the viral evolution and to improve the biosafety of antiviral therapy and PPV-based vectors.

Mitochondrial Genomics: structure, inheritance and phylogenetic utility of the Mitochondrial genome in Bivalvia and Insecta

This PhD Thesis is the result of my research activity in the last three years. My main research interest was centered on the evolution of mitochondrial genome (mtDNA), and on its usefulness as a phylogeographic and phylogenetic marker at different taxonomic levels in different taxa of Metazoa. From a methodological standpoint, my main effort was dedicated to the sequencing of complete mitochondrial genomes, and the approach to whole-genome sequencing was based on the application of Long-PCR and shotgun sequences. Moreover, this research project is a part of a bigger sequencing project of mtDNAs in many different Metazoans’ taxa, and I mostly dedicated myself to sequence and analyze mtDNAs in selected taxa of bivalves and hexapods (Insecta). Sequences of bivalve mtDNAs are particularly limited, and my study contributed to extend the sampling. Moreover, I used the bivalve Musculista senhousia as model taxon to investigate the molecular mechanisms and the evolutionary significance of their aberrant mode of mitochondrial inheritance (Doubly Uniparental Inheritance, see below). In Insects, I focused my attention on the Genus Bacillus (Insecta Phasmida). A detailed phylogenetic analysis was performed in order to assess phylogenetic relationships within the genus, and to investigate the placement of Phasmida in the phylogenetic tree of Insecta. The main goal of this part of my study was to add to the taxonomic coverage of sequenced mtDNAs in basal insects, which were only partially analyzed.

Study of two novel streptococcus pneumoniae operons: pilus islet 2 and the lantibiotic operon

Analysis of publicly available genomes of Streptococcus pneumoniae has led to the identification of a new genomic element resembling gram-positive pilus islets (PIs). Here, we demonstrate that this genomic region, herein referred to as PI-2 (containing the genes pitA, sipA, pitB, srtG1, and srtG2) codes for a novel functional pilus in pneumococcus. Therefore, there are two pilus islets identified so far in this pathogen (PI-1 and PI-2). Polymerization of the PI-2 pilus requires the backbone protein PitB as well as the sortase SrtG1 and the signal peptidase-like protein SipA. PI-2 is associated with serotypes 1, 2, 7F, 19A, and 19F, considered to be emerging in both industrialized and developing countries. Interestingly, strains belonging to clonal complex 271 (CC271) contain both PI-1 and PI-2, as revealed by genome analyses. In these strains both pili are surface exposed and independently assembled. Furthermore, in vitro experiments provide evidence that the pilus encoded by PI-2 of S. pneumoniae is involved in adherence. Thus, pneumococci encode at least two types of pili that may play a role in the initial host cell contact to the respiratory tract. In addition, the pilus proteins are potential antigens for inclusion in a new generation of pneumococcal vaccines. Adherence by pili could represent important factor in bacterial community formation, since it has been demonstrated that bacterial community formation plays an important role in pneumococcal otitis media. In vitro quantification of bacterial community formation by S. pneumoniae was performed in order to investigate the possible role of pneumococcal pili to form communities. By using different growth media we were not able to see clear association between pili and community formation. But our findings revealed that strains belonging to MLST clonal complex CC15 efficiently form bacterial communities in vitro in a glucose dependent manner. We compared the genome of forty-four pneumococcal isolates discovering four open reading frames specifically associated with CC15. These four genes are annotated as members of an operon responsible for the biosynthesis of a putative lanctibiotic peptide, described to be involved in bacterial community formation. Our experiments show that the lanctibiotic operon deletion affects glucose mediated community formation in CC 15 strain INV200. Moreover, since glucose consumption during bacterial growth produce an acidic environment, we tested bacterial community formation at different pH and we showed that the lanctibiotic operon deletion affected pH mediated community formation in CC 15 strain INV200. In conclusion, these data demonstrate that the putative lanctibiotic operon is associated with pneumococcal CC 15 strains in vitro bacterial community formation.

New insight on a Group A Streptococcus protective antigen contributing to bacterial cell division

Bioinformatic analysis of Group A Streptococcus (GAS) genomes aiming at the identification of new vaccine antigens, revealed the presence of a gene coding for a putative surface-associated protein, named GAS40, inducing protective antibodies in an animal model of sepsis. The aim of our study was to unravel the involvement of GAS40 in cell division processes and to identify the putative interactor. Firstly, bioinformatic analysis showed that gas40 shares homology with ezrA, a gene coding for a negative regulator of Z-ring formation during cell division process. Both scanning and transmission electron microscopy indicated morphological differences between wild-type and the GAS40 knock-out mutant strain, with the latter showing an impaired capacity to divide resulting in the formation of very long chains. Moreover, when the localization of the antigen on the bacterial surface was analyzed, we found that in bacteria grown at exponential phase GAS40 specifically localized at septum, indicating a possible role in cell division. Furthermore, by ELISA and co-sedimentation assays, we found that GAS40 is able to interact with FtsZ, a protein involved in Z-ring formation during cell division process. These data together with the co-localization of GAS40/FtsZ at bacterial septum demonstrated by by confocal microscopy, strongly support the hypothesis for a key role of GAS40 in bacterial cell division.

Structural and functional analysis of centromeric chromatin

Animal neocentromeres are defined as ectopic centromeres that have formed in non-centromeric locations and avoid some of the features, like the DNA satellite sequence, that normally characterize canonical centromeres. Despite this, they are stable functional centromeres inherited through generations. The only existence of neocentromeres provide convincing evidence that centromere specification is determined by epigenetic rather than sequence-specific mechanisms. For all this reasons, we used them as simplified models to investigate the molecular mechanisms that underlay the formation and the maintenance of functional centromeres. We collected human cell lines carrying neocentromeres in different positions. To investigate the region involved in the process at the DNA sequence level we applied a recent technology that integrates Chromatin Immuno-Precipitation and DNA microarrays (ChIP-on-chip) using rabbit polyclonal antibodies directed against CENP-A or CENP-C human centromeric proteins. These DNA binding-proteins are required for kinetochore function and are exclusively targeted to functional centromeres. Thus, the immunoprecipitation of DNA bound by these proteins allows the isolation of centromeric sequences, including those of the neocentromeres. Neocentromeres arise even in protein-coding genes region. We further analyzed if the increased scaffold attachment sites and the corresponding tighter chromatin of the region involved in the neocentromerization process still were permissive or not to transcription of within encoded genes. Centromere repositioning is a phenomenon in which a neocentromere arisen without altering the gene order, followed by the inactivation of the canonical centromere, becomes fixed in population. It is a process of chromosome rearrangement fundamental in evolution, at the bases of speciation. The repeat-free region where the neocentromere initially forms, progressively acquires extended arrays of satellite tandem repeats that may contribute to its functional stability. In this view our attention focalized to the repositioned horse ECA11 centromere. ChIP-on-chip analysis was used to define the region involved and SNPs studies, mapping within the region involved into neocentromerization, were carried on. We have been able to describe the structural polymorphism of the chromosome 11 centromeric domain of Caballus population. That polymorphism was seen even between homologues chromosome of the same cells. That discovery was the first described ever. Genomic plasticity had a fundamental role in evolution. Centromeres are not static packaged region of genomes. The key question that fascinates biologists is to understand how that centromere plasticity could be combined to the stability and maintenance of centromeric function. Starting from the epigenetic point of view that underlies centromere formation, we decided to analyze the RNA content of centromeric chromatin. RNA, as well as secondary chemically modifications that involve both histones and DNA, represents a good candidate to guide somehow the centromere formation and maintenance. Many observations suggest that transcription of centromeric DNA or of other non-coding RNAs could affect centromere formation. To date has been no thorough investigation addressing the identity of the chromatin-associated RNAs (CARs) on a global scale. This prompted us to develop techniques to identify CARs in a genome-wide approach using high-throughput genomic platforms. The future goal of this study will be to focalize the attention on what strictly happens specifically inside centromere chromatin.

Mappaggio fine di un QTL principale per la resa in granella sul cromosoma 3B di frumento duro

In durum wheat, two major QTL for grain yield (Qyld.idw-2B and Qyld.idw-3B) and related traits were identified in a recombinant population derived from Kofa and Svevo (Maccaferri et al. 2008). To further investigate the genetic and physiological basis of allelic variation for this important trait, the fine mapping of Qyld.idw-2B e Qyld.idw-3B was done during the PhD. In this regard, new molecular markers were added to increase the map resolution in the target interval. For Qyld.idw-2B region COS markers derived from the synteny between wheat and rice/ sorghum /brachypodiu genomes were screened. While for Qyld.idw-3B region SSR, ISBP and COS markers obtained from BAC end-sequences and BAC sequences generated during the construction of the 3B physical map (Paux et al., 2008) were screened. In the RIL population a final map resolution of 2,8 markers/cM for Qyld.idw-2B and 0,6 markers/cM for Qyld.idw-3B were obtained. Eighteen pairs of near-isogenic lines (NILs) for Qyld.idw-3B were obtained from F4:5 heterogeneous inbred families. In order to confirm the phenotypic effect of the QTL all pairs were evaluated in field trials (2010 and 2011) for all traits. Three pairs of NILs, with contrasted haplotypes at the target region, were crossed to produce a large F2 population (ca. 7,500 plants in total) that was screened for the identification of recombinants. A total of 233 homozygous F4:5 segmental isolines were obtained and the phenotypic and genotypic characterization of these materials were done. A fine mapping for Qyld.idw-3B was obtained and the QTL peak was identified in a interval of 0,4 cM. All markers were anchored to the Chinese Spring physical map of chr. 3B, which allowed us to identify the BAC Contigs spanning the QTL region and to assign the QTL peak to Contig 954. Sequencing of this contig has revealed the presence of 42 genes.

Virus oncogeni a dna in tumori umani: studio delle infezioni da papillomavirus e poliomavirus in neoplasie di diversa origine.

I virus tumorali inducono oncogenesi nel loro ospite naturale o in sistemi animali sperimentali, manipolando diverse vie cellulari. Ad oggi, sono stati identificati sette virus capaci di causare specifici tumori umani. Inoltre HPV, JCV ed SV40, sono stati associati con un grande numero di tumori umani in sedi corporee non convenzionali, ma, nonostante molti anni di ricerca, nessuna eziologia virale è stata ancora confermata. Lo scopo di questo studio è stato di valutare la presenza ed il significato sia di JCV ed SV40 in tumori ossei umani, e di HPV nel carcinoma della mammella (BC), galattoforectomie (GF), secrezioni mammarie patologiche (ND) e glioblastoma multiforme (GBM). Tecniche di biologia molecolare sono state impiegate per esaminare campioni di tessuto tumorale di 70 tumori ossei (20 osteosarcomi [OS], 20 tumori a cellule giganti [TCG], 30 condrosarcomi [CS]), 168 BCs , 30 GFs, 59 GBM e 30 campioni di ND. Il genoma di SV40 e JCV è stato trovato nel 70% dei CS + 20% degli OS, e nel 13% dei CS +10% dei TCG, rispettivamente. Il DNA di HPV è stato rilevato nel 30% dei pazienti con BC, nel 27% dei campioni GF e nel 13% dei NDs. HPV16 è stato il genotipo maggiormente osservato in tutti questi campioni, seguito da HPV18 e HPV35. Inoltre, il DNA di HPV è stato trovato nel 22% dei pazienti con GBM, in questo tumore HPV6 era il tipo più frequentemente rilevato, seguito da HPV16. L’ ISH ha mostrato che il DNA di HPV è situato all’interno di cellule tumorali mammarie e di GBM. I nostri risultati suggeriscono un possibile ruolo di JCV, SV40 e HPV in questi tumori, se non come induttori come promotori del processo neoplastico, tuttavia diversi criteri devono ancora essere soddisfatti prima di chiarirne il ruolo.

Comparing genetic and linguistic diversity in African populations with a focus on the Khoisan of southern Africa

The interaction between disciplines in the study of human population history is of primary importance, profiting from the biological and cultural characteristics of humankind. In fact, data from genetics, linguistics, archaeology and cultural anthropology can be combined to allow for a broader research perspective. This multidisciplinary approach is here applied to the study of the prehistory of sub-Saharan African populations: in this continent, where Homo sapiens originally started his evolution and diversification, the understanding of the patterns of human variation has a crucial relevance. For this dissertation, molecular data is interpreted and complemented with a major contribution from linguistics: linguistic data are compared to the genetic data and the research questions are contextualized within a linguistic perspective. In the four articles proposed, we analyze Y chromosome SNPs and STRs profiles and full mtDNA genomes on a representative number of samples to investigate key questions of African human variability. Some of these questions address i) the amount of genetic variation on a continental scale and the effects of the widespread migration of Bantu speakers, ii) the extent of ancient population structure, which has been lost in present day populations, iii) the colonization of the southern edge of the continent together with the degree of population contact/replacement, and iv) the prehistory of the diverse Khoisan ethnolinguistic groups, who were traditionally understudied despite representing one of the most ancient divergences of modern human phylogeny. Our results uncover a deep level of genetic structure within the continent and a multilayered pattern of contact between populations. These case studies represent a valuable contribution to the debate on our prehistory and open up further research threads.

Functional and Genetic characterization of new genomic islands from an E. coli strain associated with neonatal meningitis.

Enterobacteriaceae genomes evolve through mutations, rearrangements and horizontal gene transfer (HGT). The latter evolutionary pathway works through the acquisition DNA (GEI) modules of foreign origin that enhances fitness of the host to a given environment. The genome of E. coli IHE3034, a strain isolated from a case of neonatal meningitis, has recently been sequenced and its subsequent sequence analysis has predicted 18 possible GEIs, of which: 8 have not been previously described, 5 fully meet the pathogenic island definition and at least 10 that seem to be of prophagic origin. In order to study the GEI distribution of our reference strain, we screened for the presence 18 GEIs a panel of 132 strains, representative of E. coli diversity. Also, using an inverse nested PCR approach we identified 9 GEI that can form an extrachromosomal circular intermediate (CI) and their respective attachment sites (att). Further, we set up a qPCR approach that allowed us to determine the excision rates of 5 genomic islands in different growth conditions. Four islands, specific for strains appertaining to the sequence type complex 95 (STC95), have been deleted in order to assess their function in a Dictyostelium discoideum grazing assays. Overall, the distribution data presented here indicate that 16 IHE3034 GEIs are more associated to the STC95 strains. Also the functional and genetic characterization has uncovered that GEI 13, 17 and 19 are involved in the resistance to phagocitation by Dictyostelium d thus suggesting a possible role in the adaptation of the pathogen during certain stages of infection.