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.

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.

Down Syndrome: Neuropsychological phenotype and mitochondrial DNA

Introduction. Down Syndrome (DS) is the most known autosomal trisomy, due to the presence in three copies of chromosome 21. Many studies were designed to identify phenotypic and clinical consequences related to the triple gene dosage. However, the general conclusion is a senescent phenotype; in particular, the most features of physiological aging, such as skin and hair changes, vision and hearing impairments, thyroid dysfunction, Alzheimer-like dementia, congenital heart defects, gastrointestinal malformations, immune system changes, appear in DS earlier than in normal age-matched subjects. The only established risk factor for the DS is advanced maternal age, responsible for changes in the meiosis of oocytes, in particular the meiotic nondisjunction of chromosome 21. In this process mitochondria play an important role since mitochondrial dysfunction, due to a variety of extrinsic and intrinsic influences, can profoundly influence the level of ATP generation in oocytes, required for a correct chromosomal segregation. Aim. The aim of this study is to investigate an integrated set of molecular genetic parameters (sequencing of complete mtDNA, heteroplasmy of the mtDNA control region, genotypes of APOE gene) in order to identify a possible association with the early neurocognitive decline observed in DS. Results. MtDNA point mutations do not accumulate with age in our study sample and do not correlate with early neurocognitive decline of DS subjects. It seems that D-loop heteroplasmy is largely not inherited and tends to accumulate somatically. Furthermore, in our study sample no association of cognitive impairment and ApoE genotype is found. Conclusions. Overall, our data cast some doubts on the involvement of these mutations in the decline of cognitive functions observed in DS.

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.

HIF1α regulates Mitochondrial biogenesis and Cellular senescence induced by Gamma Radiation

Cellular response to γ-rays is mediated by ATM-p53 axis. When p53 is phosphorylated, it can transactivate several genes to induce permanent cell cycle arrest (senescence) or apoptosis. Epithelial and mesenchymal cells are more resistant to radiation-induced apoptosis and respond mainly by activating senescence. Hence, tumor cells in a senescent state might remain as “dormant” malignant in fact through disruption of p53 function, cells may overcome growth arrest. Oncocytic features were acquired in the recurring neoplasia after radiation therapy in patient with colonrectal cancer. Oncocytic tumors are characterized by aberrant biogenesis and are mainly non-aggressive neoplasms. Their low proliferation degree can be explained by chronic destabilization of HIF1α, which presides to adaptation to hypoxia and also plays a pivotal role in hypoxia-related radio-resistance. The aim of the present thesis was to verify whether mitochondrial biogenesis can be induced following radiation treatment, in relation of HIF1α status and whether is predictive of a senescence response. In this study was demonstrate that mitochondrial biogenesis parameters like mitochondrial DNA copy number could be used for the prediction of hypoxic status of tissue after radiation treatment. γ-rays induce an increase of mitochondrial mass and function, in response to a genotoxic stress that pushes cells into senescence. Mitochondrial biogenesis is only indirectly regulated by p53, whose activation triggers a MDM2-mediated HIF1α degradation, leading to the release of PGC-1β inhibition by HIF1α. On the other hand, this protein blunts the mitochondrial response to γ-rays as well as the induction of p21-mediated cell senescence, indicating prevalence of the hypoxic over the genotoxic response. Finally in vivo, post-radiotherapy mtDNA copy number increase well correlates with lack of HIF1α increase in the tissue, concluding this may be a useful molecular tool to infer the trigger of a hypoxic response during radiotherapy, which may lead to failure of activation of senescence.

Mitochondrial respiratory supercomplex association limits production of reactive oxygen species from Complex I

Evidence accumulated in the last ten years has demonstrated that a large proportion of the mitochondrial respiratory chain complexes in a variety of organisms is arranged in supramolecular assemblies called supercomplexes or respirasomes. Besides conferring a kinetic advantage (substrate channeling) and being required for the assembly and stability of Complex I, indirect considerations support the view that supercomplexes may also prevent excessive formation of reactive oxygen species (ROS) from the respiratory chain. Following this line of thought we have decided to directly investigate ROS production by Complex I under conditions in which the complex is arranged as a component of the supercomplex I1III2 or it is dissociated as an individual enzyme. The study has been addressed both in bovine heart mitochondrial membranes and in reconstituted proteoliposomes composed of complexes I and III in which the supramolecular organization of the respiratory assemblies is impaired by: (i) treatment either of bovine heart mitochondria or liposome-reconstituted supercomplex I-III with dodecyl maltoside; (ii) reconstitution of Complexes I and III at high phospholipids to protein ratio. The results of this investigation provide experimental evidence that the production of ROS is strongly increased in either model; supporting the view that disruption or prevention of the association between Complex I and Complex III by different means enhances the generation of superoxide from Complex I . This is the first demonstration that dissociation of the supercomplex I1III2 in the mitochondrial membrane is a cause of oxidative stress from Complex I. Previous work in our laboratory demonstrated that lipid peroxidation can dissociate the supramolecular assemblies; thus, here we confirm that preliminary conclusion that primary causes of oxidative stress may perpetuate reactive oxygen species (ROS) generation by a vicious circle involving supercomplex dissociation as a major determinant.

Role of Magnesium and its mitochondrial transporter MRS2 in the modulation of drug-induced apoptosis leading to multidrug resistance phenotype

Magnesium is an essential element for many biological processes crucial for cell life and proliferation. Growing evidences point out a role for this cation in the apoptotic process and in developing multi drug resistance (MDR) phenotype. The first part of this study aimed to highlight the involvement of the mitochondrial magnesium channel MRS2 in modulating drug-induced apoptosis. We generated an appropriate transgenic cellular system to regulate expression of MRS2 protein. The cells were then exposed to two different apoptotic agents commonly used in chemotherapy. The obtained results showed that cells overexpressing MRS2 channel are less responsiveness to pharmacological insults, looking more resistant to the induced apoptosis. Moreover, in normal condition, MRS2 overexpression induces higher magnesium uptake into isolated mitochondria respect to control cells correlating with an increment of total intracellular magnesium concentration. In the second part of this research we investigated whether magnesium intracellular content and compartmentalization could be used as a signature to discriminate MDR tumour cells from their sensitive counterparts. As MDR model we choose colon carcinoma cell line sensitive and resistant to doxorubicin. We exploited a standard-less approach providing a complete characterization of whole single-cells by combining X-Ray Fluorescence Microscopy , Atomic Force Microscopy and Scanning Transmission X-ray Microscopy. This method allows the quantification of the intracellular spatial distribution and total concentration of magnesium in whole dehydrated cells. The measurements, carried out in 27 single cells, revealed a different magnesium pattern for both concentration and distribution of the element in the two cellular strains. These results were then confirmed by quantifying the total amount of intracellular magnesium in a large populations of cells by using DCHQ5 probe and traditional fluorimetric technique.

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.

Shikonin directly targets mitochondria and causes mitochondrial dysfunction in cancer cells

Chemotherapy is a mainstay of cancer treatment. Due to increased drug resistance and the severe side effects of currently used therapeutics, new candidate compounds are required for improvement of therapy success. Shikonin, a natural naphthoquinone, was used in traditional Chinese medicine for the treatment of different inflammatory diseases and recent studies revealed the anticancer activities of shikonin. We found that shikonin has strong cytotoxic effects on 15 cancer cell lines, including multidrug-resistant cell lines. Transcriptome-wide mRNA expression studies showed that shikonin induced genetic pathways regulating cell cycle, mitochondrial function, levels of reactive oxygen species, and cytoskeletal formation. Taking advantage of the inherent fluorescence of shikonin, we analyzed its uptake and distribution in live cells with high spatial and temporal resolution using flow cytometry and confocal microscopy. Shikonin was specifically accumulated in the mitochondria, and this accumulation was associated with a shikonin-dependent deregulation of cellular Ca(2+) and ROS levels. This deregulation led to a breakdown of the mitochondrial membrane potential, dysfunction of microtubules, cell-cycle arrest, and ultimately induction of apoptosis. Seeing as both the metabolism and the structure of mitochondria show marked differences between cancer cells and normal cells, shikonin is a promising candidate for the next generation of chemotherapy.

OPA1 isoforms and protein domains in the rescue of mitochondrial dysfunctions

Mutations in OPA1 gene have been identified in the majority of patients with Dominant Optic Atrophy (DOA), a blinding disease, and the syndromic form DOA-plus. OPA1 protein is a mitochondrial GTPase involved in various mitochondrial functions, present in humans in eight isoforms, resulting from alternative splicing and proteolytic processing. In this study we have investigated the specific role of each isoform through expression in OPA-/- MEFs, by evaluating their ability to improve the defective mitochondrial phenotypes. All isoforms were able to rescue the energetic efficiency, mitochondrial DNA (mtDNA) content and cristae integrity, but only the presence of both long and short forms could recover the mitochondrial morphology. In order to identify the OPA1 protein domains crucial for its functions, we selected and modified the isoform 1, shown to be one of the most efficient in preserving mitochondrial phenotype, to express three specific OPA1 variants, namely: one with a different N-terminus portion, one unable to generate short form owing to deletion of S1 cleavage site and one with a defective GTPase domain. We demonstrated that the simultaneous presence of the N- and C-terminus of OPA1 was essential for the mtDNA maintenance; a cleavable isoform generating s-forms was necessary to completely rescue the energetic competence and the presence of the C-terminus was sufficient to partially recover the cristae ultrastructure. Lastly, several pathogenic OPA1 mutations were inserted in MEF clones and the biochemical features investigated, to correlate the defective phenotypes with the clinical severity of patients. Our results clearly indicate that this cell model reflects very well the clinical characteristics of the patients, and therefore can be proposed as an useful tool to shed light on the pathomechanism underlying DOA.

Inference of Raja miraletus population structure using mitochondrial and nuclear DNA: comparing the resolution power of molecular markers in exploring species and population boundaries

This study focused on the role of oceanographic discontinuities and the presence of transitional areas in shaping the population structure and the phylogeography of the Raja miraletus species complex, coupled with the test of the effective occurrence of past speciation events. The comparisons between the Atlantic African and the North-Eastern Atlantic-Mediterranean geographic populations were unravelled using both Cytochrome Oxidase I and eight microsatellite loci. This approach guaranteed a robust dataset for the identification of a speciation event between the Atlantic African clade, corresponding to the ex Raja ocellifera nominal species, and the NE Atlantic-Mediterranean R. miraletus clade. As a matter of fact, the origin of the Atlantic Africa and the NE Atlantic-Mediterranean deep split dated about 11.74MYA and was likely due to the synergic influence currents and two upwelling areas crossing the Western African Waters. Within the Mediterranean Sea, particular attention was also paid to the transitional area represented by Adventura and Maltese Bank, that might have contributed in sustaining the connectivity of the Western and the Eastern Mediterranean geographical populations. Furthermore, the geology of the easternmost part of Sicily and the geo-morphological depression of the Calabrian Arc could have driven the differentiation of the Eastern Mediterranean Sea. Although bathymetric and oceanographic discontinuity could represent barriers to dispersal and migration between Eastern and Western Mediterranean samples, a clear and complete genetic separation among them was not detected. Results produced by this work identified a speciation event defining Raja ocellifera and R. miraletus as two different species, and describing the R. miraletus species complex as the most ancient cryptic speciation event in the family Rajidae, representing another example of how strictly connected the environment, the behavioural habits and the evolutionary and ecologic drivers are.

First assessment on genetic structure and phylogeography of Mediterranean blue shark (prionace glauca, L. 1758) population using mitochondrial gene variation: a comparison with the Atlantic.

The blue shark, Prionace glauca, is one of the most vagile shark species worldwide distributed. The particular body shape allows blue sharks make transoceanic movements, leading to a circumglobal distribution. Due to its reproductive cycle, an extraordinarily high number of specimens is globally registered but, even if it is still a major bycatch of longline fishery rather than a commercial target, it is characterized by a high vulnerability. In this perspective it is important to increase the amount of informations regarding its population extent in the different worldwide areas, evaluating the possible phylogeographic patterns between different locations. This study, included in the "MedBlueSGen" European project, aims exactly at filling a gap in knowledges regarding the genetic population structure of the Mediterranean blue sharks, which has never been investigated before, with a comparison with the North-Eastern Atlantic blue shark population. To reach this objective, we used a dataset of samples from different Mediterranean areas implementing it with some samples from North-Eastern Atlantic. Analyzing the variability of the two mitochondrial markers control region and cytochrome b, with the design of new species-specific primer pairs, we assessed the mitochondrial genetic structure of Mediterranean and North-Eastern Atlantic samples, focusing on the analysis of their possible connectivity, and we tried to reconstruct their demographic history and population size. Data analyses highlighted the absence of a genetic structuring within the Mediterranean and among it and North-Eastern Atlantic, suggesting that the Strait of Gibraltar doesn't represent a phylogeographic barrier. These results are coherent to what has been found in similar investigations on other worldwide blue shark populations. Analysis of the historical demographic trend revealed a general stable pattern for the cytochrome-b and a slightly population expansion for the control region marker.

Mitochondrial tRNA(Leu(UUR)) mutation m.3302A > G presenting as childhood-onset severe myopathy: threshold determination through segregation study

Mitochondrial tRNA(Leu(UUR)) mutation m.3302A > G is associated with respiratory chain complex I deficiency and has been described as a rare cause of mostly adult-onset slowly progressive myopathy. Five families with 11 patients have been described so far; 5 of them died young due to cardiorespiratory failure. Here, we report on a segregation study in a family with an index patient who already presented at the age of 18 months with proximal muscular hypotonia, abnormal fatigability, and lactic acidosis. This early-onset myopathy was rapidly progressive. At 8 years, the patient is wheel-chair bound, requires nocturnal assisted ventilation, and suffers from recurrent respiratory infections. Severe complex I deficiency and nearly homoplasmy for m.3302A > G were found in muscle. We collected blood, hair, buccal swabs and muscle biopsies from asymptomatic adults in this pedigree and determined heteroplasmy levels in these tissues as well as OXPHOS activities in muscle. All participating asymptomatic adults had normal OXPHOS activities. In contrast to earlier reports, we found surprisingly little variation of heteroplasmy levels in different tissues of the same individual. Up to 45% mutation load in muscle and up to 38% mutation load in other tissues were found in non-affected adults. The phenotypic spectrum of tRNA(Leu(UUR)) m.3302A > G mutation seems to be wider than previously described. A threshold of more than 45% heteroplasmy in muscle seems to be necessary to alter complex I activity leading to clinical manifestation. The presented data may be helpful for prognostic considerations and counseling in affected families.