Efficient isolation of polymorphic microsatellites from high-throughput sequence data based on number of repeats.

Transcriptome data are a good resource to develop microsatellites due to their potential in targeting candidate genes. However, developing microsatellites can be a time-consuming enterprise due to the numerous primer pairs to be tested. Therefore, the use of methodologies that make it efficient to identify polymorphic microsatellites is desirable. Here we used a 62,038 contigs transcriptome assembly, obtained from pyrosequencing a peacock blenny (Salaria pavo) multi-tissue cDNA library, to mine for microsatellites and in silico evaluation of their polymorphism. A total of 4190 microsatellites were identified in 3670 unique unigenes, and from these microsatellites, in silico polymorphism was detected in 733. We selected microsatellites based either on their in silico polymorphism and annotation results or based only on their number of repeats. Using these two approaches, 28 microsatellites were successfully amplified in twenty-six individuals, and all but 2 were found to be polymorphic, being the first genetic markers for this species. Our results showed that the strategy of selection based on number of repeats is more efficient in obtaining polymorphic microsatellites than the strategy of in silico polymorphism (allelic richness was 8.2±3.85 and 4.56±2.45 respectively). This study demonstrates that combining the knowledge of number of repeats with other predictors of variability, for example in silico microsatellite polymorphism, improves the rates of polymorphism, yielding microsatellites with higher allelic richness, and decreases the number of monomorphic microsatellites obtained.

Unexpected high genetic diversity at the extreme northern geographic limit of Taurulus bubalis (Euphrasen, 1786).

The longspined bullhead (Taurulus bubalis, Euphrasen 1786) belongs to the family Cottidae and is a rocky shore species that inhabits the intertidal zones of the Eastern Atlantic since Iceland, southward to Portugal and also the North Sea and Baltic, northward to the Gulf of Finland, with some occurrences in the northern Mediterranean coasts eastward to the Gulf of Genoa. We analysed the phylogeographic patterns of this species using mitochondrial and nuclear markers in populations throughout most of its distributional range in west Europe. We found that T. bubalis has a relatively shallow genealogy with some differentiation between Atlantic and North Sea. Genetic diversity was homogeneous across all populations studied. The possibility of a glacial refugium near the North Sea is discussed. In many, but not all, marine temperate organisms, patterns of diversity are similar across the species range. If this phenomenon proves to be most common in cold adapted species, it may reflect the availability of glacial refugia not far from their present-day northern limits.

Northern refugia and recent expansion in the North Sea: the case of the wrasse Symphodus melops (Linnaeus, 1758).

Pleistocene climate changes have imposed extreme conditions to intertidal rocky marine communities, forcing many species to significant range shifts in their geographical distributions. Phylogeographic analyses based on both mitochondrial and nuclear genetic markers provide a useful approach to unravel phylogeographic patterns and processes of species after this time period, to gain general knowledge of how climatic changes affect shifts in species distributions. We analyzed these patterns on the corkwing wrasse (Symphodus melops, Labridae), a rocky shore species inhabiting North Sea waters and temperate northeastern Atlantic Ocean from Norway to Morocco including the Azores, using a fragment of the mitochondrial control region and the first intron of the nuclear S7 ribosomal protein gene. We found that S. melops shows a clear differentiation between the Atlantic and the Scandinavian populations and a sharp contrast in the genetic diversity, high in the south and low in the north. Within each of these main geographic areas there is little or no genetic differentiation. The species may have persisted throughout the last glacial maximum in the southern areas as paleotemperatures were not lower than they are today in North Scandinavia. The North Sea recolonization most likely took place during the current interglacial and is dominated by a haplotype absent from the south of the study area, but present in Plymouth and Belfast. The possibility of a glacial refugium in or near the English Channel is discussed.

Lusitanian toadfish song reflects male quality.

Lusitanian toadfish males that provide parental care rely on acoustic signals (the boatwhistle) to attract females to their nest. We test the hypothesis that male quality, namely male size and condition that are relevant for parental success, is reflected in vocal activity and boatwhistle characteristics and thus advertised to females. We recorded 22 males over a week during the peak of the breeding season. Calling rate and calling effort (percentage of time spent calling) strongly reflected male condition (lipid content of somatic muscles) and to a smaller extent sonic muscle hypertrophy and larger gonads. Males in better condition (increased body lipid and relative higher liver mass) also contracted the sonic muscles at faster rate as shown by the shorter boatwhistle pulse periods. Amplitude modulation reflected the degree of sonic muscle hypertrophy. None of the measured male quality parameters were good predictors of boatwhistle duration and dominant frequency. Altogether this study strongly suggests that Lusitanian toadfish males advertise their quality to females primarily with boatwhistle calling rate and calling effort, which mainly reflect male condition. Because pulse period had low variability, consistent with the existence of a vocal central pattern generator, we suggest that males that sustain sonic muscles contraction at a very fast rate close to their physiological limit may be honestly advertising their quality (condition). Similarly, males that produce boatwhistles with higher amplitude modulation, a feature that seems dependent on sonic muscle hypertrophy, could be more attractive to females.

Genetic divergence in the Atlantic-Mediterranean Montagu's blenny, Coryphoblennius galerita (Linnaeus 1758) revealed by molecular and morphological characters.

Coryphoblennius galerita is a small intertidal fish with a wide distribution and limited dispersal ability, occurring in the northeastern Atlantic and Mediterranean. In this study, we examined Atlantic and Mediterranean populations of C. galerita to assess levels of genetic divergence across populations and to elucidate historical and contemporary factors underlying the distribution of the genetic variability. We analyse three mitochondrial and one nuclear marker and 18 morphological measurements. The combined dataset clearly supports the existence of two groups of C. galerita: one in the Mediterranean and another in the northeastern Atlantic. The latter group is subdivided in two subgroups: Azores and the remaining northeastern Atlantic locations. Divergence between the Atlantic and the Mediterranean can be the result of historical isolation between the populations of the two basins during the Pleistocene glaciations. Present-day barriers such as the Gibraltar Strait or the 'Almeria-Oran jet' are also suggested as responsible for this isolation. Our results show no signs of local extinctions during the Pleistocene glaciations, namely at the Azores, and contrast with the biogeographical pattern that has been observed for Atlantic-Mediterranean warm-water species, in which two groups of populations exist, one including the Mediterranean and the Atlantic coast of western Europe, and another encompassing the western tropical coast of Africa and the Atlantic islands of the Azores, Madeira and Canaries. Species like C. galerita that tolerate cooler waters, may have persisted during the Pleistocene glaciations in moderately affected locations, thus being able to accumulate genetic differences in the more isolated locations such as the Azores and the Mediterranean. This study is one of the first to combine morphological and molecular markers (mitochondrial and nuclear) with variable rates of molecular evolution to the study of the relationships of the Atlantic and Mediterranean populations of a cool-water species.

Historical population dynamics and demography of the Eastern Atlantic pomacentrid Chromis limbata (Valenciennes, 1833).

Recent studies have focused on the relationship between the marine fauna of the Eastern Atlantic and the Mediterranean Sea, but within the Atlantic, little is known about genetic relationships between populations of the Macaronesian islands. In this study, we tested whether the paleo-climatology and paleo-oceanography of the region could predict the genetic relationships among three Eastern Atlantic populations (Azores, Madeira, and Canaries) of a damselfish, Chromis limbata, and compared our results with its Mediterranean and adjacent Atlantic sister species, Chromis chromis. We combined phylogeographic and coalescent approaches using the fast evolving mitochondrial control region gene. No population structure was found for the three archipelagos. The coalescence time estimated for C. limbata (0.857-1.17 Mya) was much greater than that estimated for C. chromis. We propose that this difference reflects differences in glaciating extents in the Northeastern Atlantic and the Mediterranean. Diversity indexes (Hd and genetic distances) together with historical demographic parameters of C. limbata (Theta and g) revealed a more stable population history when compared to C. chromis. Our results suggest that the Macaronesian populations of C. limbata have probably been less affected by the last glaciation than the Mediterranean populations of C. chromis. Migration across the three archipelagos was estimated and a prevailing northwest trend was detected. This result supports the idea of a colonization of the Azores by warm water fish from Madeira or the westernmost Canary islands which acted as major glacial refugia for the tropical and subtropical marine fauna during the glaciations.

Historical colonization and demography of the Mediterranean damselfish, Chromis chromis.

The desiccation of the Mediterranean Sea during the Messinian Salinity Crisis 6.0-5.3 million years ago (Ma), caused a major extinction of the marine ichthyofauna of the Mediterranean. This was followed by an abrupt replenishment of the Mediterranean from the Atlantic after the opening of the Strait of Gibraltar. In this study, we combined demographic and phylogeographic approaches using mitochondrial and nuclear DNA markers to test the alternative hypotheses of where (Atlantic or Mediterranean) and when (before or after the Messinian Salinity Crisis) speciation occurred in the Mediterranean damselfish, Chromis chromis. The closely related geminate transisthmian pair Chromis multilineata and Chromis atrilobata was used as a way of obtaining an internally calibrated molecular clock. We estimated C. chromis speciation timing both by determining the time of divergence between C. chromis and its Atlantic sister species Chromis limbata (0.93-3.26 Ma depending on the molecular marker used, e.g. 1.23-1.39 Ma for the control region), and by determining the time of coalescence for C. chromis based on mitochondrial control region sequences (0.14-0.21 Ma). The time of speciation of C. chromis was always posterior to the replenishment of the Mediterranean basin, after the Messinian Salinity Crisis. Within the Mediterranean, C. chromis population structure and demographic characteristics revealed a genetic break at the Peloponnese, Greece, with directional and eastbound gene flow between western and eastern groups. The eastern group was found to be more recent and with a faster growing population (coalescent time = 0.09-0.13 Ma, growth = 485.3) than the western group (coalescent time = 0.13-0.20 Ma, growth = 325.6). Our data thus suggested a western origin of C. chromis, most likely within the Mediterranean. Low sea water levels during the glacial periods, the hydrographic regime of the Mediterranean and dispersal restriction during the short pelagic larval phase of C. chromis (18-19 days) have probably played an important role in C. chromis historical colonization.

Heterozygous indels as useful tools in the reconstruction of DNA sequences and in the assessment of ploidy level and genomic constitution of hybrid organisms.

In this paper we describe a simple approach using double peaks in chromatograms generated as artefacts in the vicinity of heterozygous indels, to identify the specific sequences present in individual strands of a given DNA fragment. This method is useful to assign bases in individuals that are heterozygous at multiple sites. In addition, the relative sizes of the double peaks help to determine the ploidy level and the relative contribution of the parental genomes in hybrids. Our interpretation was confirmed with the analysis of artificial mixtures of DNA of two different species. Results were robust with varying PCR and sequencing conditions. The applicability of this method was demonstrated in hybrids of the Squalius alburnoides complex and in heterozygotes of Chondrostoma oligolepis. Far from being limited to these fish models and the gene where it was tested (beta-actin), this sequence reconstruction methodology is expected to have a broader application.