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.

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.

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.

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.

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.

Molecular diversity of Primula apennina and phylogeny of Primula Subsection Euauricula

Primula apennina Widmer is endemic to the North Apennines (Italy). ISSR were used to detect the genetic diversity within and among six populations representative of the species distribution range. High levels of genetic diversity were revealed both at population (PPB = 75.92%, HS = 0.204, Hpop = 0.319) and at species level (PPB = 96.95%, HT = 0.242, Hsp = 0.381). Nei gene diversity statistics (15.7%), Shannon diversity index (16.3%) and AMOVA (14%) detected a moderate level of interpopulation diversity. Principal coordinate and bayesian analyses clustered the populations in three major groups along a geographic gradient. The correlation between genetic and geographic distances was positive (Mantel test, r = 0.232). All together, these analyses revealed a weak but significant spatial genetic structure in P. apennina, with gene flow acting as a homogenizing force that prevents a stronger differentiation of populations. Conservation measures are suggested based on the observed pattern of genetic variability. P. apennina belongs to Primula subsect. Euauricula which includes 15 species distributed on the whole Alps and Apennines. A phylogenetic analysis was carried out using AFLP markers in order both to clarify the relationships among the species of subsection Euauricula that remained unresolved in previous works and to make some hypoteses on their evolutive dynamics. NJ, PCO and BAPS analyses strongly confirmed the monophyly of P. subsect. Euauricula and all the species form strongly supported clades. NJ tree topology suggested a simultaneous fragmentations of ancestral species in a large number of isolated populations that survived in refugia along the unglaciated margins of the Alps in response to the Pleistocene climatic oscillations.

Phylogeny and population dynamics in European Reticulitermes and Kalotermes genera (Insecta, Isoptera)

The aim of this thesis is to detect the phylogeny and the population dynamics of the European termites of the genera Reticulitermes and Kalotermes, by the use of different mitochondrial (16S, COI/tRNA/COII, CR) and nuclear (microsatellites and Inter-SINE) molecular markers. In the phylogenetic analyses, the obtained results have well defined the cladogenetic events that generated the nowadays species biodiversity of the genus Reticulitermes, while the analysis of the Kalotermes flavicollis taxon showed the presence of at least four genetic clades, defined on the basis of the geographical distance. The second part of the thesis is centred on the population dynamics of two species: Reticulitermes urbis and Kalotermes flavicollis. The first species, native of the Balkans, is known to be present in some cities of Italy and France. I’ve analyzed the colony genetic structure of the introduced population of Bagnacavallo (RA, Italy), using nine microsatellite loci. The obtained results are in accordance with those obtained from another population in France: this species in fact confirms its invasive and infestation capacities. The analysis of the natural population of K. flavicollis, performed with a combination of mitochondrial (control region) and nuclear (microsatellites and I-SINE) markers, clearly evidenced the presence of two genetic lineages that coexist in the same area. Moreover, results clearly indicate that the cross-breeding is allowed. Finally, the whole results are discussed in a comparative view to better understand the differences in ecology, evolutionary dynamics and colony social structure between these two genera.

Protostrongylids nematodes as a tool to study the biogeography of wild mammals

The European brown hare (Lepus europaeus Pallas, 1778) is an important game species, distributed across Europe and introduced in other regions. Recently, a geographically isolated population, closely related to an ancestral lineage of Lepus europaeus meridiei, was found on Pianosa Island, off the coast of Tuscany, Italy (Mengoni et al., 2018). Thus, the unique opportunity to explore the evolution and genetic structure of its helminth parasites was added to its exceptional isolation condition. Various lungworm species within the genus Protostrongylus (Nematoda: Protostrongylidae) are described in European brown hares. Our aim was to analyze the parasite population through morphological and molecular approaches in order to study the biogeography of the European brown hares (L. e. meridiei) population from Pianosa Island. Moreover, we investigated the morphology of a monospecific genus, i.e. Orthostrongylus, considering its quite intrigant descriptive history and its still unclear and debated classification. Nuclear and mitochondrial markers were used based on their resolution power and expected polymorphism; the whole Internal Transcribed Spacer 1 and 2 (ITS), including the 5.8S rRNA sequence and the Large Subunit (28S) were used, as nuclear genes, for confirmation of the species identification. Conversely, the cytochrome oxidase c subunit I (COI) was used, as mithocondrial genes, to assess interspecific genetic relationships. Molecular analysis corroborated the morphological identification since all the generated ITS and LSU sequences were 100% consistent with the species Protostrongylus oryctolagi and Orthostrongylus macrotis. The paucity of molecular data existent about this genus of parasites underlines the need for more insight’s studies. An in-depth analysis of broncho-pulmonary parasites and the host-parasites relationships along with the improvement of the use of mitochondrial genes, as well as the assessment of new polymorphic markers could contribute to an extensive understanding of parasites fauna and taxonomy, as well as their relationship with wild mammals’ hosts.

Analyzing the dependence structure of microarray data: a copula–based approach

The main aim of this Ph.D. dissertation is the study of clustering dependent data by means of copula functions with particular emphasis on microarray data. Copula functions are a popular multivariate modeling tool in each field where the multivariate dependence is of great interest and their use in clustering has not been still investigated. The first part of this work contains the review of the literature of clustering methods, copula functions and microarray experiments. The attention focuses on the K–means (Hartigan, 1975; Hartigan and Wong, 1979), the hierarchical (Everitt, 1974) and the model–based (Fraley and Raftery, 1998, 1999, 2000, 2007) clustering techniques because their performance is compared. Then, the probabilistic interpretation of the Sklar’s theorem (Sklar’s, 1959), the estimation methods for copulas like the Inference for Margins (Joe and Xu, 1996) and the Archimedean and Elliptical copula families are presented. In the end, applications of clustering methods and copulas to the genetic and microarray experiments are highlighted. The second part contains the original contribution proposed. A simulation study is performed in order to evaluate the performance of the K–means and the hierarchical bottom–up clustering methods in identifying clusters according to the dependence structure of the data generating process. Different simulations are performed by varying different conditions (e.g., the kind of margins (distinct, overlapping and nested) and the value of the dependence parameter ) and the results are evaluated by means of different measures of performance. In light of the simulation results and of the limits of the two investigated clustering methods, a new clustering algorithm based on copula functions (‘CoClust’ in brief) is proposed. The basic idea, the iterative procedure of the CoClust and the description of the written R functions with their output are given. The CoClust algorithm is tested on simulated data (by varying the number of clusters, the copula models, the dependence parameter value and the degree of overlap of margins) and is compared with the performance of model–based clustering by using different measures of performance, like the percentage of well–identified number of clusters and the not rejection percentage of H0 on . It is shown that the CoClust algorithm allows to overcome all observed limits of the other investigated clustering techniques and is able to identify clusters according to the dependence structure of the data independently of the degree of overlap of margins and the strength of the dependence. The CoClust uses a criterion based on the maximized log–likelihood function of the copula and can virtually account for any possible dependence relationship between observations. Many peculiar characteristics are shown for the CoClust, e.g. its capability of identifying the true number of clusters and the fact that it does not require a starting classification. Finally, the CoClust algorithm is applied to the real microarray data of Hedenfalk et al. (2001) both to the gene expressions observed in three different cancer samples and to the columns (tumor samples) of the whole data matrix.

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.

Genome characterization through a mathematical model of the genetic code: an analysis of the whole chromosome 1 of A. thaliana

The objective of this work is to characterize the genome of the chromosome 1 of A.thaliana, a small flowering plants used as a model organism in studies of biology and genetics, on the basis of a recent mathematical model of the genetic code. I analyze and compare different portions of the genome: genes, exons, coding sequences (CDS), introns, long introns, intergenes, untranslated regions (UTR) and regulatory sequences. In order to accomplish the task, I transformed nucleotide sequences into binary sequences based on the definition of the three different dichotomic classes. The descriptive analysis of binary strings indicate the presence of regularities in each portion of the genome considered. In particular, there are remarkable differences between coding sequences (CDS and exons) and non-coding sequences, suggesting that the frame is important only for coding sequences and that dichotomic classes can be useful to recognize them. Then, I assessed the existence of short-range dependence between binary sequences computed on the basis of the different dichotomic classes. I used three different measures of dependence: the well-known chi-squared test and two indices derived from the concept of entropy i.e. Mutual Information (MI) and Sρ, a normalized version of the “Bhattacharya Hellinger Matusita distance”. The results show that there is a significant short-range dependence structure only for the coding sequences whose existence is a clue of an underlying error detection and correction mechanism. No doubt, further studies are needed in order to assess how the information carried by dichotomic classes could discriminate between coding and noncoding sequence and, therefore, contribute to unveil the role of the mathematical structure in error detection and correction mechanisms. Still, I have shown the potential of the approach presented for understanding the management of genetic information.

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.