Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

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.

Study of variability and genetic structure of European populations of Myotis emarginatus and Myotis capaccinii (Chiroptera, Vespertilionidae)

The Geoffroy’s bat Myotis emarginatus is mainly present in southern, south-eastern and central Europe (Červerný, 1999) and is often recorded from northern Spain (Quetglas, 2002; Flaquer et al., 2004). It has demonstrated the species’ preference for forest. Myotis capaccinii, confined to the Mediterranean (Guille´n, 1999), is classified as ‘vulnerable’ on a global scale (Hutson, Mickleburgh & Racey, 2001). In general, the species preferred calm waters bordered by well-developed riparian vegetation and large (> 5 m) inter-bank distances (Biscardi et al. 2007). In this study we present the first results about population genetic structure of these two species of genus Myotis. We used two methods of sampling: invasive and non-invasive techniques. A total of 323 invasive samples and a total of 107 non-invasive samples were collected and analyzed. For Myotis emarginatus we have individuated for the first time a set of 7 microsatellites, which can work on this species, started from a set developed on Myotis myotis (Castella et al. 2000). We developed also a method for analysis of non-invasive samples, that given a good percentage of positive analyzed samples. The results have highlighted for the species Myotis emarginatus the presence on the European territory of two big groups, discovered by using the microsatellites tracers. On this species, 33 haplotypes of Dloop have been identified, some of them are presented only in some colonies. We identified respectively 33 haplotypes of Dloop and 10 of cytB for Myotis emarginatus and 25 of dloop and 15 of cytB for Myotis capaccinii. Myotis emarginatus’ results, both microsatellites and mtDNA, show that there is a strong genetic flow between different colonies across Europe. The results achieved on Myotis capaccinii are very interesting, in this case either for the microsatellites or the mitochondrial DNA sequences, and it has been highlighted a big difference between different colonies.

Comparing genetic and linguistic diversity in African populations with a focus on the Khoisan of southern Africa

The interaction between disciplines in the study of human population history is of primary importance, profiting from the biological and cultural characteristics of humankind. In fact, data from genetics, linguistics, archaeology and cultural anthropology can be combined to allow for a broader research perspective. This multidisciplinary approach is here applied to the study of the prehistory of sub-Saharan African populations: in this continent, where Homo sapiens originally started his evolution and diversification, the understanding of the patterns of human variation has a crucial relevance. For this dissertation, molecular data is interpreted and complemented with a major contribution from linguistics: linguistic data are compared to the genetic data and the research questions are contextualized within a linguistic perspective. In the four articles proposed, we analyze Y chromosome SNPs and STRs profiles and full mtDNA genomes on a representative number of samples to investigate key questions of African human variability. Some of these questions address i) the amount of genetic variation on a continental scale and the effects of the widespread migration of Bantu speakers, ii) the extent of ancient population structure, which has been lost in present day populations, iii) the colonization of the southern edge of the continent together with the degree of population contact/replacement, and iv) the prehistory of the diverse Khoisan ethnolinguistic groups, who were traditionally understudied despite representing one of the most ancient divergences of modern human phylogeny. Our results uncover a deep level of genetic structure within the continent and a multilayered pattern of contact between populations. These case studies represent a valuable contribution to the debate on our prehistory and open up further research threads.

Applied genomics to monitor the genetic diversity and integrity of honey bee (Apis mellifera) populations and define sustainable strategies to benefit the beekeeping sector

Honey bees are considered keystone species in ecosystem, the effect of harmful pesticides for the honey bees, the action of extreme climatic waves and their consequence on honey bees health can cause the loss of many colonies which could contribute to the reduction of the effective population size and incentive the use of non-autochthonous queens to replace dead colonies. Over the last decades, the use of non-ligustica bee subspecies in Italy has increased and together with the mentioned phenomena exposed native honey bees to hybridization, laeding to a dramatic loss of genetic erosion and admixture. Healthy genetic diversity within honey bee populations is critical to provide tolerance and resistance to current and future threatening. Nowadays it is urgent to design strategies for the conservation of local subspecies and their valorisation on a productive scale. In this Thesis we applied genomics tool for the analysis of the genetic diversity and the genomic integrity of honey bee populations in Italy are described. In this work mtDNA based methods are presented using honey bee DNA or honey eDNA as source of information of the genetic diversity of A. mellifera at different level. Taken together, the results derived from these studies should enlarge the knowledge of the genetic diversity and integrity of the honey bee populations in Italy, filling the gap of information necessary to design efficient conservation programmes. Furthermore, the methods presented in these works will provide a tool for the honey authentication to sustain and valorise beekeeping products and sector against frauds.

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.

Phylogeny and genetic diversity of Italian species of hares (genus Lepus)

In this study we have analysed the genetic variability in ca. 700 samples belonging to six species of genus Lepus using maternal and biparental molecular markers (mitochondrial DNA, microsatellites, Single Nucleotide Polimorphisms). We aimed to reconstruct the phylogenetic relationships of species of hares living in Europe, and assess the occurrence of hybridization between the European hare Lepus europaeus and the Italian hare Lepus corsicanus. Results showed a deep genetic differentiation and absence of hybridization between L. corsicanus and L. europaeus, confirming that they are distinct and distantly related biological species. In contrast, we showed small genetic distances and a close phylogenetic relationship between the Italian hare and Cantabrian hare L. castroviejoi, which suggest a deeper evaluation of their taxonomic status. Populations of L. corsicanus are geographically differentiated. In particular, the peninsular and Sicilian populations of Italian hares are sharply genetically distinct, which calls for avoiding any translocation between Italy and Sicily. Information on genetic variability and population structure is being used to implement the Italian Action Plan for L. corsicanus.

The genetic signature of Neolithic in Greece

The Neolithic is characterized by the transition from a subsistence economy, based on hunting and gathering, to one based on food producing. This important change was paralleled by one of the most significant demographic increase in the recent history of European populations. The earliest Neolithic sites in Europe are located in Greece. However, the debate regarding the colonization route followed by the Middle-eastern farmers is still open. Based on archaeological, archaeobotanical, craniometric and genetic data, two main hypotheses have been proposed. The first implies the maritime colonization of North-eastern Peloponnesus from Crete, whereas the second points to an island hopping route that finally brought migrants to Central Greece. To test these hypotheses using a genetic approach, 206 samples were collected from the two Greek regions proposed as the arrival point of the two routes (Korinthian district and Euboea). Expectations for each hypothesis were compared with empirical observations based on the analysis of 60 SNPs and 26 microsatellite loci of Y-chromosome and mitochondrial DNA hypervariable region I. The analysis of Y-chromosome haplogroups revealed a strong genetic affinity of Euboea with Anatolian and Middle-eastern populations. The inferences of the time since population expansion suggests an earlier usage of agriculture in Euboea. Moreover, the haplogroup J2a-M410, supposed to be associated with the Neolithic transition, was observed at higher frequency and variance in Euboea showing, for both these parameters, a decreasing gradient moving from this area. The time since expansion estimates for J2a-M410 was found to be compatible with the Neolithic and slightly older in Euboea. The analysis of mtDNA resulted less informative. However, a higher genetic affinity of Euboea with Anatolian and Middle-eastern populations was confirmed. These results taken as a whole suggests that the most probable route followed by Neolithic farmers during the colonization of Greece was the island hopping route.

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.

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.

Evolutionary epigenetics of modern human populations

Epigenetic variability is a new mechanism for the study of human microevolution, because it creates both phenotypic diversity within an individual and within population. This mechanism constitutes an important reservoir for adaptation in response to new stimuli and recent studies have demonstrated that selective pressures shape not only the genetic code but also DNA methylation profiles. The aim of this thesis is the study of the role of DNA methylation changes in human adaptive processes, considering the Italian peninsula and macro-geographical areas. A whole-genome analysis of DNA methylation profile across the Italian penisula identified some genes whose methylation levels differ between individuals of different Italian districts (South, Centre and North of Italy). These genes are involved in nitrogen compound metabolism and genes involved in pathogens response. Considering individuals with different macro-geographical origins (individuals of Asians, European and African ancestry) more significant DMRs (differentially methylated regions) were identified and are located in genes involved in glucoronidation, in immune response as well as in cell comunication processes. A "profile" of each ancestry (African, Asian and European) was described. Moreover a deepen analysis of three candidate genes (KRTCAP3, MAD1L and BRSK2) in a cohort of individuals of different countries (Morocco, Nigeria, China and Philippines) living in Bologna, was performed in order to explore genetic and epigenetic diversity. Moreover this thesis have paved the way for the application of DNA methylation for the study of hystorical remains and in particular for the age-estimation of individuals starting from biological samples (such as teeth or blood). Noteworthy, a mathematical model that considered methylation values of DNA extracted from cementum and pulp of living individuals can estimate chronological age with high accuracy (median absolute difference between age estimated from DNA methylation and chronological age was 1.2 years).

Assessment of the population structure and temporal changes in spatial dynamics and genetic characteristics of the Atlantic bluefin tuna under a fishery independent framework.

As a large and long-lived species with high economic value, restricted spawning areas and short spawning periods, the Atlantic bluefin tuna (BFT; Thunnus thynnus) is particularly susceptible to over-exploitation. Although BFT have been targeted by fisheries in the Mediterranean Sea for thousands of years, it has only been in these last decades that the exploitation rate has reached far beyond sustainable levels. An understanding of the population structure, spatial dynamics, exploitation rates and the environmental variables that affect BFT is crucial for the conservation of the species. The aims of this PhD project were 1) to assess the accuracy of larval identification methods, 2) determine the genetic structure of modern BFT populations, 3) assess the self-recruitment rate in the Gulf of Mexico and Mediterranean spawning areas, 4) estimate the immigration rate of BFT to feeding aggregations from the various spawning areas, and 5) develop tools capable of investigating the temporal stability of population structuring in the Mediterranean Sea. Several weaknesses in modern morphology-based taxonomy including demographic decline of expert taxonomists, flawed identification keys, reluctance of the taxonomic community to embrace advances in digital communications and a general scarcity of modern user-friendly materials are reviewed. Barcoding of scombrid larvae revealed important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. Using a Genotyping-by-Sequencing a panel of 95 SNPs was developed and used to characterize the population structuring of BFT and composition of adult feeding aggregations. Using novel molecular techniques, DNA was extracted from bluefin tuna vertebrae excavated from late iron age, ancient roman settlements Byzantine-era Constantinople and a 20th century collection. A second panel of 96 SNPs was developed to genotype historical and modern samples in order to elucidate changes in population structuring and allele frequencies of loci associated with selective traits.

The genetic prehistory of Northern Calabria. Archeogenomic investigation of the human remains from Grotta di Pietra Sant’Angelo and Grotta della Monaca.

Given its strategic position at the center of the Mediterranean Basin, and its unique history of contacts and migrations, Calabria is an ideal region to decipher the genetic traces of at least some of the numerous demographic prehistoric events in Southern Italy. This thesis focuses on the genetic and social changes of ancient inhabitants of Calabria, covering a timeline of approximately 7000 to 3300 years ago, ranging from the Middle Neolithic to the Middle Bronze Age. We generated the first genome-wide data from Calabria, by focusing on the single inhumation of “Grotta di Pietra Sant’Angelo” (San Lorenzo Bellizzi, Cosenza) and on the vast community found buried in “Grotta della Monaca” (Sant’Agata di Esaro, Cosenza). Supported by archaeological evidence, the primary objective of this research was to employ paleogenomic evidence to decipher funerary customs, social organization, family ties, and demographic shifts in Southern Italy over a period extending beyond three millennia. The possibility of gender-related burial practices and kinship ties among the deceased was also explored. Subsequently, the biogeographical origin and ancestry of prehistoric people of Calabria was contextualized within the broad landscape of existing data on Mediterranean populations. By generating the first genomic evidence from prehistoric Calabria, unresolved questions were addressed, related to the appearance and persistence of distinct genomic components, such as the Iran-related and the Steppe-related ancestry, whose impact on ancient Southern Italian genomes remains uncharted.