Cryptic lineages and Pleistocene population expansion in a Brazilian Cerrado frog

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Diversity and genetic population structure of the Brazilian sharpnose shark Rhizoprionodon lalandii

Similar to many small, range-restricted elasmobranchs, the Brazilian sharpnose shark (Rhizoprionodon lalandii) is listed as 'data deficient' by the International Union for the Conservation of Nature (IUCN). Data on stock assessment and sustainability are scarce, and there is no information on population structure. This constitutes a management problem because this shark comprises approximately 50% of the catch of small coastal sharks in Brazil. In this study, populations of R. lalandii distributed from the Caribbean to southern Brazil were investigated using sequences from the mitochondrial DNA control region. Analysis of molecular variance revealed strong structuring between population samples from the Caribbean and those from the Brazilian coast (F{cyrillic}ST=0.254, P<0.0001). Significant differences in the rates of genetic diversity between these major areas were also detected. The observed levels of population structuring are likely to be driven by female phylopatry. Therefore, the identification of both mating and nursery areas with parallel ban/restriction of fishing in these areas may be critical for the long-term sustainability of these populations. © 2013 John Wiley & Sons, Ltd.

High Levels of Diversity Uncovered in a Widespread Nominal Taxon: Continental Phylogeography of the Neotropical Tree Frog Dendropsophus minutus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phylogeographic analysis in populations of Protonectarina sylveirae (Saussure, 1854) (Hymenoptera, Vespidae, Polistinae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Skin lesions on blue whales off southern Chile: Possible conservation implications?

We report on three types of skin lesions in a population of blue whales, Balaenoptera musculus, off the northwestern coast of Isla Grande de Chiloe, Chile. These lesions were: (1) cookie-cutter shark, Isistius brasilensis, bites, (2) vesicular or blister lesions, and (3) a tattoo-like skin disease. The presence of these lesions was determined by the examining photos collected in 2006 and 2007 for a blue whale photo-identification project. We examined 289 photographs of 68 individuals for lesions. The cookie-cutter shark lesions are common on these blue whales and similar to those reported from other species of cetaceans. Skin peeling or shedding was observed on some whales and is believed to be a normal condition. Based on the photographs examined to date the vesicular lesions are more common than the tattoo-like lesions. The tattoo-like skin lesions was observed just on a single whale in 2007. The blister lesions were common on whales in both 2006 and 2007. The presence of blister lesions in both years may indicate that this “disease” will be present in the population for a long time. It is unknown if these lesions contribute to mortality of blue whales frequenting Chilean waters, but the tattoo-like skin lesions if shown to be a pox virus could cause neonatal and calf mortality. Additional investigations are needed that, as a minimum, must include the histological and genetic examination of the two types of disease from live or dead whales, especially the tattoo-like skin lesions. Until this work is undertaken, it will be impossible to determine if these lesions pose a conservation risk to the blue whales off Chile.

Phylogeographic analysis in populations of Protonectarina sylveirae (Saussure, 1854) (Hymenoptera, Vespidae, Polistinae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phylogeography of an Atlantic forest passerine reveals demographic stability through the last glacial maximum

In this study we analyzed the phylogeographic pattern and historical demography of an endemic Atlantic forest (AF) bird, Basileuterus leucoblepharus, and test the influence of the last glacial maximum (LGM) on its population effective size using coalescent simulations. We address two main questions: (i) Does B. leucoblepharus present population genetic structure congruent with the patterns observed for other AF organisms? (ii) How did the LGM affect the effective population size of B. leucoblepharus? We sequenced 914 bp of the mitochondrial gene cytochrome b and 512 bp of the nuclear intron 5 of beta-fibrinogen of 62 individuals from 15 localities along the AF. Both molecular markers revealed no genetic structure in B. leucoblepharus. Neutrality tests based on both loci showed significant demographic expansion. The extended Bayesian skyline plot showed that the species seems to have experienced demographic expansion starting around 300,000 years ago, during the late Pleistocene. This date does not coincide with the LGM and the dynamics of population size showed stability during the LGM. To further test the effect of the LGM on this species, we simulated seven demographic scenarios to explore whether populations suffered specific bottlenecks. The scenarios most congruent with our data were population stability during the LGM with bottlenecks older than this period. This is the first example of an AF organism that does not show phylogeographic breaks caused by vicariant events associated to climate change and geotectonic activities in the Quaternary. Differential ecological, environmental tolerances and habitat requirements are possibly influencing the different evolutionary histories of these organisms. Our results show that the history of organism diversification in this megadiverse Neotropical forest is complex. Crown Copyright (c) 2012 Published by Elsevier Inc. All rights reserved.

DNA barcoding taxonomy and diversity of the Mediterranean sharks and rays

Elasmobranchs are an important by-catch of commercial fisheries targeting bony fishes. Fisheries targeting sharks are rare, but usually almost all specimen bycatched are marketed. They risk extinction if current fishing pressure continues (Ferretti et al., 2008). Accurate species identification is critical for the design of sustainable fisheries and appropriate management plans, especially since not all species are equally sensitive to fishing pressure (Walker & Hislop 1998). The identification of species constitutes the first basic step for biodiversity monitoring and conservation (Dayrat B et al., 2005). More recently, mtDNA sequencing has also been used for species identification and its use has become widespread under the DNA Barcode initiative (e.g. Hebert et al. 2003a, 2003b; Ward et al. 2005, 2008a; Moura et al 2008; Steinke et al. 2009). The aims of this work were: 1) identify sharks and skates species using DNA barcode; 2) compare species of different provenance; 3) use DNA barcode for misidentified species. Using DNA barcode 15 species of sharks (Alopias vulpinus, Centrophorus granulosus, Cetorhinus maximus, Dalatias licha, Etmopterus spinax, Galeorhinus galeus, Galeus melastomus, Heptranchias perlo, Hexanchus griseus, Mustelus mustelus, Mustelus punctulatus, Oxynotus centrina, Scyliorhinus canicula Squalus acanthias, Squalus blainville), 1 species of chimaera (Chimaera monstrosa) and 21 species of rays/skayes (Dasyatis centroura, Dasyatis pastinaca, Dasyatis sp., Dipturus nidarosiensis, Dipturus oxyrinchus, Leucoraja circularis, Leucoraja melitensis, Myliobatis aquila, Pteromylaeus bovinus, Pteroplatytrygon violacea, Raja asterias, Raja brachyura, Raja clavata, Raja miraletus, Raja montagui, Raja radula, Raja polystigma, Raja undulata, Rostroraja alba, Torpedo marmorata, Torpedo nobiliana, Torpedo torpedo) was identified.

Molecular phylogeography of the European coastal plants Crithmum maritimum L., Halimione portulacoides (L.) Aellen, Salsola kali L. and Calystegia soldanella (L.) R. Br.

The intraspecific phylogeography of four European coastal plants, Crithmum maritimum, Halimione portulacoides, Salsola kali and Calystegia soldanella, was inferred from AFLP and ITS data. Only in C. maritimum, H. portulacoides and S. kali, a spatial genetic structure was revealed. The phylogeographic similarities and dissimilarities of these species include: (1) All three have distinct Black/Aegean and Adriatic Sea clusters. (2) Salsola kali and H. portulacoides show a distinct Atlantic/North Sea/Baltic Sea cluster, while Atlantic and eastern Spanish material of C. maritimum clustered together. (3) In the west Mediterranean, only S. kali forms a single cluster, while both H. portulacoides and C. maritimum display a phylogeographic break in the vicinity of the southern French coast. For S. kali, AFLP and ITS data concur in identifying separate Atlantic, east and west Mediterranean clades. All these patterns are postulated to result from both temperature changes during the last glacial and contemporary sea currents. No geographic AFLP structure was revealed in C. soldanella, both at the range-wide and population level. This was attributed to the remarkable seed dispersal ability of this species and possibly its longevity and clonal growth, preserving a random pattern of genetic variation generated by long-distance seed dispersal over long time periods.

European phylogeography of the coastal plants Cakile maritima Scop. (Brassicaceae) and Eryngium maritimum L. (Apiaceae)

Linear dispersal systems, such as coastal habitats, are well suited for phylogeographic studies because of their low spatial complexity compared to three dimensional habitats. Widely distributed coastal plant species additionally show azonal and often essentially continuous distributions. These properties, firstly, make it easier to reconstruct historical distributions of coastal plants and, secondly, make it more likely that present distributions contain both Quaternary refugia and recently colonized areas. Taken together this makes it easier to formulate phylogeographic hypotheses. This work investigated the phylogeography of Cakile maritima and Eryngium maritimum, two species growing in sandy habitats along the north Atlantic Ocean and the Mediterranean Sea coasts on two different spatial scales using AFLP data. The genetic structure of these species was investigated by sampling single individuals along most of their distributions from Turkey to south Sweden. On a regional scale the population genetic structure of both species was also studied in detail in the Bosporus and Dardanelles straits, the Strait of Gibraltar and along a continuous stretch of dunes in western France. Additionally, populations of C. maritima were investigated in the Baltic Sea/Kattegat/North Sea area. Over the complete sampling range the species show both differences and similarities in their genetic structure. In the Mediterranean Sea, both species contain Aegean Sea/Black Sea and west Mediterranean clusters. Cakile maritima additionally shows a clustering of Ionian Sea/Adriatic Sea collections. Further, both species show a subdivision of Atlantic Ocean/North Sea/Baltic Sea material from Mediterranean. Within the Atlantic Ocean group, C. maritima from the Baltic Sea and the most northern Atlantic localities form an additional cluster while no such substructure was found in E. maritimum. In all three instances where population genetic investigations of both species were performed in the same area, the results showed almost complete congruency of spatial genetic patterns. In the Aegean/Black Sea/Marmara region a subdivision of populations into a Black Sea, a Sea of Marmara and an Aegean Sea group is shared by both species. In addition the Sea of Marmara populations are more close to the Aegean Sea populations than they are to the Black Sea populations in both cases. Populations from the Atlantic side of the Strait of Gibraltar are differentiated from those on the Mediterranean side in both species, a pattern that confirms the results of the wide scale study. Along the dunes of West France no clear genetic structure could be detected in any of the species. Additionally, the results from the Baltic Sea/North Sea populations of C. maritima did not reveal any geographical genetic pattern. It is postulated that the many congruencies between the species are mainly due to a predominantly sea water mediated seed dispersal in both species and their shared sandy habitat. The results are compared to hypothetical distributions for the last glacial maximum based on species specific temperature requirements. It is argued that in both species the geographical borders of the clusters in the Mediterranean area were not affected by quaternary temperature changes and that the Aegean/Black Sea/Marmara cluster, and possibly the Ionian Sea/Adriatic Sea cluster in C. maritima, is the result of sea currents that isolate these basins from the rest of the sampled areas. The genetic gap in the Strait of Gibraltar between Atlantic Ocean and Mediterranean Sea populations in both species is also explained in terms of sea currents. The existence of three subgroups corresponding to the Aegean Sea, Black Sea and Sea of Marmara basins is suggested to have arisen due to geographical isolation during periods of global sea regressions in the glacials. The population genetic evidence was inconclusive regarding the Baltic Sea cluster of C. Maritima from the wide scale study. The results of this study are very similar to those of an investigation of three other coastal plant species over a similar range. This suggests that the phylo-geographic patterns of widespread coastal plants may be more predictable than those of other terrestrial plants.

Molecular phylogeography of the Woodland Ringlet (Erebia medusa [Denis and Schiffermüller] 1775) in Europe

Phylogeography is a recent field of biological research that links phylogenetics to biogeography through deciphering the imprint that evolutionary history has left on the genetic structure of extant populations. During the cold phases of the successive ice ages, which drastically shaped species’ distributions since the Pliocene, populations of numerous species were isolated in refugia where many of them evolved into different genetic lineages. My dissertation deals with the phylogeography of the Woodland Ringlet (Erebia medusa [Denis and Schiffermüller] 1775) in Central and Eastern Europe. This Palaearctic butterfly species is currently distributed from central France and south eastern Belgium over large parts of Central Europe and southern Siberia to the Pacific. It is absent from those parts of Europe with mediterranean, oceanic and boreal climates. It was supposed to be a Siberian faunal element with a rather homogeneous population structure in Central Europe due to its postglacial expansion out of a single eastern refugium. An already existing evolutionary scenario for the Woodland Ringlet in Central and Eastern Europe is based on nuclear data (allozymes). To know if this is corroborated by organelle evolutionary history, I sequenced two mitochondrial markers (part of the cytochrome oxydase subunit one and the control region) for populations sampled over the same area. Phylogeography largely relies on the construction of networks of uniparentally inherited haplotypes that are compared to geographic haplotype distribution thanks to recent developed methods such as nested clade phylogeographic analysis (NCPA). Several ring-shaped ambiguities (loops) emerged from both haplotype networks in E. medusa. They can be attributed to recombination and homoplasy. Such loops usually avert the straightforward extraction of the phylogeographic signal contained in a gene tree. I developed several new approaches to extract phylogeographic information in the presence of loops, considering either homoplasy or recombination. This allowed me to deduce a consistent evolutionary history for the species from the mitochondrial data and also adds plausibility for the occurrence of recombination in E. medusa mitochondria. Despite the fact that the control region is assumed to have a lack of resolving power in other species, I found a considerable genetic variation of this marker in E. medusa which makes it a useful tool for phylogeographic studies. In combination with the allozyme data, the mitochondrial genome supports the following phylogeographic scenario for E. medusa in Europe: (i) a first vicariance, due to the onset of the Würm glaciation, led to the formation of several major lineages, and is mirrored in the NCPA by restricted gene flow, (ii) later on further vicariances led to the formation of two sub-lineages in the Western lineage and two sub-lineages in the Eastern lineage during the Last Glacial Maximum or Older Dryas; additionally the NCPA supports a restriction of gene flow with isolation by distance, (iii) finally, vicariance resulted in two secondary sub-lineages in the area of Germany and, maybe, to two other secondary sub-lineages in the Czech Republic. The last postglacial warming was accompanied by strong range expansions in most of the genetic lineages. The scenario expected for a presumably Siberian faunal element such as E. medusa is a continuous loss of genetic diversity during postglacial westward expansion. Hence, the pattern found in this thesis contradicts a typical Siberian origin of E. medusa. In contrast, it corroboratess the importance of multiple extra-Mediterranean refugia for European fauna as it was recently assumed for other continental species.

Genetic structure and connectivity between populations of two common Mediterranean sessile invertebrates

Population genetic and phylogeography of two common mediterranean species were studied in 10 localities located on the coasts of Toscana, Puglia and Calabria. The aim of the study was to verify the extent of genetic breaks, in areas recognized as boundaries between Mediterranean biogeographic sectors. From about 100 sequences obtained from the mitochondrial Cytochrome Oxidase subunit I (COI) gene of Halocynthia papillosa and Hexaplex trunculus genetic diversity, genetic structure at small and large distances and demographic history of both specieswere analyzed. No evidences of genetic breaks were found for the two species in Toscana and Puglia. The genetic structure of H. trunculus evidences the extent of a barrier to gene flow localized in Calabria, which could be represented by the Siculo-Tunisian Strait and the Strait of Messina. The observed patterns showed similar level of gene flow at small distances in both species, although the two species have different larval ecology. These results suggest that other factors, such as currents, local dynamics and seasonal temperatures, influence the connectivity along the Italian peninsula. The geographic distribution of the haplotypes shows that H. papillosacould represent a single genetic pool in expansion, whereas H. trunculus has two distinct genetic pools in expansion. The demographic pattern of the two species suggests that Pleistocene sea level oscillations, in particular of the LGM, may have played a key role in shaping genetic structure of the two species. This knowledge provides basic information, useful for the definition of management plans, or for the design of a network of marine protected areas along the Italian peninsula.

Genetic characterisation and global comparison of ancient Mediterranean Great White Sharks population (Carcharodon carcharias, L.1758)

This study poses as its objective the genetic characterization of the ancient population of the Great White shark, Carcharodon carcharias, L.1758, present in the Mediterranean Sea. Using historical evidence, for the most part buccal arches but also whole, stuffed examples from various national museums, research institutes and private collections, a dataset of 18 examples coming from the Mediterranean Sea has been created, in order to increase the informations regarding this species in the Mediterranean. The importance of the Mediterranean provenance derives from the fact that a genetic characterization of this species' population does not exist, and this creates gaps in the knowledge of this species in the Mediterranean. The genetic characterization of the individuals will initially take place by the extraction of the ancient DNA and the analysis of the variations in the sequence markers of the mitochondrial DNA. This approach has allowed the genetic comparison between ancient populations of the Mediterranean and contemporary populations of the same geographical area. In addition, the genetic characterization of the population of white sharks of the Mediterranean, has allowed a genetic comparison with populations from global "hot spots", using published sequences in online databases (NCBI, GenBank). Analyzing the variability of the dataset, both in terms space and time, I assessed the evolutionary relationships of the Mediterranean population of Great Whites with the global populations (Australia/New Zealand, South Africa, Pacific USA, West Atlantic), and the temporal trend of the Mediterranean population variability. This method based on the sequencing of two portions of mitochondrial DNA genes, markers showed us how the population of Great White Sharks in the Mediterranean, is genetically more similar to the populations of the Australia Pacific ocean, American Pacific Ocean, rather than the population of South Africa, and showing also how the population of South Africa is abnormally distant from all other clusters. Interestingly, these results are inconsistent with the results from tagging of this species. In addition, there is evidence of differences between the ancient population of the Mediterranean with the modern one. This differentiation between the ancient and modern population of white shark can be the result of events impacting on this species occurred over the last two centuries.

Comparative phylogeography of two co-distributed arctic-alpine freshwater insect species in Europe

The distribution pattern of European arctic-alpine disjunct species is of growing interest among biogeographers due to the arising variety of inferred demographic histories. In this thesis I used the co-distributed mayfly Ameletus inopinatus and the stonefly Arcynopteryx compacta as model species to investigate the European Pleistocene and Holocene history of stream-inhabiting arctic-alpine aquatic insects. I used last glacial maximum (LGM) species distribution models (SDM) to derive hypotheses on the glacial survival during the LGM and the recolonization of Fennoscandia: 1) both species potentially survived glacial cycles in periglacial, extra Mediterranean refugia, and 2) postglacial recolonization of Fennoscandia originated from these refugia. I tested these hypotheses using mitochondrial sequence (mtCOI) and species specific microsatellite data. Additionally, I used future SDM to predict the impact of climate change induced range shifts and habitat loss on the overall genetic diversity of the endangered mayfly A. inopinatus.rnI observed old lineages, deep splits, and almost complete lineage sorting of mtCOI sequences between mountain ranges. These results support the hypothesis that both species persisted in multiple periglacial extra-Mediterranean refugia in Central Europe during the LGM. However, the recolonization of Fennoscandia was very different between the two study species. For the mayfly A. inopinatus I found strong differentiation between the Fennoscandian and all other populations in sequence and microsatellite data, indicating that Fennoscandia was recolonized from an extra European refugium. High mtCOI genetic structure within Fennoscandia supports a recolonization of multiple lineages from independent refugia. However, this structure was not apparent in the microsatellite data, consistent with secondary contact without sexual incompability. In contrast, the stonefly A. compacta exhibited low genetic structure and shared mtCOI haplotypes among Fennoscandia and the Black Forest, suggesting a shared Pleistocene refugium in the periglacial tundrabelt. Again, there is incongruence with the microsatellite data, which could be explained with ancestral polymorphism or female-biased dispersal. Future SDM projects major regional habitat loss for the mayfly A. inopinatus, particularly in Central European mountain ranges. By relating these range shifts to my population genetic results, I identified conservation units primarily in Eastern Europe, that if preserved would maintain high levels of the present-day genetic diversity of A. inopinatus and continue to provide long-term suitable habitat under future climate warming scenarios.rnIn this thesis I show that despite similar present day distributions the underlying demographic histories of the study species are vastly different, which might be due to differing dispersal capabilities and niche plasticity. I present genetic, climatic, and ecological data that can be used to prioritize conservation efforts for cold-adapted freshwater insects in light of future climate change. Overall, this thesis provides a next step in filling the knowledge gap regarding molecular studies of the arctic-alpine invertebrate fauna. However, there is continued need to explore the phenomenon of arctic-alpine disjunctions to help understand the processes of range expansion, regression, and lineage diversification in Europe’s high latitude and high altitude biota.

Temporal variation of genetic diversity and differentiation of albacore, Thunnus alalunga, in the Mediterranean Sea

This study is on albacore (Thunnus alalunga, Bonnaterre 1788), an epi- and mesopelagic oceanic tuna species cosmopolitan in the tropical and temperate waters of all oceans including the Mediterranean Sea, extending in a broad band between 40°N and 40°S. What it’s known about albacore population structure is based on different studies that used fisheries data, RFLP, mtDNA control region and nuDNA markers, blood lectins analysis, individual tags and microsatellite. At the moment, for T. alalunga six management units are recognized: the North Pacific, South Pacific, Indian, North Atlantic, South Atlantic and Mediterranean stocks. In this study I have done a temporal and spatial comparison of genetic variability between different Mediterranean populations of Thunnus alalunga matching an historical dataset ca. from 1920s composed of 43 individuals divided in 3 populations (NADR, SPAIN and CMED) with a modern dataset composed of 254 individuals and 7 populations (BAL, CYP, LIG, TYR, TUR, ADR, ALB). The investigation was possible using a panel of 94 nuclear SNPs, built specifically for the target species at the University of Basque Country UPV/EHU. First analysis done was the Hardy-Weinberg, then the number of clusters (K) was determined using STRUCTURE and to assess the genetic variability, allele frequencies, the average number of alleles per locus, expected (He) and observed (Ho) heterozygosis, and the index of polymorphism (P) was used the software Genetix. Historical and modern samples gives different results, showing a clear loss of genetic diversity over time leading to a single cluster in modern albacore instead of the two found in historical samples. What this study reveals is very important for conservation concerns, and additional research endeavours are needed.