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.

Organismi acquatici e ambiente: meccanismi biochimici di interazione, risposta e adattamento

The research is focused on the relationship between some Mg2+-dependent ATPase activities of plasma- and mitochondrial membranes from tissues of cultured marine bivalve molluscs and potentially stressful environmental conditions, such as the exposure to contaminants both of natural origin (ammonia nitrogen, the main contaminant of aquaculture plants) and of anthropic source (alkyltins). The two filter-feeding bivalve species selected colonize different habitats: the common mussel Mytilus galloprovincialis binds to hard substrates and the Philippine clam Tapes philippinarum burrows into sea bottom sandy beds. The choice of typical species of coastal waters, extremely suitable for environmental studies due to their features of poor motility, resistance to transport and great filtering efficiency, may constitute a model to evaluate responses to contaminants of membrane-bound enzyme activities involved in key biochemical mechanisms, namely cell ionic regulation and mitochondrial energy production. In vitro and in vitro approaches have been pursued. In vitro assays were carried out by adding the contaminants (NH4Cl and alkyltins) directly to the ATPase reaction media. In vivo experiments were carried out by exposing mussels to various tributyl tin (TBT) concentrations under controlled conditions in aquaria. ATPase activities were determined spectrophotometrically according to the principles of the method of Fiske and Subbarow (1925). The main results obtained are detailed below. In Tapes philippinarum the interaction of NH4 +, the main form of ammonia nitrogen at physiological and seawater pHs, with the Na,K-ATPase and the ouabaininsensitive Na-ATPase was investigated in vitro on gill and mantle microsomal membranes. The proven replacement by NH4 +of K+ in the activation of the Na,KATPase and of Na+ in the activation of the ouabain-insensitive ATPase displayed similar enzyme affinity for the substituted cation. on the one hand this finding may represent one of the possible mechanisms of ammonia toxicity and, on the other, it supports the hypothesis that NH4 + can be transported across the plasma membrane through the two ATPases. In this case both microsomal ATPases may be involved and co-operate, at least under peculiar circumstances, to nitrogen excretion and ammonia detoxification mechanisms in bivalve molluscs. The two ATPase activities stimulated by NH4 + maintained their typical response to the glycoside ouabain, specific inhibitor of the Na,K-ATPase, being the Na++ NH4 +-activated ATPase even more susceptive to the inhibitor and the ouabain-insensitive ATPase activity activated indifferently by Na+ or NH4 + unaffected by up to 10-2 M ouabain. In vitro assays were carried out to evaluate the response of the two Na-dependent ATPases to organotins in clams and mussels and to investigate the interaction of TBT with mussel mitochondrial oligomycin-sensitive Mg-ATPase. Since no literature data were available, the optimal assay conditions and oligomycin sensitivity of mussel mitochondrial MgATPase were determined. In T. philippinarum the ouabain-insensitive Na-ATPase was found to be refractory to TBT both in the gills and in the mantle, whereas the Na,K-ATPase was progressively inhibited by increasing TBT doses; the enzyme inhibition was more pronounced in the gills than in the mantle. In both tissues of M. galloprovincialis the Na,K-ATPase inhibition by alkyltins decreased in the order TBT>DBT(dibutyltin)>>MBT(monobutyltin)=TeET(tetraethyltin) (no effect). Mussel Na-ATPase confirmed its refractorimess to TBT and derivatives both in the gills and in the mantle. These results indicate that the Na,K-ATPase inhibition decreases as the number of alkyl chains bound to tin decreases; however a certain polarity of the organotin molecule is required to yield Na,K-ATPase inhibition, since no enzyme inhibition occurred in the presence of tetraalkyl-substituted derivatives such as TeET . Assays carried out in the presence of the dithioerythritol (DTE) pointed out that the sulphhydrylic agent is capable to prevent the Na,K-ATPase inhibition by TBT, thus suggesting that the inhibitor may link to -SH groups of the enzyme complex.. Finally, the different effect of alkyltins on the two Na-dependent ATPases may constitute a further tool to differentiate between the two enzyme activities. These results add to the wealth of literature data describing different responses of the two enzyme activities to endogenous and exogenous modulators . Mussel mitochondrial Mg-ATPase was also found to be in vitro inhibited by TBT both in the gills and in the mantle: the enzyme inhibition followed non competitive kinetics. The failed effect of DTE pointed out that in this case the interaction of TBT with the enzyme complex is probably different from that with the Na,K-ATPase. The results are consistent with literature data showing that alkyltin may interact with enzyme structures with different mechanisms. Mussel exposure to different TBT sublethal doses in aquaria was carried out for 120 hours. Two samplings (after 24 and 120 hrs) were performed in order to evaluate a short-term response of gill and mantle Na,K-ATPase, ouabain-insensitive Na-ATPase and Mg-ATPase activities. The in vivo response to the contaminants of the enzyme activities under study was shown to be partially different from that pointed out in the in vitro assays. Mitochondrial Mg-ATPase activity appeared to be activated in TBTexposed mussels with respect to control ones, thus confirming the complexity of evaluating in vivo responses of the enzyme activities to contaminants, due to possible interactions of toxicants with molluscan metabolism. Concluding, the whole of data point out that microsomal and mitochondrial ATPase activities of bivalve molluscs are generally responsive to environmental contaminants and suggest that in some cases membrane-bound enzyme activities may represent the molecular target of their toxicity. Since the Na,K-ATPase, the Na-ATPase and the Mg-ATPase activities are poorly studied in marine bivalves, this research may contribute to enlarge knowledge in this quite unexplored field.

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.

Mitochondrial inheritance and germ line determination in a bivalve species with Doubly Uniparental Inheritance (DUI)

Mitochondria are inherited maternally in most metazoans. However, in some bivalves, two mitochondrial lineages are present: one transmitted through eggs (F), the other through sperm (M). This is called Doubly Uniparental Inheritance (DUI). During male embryo development, spermatozoon mitochondria aggregate and end up in the primordial germ cells, while they are dispersed in female embryos. The molecular mechanisms of segregation patterns are still unknown. In the DUI species Ruditapes philippinarum, I examined sperm mitochondria distribution by MitoTracker, microtubule staining and TEM, and I localized germ line determinants with immunocytochemical analysis. I also analyzed the gonad transcriptome, searching for genes involved in reproduction and sex determination. Moreover, I analyzed an M-type specific open reading frame that could be responsible for maintenance/degradation of M mitochondria during embryo development. These transcripts were also localized in tissues using in situ hybridization. As in Mytilus, two distribution patterns of M mitochondria were detected in R. philippinarum, supporting that they are related to DUI. Moreover, the first division midbody concurs in positioning aggregated M mitochondria on the animal-vegetal axis of the male embryo: in organisms with spiral segmentation this zone is not involved in further cleavages, so aggregation is maintained. Moreover, sperm mitochondria reach the same embryonic area where germ plasm is transferred, suggesting their contribution in male germ line formation. The finding of reproduction and ubiquitination transcripts led to formulate a model in which ubiquitination genes stored in female oocytes during gametogenesis would activate sex-gene expression in the early embryonic developmental stages (preformation). Only gametogenetic cells were labeled by in situ hybridization, proving their specific transcription in developing gametes. Other than having a role in sex determination, some ubiquination factors could also be involved in mitochondrial inheritance, and their differential expression could be responsible for the different fate of sperm mitochondria in the two sexes.

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.

Molecular mechanisms controlling physiological plasticity in marine mussels under the influence of natural and anthropogenic stress factors

Marine mussels are exceptionally well-adapted to live in transitional habitats where they are exposed to fluctuating environmental parameters and elevated levels of natural and anthropogenic stressors throughout their lifecycle. However, there is a dearth of information about the molecular mechanisms that assist in dealing with environmental changes. This project aims to investigate the molecular mechanisms governing acclimatory and stress responses of the Mediterranean mussel (Mytilus galloprovincialis) by addressing relevant life stages and environmental stressors of emerging concern. The experimental approach consisted of two phases to explore (i) the physiological processes at early life history and the consequences of plastic pollution and (ii) the adult physiology processes under natural habitats. As the first phase, I employed a plastic leachate (styrene monomer), and polystyrene microplastics to understand the modulation of cytoprotective mechanisms during the early embryo stages. Results revealed the onset of transcriptional impairments of genes involved in MXR-related transporters and other physiological processes induced by styrene and PS-MPs. In the second phase, as a preliminary analysis, microbiota profile of adult mussels at the tissue scale and its surrounding water was explored to understand microbiota structures that may reflect peculiar adaptations to the respective tissue functions. The broader experiment has been implemented to understand the variability of transcriptional profiles in the mussel digestive glands in the natural setting. All the genes employed in this study have shown possibilities to use as molecular biomarker responses throughout the year for monitoring the physiology of mussels living in a particular environment and, in turn, more properly detecting changes in the environment. As a whole, my studies provide insights into the interactions between environmental parameters, and intrinsic characters, and physiology of marine bivalves, and it could help to interpretation of responses correctly under stress conditions and climate change scenarios.