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.

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.

Evaluation of Microbial Contamination in Bivalve Mollusks: Epidemiology and Diagnosis

Shellfish are filter-feeding organisms that can accumulate many bacteria and viruses. Considering that depuration procedures are not effective in removal of certain microorganisms, shellfish-borne diseases are frequent in many parts of the world, and their control must rely primarily on investigation of prevalence of human pathogens in shellfish and water environment. However, the diffusion of enteric viruses and Vibrio bacteria is not known in many geographical areas, for example in Sardinia, Italy. A survey aimed at investigating the prevalence of Norovirus (NoV), hepatitis A virus (HAV), V. parahaemolyticus, V. cholerae and V. vulnificus was carried out, analyzing both local and imported purified, non-purified and retail shellfish from North Italy and Sardinia. Shellfish from both areas were found contaminated by NoVs, HAV and Vibrio, including retail and purified animals. Molecular analysis evidenced different NoV genogroups and genotypes, including bovine NoVs, as well as pathogenic Vibrio strains, underlining the risk for shellfish consumers. However, also other approaches are needed to control the diffusion of shellfish-borne diseases. It was originally thought that enteric viruses are passively accumulated by shellfish. Recently, it was proven that NoVs bind to specific carbohydrate ligands in oysters, and various NoV strains are characterized by a different bioaccumulation pattern. To deepen the knowledge on this argument, a study was carried out, analyzing bioaccumulation of up to 8 different NoV strains in four different species of shellfish. Different bioaccumulation patterns were observed for each shellfish species and NoV strain used, potentially important in setting up effective shellfish purification protocols. Finally, a novel study of evaluation of viral contamination in shellfish from the French Atlantic coast was carried out following the passage of Xynthia tempest over Western Europe which caused massive destruction. Different enteric viruses were found over a one month period, evidencing the potential of these events of contaminating shellfish.

High-resolution reconstruction of Holocene climate variability and environmental changes in the North Pacific using bivalve shells

Bivalve mollusk shells are useful tools for multi-species and multi-proxy paleoenvironmental reconstructions with a high temporal and spatial resolution. Past environmental conditions can be reconstructed from shell growth and stable oxygen and carbon isotope ratios, which present an archive for temperature, freshwater fluxes and primary productivity. The purpose of this thesis is the reconstruction of Holocene climate and environmental variations in the North Pacific with a high spatial and temporal resolution using marine bivalve shells. This thesis focuses on several different Holocene time periods and multiple regions in the North Pacific, including: Japan, Alaska (AK), British Columbia (BC) and Washington State, which are affected by the monsoon, Pacific Decadal Oscillation (PDO) and El Niño/Southern Oscillation (ENSO). Such high-resolution proxy data from the marine realm of mid- and high-latitudes are still rare. Therefore, this study contributes to the optimization and verification of climate models. However, before using bivalves for environmental reconstructions and seasonality studies, life history traits must be well studied to temporally align and interpret the geochemical record. These calibration studies are essential to ascertain the usefulness of selected bivalve species as paleoclimate proxy archives. This work focuses on two bivalve species, the short-lived Saxidomus gigantea and the long-lived Panopea abrupta. Sclerochronology and oxygen isotope ratios of different shell layers of P. abrupta were studied in order to test the reliability of this species as a climate archive. The annual increments are clearly discernable in umbonal shell portions and the increments widths should be measured in these shell portions. A reliable reconstruction of paleotemperatures may only be achieved by exclusively sampling the outer shell layer of multiple contemporaneous specimens. Life history traits (e.g., timing of growth line formation, duration of the growing season and growth rates) and stable isotope ratios of recent S. gigantea from AK and BC were analyzed in detail. Furthermore, a growth-temperature model based on S. gigantea shells from Alaska was established, which provides a better understanding of the hydrological changes related to the Alaska Coastal Current (ACC). This approach allows the independent measurement of water temperature and salinity from variations in the width of lunar daily growth increments of S. gigantea. Temperature explains 70% of the variability in shell growth. The model was calibrated and tested with modern shells and then applied to archaeological specimens. The time period between 988 and 1447 cal yrs BP was characterized by colder (~1-2°C) and much drier (2-5 PSU) summers, and a likely much slower flowing ACC than at present. In contrast, the summers during the time interval of 599-1014 cal yrs BP were colder (up to 3°C) and fresher (1-2 PSU) than today. The Aleutian Low may have been stronger and the ACC was probably flowing faster during this time.

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.

Calibration study of bivalve shells : implications for paleoenvironmental reconstruction

The accretionary shells of bivalve mollusks can provide environmental information, such as water temperature, precipitation, freshwater fluxes, primary productivity and anthropogenic activities in the form of variable growth rates and variable geochemical properties, such as stable oxygen and carbon isotopes. However, paleoenvironmental reconstructions are constrained by uncertainties about isotopic equilibrium fractionation during shell formation, which is generally acknowledged as a reasonable assumption for bivalves, but it has been disputed in several species. Furthermore, the variation in shell growth rates is accepted to rely on multiple environmental variables, such as temperature, food availability and salinity, but can differ from species to species. Therefore, it is necessary to perform species-specific calibration studies for both isotope proxies and shell growth rates before they can be used with confidence for environmental interpretations of the past. Accordingly, the principal objective of this Ph.D research is to examine the reliability of selected bivalve species, the long-lived Eurhomalea exalbida (Dillwyn), the short-lived and fast growing species Paphia undulata (Born 1778), and the freshwater mussel Margaritifera falcata (Gould 1850), as paleoenvironmental proxy archives.rnThe first part is focused on δ18Oshell and shell growth history of live-collected E. exalbida from the Falkland Islands. The most remarkable finding, however, is that E. exalbida formed its shell with an offset of -0.48‰ to -1.91‰ from the expected oxygen isotopic equilibrium with the ambient water. If this remained unnoticed, paleotemperature estimates would overestimate actual water temperatures by 2.1-8.3°C. With increasing ontogenetic age, the discrepancy between measured and reconstructed temperatures increased exponentially, irrespective of the seasonally varying shell growth rates. This study clearly demonstrates that, when the disequilibrium fractionation effect is taken into account, E. exalbida can serve as a high-resolution paleoclimate archive for the southern South America. The species therefore provides quantifiable temperature estimates, which yields new insights into long-term paleoclimate dynamics for mid to high latitudes on the southern hemisphere.rnThe stable carbon isotope of biogenic carbonates is generally considered to be useful for reconstruction of seawater dissolved inorganic carbon. The δ13Cshell composition of E. exalbida was therefore, investigated in the second part of this study. This chapter focuses on inter-annual and intra-annual variations in δ13Cshell. Environmental records in δ13Cshell are found to be strongly obscured by changes in shell growth rates, even if removing the ontogenetic decreasing trend. This suggests that δ13Cshell in E. exalbida may not be useful as an environmental proxy, but a potential tool for ecological investigations. rnIn addition to long-lived bivalve species, short-lived species that secrete their shells extremely fast, can also be useful for environmental reconstructions, especially as a high-resolution recorder. Therefore, P. undulata from Daya Bay, South China Sea was utilized in Chapter 4 to evaluate and establish a potential proxy archive for past variations of the East Asian monsoon on shorter time-scales. The δ18Oshell can provide qualitative estimates of the amount of monsoonal rain and terrestrial runoff and the δ13Cshell likely reflect the relative amount of isotopically light terrestrial carbon that reaches the ocean during the summer monsoon season. Therefore, shells of P. undulata can provide serviceable proxy archives to reconstruct the frequency of exceptional summer monsoons in the past. The relative strength of monsoon-related precipitation and associated changes in ocean salinity and the δ13C ratios of the dissolved inorganic carbon signature (δ13CDIC) can be estimated from the δ18Oshell and δ13Cshell values as well as shell growth patterns. rnIn the final part, the freshwater pearl shell M. falcata from four rivers in British Columbia, Canada was preliminarily studied concerning the lifespans and the shell growth rates. Two groups separated by the Georgia Strait can be clearly distinguished. Specimens from the western group exhibit a shorter lifespan, while the eastern group live longer. Moreover, the average lifespan seems to decrease from south to north. The computed growth equations from the eastern and western groups differ as well. The western group exhibits a lower growth rate, while bivalves from the eastern group grow faster. The land use history seems to be responsible for the differences in lifespans of the specimens from the two groups. Differences in growth rate may be induced by differences in water temperature or nutrient input also related to the land use activities.

Thermal plasticity of mitochondria: a latitudinal comparison between Southern Ocean molluscs

Mitochondrial volume density (Vv((mt,f))), cristae surface density (Sv((im,mt))), cristae surface area (Sv((im,f))) and citrate synthase (CS) activity were analysed as indicators of thermal acclimation in foot muscle of the limpet, Nacella concinna, and the clam, Laternula elliptica, collected from 4 locations within the Southern Ocean, South Georgia (54 degrees S, N. concinna only), Signy (60 degrees S), Jubany (L. elliptica only -62 degrees S) and Rothera (67 degrees S). Animals were acclimated to 0.0 degrees C whilst a sub-set of N. concinna (South Georgia, Signy and Rothera) and L. elliptica (Rothera) were acclimated to 3.0 degrees C. At 0.0 degrees C N. concinna had higher Vv((mt,f)), Sv((im,mt)), Sv((im,f)) and muscle fibre specific CS activity than L. elliptica which correlated with the more active life style of N. concinna. However, mitochondrial density was very low, 1-2% in both species, suggesting that low temperature compensation of mitochondrial density is not a universal evolutionary response of Antarctic marine ectotherms. Both Sv((im,mt)) and Sv((im,f)) were reduced by warm acclimation of N. concinna. South Georgia N. concinna maintained muscle fibre specific CS activity after acclimation, in contrast to N. concinna from Rothera and Signy and L. elliptica from Rothera, indicating that they have the physiological plasticity to respond to their warmer, more variable thermal environment.

Sediment core fossils in ancient Lake Ohrid: testing for faunal change since the Last Interglacial

Ancient Lake Ohrid is probably of early Pleistocene or Pliocene origin and amongst the few lakes in the world harbouring an outstanding degree of endemic biodiversity. Although there is a long history of evolutionary research in Lake Ohrid, particularly on molluscs, a mollusc fossil record has been missing up to date. For the first time, gastropod and bivalve fossils are reported from the basal, calcareous part of a 2.6 m long sediment succession (core Co1200) from the north-eastern part of Lake Ohrid. Electron spin resonance (ESR) dating of mollusc shells from the same stratigraphic level yielded an age of 130 ± 28 ka. Lithofacies III sediments, i.e. a stratigraphic subdivision comprising the basal succession of core Co1200 between 181.5–263 cm, appeared solid, greyish-white, and consisted almost entirely of silt-sized endogenic calcite (CaCO3>70%) and intact and broken mollusc shells. Here we compare the faunal composition of the thanatocoenosis with recent mollusc associations in Lake Ohrid. A total of 13 mollusc species (9 gastropod and 4 bivalve species) could be identified within Lithofacies III sediments. The value of sediment core fossils for reconstructing palaeoenvironmental settings was evaluated and the agreement between sediment and palaeontological proxies was tested. The study also aims at investigating major faunal changes since the Last Interglacial and searching for signs of extinction events. The combined findings of the ecological study and the sediment characteristics suggest deposition in a shallow water environment during the Last Interglacial. The fossil fauna exclusively included species also found in the present fauna, i.e. no extinction events are evident for this site since the Last Interglacial. The thanatocoenosis showed the highest similarity with recent Intermediate Layer (5–25 m water depth) mollusc assemblages. The demonstrated existence of a mollusc fossil record in Lake Ohrid sediment cores also has great significance for future deep drilling projects. It can be hoped that a more far reaching mollusc fossil record will then be obtained, enabling insight into the early evolutionary history of Lake Ohrid.

An Aquaculture-Based Method for Calibrated Bivalve Isotope Paleothermometry

To quantify species- specific relationships between bivalve carbonate isotope geochemistry ( delta O-18(c)) and water conditions ( temperature and salinity, related to water isotopic composition [delta O-18(w)]), an aquaculture-based methodology was developed and applied to Mytilus edulis ( blue mussel). The four- by- three factorial design consisted of four circulating temperature baths ( 7, 11, 15, and 19 degrees C) and three salinity ranges ( 23, 28, and 32 parts per thousand ( ppt); monitored for delta O-18(w) weekly). In mid- July of 2003, 4800 juvenile mussels were collected in Salt Bay, Damariscotta, Maine, and were placed in each configuration. The size distribution of harvested mussels, based on 105 specimens, ranged from 10.9 mm to 29.5 mm with a mean size of 19.8 mm. The mussels were grown in controlled conditions for up to 8.5 months, and a paleotemperature relationship based on juvenile M. edulis from Maine was developed from animals harvested at months 4, 5, and 8.5. This relationship [ T degrees C = 16.19 (+/- 0.14) - 4.69 (+/- 0.21) {delta O-18(c) VPBD - delta O-18(w) VSMOW} + 0.17 (+/- 0.13) {delta O-18(c) VPBD - delta O-18(w) VSMOW}(2); r(2) = 0.99; N = 105; P < 0.0001] is nearly identical to the Kim and O'Neil ( 1997) abiogenic calcite equation over the entire temperature range ( 7 - 19 degrees C), and it closely resembles the commonly used paleotemperature equations of Epstein et al. ( 1953) and Horibe and Oba ( 1972). Further, the comparison of the M. edulis paleotemperature equation with the Kim and O'Neil ( 1997) equilibrium- based equation indicates that M. edulis specimens used in this study precipitated their shell in isotopic equilibrium with ambient water within the experimental uncertainties of both studies. The aquaculture- based methodology described here allows similar species- specific isotope paleothermometer calibrations to be performed with other bivalve species and thus provides improved quantitative paleoenvironmental reconstructions.

Stable isotopic ratios, shell lengths and organochlorine compounds in four bivalve species collected in July 2009 around Svalbard

Organochlorine compounds (OC) were determined in Arctic bivalves (Mya truncata, Serripes groenlan-dicus, Hiatella arctica and Chlamys islandica) from Svalbard with regard to differences in geographic location, species and variations related to their size and age. Higher chlorinated polychlorinated biphenyls (PCB 101-PCB 194), chlordanes and alpha-hexachlorocyclohexane (alpha-HCH) were consistently detected in the bivalves and PCBs dominated the OC load in the organisms. OC concentrations were highest in Mya truncata and the lowest in Serripes groenlandicus. Species-specific OC levels were likely related to differences in the species' food source, as indicated by the d13C results, rather than size and age. Higher OC concentrations were observed in bivalves from Kongsfjorden compared to the northern sampling locations Liefdefjorden and Sjuoyane. The spatial differences might be related to different water masses influencing Kongsfjorden (Atlantic) and the northern locations (Arctic), with differing phytoplankton bloom situations.