Examining the role of parasites in limiting unidirectional gene flow between lake and river sticklebacks

Parasites are important selective agents with the potential to limit gene flow between host populations by shaping local host immunocompetence. We report on a contact zone between lake and river three-spined sticklebacks (Gasterosteus aculeatus) that offers the ideal biogeographic setting to explore the role of parasite-mediated selection on reproductive isolation. A waterfall acts as a natural barrier and enforces unidirectional migration from the upstream river stickleback population to the downstream river and lake populations. We assessed population genetic structure and parasite communities over four years. In a set of controlled experimental infections, we compared parasite susceptibility of upstream and downstream fish by exposing laboratory-bred upstream river and lake fish, as well as hybrids, to two common lake parasite species: a generalist trematode parasite, Diplostomum pseudospathaceum, and a host-specific cestode, Schistocephalus solidus. We found consistent genetic differentiation between upstream and downstream populations across four sampling years, even though the downstream river consisted of ~10% first-generation migrants from the upstream population as detected by parentage analysis. Fish in the upstream population had lower genetic diversity and were strikingly devoid of macroparasites. Through experimental infections, we demonstrated that upstream fish and their hybrids had higher susceptibility to parasite infections than downstream fish. Despite this, naturally sampled upstream migrants were less infected than downstream residents. Thus, migrants coming from a parasite-free environment may enjoy an initial fitness advantage, but their descendants seem likely to suffer from higher parasite loads. Our results suggest that adaptation to distinct parasite communities can influence stickleback invasion success and may represent a barrier to gene flow, even between close and connected populations.

Cryptic haplotype-specific gamete selection yields offspring with optimal MHC immune genes

Females choose specific mates in order to produce fitter offspring. However, several factors interfere with females' control over fertilization of their eggs, including sneaker males and phenotypically unpredictable allele segregation during meiosis. Mate choice at the individual level thus provides only a poor approximation for obtaining the best genetic match. Consequently, postcopulatory sperm selection by female oocytes has been proposed as a mechanism to achieve complementary combinations of parental haplotypes. Here, using controlled in vitro fertilization of three‐spined stickleback eggs, we find haplotype‐specific fertilization bias toward gametes with complementary major histocompatibility complex (MHC) immunogenes. The resulting zygote (and thus offspring) genotypes exhibit an intermediate level of individual MHC diversity that was previously shown to confer highest pathogen resistance. Our finding of haplotype‐specific gamete selection thus represents an intriguing mechanism for fine‐tuned optimization of the offspring's immune gene composition and an evolutionary advantage in the Red Queen dynamics of host‐parasite coevolution.

Transcriptome profiling of immune tissues reveals habitat-specific gene expression between lake and river sticklebacks

The observation of habitat-specific phenotypes suggests the action of natural selection. The three-spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to diverse freshwater habitats across the northern hemisphere since the last glaciation, while giving rise to recurring phenotypes associated with specific habitats. Parapatric lake and river populations of sticklebacks harbour distinct parasite communities, a factor proposed to contribute to adaptive differentiation between these ecotypes. However, little is known about the transcriptional response to the distinct parasite pressure of those fish in a natural setting. Here, we sampled wild-caught sticklebacks across four geographical locations from lake and river habitats differing in their parasite load. We compared gene expression profiles between lake and river populations using 77 whole-transcriptome libraries from two immune-relevant tissues, the head kidney and the spleen. Differential expression analyses revealed 139 genes with habitat-specific expression patterns across the sampled population pairs. Among the 139 differentially expressed genes, eight are annotated with an immune function and 42 have been identified as differentially expressed in previous experimental studies in which fish have been immune challenged. Together, these findings reinforce the hypothesis that parasites contribute to adaptation of sticklebacks in lake and river habitats.

Genomics of divergence along a continuum of parapatric population differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.

A variety of evolutionary forces influence the genomic landscape of divergence during ecological speciation. Here we characterize the evolution of genomic divergence patterns based on 60 fully sequenced three-spined stickleback genomes, contrasting lake and river populations that differ in parasite abundance. Our comparison of the size and abundance of divergent regions in the genomes across a continuum of population differentiation suggests that selection and the hitchhiking effect on neutral sites mainly contributes to the observed heterogeneous patterns of genomic divergence. Additional divergent regions of the genome can be explained by a local reduction of gene flow. Our description of genomic divergence patterns across a continuum of population differentiation combined with an analysis of molecular signatures of evolution highlights how adaptation shapes the differentiation of sticklebacks in freshwater habitats.

Genome-wide patterns of standing genetic variation in a marine population of three-spined sticklebacks

Since the end of the Pleistocene, the three-spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to various freshwater habitats probably originating from ancestral marine populations. Standing genetic variation and the underlying genomic architecture both have been speculated to contribute to recent adaptive radiations of sticklebacks. Here, we expand on the current genomic resources of this fish by providing extensive genome-wide variation data from six individuals from a marine (North Sea) stickleback population. Using next-generation sequencing and a combination of paired-end and mate-pair libraries, we detected a wide size range of genetic variation. Among the six individuals, we found more than 7% of the genome is polymorphic, consisting of 2599111 SNPs, 233464 indels and structural variation (SV) (>50 bp) such as 1054 copy-number variable regions (deletions and duplications) and 48 inversions. Many of these polymorphisms affect gene and coding sequences. Based on SNP diversity, we determined outlier regions concordant with signatures expected under adaptive evolution. As some of these outliers overlap with pronounced regions of copy-number variation, we propose the consideration of such SV when analysing SNP data from re-sequencing approaches. We further discuss the value of this resource on genome-wide variation for further investigation upon the relative contribution of standing variation on the parallel evolution of sticklebacks and the importance of the genomic architecture in adaptive radiation.

Gene flow between species of Lake Victoria haplochromine fishes

The haplochromine cichlid fishes of Lake Victoria (LV), East Africa, are a textbook example of adaptive radiation-a rapid divergence of multiple morphologically distinguishable forms from a few founding lineages. The forms are generally believed to constitute a "flock" of several hundred reproductively isolated species in a dozen or so genera. This belief has, until now, not been subjected to a test, however. Here, we compare genetic variation at 11 loci in 10 haplochromine populations of 6 different species. Although the genetic diversity in the populations is quite high, using a variety of statistical tests, we find no evidence of genetic differentiation among the populations of LV haplochromines. On genetic distance trees, populations of the same species intermingle with those of different species. At the molecular level, the species are indistinguishable from one another. Genetic comparisons with closely related species in 2 crater lakes indicate that the species within LV continue exchanging genes. These observations have important implications for phylogenetic reconstruction. The approach used in this study is applicable to other instances of adaptive radiation.

Phylogeny of Antigen-Processing Enzymes: Cathepsins of a Cephalochordate, an Agnathan and a Bony Fish

Cathepsins are enzymes that have been cleaving peptide bonds of lysosomal proteins probably since lysosomes appeared in early eucaryotes. When the adaptive system emerged in gnathostomes, cathepsins were recruited to produce peptides for loading onto the major histocompatibility complex class II molecules and for degrading the class II-associated invariant chain just before the loading. The circumstances under which this recruitment took place are unclear because the knowledge about vertebrate cathepsins is limited largely to mammals. To shed light on the recruitment, 10 amphioxus, one lamprey and one cichlid fish cathepsin cDNA clone were characterized and analysed phylogenetically. Disregarding cathepsin O, whose phylogenetic position is uncertain, the analysis confirms the existence of two old lines of descent, the B and the L lineages of cathepsins, which diverged from each other early in the evolution of eucaryotes. The B lineage encompasses cathepsins B, C and Z (X). The L lineage splits off sublineages encompassing cathepsins F and W before the plant-animal separation and cathepsin H early in the evolution of the metazoa. The remaining cathepsins belonging to the L lineage diverged from one another during the evolution of vertebrates: S, K and L before the emergence of bony fishes, and the group of rodent placentally expressed cathepsins [J (P), M, Q, R, 3, 6, 7 and 8] as well as the testis/ova-expressed cathepsins (testins) probably after the divergence of rodents from primates. The part possibly played by the adaptive immune system in some of these divergences is discussed.

Identification of chemokines and a chemokine receptor in cichlid fish, shark, and lamprey

Chemokines are small, inducible, structurally related proteins that guide cells expressing the right chemokine receptors to sites of immune response. They have been identified and studied extensively in mammals, but little is known about their presence in other vertebrate groups. Here we describe seven new chemokines in bony fish and one in a cartilaginous fish, as well as one chemokine receptor in a jawless vertebrate. All eight chemokines belong to the SCYA (CC) subfamily characterized by four conserved cysteine residues of which the first two are adjacent. The chemokine receptor is of the CXCR4 type. Phylogenetic analysis does not reveal any clear evidence of orthology of fish and human chemokines. Although the divergence of the subfamilies began before the fish-tetrapod split, much of the divergence within the subfamilies took place separately in the two vertebrate groups. The existence of a chemokine receptor in the lamprey indicates that chemokines are apparently also present in the Agnatha.

Linkage relationships of genes coding for alpha2-macroglobulin, C3 and C4 in the zebrafish: implications for the evolution of the complement and Mhc systems

The alpha2-macroglobulin (A2M) and the complement components C3 and C4 are related proteins derived from a common ancestor. Theoretically, this derivation could have occurred either by tandem duplications of their encoding genes or by polyploidization involving chromosomal segments, a chromosome or the whole genome. In tetrapods the A2M-, C3- and C4-encoding genes are generally each located on a different chromosome. This observation has been interpreted as supporting their origin by polyploidization. We identified and mapped (with the help of a radiation hybrid panel of cell lines) the A2M, C3 and C4 loci in the zebrafish, Danio rerio. Each of the three types of loci is present in the zebrafish in multiple copies, but all of the identified copies of a given type map to the same region in linkage groups 1 (C3) and 15 (A2M, C4). The A2M and C4 loci are mapped in the same region not linked to any of the class I or class II major histocompatibility complex (Mhc) loci. These observations are interpreted as supporting the origin of the A2M family of genes by tandem duplications, followed by the dispersal of the copies to different chromosomes. It is also argued that the association of C4 with the class I/II loci in tetrapods is accidental and without functional significance.

Classification and phylogenetic relationships of African tilapiine fishes inferred from mitochondrial DNA sequences

African cichlid fishes are composed of two major lineages, the haplochromines and the tilapiines. Whereas the phylogenetic relationships of the haplochromines have been studied extensively, primarily because of their spectacular adaptive radiations in the Great Lakes of East Africa, little is known about the relationships among the tilapiine species, despite the fact that they have become an important component of African, indeed world, aquaculture. To remedy this situation, molecular phylogenetic analysis of tilapiine fishes was undertaken. A segment of mitochondrial DNA encompassing the terminal part of the tRNA(Pro) gene and the most variable part of the control region was amplified by the polymerase chain reaction with DNA samples isolated from 42 tilapiine species, and the amplification products were subjected to heteroduplex analysis and sequencing. Phylogenetic trees based on 68 sequences revealed the existence of 11 sequence groups and 11 single-sequence branches. The groups, designated Ti1 through Ti11, were distinguished by specific combinations of diagnostic substitutions, formation of monophyletic clusters, and separation by genetic distances in excess of 0.04. Although the relationships among the groups could not be resolved, the sequences separated Oreochromis and Sarotherodon from Tilapia, as defined by Trewavas. The Oreochromis sequences clustered with the Sarotherodon sequences and thus supported the hypothesis that the mouthbrooding behavior of the tilapiine fishes evolved only once from the substrate-spawning behavior. Since on phylogenetic trees the O. alcalicus (sub)species were always separated from O. amphimelas by other Oreochromis species, it was concluded that the adaptation to life in water with a high salt concentration and high pH values evolved independently at least twice in the tilapiine fishes. The tilapiines diverged from the haplochromines more than 8 million years ago; most of the intragroup divergences among the tilapiines took place an estimated 1.1 to 6 million years ago.

Molecular phylogeny of Philippine freshwater sardines based on mitochondrial DNA analysis

The commercially important Sardinella species (family Clupeidae or herrings) usually thrive in marine environments. An exception is Sardinella tawilis of Taal Lake, Batangas, Philippines, the only known freshwater sardine. This species is believed to have immigrated from Balayan Bay to the lake when it was formed in the course of volcanic eruptions some 240 years ago. To determine the relationship of S. tawilis to the marine species S. albella, S. fimbriata, and S. longiceps from the Balayan Bay we sequenced 358 bp of the cytochrome b gene and the mitochondrial control region. The cytochrome b gene was highly conserved and contained little phylogenetic information. The control region sequences, however, demonstrated two highly diversified main haplotypes grouping S. tawilis with S. albella, as shown by maximum parsimony and neighbor-joining analysis. The haplotypes are characterized by the presence of an 81 bp indel and up to eight 35 bp tandem repeat elements. The repeat copy number varied within individuals of S. tawilis and S. albella, thus showing heteroplasmy in these two species only. The analysis of two subpopulations of S. tawilis revealed restricted substitutions that may indicate the beginning of genetic differentiation of the two subpopulations.