TNFRSF13B Genetic variability an anthropological - evolutionary approach to Biomedical Research

In the recent years TNFRSF13B coding variants have been implicated by clinical genetics studies in Common Variable Immunodeficiency (CVID), the most common clinically relevant primary immunodeficiency in individuals of European ancestry, but their functional effects in relation to the development of the disease have not been entirely established. To examine the potential contribution of such variants to CVID, the more comprehensive perspective of an evolutionary approach was applied in this study, underling the belief that evolutionary genetics methods can play a role in dissecting the origin, causes and diffusion of human diseases, representing a powerful tool also in human health research. For this purpose, TNFRSF13B coding region was sequenced in 451 healthy individuals belonging to 26 worldwide populations, in addition to 96 control, 77 CVID and 38 Selective IgA Deficiency (IgAD) individuals from Italy, leading to the first achievement of a global picture of TNFRSF13B nucleotide diversity and haplotype structure and making suggestion of its evolutionary history possible. A slow rate of evolution, within our species and when compared to the chimpanzee, low levels of genetic diversity geographical structure and the absence of recent population specific selective pressures were observed for the examined genomic region, suggesting that geographical distribution of its variability is more plausibly related to its involvement also in innate immunity rather than in adaptive immunity only. This, together with the extremely subtle disease/healthy samples differences observed, suggests that CVID might be more likely related to still unknown environmental and genetic factors, rather than to the nature of TNFRSF13B variants only.

Evolutionary genetics of lactase persistence in Eurasian human populations

Although ability to digest lactose generally declines after weaning in all mammals, in some human populations it persists also in adult individuals, a condition named lactase persistence (LP). Studies on the prevalence of the LP phenotype in worldwide human populations have shown that the frequency of this trait is highly variable in different ethnic groups, appearing to be positively correlated with the importance of milk in the diet. In particular, several single-nucleotide polymorphisms (SNPs) in the proximity of the LCT gene have been proved to be associated with LP. Nevertheless, few studies have till now analyzed genetic variation underlying LP in a wide set of Eurasian populations and, especially, in the Italian one. In the present study, we thus typed 40 SNPs surrounding the LCT gene in more than 1,000 samples from Italian and Arabic peninsulas to investigate patterns of LP-related genetic diversity in two regions which have played a pivotal role in the recent human evolutionary history according to their geographical position and historical/archaeological records. Our results underline a high and complex variability of the explored genomic region in both studied populations. In particular, a clear diversification of Northern Italian groups from the rest of the peninsula, was observed, with the formers being genetically more similar to Northern European populations than to Southern Italians. These observation are consistent with known decreasing pattern of LP from Northern to Southern Italy and suggest the possibility of an independent evolution of LP-associated genotypes in Northern Italy. A similar scenario was observed in the Arabian peninsula, with Dhofari Arabs from Southern Oman and Yemeni clustering together with respect to Arabs from Northern Oman and the subgroup of Omanis of Asian origin which appeared instead to be genetically closer to Europeans than to the rest of Arabic groups.

TiMMing: developing an innovative suite of bioinformatic tools to harmonize and track the origin of copy number alterations in the evolutive history of multiple myeloma

Multiple Myeloma (MM) is a hematologic cancer with heterogeneous and complex genomic landscape, where Copy Number Alterations (CNAs) play a key role in the disease's pathogenesis and prognosis. It is of biological and clinical interest to study the temporal occurrence of early alterations, as they play a disease "driver" function by deregulating key tumor pathways. This study presents an innovative bioinformatic tools suite created for harmonizing and tracing the origin of CNAs throughout the evolutionary history of MM. To this aim, large cohorts of newly-diagnosed MM (NDMM, N=1582) and Smoldering-MM (SMM, N=282) were aggregated. The tools developed in this study enable the harmonization of CNAs as obtained from different genomic platforms in such a way that a high statistical power can be obtained. By doing so, the high numerosity of those cohorts was harnessed for the identification of novel genes characterized as "driver" (NFKB2, NOTCH2, MAX, EVI5 and MYC-ME2-enhancer), and the generation of an innovative timing model, implemented with a statistical method to introduce confidence intervals in the CNAs-calls. By applying this model on both NDMM and SMM cohorts, it was possible to identify specific CNAs (1q(CKS1B)amp, 13q(RB1)del, 11q(CCND1)amp and 14q(MAX)del) and categorize them as "early"/ "driver" events. A high level of precision was guaranteed by the narrow confidence intervals in the timing estimates. These CNAs were proposed as critical MM alterations, which play a foundational role in the evolutionary history of both SMM and NDMM. Finally, a multivariate survival model was able to identify the independent genomic alterations with the greatest effect on patients’ survival, including RB1-del, CKS1B-amp, MYC-amp, NOTCH2-amp and TRAF3-del/mut. In conclusion, the alterations that were identified as both "early-drivers” and correlated with patients’ survival were proposed as biomarkers that, if included in wider survival models, could provide a better disease stratification and an improved prognosis definition.

Image analysis in the morphological and morphometric study of teeth

The subject of this doctoral dissertation concerns the definition of a new methodology for the morphological and morphometric study of fossilized human teeth, and therefore strives to provide a contribution to the reconstruction of human evolutionary history that proposes to extend to the different species of hominid fossils. Standardized investigative methodologies are lacking both regarding the orientation of teeth subject to study and in the analysis that can be carried out on these teeth once they are oriented. The opportunity to standardize a primary analysis methodology is furnished by the study of certain early Neanderthal and preneanderthal molars recovered in two caves in southern Italy [Grotta Taddeo (Taddeo Cave) and Grotta del Poggio (Poggio Cave), near Marina di Camerata, Campania]. To these we can add other molars of Neanderthal and modern man of the upper Paleolithic era, specifically scanned in the paleoanthropology laboratory of the University of Arkansas (Fayetteville, Arkansas, USA), in order to increase the paleoanthropological sample data and thereby make the final results of the analyses more significant. The new analysis methodology is rendered as follows: 1. Standardization of an orientation system for primary molars (superior and inferior), starting from a scan of a sample of 30 molars belonging to modern man (15 M1 inferior and 15 M1 superior), the definition of landmarks, the comparison of various systems and the choice of a system of orientation for each of the two dental typologies. 2. The definition of an analysis procedure that considers only the first 4 millimeters of the dental crown starting from the collar: 5 sections parallel to the plane according to which the tooth has been oriented are carried out, spaced 1 millimeter between them. The intention is to determine a method that allows for the differentiation of fossilized species even in the presence of worn teeth. 3. Results and Conclusions. The new approach to the study of teeth provides a considerable quantity of information that can better be evaluated by increasing the fossil sample data. It has been demonstrated to be a valid tool in evolutionary classification that has allowed (us) to differentiate the Neanderthal sample from that of modern man. In a particular sense the molars of Grotta Taddeo, which up until this point it has not been possible to determine with exactness their species of origin, through the present research they are classified as Neanderthal.

Metodi molecolari nello studio della variabilità genetica di Tuber borchii

Tuber borchii (Ascomycota, order Pezizales) is highly valued truffle sold in local markets in Italy. Despite its economic importance, knowledge on its distribution and population variation is scarce. The objective of this work was to investigate the evolutionary forces shaping the genetic structure of this fungus using coalescent and phylogenetic methods to reconstruct the evolutionary history of populations in Italy. To assess population structure, 61 specimens were collected from 11 different Provinces of Italy. Sampling was stratified across hosts and habitats to maximize coverage in native oak and pine stands and both mychorrizae and fruiting bodies were collected. Samples were identified considering anatomo-morphological characters. DNA was extracted and both multilocus (AFLP) and single-locus (18 loci from rDNA, nDNA, and mtDNA) approaches were used to look for polymorphisms. Screening AFLP profiles, both Jaccard and Dice coefficients of similarity were utilized to transform binary matrix into a distance matrix and then to desume Neighbour-Joining trees. Though these are only preliminary examinations, phylogenetic trees were totally concordant with those deriving from single locus analyses. Phylogenetic analyses of the nuclear loci were performed using maximum likelihood with PAUP and a combined phylogenetic inference, using Bayesian estimation with all nuclear gene regions, was carried out. To reconstruct the evolutionary history, we estimated recurrent migration, migration across the history of the sample, and estimated the mutation and approximate age of mutations in each tree using SNAP Workbench. The combined phylogenetic tree using Bayesian estimation suggests that there are two main haplotypes that are difficult to be differentiated on the basis of morphology, of ecological parameters and symbiontic tree. Between these two lineages, that occur in sympatry within T. borchii populations, there is no evidence of recurrent migration. However, migration over the history of the sample was asymmetrical suggesting that isolation was a result of interrupted gene flow followed by range expansion. Low levels of divergence between the haplotypes indicate that there are likely to be two cryptic species within the T. borchii population sampled. Our results suggest that isolation between populations of T. borchii could have led to reproductive isolation between two lineages. This isolation is likely due to sympatric speciation caused by a multiple colonization from different refugia or a recent isolation. In attempting to determinate whether these haplotypes represent separate species or a partition of the same species we applied Biological and Mechanistic species Concepts. Notwithstanding, further analyses are necessary to evaluate if selection favoured premating or post-mating isolation.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

Environmental influence on the phenotype: Morphological variation in human dentition

This work is about the role that environment plays in the production of evolutionary significant variations. It starts with an historical introduction about the concept of variation and the role of environment in its production. Then, I show how a lack of attention to these topics may lead to serious mistakes in data interpretation. A statistical re-analysis of published data on the effects of malnutrition on dental eruption, shows that what has been interpreted as an increase in the mean value, is actually linked to increase of variability. In Chapter 3 I present the topic of development as a link between variability and environmental influence, giving a review of the possible mechanisms by which development influences evolutionary dynamics. Chapter 4 is the core chapter of the thesis; I investigated the role of environment in the development of dental morphology. I used dental hypoplasia as a marker of stress, characterizing two groups. Comparing the morphology of upper molars in the two groups, three major results came out: (i) there is a significant effect of environmental stressors on the overall morphology of upper molars; (ii) the developmental response increases morphological variability of the stressed population; (iii) increase of variability is directional: stressed individuals have increased cusps dimensions and number. I also hypothesized the molecular mechanisms that could be responsible of the observed effects. In Chapter 5, I present future perspectives for developing this research. The direction of dental development response is the same direction of the trend in mammalian dental evolution. Since malnutrition triggers the developmental response, and this particular kind of stressor must have been very common in our class evolutionary history, I propose the possibility that environmental stress actively influenced mammals evolution. Moreover, I discuss the possibility of reconsidering the role of natural selection in the evolution of dental morphology.

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.

Coastal Forests of Kenya-Ecology, Biodiversity & Conservation

Aims: the broad objective of this study is to investigate the ecological, biodiversity and conservation status of the coastal forests of Kenya fragments. The specific aims of the study are: (1) to investigate current quantitative trends in plant diversity; (2) develop a spatial and standardised vegetation database for the coastal forests Kenya; (3) investigate forest structure, species diversity and composition across the forests; (4) investigate the effect of forest fragment area on plant species diversity; (5) investigate phylogenetic diversity across these coastal remnants (6) assess vulnerability and provide conservation perspectives to concrete policy issues; (7) investigate plant and butterfly diversity correlation. Methods: I performed various analytical methods including species diversity metrics; multiple regression models for species-area relationship and small island effect; non-metric multidimensional scaling; ANOSIM; PERMANOVA; multiplicative beta diversity partitioning; species accumulation curve and species indicator analysis; statistical tests, rarefaction of species richness; phylogenetic diversity metrics of Phylogenetic diversity index, mean pairwise distance, mean nearest taxon distance, and their null-models: and Co-correspondence analysis. Results: developed the first large standardised, spatial and geo-referenced vegetation database for coastal forests of Kenya consisting of 600 plant species, across 25 forest fragments using 158 plots subdivided into 3160 subplots, 18 sacred forests and seven forest reserves; species diversity, composition and forest structure was significantly different across forest sites and between forest reserves and sacred forests, higher beta diversity, species-area relationship explained significant variability of plant diversity, small Island effect was not evident; sacred forests exhibited higher phylogenetic diversity compared to forest reserves; the threatened Red List species contributed higher evolutionary history; a strong correlation between plants and butterfly diversity. Conclusions: This study provides for the first time a standardized and large vegetation data. Results emphasizes need to improve sacred forests protection status and enhance forest connectivity across forest reserves and sacred forests.

Multiscale comparative analysis of marine biominerals and their ecological implications

Marine biomineralizing organisms provide a fundamental link between biology and environment. Calcified structure are important archives that can provide us main means of understanding organism adaptation, habits, environmental characteristics, and to look back in time and explore the past climate and their evolutionary history. In fact, biomineralized structures retain an unparalleled record of current and past ocean conditions through the investigation of their microchemistry and isotopes. This thesis considers aspects of two different biomineralization systems: fish otolith and coral skeletons at macro-, micro- and nanoscale, with the aim to understand how their morphology, structural characteristics and compositions can provide information of their functionality, and the environmental, behavioural, and evolutionary context in which organisms are framed. To this end, I applied a multidisciplinary approach in the scope to investigate calcified structures as “information recorders” and as models to study the phenotypic plasticity.

Genes involved in germline regulative pathways in metazoa: an integrative approach to investigate a key feature of animal biology

In Metazoa, the germline represents the cell lineage devoted to transmission of genetic heredity across generations. Its functions intuitively evoke the crucial roles that it plays in the development of a new organism and in the evolution of the species. Germline establishment is tightly tied to animal multicellularity itself, in which the complex differentiation of cell lineages is favoured by the confinement of totipotency in specific cell populations. In the present thesis, I addressed the subject of germline characterization in animals through different approaches, in an attempt to cover different sides and scales. First, I investigated the extent and nature of shared differentially transcribed molecular factors in 10 different species germline-related lineages. I observed that newly evolved genes are less likely to be involved in germline-related mechanisms and that the mostly shared transcriptional signal across the species considered was the upregulation of genes associated to proper DNA replication, instead of the expected transcriptional and post-transcriptional regulation, that apparently have a higher level of lineage-specificity. I then focused on the evolutionary history of Tudor domain containing proteins, a gene family that underwent germline-associated expansions in animals. Using data from 24 holozoan phyla, I could confirm the previously proposed evolution of the Tudor domain secondary structure. Also, I associated lineage-specific family reductions and expansions to peculiar genomic dynamics and to the evolution of germline-associated piRNA pathway of retrotransposon silencing. Lastly, I characterized and investigated the expression of the Tudor protein TDRD7 in the clam Ruditapes philippinarum. Through immunolocalization, I could compare its expression profiles in gametogenic specimens to the previously characterized germline marker vasa. Combining results with literature, I proposed that, in this species, TDRD7 is involved in the assembly of germ granules, i.e. cytoplasmic structures associated to germline differentiation in virtually all animals, but whose assemblers can be taxon specific.

Multidisciplinary approach to horse welfare in different management systems

From its domestication until nowadays, the horse has assumed multiple roles in human society. Over time, this condition and the lack of specific regulation have led to the development of different kinds of management systems for this species. This Ph.D. research project aims to investigate horses' welfare in different management practices and housing systems, considering a multidisciplinary approach, taking into account biological function, naturalness, and affective dimension. The results are presented in five articles that evidence risk factors that can mine horse welfare, and examine tools and parameters that can be employed for its assessment. Our research shows the importance of considering the evolutionary history and the species-specific and behavioural needs of horses in their management and housing. Sociality, free movement, diet composition and foraging routine, and the workload that these animals undergo are important factors that should be taken into account. Furthermore, this research has evidenced the importance of employing different parameters (e.g., behaviour, endocrinological parameters, and immune activity) in welfare assessment and proposes the use of horsehair DHEA (dehydroepiandrosterone) as a possible useful additional non-invasive measure for the investigation of long-term stress conditions. Finally, our results underline the importance of considering the affective dimension in welfare research. Recently, Judgement Bias Tests (JBT), which are based on the influence of affective states on the decision-making process, have been widely employed in animal welfare research. However, our studies show that the use of spatial JBT in horses can have some limitations. Still today several management systems do not fulfill species-specific needs of horses, thus the implementation of specific regulations could ameliorate horse welfare. A multidisciplinary approach to welfare assessment is fundamental, but it should be always remembered the individual and its own characteristics, which can influence not only physiological, immunological, and behavioural responses but also emotional and cognitive dimensions.

Exploring the internal kinematics of Galactic Globular Cluster cores

Globular clusters (GCs) are traditionally described as simple quasi-relaxed non-rotating stellar systems, characterized by spherical symmetry and isotropy in velocity space. However, recent studies have shown deviations from isotropic velocity distributions and significant internal rotation in many GCs, suggesting that their internal structure and kinematics are more complex than previously thought. The aim of this thesis is to investigate the internal kinematics of Galactic Globular Clusters (GGCs) as part of the Multi-Instrument Kinematic Survey (MIKiS), which exploits the capabilities of different ESO-VLT spectrographs to obtain comprehensive velocity dispersion (VD) and rotation profiles of GGCs. Moreover, this thesis has the particular goal of unraveling the kinematics of GC cores, which are still largely unexplored, by taking advantage of the exceptional spatial resolution of the adaptive-optics assisted integral-field spectrograph MUSE/NFM. The thesis presents a thorough kinematic study of three GGCs NGC 1904, NGC 6440, and NGC 6569. By combining the data sets acquired with four different spectrographs, we obtained the radial velocity (RV) of more than 1000 individual stars in each cluster, sampling from the innermost to the outermost regions. This allowed us to obtain the entire VD profile of each cluster and exclude the presence of an intermediate-mass black hole in the core of NGC 1904, at odds with previous findings obtained from integrated-light spectra. The studies also revealed signatures of internal rotation in each of the GCs studied. These results, supported by those of N-body simulations, prove that GCs were born with a significant initial rotation that they gradually lost through internal two-body relaxation and angular momentum loss carried away by escaping stars. Furthermore, we derived the structural parameters of NGC 6440 and NGC 6569, obtaining a comprehensive overview of the internal kinematics and structure of these GCs, which is necessary to properly reconstruct the evolutionary history of these systems.

Genetic characterization, population history and evolutionary medicine perspective in two native south american populations: Yanesha and Wichi

Two Amerindian populations from the Peruvian Amazon (Yanesha) and from rural lowlands of the Argentinean Gran Chaco (Wichi) were analyzed. They represent two case study of the South American genetic variability. The Yanesha represent a model of population isolated for long-time in the Amazon rainforest, characterized by environmental and altitudinal stratifications. The Wichi represent a model of population living in an area recently colonized by European populations (the Criollos are the population of the admixed descendents), whose aim is to depict the native ancestral gene pool and the degree of admixture, in relation to the very high prevalence of Chagas disease. The methods used for the genotyping are common, concerning the Y chromosome markers (male lineage) and the mitochondrial markers (maternal lineage). The determination of the phylogeographic diagnostic polymorphisms was carried out by the classical techniques of PCR, restriction enzymes, sequencing and specific mini-sequencing. New method for the detection of the protozoa Trypanosoma cruzi was developed by means of the nested PCR. The main results show patterns of genetic stratification in Yanesha forest communities, referable to different migrations at different times, estimated by Bayesian analyses. In particular Yanesha were considered as a population of transition between the Amazon basin and the Andean Cordillera, evaluating the potential migration routes and the separation of clusters of community in relation to different genetic bio-ancestry. As the Wichi, the gene pool analyzed appears clearly differentiated by the admixed sympatric Criollos, due to strict social practices (deeply analyzed with the support of cultural anthropological tools) that have preserved the native identity at a diachronic level. A pattern of distribution of the seropositivity in relation to the different phylogenetic lineages (the adaptation in evolutionary terms) does not appear, neither Amerindian nor European, but in relation to environmental and living conditions of the two distinct subpopulations.

Unlocking ecological history using fish remains: Eco-evolutionary consequences of exploitation in the Atlantic bluefin tuna

During recent decades, the health of ocean ecosystems and fish populations has been threatened by overexploitation, pollution, and anthropogenic-driven climate change. Due to a lack of long-term data, we have a poor understanding of when intensive exploitation began and what impact anthropogenic activities have had on the ecology and evolution of fishes. Such information is crucial to recover degraded and depleted marine ecosystems and fish populations, maximise their productivity in-line with historical levels, and predict their future dynamics. In this thesis, I evaluate anthropogenic impacts on the iconic Atlantic bluefin tuna (Thunnus thynnus; BFT), one of the longest and recently most intensely exploited marine fishes, with a tremendous cultural and economic importance. Using a long-time series of archaeological and archived faunal remains (bones) dating back to approximately two millennia ago, I apply morphological, isotopic, and genomic techniques to perform the first studies on long-term BFT size and growth, diet and habitat use, and demography and adaptation, and produce the first genome-wide data on this species. My findings suggest that exploitation had impacted BFT foraging behaviour by the ~16th century when coastal ecosystem degradation induced a pelagic shift in diet and habitat use. I reveal that BFT biomass began to decline much earlier than hitherto documented, by the 19th century, consistent with intensive tuna trap catches during this period and catch-at-size increasing. I find that BFT juvenile growth had increased by the early 1900s (and more dramatically by the 21st century) which may reflect an evolutionary response to size selective harvest–which I find putative genomic signatures of. Further, I observed that BFT foraging behaviours have been modified following overexploitation during the 20th century, which previously included a isotopically distinct, Black Sea niche. Finally, I show that despite biomass declining from centuries ago, BFT has retained genomic diversity.