Ciro Rico

470 million years of conservation of microsatellite loci among fish species

Primers for 18 microsatellite loci originally isolated from whiting (6), stickleback (Gasterosteus aculeatus, n — 5) and cod (Gadus morhua, n = 7) were tested across a panel of diverse fish species, representing the three principal superclasses and most principal superorders of fish, to examine conservation of microsatellite regions across distantly related taxa. Three methods were used. First, amplified fragments were analysed by Southern blotting using the relevant microsatellite motif probes. A total of 17 of the tested primer pairs gave a product in the expected size range in at least four of 11 tested species. Second, for two study loci the amplified polymerase chain reaction products were cloned and sequenced in five fish species to reveal a high level of conservation of the flanking and microsatellite sequences. Finally, the 17 loci successfully amplified in non-source species were tested for polymorphism in groups of unrelated individuals from nine species, in several cases revealing extensive polymorphism. Levels of polymorphism were generally high in species from which the loci were derived or among closely related species. The conservation of flanking sequences for particular microsatellite motifs over the span of fish evolution represented in the test species (470 million years) far exceeds that hitherto reported and lends support to the suggestion (derived from studies of whales and marine turtles) that the rate of base substitution in nuclear and mitochondrial sequences is lower in aquatic than terrestrial organisms. A further explanation could be that these sequences, although generally considered neutral, may play an important role in eukaryotic genomes, and may be under strong selective constraints. The study suggests that heterologous primers will be a ready source of polymorphic markers among fish species, but also indicates that caution should be used in cross-species comparisons of variability.

Unusually fine–scale genetic structuring found in rapidly speciating Malawi cichlid fishes

Mechanisms behind the explosive radiation of over 500 cichlid fish species from a single founding population in Lake Malawi during the last 700 000 years are poorly understood. Recent studies have suggested that the degree of population subdivision among the habitat patches within the lake may be responsible, but the evidence has been circumstantial: lack of a dispersal stage in haplochromine cichlids; genetic and colour variation among populations separated by large–scale geographical barriers; and fluctuating lake levels. One reason for the rapidity of speciation in these fishes may be that population subdivision is on a much finer scale than previously thought. Here we quantify the level of population subdivision and estimate migration at a scale of 700–1400 m, in order to investigate whether cichlid populations are sufficiently isolated from each other for allopatric divergence and perhaps speciation to take place. Using six microsatellite loci, we demonstrate the existence of highly significant genetic differentiation between subpopulations on adjacent headlands in each of four rock–dwelling haplochromine cichlid species. Our results suggest that these fish populations are divided into thousands of subunits among which genetic divergence is currently occurring, and that this may provide unprecedented opportunities for allopatric speciation.

Assortative mating among rock-dwelling cichlid fishes supports high estimates of species richness from Lake Malawi

It has been estimated that Lake Malawi, Africa, contains 500–650 endemic species of cichlid fishes, the largest number of vertebrate species endemic to any comparable sized area on the planet. As many of these putative species cannot be distinguished anatomically, these estimates of species richness depend to a great extent on the assumption that sympatrically occurring male colour morphs represent biological species. We have tested this assumption using a combination of behavioural observations of courtship and microsatellite DNA analysis for six putative species of the Pseudotropheus (Tropheops) complex and three of the Pseudotropheus (Maylandia) complex occurring sympatrically at Nkhata Bay. We were unable to demonstrate assortative courtship for the species pairs Pseudotropheus (Maylandia) zebra/P. ‘gold zebra’ or P. (Tropheops) ‘band’/P. (T.) ‘rust’ because we were unable to distinguish between the females of these taxa. All other taxa showed clear assortative courtship, except for P. (T.) ‘deep’, a deep‐water species which was rarely observed. Fixation indices (θST for the infinite allele model, and RST for the stepwise mutation model) calculated from six microsatellite DNA loci demonstrated significant deviations from panmixia in all pairwise comparisons of putative species, indicating little or no gene flow between populations. All taxa showed high levels of allelic diversity providing evidence that genetic bottlenecking may have been of limited importance in the speciation process. Assortative mating among taxa differing only in male colouration is supportive of theories that speciation in these fishes has been driven by sexual selection by female choice.

Disentangling vector-borne transmission networks: a universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals.

Emerging infectious diseases represent a challenge for global economies and public health. About one fourth of the last pandemics have been originated by the spread of vector-borne pathogens. In this sense, the advent of modern molecular techniques has enhanced our capabilities to understand vector-host interactions and disease ecology. However, host identification protocols have poorly profited of international DNA barcoding initiatives and/or have focused exclusively on a limited array of vector species. Therefore, ascertaining the potential afforded by DNA barcoding tools in other vector-host systems of human and veterinary importance would represent a major advance in tracking pathogen life cycles and hosts. Here, we show the applicability of a novel and efficient molecular method for the identification of the vertebrate hosts DNA contained in the midgut of blood-feeding arthropods. To this end, we designed a eukaryote-universal forward primer and a vertebrate-specific reverse primer to selectively amplify 758 base pairs (bp) of the vertebrate mitochondrial Cytochrome c Oxidase Subunit I (COI) gene. Our method was validated using both extensive sequence surveys from the public domain and Polymerase Chain Reaction (PCR) experiments carried out over specimens from different Classes of vertebrates (Mammalia, Aves, Reptilia and Amphibia) and invertebrate ectoparasites (Arachnida and Insecta). The analysis of mosquito, culicoid, phlebotomie, sucking bugs, and tick bloodmeals revealed up to 40 vertebrate hosts, including 23 avian, 16 mammalian and one reptilian species. Importantly, the inspection and analysis of direct sequencing electropherograms also assisted the resolving of mixed bloodmeals. We therefore provide a universal and high-throughput diagnostic tool for the study of the ecology of haematophagous invertebrates in relation to their vertebrate hosts. Such information is crucial to support the efficient management of initiatives aimed at reducing epidemiologic risks of arthropod vector-borne pathogens, a priority for public health.

Macrogeographical population differentiation in oceanic environments: a case study of European hake (Merluccius merluccius), a commercially important fish

Polymorphism at six microsatellite loci was used to study genetic variability and population structure in six geographically distant natural populations of European hake (Merluccius merluccius L.). Four hundred and eighty‐three individuals were sampled from Trondheimsfjord in Norway, the Celtic Sea, the southern Bay of Biscay, Faro off Portugal, the Mediterranean Sea north of the coast of Tunisia and the Adriatic Sea. Population subdivision was found between Mediterranean and Atlantic samples, θ = 0.029 (P < 0.001). No substructuring was found between samples within the Mediterranean Sea, θ = 0.003 and RST = 0.007 (P > 0.05). The Atlantic population structure appears to be more complex than previously suggested by the placement of stock boundaries by the International Council for the Exploration of the Seas (ICES). Analyses based on various models of microsatellite evolution all suggest that differentiation exists between Bay of Biscay and Portugese samples, θ = 0.013 (P < 0.001), RST = 0.036 (P < 0.001) which are currently managed as one stock. By contrast, fixation indices indicated no differentiation between southern Bay of Biscay samples and Celtic Sea samples, θ = 0.003 (P = 0.02), φST = 0.007 (P = 0.10) which are managed as separate stocks. These results suggest that if the observed trends are stable through time, current management policy of European hake may need revision.

Temporal and spatial genetic variation in spawning grounds of European hake (Merluccius merluccius) in the Bay of Biscay.

Polymorphism at five microsatellite loci were screened to determine the genetic variability and the temporal stability of population structure in natural populations of European hake (Merluccius merluccius, L.) within the Bay of Biscay. In addition, the control region (900 bp) and two protein coding genes (ATPase, subunits 6 and 8, 842 bp and a partial sequence of the ND1, 800 bp) of the mitochondrial DNA (mtDNA) were sequenced from geographically distant populations from the extremes of the species range. One hundred individuals from either side of a supposed stock boundary within the bay were collected in autumn 1997. This sampling strategy was repeated during hake spawning seasons in late spring of 1998 and 1999. Low levels of population subdivision were found between putative populations within years. Similarly, low levels of differentiation were found between autumn 1997 northern samples and spring 1998 southern samples which were collected 7 months later on spawning grounds. These results are discussed in relation to ecological, behavioural and oceanographic information. Sampling effects, which may influence these results, are also discussed. Theta (θ) estimates were significantly different from zero in every other pairwise comparison between geographical areas (north and south of the Bay of Biscay) and between years within the same area (P < 0.05). Hierarchical analysis of molecular variance (amova) does not confirm the temporal persistence of population structure. These results are discussed in relation to variance in reproductive success, and temporal spawning patterns, which may exist within the bay. mtDNA variability was very low between geographically distant samples from Norway and the Mediterranean Sea with only 10 variable sites found in a total of 2542 bp of mtDNA, these differences being exclusively in the d‐loop.

MHC adaptive divergence between closely related and sympatric African cichlids.

Background The haplochromine cichlid species assemblages of Lake Malawi and Victoria represent some of the most important study systems in evolutionary biology. Identifying adaptive divergence between closely-related species can provide important insights into the processes that may have contributed to these spectacular radiations. Here, we studied a pair of sympatric Lake Malawi species, Pseudotropheus fainzilberi and P. emmiltos, whose reproductive isolation depends on olfactory communication. We tested the hypothesis that these species have undergone divergent selection at MHC class II genes, which are known to contribute to olfactory-based mate choice in other taxa. Methodology/Principal Findings Divergent selection on functional alleles was inferred from the higher genetic divergence at putative antigen binding sites (ABS) amino acid sequences than at putatively neutrally evolving sites at intron 1, exon 2 synonymous sequences and exon 2 amino acid residues outside the putative ABS. In addition, sympatric populations of these fish species differed significantly in communities of eukaryotic parasites. Conclusions/Significance We propose that local host-parasite coevolutionary dynamics may have driven adaptive divergence in MHC alleles, influencing odor-mediated mate choice and leading to reproductive isolation. These results provide the first evidence for a novel mechanism of adaptive speciation and the first evidence of adaptive divergence at the MHC in closely related African cichlid fishes.

Evidence for male-biased dispersal in lake malawi cichlids from microsatellites

This study addressed within‐population dispersal patterns among the species‐rich Lake Malawi cichlids, specifically among the rock‐dwelling mbuna group. Relatedness values were calculated for 160 individuals belonging to two species from known locations in the field by screening six highly polymorphic microsatellite loci. The results from both species indicate that spatially adjacent females have higher average relatedness values than those separated by larger distances, but that this pattern is reversed in males. This therefore provides firm evidence for male‐biased dispersal within the Malawian cichlid flock.

Phylogeography and local endemism of the native Mediterranean brine shrimp Artemia salina (Branchiopoda: Anostraca)

There has been a recent appreciation of the ecological impacts of zooplanktonic species invasions. The North American brine shrimp Artemia franciscana is one such alien invader in hyper‐saline water ecosystems at a global scale. It has been shown to outcompete native Artemia species, leading to their local extinction. We used partial sequences of the mitochondrial Cytochrome c Oxidase Subunit 1 (COI or cox1) gene to investigate the genetic diversity and phylogeography of A. salina, an extreme halophilic sexual brine shrimp, over its known distribution range (Mediterranean Basin and South Africa) and to assess the extent of local endemism, the degree of population structure and the potential impact of traditional human saltpan management on this species. We also examined the phylogenetic relationships in the genus Artemia using COI sequences. Our results show extensive regional endemism and indicate an early Pleistocene expansion of A. salina in the Mediterranean Basin. Subsequent population isolation in a mosaic of Pleistocene refugia is suggested, with two or three refugia located in the Iberian Peninsula. Two instances of long‐distance colonization were also observed. Surprisingly, given its strong phylogeographical structure, A. salina showed a signature of correlation between geographical and genetic distance. Owing to strong ‘priority effects’, extensive population differentiation is retained, despite dispersal via migrant birds and human management of saltpans. The foreseeable expansion of A. franciscana is likely to be followed by substantial loss of genetic diversity in Mediterranean A. salina. Large genetic divergences between Mediterranean and South African A. salina suggest that the latter deserves species status.

Genetic mosaic in a marine species flock.

We used molecular approaches to study the status of speciation in coral reef fishes known as hamlets (Serranidae: Hypoplectrus). Several hamlet morphospecies coexist on Caribbean reefs, and mate assortatively with respect to their strikingly distinct colour patterns. We provide evidence that, genetically, the hamlets display characteristics common in species flocks on land and in freshwaters. Substitutions within two mitochondrial DNA (mtDNA) protein‐coding genes place hamlets within a monophyletic group relative to members of two related genera (Serranus and Diplectrum), and establish that the hamlet radiation must have been very recent. mtDNA distances separating hamlet morphospecies were slight (0.6 ± 0.04%), yielding a coalescent estimate for the age of the hamlet flock of approximately 430 000 years. Morphospecies did not sort into distinct mtDNA haplotype phylogroups, and alleles at five hypervariable microsatellite loci were shared broadly across species boundaries. None the less, molecular variation was not distributed at random. Analyses of mtDNA haplotype frequencies and nested clades in haplotype networks revealed significant genetic differences between geographical regions and among colour morphospecies. We also observed significant microsatellite differentiation between geographical regions and in Puerto Rico, among colour morphospecies; the latter providing evidence for reproductive isolation between colour morphospecies at this locale. In our Panama collection, however, colour morphospecies were mostly genetically indistinguishable. This mosaic pattern of DNA differentiation implies a complex interaction between population history, mating behaviour and geography and suggests that porous boundaries separate species in this flock of brilliantly coloured coral reef fishes.