Mitochondrial Inheritance, Mito-nuclear Coevolution, and Sex-associated Genes in Bivalve Molluscs: Transcriptomics and Molecular Evolution

Bivalvia represents an ancient taxon including around 25,000 living species that have adapted to a wide range of environmental conditions, and show a great diversity in body size, shell shapes, and anatomic structure. Bivalves are characterized by highly variable genome sizes and extremely high levels of heterozygosity, which obstacle complete and accurate genome assemblies and hinder further genomic studies. Moreover, some bivalve species presented a stable evolutionary exception to the strictly maternal inheritance of mitochondria, namely doubly uniparental inheritance (DUI), making these species a precious model to study mitochondrial biology. During my PhD, I focused on a DUI species, the Manila clam Ruditapes philippinarum, and my work was two-folded. First, taking advantage of a newly assembled draft genome and a large RNA-seq dataset from different tissues of both sexes, I investigated 1) the role of gene expression and alternative splicing in tissue differentiation; 2) the relationship across tissue specificity, regulatory network connectivity, and sequence evolution; 3) sexual contrasting genetic markers potentially associated with sexual differentiation. The detailed information for this part is in Chapter 2. Second, using the same RNA-seq data, I investigated how nuclear oxidative phosphorylation (OXPHOS) genes coordinate with two divergent mitochondrial genomes in DUI species (mito-nuclear coordination and coevolution). To address this question, I compared transcription, polymorphism, and synonymous codon usage in the mitochondrial and nuclear OXPHOS genes of R. philippinarum in Chapter 3. To my knowledge, this thesis represents the first study exploring the role of alternative splicing in tissue differentiation, and the first study analyzing both transcriptional regulation and sequence evolution to investigate the coordination of OXPHOS genes in bivalves.

A molecular phylogeny of bivalve mollusks: ancient radiations and divergences as revealed by mitochondrial genes

The main scope of my PhD is the reconstruction of the large-scale bivalve phylogeny on the basis of four mitochondrial genes, with samples taken from all major groups of the class. To my knowledge, it is the first attempt of such a breadth in Bivalvia. I decided to focus on both ribosomal and protein coding DNA sequences (two ribosomal encoding genes -12s and 16s -, and two protein coding ones - cytochrome c oxidase I and cytochrome b), since either bibliography and my preliminary results confirmed the importance of combined gene signals in improving evolutionary pathways of the group. Moreover, I wanted to propose a methodological pipeline that proved to be useful to obtain robust results in bivalves phylogeny. Actually, best-performing taxon sampling and alignment strategies were tested, and several data partitioning and molecular evolution models were analyzed, thus demonstrating the importance of molding and implementing non-trivial evolutionary models. In the line of a more rigorous approach to data analysis, I also proposed a new method to assess taxon sampling, by developing Clarke and Warwick statistics: taxon sampling is a major concern in phylogenetic studies, and incomplete, biased, or improper taxon assemblies can lead to misleading results in reconstructing evolutionary trees. Theoretical methods are already available to optimize taxon choice in phylogenetic analyses, but most involve some knowledge about genetic relationships of the group of interest, or even a well-established phylogeny itself; these data are not always available in general phylogenetic applications. The method I proposed measures the "phylogenetic representativeness" of a given sample or set of samples and it is based entirely on the pre-existing available taxonomy of the ingroup, which is commonly known to investigators. Moreover, it also accounts for instability and discordance in taxonomies. A Python-based script suite, called PhyRe, has been developed to implement all analyses.

Caratterizzazione genomica di virus influenzali suini in Italia mediante next generation sequencing

I sottotipi H1N1, H1N2 e H3N2 di influenza A virus sono largamente diffusi nella popolazione suina di tutto il mondo. Nel presente lavoro è stato sviluppato un protocollo di sequenziamento di c.d. nuova generazione, su piattaforma Ion Torrent PGM, idoneo per l’analisi di tutti i virus influenzali suini (SIV). Per valutare l’evoluzione molecolare dei SIV italiani, sono stati sequenziati ed analizzati mediante analisi genomica e filogenetica un totale di sessantadue ceppi di SIV appartenenti ai sottotipi H1N1, H1N2 e H3N2, isolati in Italia dal 1998 al 2014. Sono stati evidenziati in sei campioni due fenomeni di riassortimento: tutti i SIV H1N2 esaminati presentavano una neuraminidasi di derivazione umana, diversa da quella dei SIV H1N2 circolanti in Europa, inoltre l’emoagglutinina (HA) di due isolati H1N2 era originata dal riassortimento con un SIV H1N1 avian-like. L’analisi molecolare dell’HA ha permesso di rivelare un’inserzione di due amminoacidi in quattro SIV H1N1 pandemici e una delezione di due aminoacidi in quattro SIV H1N2, entrambe a livello del sito di legame con il recettore cellulare. E’ stata inoltre evidenziata un’elevata omologia di un SIV H1N1 con ceppi europei isolati negli anni ’80, suggerendo la possibile origine vaccinale di questo virus. E’ stato possibile, in aggiunta, applicare il nuovo protocollo sviluppato per sequenziare un virus influenzale aviare altamente patogeno trasmesso all’uomo, direttamente da campione biologico. La diversità genetica nei SIV esaminati in questo studio conferma l’importanza di un continuo monitoraggio della costellazione genomica dei virus influenzali nella popolazione suina.

Molecular Systematics and Traits Evolution in Phasmatodea Leach, 1815 (Hexapoda; Insecta)

Phasmatodea Leach, 1815 (Hexapoda; Insecta) is a polyneopteran order which counts approximately 3000 described species, often known for their remarkable forms of mimicry. In this thesis, I provide a comprehensive systematic framework which includes over 180 species never considered in a phylogenetic framework: the latter can facilitate a better understanding of the processes underlying phasmids evolutionary history. The clade represents in fact an incredible testing ground to study trait evolution and its striking disparity of reproductive strategies and wing morphologies have been of great interest to the evolutionary biology community. Phasmids wings represent one of the first and most notable rejection of Dollo’s law and they played a central role in initiating a long- standing debate on the irreversibility of complex traits loss. Macroevolutionary analyses presented here confirm that wings evolution in phasmids is a reversible process even when possible biases - such as systematic uncertainty and trait-dependent diversification rates - are considered. These findings remark how complex traits can evolve in a dynamic, reversible manner and imply that their molecular groundplan can be preserved despite its phenotypical absence. This concept has been further tested with phylogenetic and transcriptomic approaches in two phasmids parthenogenetic lineages and a bisexual congeneric of the European Bacillus species complex. Leveraging a gene co-expression network approach, male gonad associated genes were retrieved in the bisexual species and then their modifications in the parthenogens were charachterized. Pleiotropy appears to constrain gene modifications associated to male reproductive structures after their loss in parthenogens, so that the lost trait molecular groundplan can be largely preserved in both transcription patterns and sequence evolution. Overall, the results presented in this thesis contribute to shape our understanding of the interplay between the phenotypic and molecular levels in trait evolution.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Effects of Global Warming on Berry Composition of cv. Sangiovese: Biochemical and Molecular Aspects and Agronomical Adaptation Approaches

Wine grape must deal with serious problems due to the unfavorable climatic conditions resulted from global warming. High temperatures result in oxidative damages to grape vines. The excessive elevated temperatures are critical for grapevine productivity and survival and contribute to degradation of grape and wine quality and yield. Elevated temperature can negatively affect anthocyanin accumulation in red grape. Particularly, cv. Sangiovese was identified to be very sensitive to such condition. The quantitative real-time PCR analysis showed that flavonoid biosynthetic genes were slightly repressed by high temperature. Also, the heat stress repressed the expression of the transcription factor “VvMYBA1” that activates the expression of UFGT. Moreover, high temperatures had repressing effects on the activity of the flavonoids biosynthetic enzymes “PAL” and “UFGT”.Anthocyanin accumulation in berry skin is due to the balance between its synthesis and oxidation. In grape cv. Sangiovese, the gene transcription and activity of peroxidases enzyme was elevated by heat stress as a defensive mechanism of ROS-scavenging. Among many isoforms of peroxidases genes, one gene (POD 1) was induced in Sangiovese under thermal stress condition. This gene was isolated and evaluated via the technique of genes transformation from grape to Petunia. Reduction in anthocyanins concentration and higher enzymatic activity of peroxidase was observed in POD 1 transformed Petunia after heat shock compared to untrasformed control. Moreover, in wine producing regions, it is inevitable for the grape growers to adopt some adaptive strategies to alleviate grape damages to abiotic stresses. Therefore, in this thesis, the technique of post veraison trimming was done to improve the coupling of phenolic and sugar ripening in Vitis vinifera L. cultivar Sangiovese. Trimming after veraison showed to be executable to slow down the rate of sugar accumulation in grape (to decrease the alcohol potential in wines) without evolution of the main berry flavonoids compounds.

The evolution of massive clumps in star forming regions

In this thesis two related arguments are investigated: - The first stages of the process of massive star formation, investigating the physical conditions and -properties of massive clumps in different evolutionary stages, and their CO depletion; - The influence that high-mass stars have on the nearby material and on the activity of star formation. I characterise the gas and dust temperature, mass and density of a sample of massive clumps, and analyse the variation of these properties from quiescent clumps, without any sign of active star formation, to clumps likely hosting a zero-age main sequence star. I briefly discuss CO depletion and recent observations of several molecular species, tracers of Hot Cores and/or shocked gas, of a subsample of these clumps. The issue of CO depletion is addressed in more detail in a larger sample consisting of the brightest sources in the ATLASGAL survey: using a radiative tranfer code I investigate how the depletion changes from dark clouds to more evolved objects, and compare its evolution to what happens in the low-mass regime. Finally, I derive the physical properties of the molecular gas in the photon-dominated region adjacent to the HII region G353.2+0.9 in the vicinity of Pismis 24, a young, massive cluster, containing some of the most massive and hottest stars known in our Galaxy. I derive the IMF of the cluster and study the star formation activity in its surroundings. Much of the data analysis is done with a Bayesian approach. Therefore, a separate chapter is dedicated to the concepts of Bayesian statistics.

Neurodegenerative diseases: molecular mechanisms and new strategies for neuroprotection

The term neurodegeneration defines numerous conditions that modify neuron’s normal functions in the human brain where is possible to observe a progressive and consistent neuronal loss. The mechanisms involved in neurodegenerative chronic and acute diseases evolution are not completely understood yet, however they share common characteristics such as misfolded proteins, oxidative stress, inflammation, excitotoxicity, and neuronal loss. Many studies have shown the frequency to develop neurodegenerative chronic diseases several years after an acute brain injury. In addition, many patients show, after a traumatic brain injury, motor and cognitive manifestations that are close to which are observed in neurodegenerative chronic patients. For this reason it is evident how is fundamental the concept of neuroprotection as a way to modulate the neurodegenerative processes evolution. Neuroinflammation, oxidative stress and the apoptotic process may be functional targets where operate to this end. Taking into account these considerations, the aim of the present study is to identify potential common pathogenetic pathways in neurodegenerative diseases using an integrated approach of preclinical studies. The goal is to delineate therapeutic strategies for the prevention of neuroinflammation, neurodegeneration and dysfunctions associated to Parkinson’s disease (PD) and cerebral ischemia. In the present study we used a murine model of PD treated with an isothiocyanate, 6-MSITC, able to quench ROS formation, restore the antioxidant GSH system, slow down the apoptotic neuronal death and counteract motor dysfunction induced by 6-OHDA. In the second study we utilized a transgenic mouse model knockout for CD36 receptor to investigate the inflammation involvement in a long term study of MCAo, which shows a better outcome after the damage induced. In conclusion, results in this study allow underlying the connection among these pathologies, and the importance of a neuroprotective strategy able to restore neurons activity where current drugs therapies have shown palliative but not healing abilities.

Unveiling the molecular mechanism of BRCA2-RAD51 interaction to tackle cancer onset

Different kinds of lesions can occur to DNA, and among them, one of the most dangerous is the double strand breaks (DSBs). Actually, DSBs can result in mutations, chromosome translocation or deletion. For this kind of lesions, depending on cell cycle phase as well as DNA-end resection, cells have developed specific repair pathways. Among these the error-free homologous recombination (HR) plays a crucial role. HR takes place during S/G2 phases, since the sister chromatids can be used as homologous templates. In this process, hRAD51 and BRCA2 are key players. hRAD51 is a recombinase of 339 amino-acids highly conserved through evolution which displays an intrinsic tendency to form oligomeric structures. BRCA2 is a very large protein of 3418 amino-acids, essential for the recruitment and accumulation of hRAD51 in the nucleus repairing-foci. BRCA2 interacts with hRAD51 through eight, so-called, BRC repeats, composed of 35-40 amino-acids. Mutations within this region have been linked to an increased risk of ovarian cancer development. In particular, several reports highlighted that missense mutations within one BRC repeat can hamper BRCA2 activity. Considering the close homology between the BRC repeats, it is striking how these mutations cannot be counterbalanced by the other non-mutated repeats preserving the function and the interactions of BRCA2 with hRAD51. To date the only interaction that has been structurally elucidated, is the one taking place amid the fourth BRC repeat and hRAD51. Only very little biophysical information is available on the interaction of the other BRC repeats with hRAD51. This thesis aims at elucidating the mechanism of hRAD51-BRCA2 interaction, by means of biophysical and structural approaches.

Spectroscopic investigation on photoreactivity, structure and polymorphism of organic molecular crystals

The study of the spectroscopic phenomena in organic solids, in combination with other techniques, is an effective tool for the understanding of the structural properties of materials based on these compounds. This Ph.D. work was dedicated to the spectroscopic investigation of some relevant processes occurring in organic molecular crystals, with the goal of expanding the knowledge on the relationship between structure, dynamics and photoreactivity of these systems. Vibrational spectroscopy has been the technique of choice, always in combination with X-ray diffraction structural studies and often the support of computational methods. The vibrational study of the molecular solid state reaches its full potential when it includes the low-wavenumber region of the lattice-phonon modes, which probe the weak intermolecular interactions and are the fingerprints of the lattice itself. Microscopy is an invaluable addition in the investigation of processes that take place in the micro-meter scale of the crystal micro-domains. In chemical and phase transitions, as well as in polymorph screening and identification, the combination of Raman microscopy and lattice-phonon detection has provided useful information. Research on the fascinating class of single-crystal-to-single-crystal photoreactions, has shown how the homogeneous mechanism of these transformations can be identified by lattice-phonon microscopy, in agreement with the continuous evolution of their XRD patterns. On describing the behavior of the photodimerization mechanism of vitamin K3, the focus was instead on the influence of its polymorphism in governing the product isomerism. Polymorphism is the additional degree of freedom of molecular functional materials, and by advancing in its control and properties, functionalities can be promoted for useful applications. Its investigation focused on thin-film phases, widely employed in organic electronics. The ambiguities in phase identification often emerging by other experimental methods were successfully solved by vibrational measurements.

Molecular and hybrid systems for solar-to-chemical energy conversion

My Ph.D. thesis was dedicated to the exploration of different paths to convert sunlight into the shape of chemical bonds, by the formation of solar fuels. During the past three years, I have focused my research on two of these, namely molecular hydrogen H2 and the reduced nicotinamide adenine dinucleotide enzyme cofactor NAD(P)H. The first could become the ideal energy carrier for a truly clean energy system; it currently represents the best chance to liberate humanity from its dependence on fossil fuels. To address this, I studied different systems which can achieve proton reduction upon light absorption. More specifically, part of my work was aimed to the development of a cost-effective and stable catalyst in combination with a well-known photochemical cycle. To this extent, I worked on transition metal oxides which, as demonstrated in this work, have been identified as promising H2 evolution catalysts, showing excellent activity, stability, and previously unreported versatility. Another branch of my work on hydrogen production dealt with the use of a new class of polymeric semiconductor materials to absorb light and convert it into H2. The second solar fuel mentioned above is a key component of the most powerful methods for chemical synthesis: enzyme catalysis. The high cost of the reduced forms prohibits large-scale utilization, so artificial photosynthetic approaches for regenerating it are being intensively studied. The first system I developed exploits the tremendous reducing properties of a scarcely known ruthenium complex which is able to reduce NAD+. Lastly, I sought to revert the classical role of the sacrificial electron donor to an active component of the system and, to boost the process, I build up an autonomous microfluidic system able to generate highly reproducible NAD(P)H amount, demonstrating the superior performance of microfluidic reactors over batch and representing another successful photochemical NAD(P)H regeneration system.