Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.

The genetic history of Italians: new-insights from uniparentally-inherited markers

Genetic differences among human groups can be ascribed both to the broad-scale extents of pre-historical and historical migrations and to the fine-scale impacts of socio-cultural and geographic heterogeneity. In this thesis, the genetic information provided by uniparental markers were exploited to address different aspects of the Italian population history, by combining macro- and micro-geographic investigations at different spatial and temporal scales. To firstly assess the overall Italian variability, Y-chromosome and mtDNA markers were deeply typed in ~900 individuals from continental Italy, Sicily and Sardinia. Sex-biased patterns and contrasting demographic histories were observed for males and females. Differential European and Mediterranean contributions were invoked to explain the paternal genetic sub-structure observed in peninsular Italy, compared to the homogeneous maternal genetic landscape. If Neolithic showed to be one principal determinant of the detected paternal structure, local insights into specific Italian regional contexts highlighted the importance of Post-Neolithic contributions. Among them, migrations from the Balkans (particularly Greece) during late Metal Ages, played a relevant role in the cultural and genetic transitions occurred in Sicily and Southern Italy. On a finer geographic and temporal perspective, the more recent layers of Italian genetic history and some aspects of the gene-culture interaction were assessed by exploring the genetic variability within two “marginal populations”: Arbereshe of Southern Italy and Partecipanza in Northern Italy. The Arbereshe are Albanian-speaking communities settled in Sicily and Calabria since the end of Middle Ages. Despite sharing common genetic and cultural backgrounds, these groups revealed diverging micro-evolutionary histories, implying different founding events and different patterns of cultural isolation and local admixture. Partecipanza is an idiosyncratic institution of Medieval origin aimed at sharing and devolving collective lands. This case-study exemplified that socio-economic stratification within the same population may induce sex-biased genetic structuring and the maintenance of otherwise hidden historical genetic traces.

Data management and data analysis in the large European projects GEHA (GEnetics of Healthy Aging) and NU-AGE (NUtrition and AGEing): a bioinformatic approach

The aging process is characterized by the progressive fitness decline experienced at all the levels of physiological organization, from single molecules up to the whole organism. Studies confirmed inflammaging, a chronic low-level inflammation, as a deeply intertwined partner of the aging process, which may provide the “common soil” upon which age-related diseases develop and flourish. Thus, albeit inflammation per se represents a physiological process, it can rapidly become detrimental if it goes out of control causing an excess of local and systemic inflammatory response, a striking risk factor for the elderly population. Developing interventions to counteract the establishment of this state is thus a top priority. Diet, among other factors, represents a good candidate to regulate inflammation. Building on top of this consideration, the EU project NU-AGE is now trying to assess if a Mediterranean diet, fortified for the elderly population needs, may help in modulating inflammaging. To do so, NU-AGE enrolled a total of 1250 subjects, half of which followed a 1-year long diet, and characterized them by mean of the most advanced –omics and non –omics analyses. The aim of this thesis was the development of a solid data management pipeline able to efficiently cope with the results of these assays, which are now flowing inside a centralized database, ready to be used to test the most disparate scientific hypotheses. At the same time, the work hereby described encompasses the data analysis of the GEHA project, which was focused on identifying the genetic determinants of longevity, with a particular focus on developing and applying a method for detecting epistatic interactions in human mtDNA. Eventually, in an effort to propel the adoption of NGS technologies in everyday pipeline, we developed a NGS variant calling pipeline devoted to solve all the sequencing-related issues of the mtDNA.

Mitochondrial DNA rearrangements during human aging: a study on liver, muscle and adipose tissue

Aging is characterized by a chronic, low-grade inflammatory state called “inflammaging”. Mitochondria are the main source of reactive oxygen species (ROS), which trigger the production of pro-inflammatory molecules. We are interested in studying the age-related modifications of the mitochondrial DNA (mtDNA), which can be affected by the lifelong exposure to ROS and are responsible of mitochondrial dysfunction. Moreover, increasing evidences show that telomere shortening, naturally occurring with aging, is involved in mtDNA damage processes and thus in the pathogenesis of age-related disorders. Thus the primary aim of this thesis was the analysis of mtDNA copy number, deletion level and integrity in different-age human biopsies from liver, vastus lateralis skeletal muscle of healthy subjects and patients with limited mobility of lower limbs (LMLL), as well as adipose tissue. The telomere length and the expression of nuclear genes related to mitobiogenesis, fusion and fission, mitophagy, mitochondrial protein quality control system, hypoxia, production and protection from ROS were also evaluated. In liver the decrease in mtDNA integrity with age is accompanied with an increase in mtDNA copy number, suggesting the existence of a “compensatory mechanism” able to maintain the functionality of this organ. Different is the case of vastus lateralis muscle, where any “compensatory pathway” is activated and mtDNA integrity and copy number decrease with age, both in healthy subjects and in patients. Interestingly, mtDNA rearrangements do not incur in adipose tissue with advancing age. Finally, in all tissues a marked gender difference appears, suggesting that aging and also gender diversely affect mtDNA rearrangements and telomere length in the three human tissues considered, likely depending on their different metabolic needs and inflammatory status.

Profili di espressione genetica ed epigenetica ad impatto prognostico e predittivo nel carcinoma squamoso e nelle lesioni potenzialmente maligne del cavo orale

Il carcinoma squamoso orale (CSO) è spesso preceduto da lesioni definite potenzialmente maligne tra cui la leucoplachia e il lichen ma una diagnosi precoce avviene ancora oggi in meno della metà dei casi. Inoltre spesso un paziente trattato per CSO svilupperà secondi tumori. Scopo del lavoro di ricerca è stato: 1) Studiare, mediante metodica di next generation sequencing, lo stato di metilazione di un gruppo di geni a partire da prelievi brushing del cavo orale al fine di identificare CSO o lesioni ad alto rischio di trasformazione maligna. 2) Valurare la relazione esistente tra sovraespressione di p16INK4A e presenza di HPV in 35 pazienti affetti da lichen 3) Valutare la presenza di marker istopatologici predittivi di comparsa di seconde manifestazioni tumorali 4) valutare la relazione clonale tra tumore primitivo e metastasi linfonodale in 8 pazienti mediante 2 metodiche di clonalità differenti: l’analisi di mtDNA e delle mutazioni del gene TP53. I risultati hanno mostrato: 1) i geni ZAP70 e GP1BB hanno presentato un alterato stato di metilazione rispettivamente nel 100% e nel 90,9% di CSO e lesioni ad alto rischio, mentre non sono risultati metilati nei controlli sani; ipotizzando un ruolo come potenziali marcatori per la diagnosi precoce nel CSO. 2)Una sovraespressione di p16INK4A è risultata in 26/35 pazienti affetti da lichen ma HPV-DNA è stato identificato in soli 4 campioni. Nessuna relazione sembra essere tra sovraespressione di p16INK4A e virus HPV. 3)L’invasione perineurale è risultato un marker predittivo della comparsa di recidiva locale e metastasi linfonodale, mentre lo stato dei margini chirurgici si è rilevato un fattore predittivo per la comparsa di secondi tumori primitivi 4) Un totale accordo nei risultati c’è stato tra analisi di mtDNA e analisi di TP53 e le due metodiche hanno identificato la presenza di 4 metastasi linfonodali non clonalmente correlate al tumore primitivo.

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.

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.