Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation.

Although population genomic studies using next generation sequencing (NGS) data are becoming increasingly common, studies focusing on phylogenetic inference using these data are in their infancy. Here, we use NGS data generated from reduced representation genomic libraries of restriction-site-associated DNA (RAD) markers to infer phylogenetic relationships among 16 species of cichlid fishes from a single rocky island community within Lake Victoria's cichlid adaptive radiation. Previous attempts at sequence-based phylogenetic analyses in Victoria cichlids have shown extensive sharing of genetic variation among species and no resolution of species or higher-level relationships. These patterns have generally been attributed to the very recent origin (<15,000 years) of the radiation, and ongoing hybridization between species. We show that as we increase the amount of sequence data used in phylogenetic analyses, we produce phylogenetic trees with unprecedented resolution for this group. In trees derived from our largest data supermatrices (3 to >5.8 million base pairs in width), species are reciprocally monophyletic with high bootstrap support, and the majority of internal branches on the tree have high support. Given the difficulty of the phylogenetic problem that the Lake Victoria cichlid adaptive radiation represents, these results are striking. The strict interpretation of the topologies we present here warrants caution because many questions remain about phylogenetic inference with very large genomic data set and because we can with the current analysis not distinguish between effects of shared ancestry and post-speciation gene flow. However, these results provide the first conclusive evidence for the monophyly of species in the Lake Victoria cichlid radiation and demonstrate the power that NGS data sets hold to resolve even the most difficult of phylogenetic challenges.

When phenotypes do not match genotypes--unexpected phenotypic diversity and potential environmental constraints in Icelandic stickleback.

Divergent lateral plate phenotypes in stickleback represent one of only a few cases known, where a single gene underlies the phenotype under divergent selection between different habitats. However, the selection pressures leading to the repeated loss of lateral plates in freshwater are still not well understood. By genotyping 838 individuals from 9 independently colonized lakes and 1 marine population in Iceland, we found 1) that only in some lakes are phenotypes associated with the expected genotype and 2) that the independent repeated occurrence of a rarely described plate phenotype is expressed in the absence of an allele that is usually associated with this phenotype. This suggests that either other genes such as modifiers might be under divergent selection between lakes or that lateral plate expression in these populations is restricted due to environmental constraints.

Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland.

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

Parallel amino acid replacements in the rhodopsins of the rockfishes (Sebastes spp.) associated with shifts in habitat depth.

Among various groups of fishes, a shift in peak wavelength sensitivity has been correlated with changes in their photic environments. The genus Sebastes is a radiation of marine fish species that inhabit a wide range of depths from intertidal to over 600 m. We examined 32 species of Sebastes for evidence of adaptive amino acid substitution at the rhodopsin gene. Fourteen amino acid positions were variable among these species. Maximum likelihood analyses identify several of these to be targets of positive selection. None of these correspond to previously identified critical amino acid sites, yet they may in fact be functionally important. The occurrence of independent parallel changes at certain amino acid positions reinforces this idea. Reconstruction of habitat depths of ancestral nodes in the phylogeny suggests that shallow habitats have been colonized independently in different lineages. The evolution of rhodopsin appears to be associated with changes in depth, with accelerated evolution in lineages that have had large changes in depth.

A method for detecting population genetic structure in diverse, high gene-flow species.

Detecting small amounts of genetic subdivision across geographic space remains a persistent challenge. Often a failure to detect genetic structure is mistaken for evidence of panmixia, when more powerful statistical tests may uncover evidence for subtle geographic differentiation. Such slight subdivision can be demographically and evolutionarily important as well as being critical for management decisions. We introduce here a method, called spatial analysis of shared alleles (SAShA), that detects geographically restricted alleles by comparing the spatial arrangement of allelic co-occurrences with the expectation under panmixia. The approach is allele-based and spatially explicit, eliminating the loss of statistical power that can occur with user-defined populations and statistical averaging within populations. Using simulated data sets generated under a stepping-stone model of gene flow, we show that this method outperforms spatial autocorrelation (SA) and Phi(ST) under common real-world conditions: at relatively high migration rates when diversity is moderate or high, especially when sampling is poor. We then use this method to show clear differences in the genetic patterns of 2 nearshore Pacific mollusks, Tegula funebralis (= Chlorostoma funebralis) and Katharina tunicata, whose overall patterns of within-species differentiation are similar according to traditional population genetics analyses. SAShA meaningfully complements Phi(ST)/F(ST), SA, and other existing geographic genetic analyses and is especially appropriate for evaluating species with high gene flow and subtle genetic differentiation.

Sampling from natural populations with RNAI reveals high outcrossing and population structure in Caenorhabditis elegans.

Despite a nearly worldwide distribution in nature, Caenorhabditis elegans exhibits low levels of genetic polymorphism, possibly as an indirect consequence of low levels of outcrossing. In the laboratory, Caenorhabditis elegans males are produced at low rates and are steadily eliminated from cultures, so that reproduction happens largely through self-fertilization in hermaphrodites. C. elegans is increasingly the focus of evolutionary research; however, natural outcrossing rates are difficult to measure because mating tests with laboratory strains are usually required to identify C. elegans. We sampled natural populations of C. elegans with an RNA interference (RNAi) assay. Heterozygosities and polymorphism patterns revealed surprisingly high levels of population structure and outcrossing (approximately 22% of individuals are estimated to be the result of outcrossing and not self-fertilization). The finding of strong local population structure, together with low levels of diversity on local and global scales, suggests a metapopulation model of frequent extinction and recolonization of local populations. The occurrence of substantial outcrossing suggests that the extinction of local populations is probably not driven by the accumulation of harmful mutations.

Population genetics of Caenorhabditis elegans: the paradox of low polymorphism in a widespread species.

Caenorhabditis elegans has become one of the most widely used model research organisms, yet we have little information on evolutionary processes and recent evolutionary history of this widespread species. We examined patterns of variation at 20 microsatellite loci in a sample of 23 natural isolates of C. elegans from various parts of the world. One-half of the loci were monomorphic among all strains, and overall genetic variation at microsatellite loci was low, relative to most other species. Some population structure was detected, but there was no association between the genetic and geographic distances among different natural isolates. Thus, despite the nearly worldwide occurrence of C. elegans, little evidence was found for local adaptation in strains derived from different parts of the world. The low levels of genetic variation within and among populations suggest that recent colonization and population expansion might have occurred. However, the patterns of variation are not consistent with population expansion. A possible explanation for the observed patterns is the action of background selection to reduce polymorphism, coupled with ongoing gene flow among populations worldwide.

Disentangling the role of phenotypic plasticity and genetic divergence in contemporary ecotype formation during a biological invasion.

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full-sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.

Evidence of adaptive evolutionary divergence during biological invasion.

Rapid phenotypic diversification during biological invasions can either arise by adaptation to alternative environments or by adaptive phenotypic plasticity. Where experimental evidence for adaptive plasticity is common, support for evolutionary diversification is rare. Here, we performed a controlled laboratory experiment using full-sib crosses between ecologically divergent threespine stickleback populations to test for a genetic basis of adaptation. Our populations are from two very different habitats, lake and stream, of a recently invaded range in Switzerland and differ in ecologically relevant morphological traits. We found that in a lake-like food treatment lake fish grow faster than stream fish, resembling the difference among wild type individuals. In contrast, in a stream-like food treatment individuals from both populations grow similarly. Our experimental data suggest that genetically determined diversification has occurred within less than 140 years after the arrival of stickleback in our studied region.

Nuclear gene variation and molecular dating of the cichlid species flock of Lake Malawi.

The cichlid fishes of Lake Malawi are famously diverse. However, evolutionary studies have been difficult because of their recent and uncertain phylogenetic history. Portions of 12 nuclear loci were sequenced in nine rock-dwelling species (mbuna) and three representatives of pelagic nonmbuna species. In contrast to the pattern of variation at mitochondrial genes, which do provide phylogenetic resolution at the level of mbuna versus nonmbuna, and among some genera, the nuclear loci were virtually devoid of phylogenetic signal. Only a small minority of variable positions were phylogenetically informative, and no phylogenetic branches are supported by more than one site. From the nuclear gene perspective the Malawian radiation appears to be a star phylogeny, as if the founding of the lake was accompanied by a partial bottleneck. The pattern is different from that found in Lake Victoria, in which nuclear loci share large amounts of ancestral variation. In the case of nuclear genes of Lake Malawi, the absence of phylogenetically informative variation suggests a relative absence of ancestral variation. Nuclear genes also differed from the mitochondria in having nearly twice the amount of divergence from Oreochromis (tilapia). An approximate splitting time between mbuna and nonmbuna lineages was estimated as 0.7 Myr. Oreochromis is estimated to have diverged from the cichlids in Lake Malawi and Lake Tanganyika about 18 MYA.

On the origin of Lake Malawi cichlid species: a population genetic analysis of divergence.

The cichlid fishes of Lake Malawi are famously diverse. However, phylogenetic and population genetic studies of their history have been difficult because of the great amount of genetic variation that is shared between species. We apply a recently developed method for fitting the "isolation with migration" divergence model to a data set of specially designed compound loci to develop portraits of cichlid species divergence. Outgroup sequences from a cichlid from Lake Tanganyika permit model parameter estimates in units of years and effective population sizes. Estimated speciation times range from 1,000 to 17,000 years for species in the genus Tropheops. These exceptionally recent dates suggest that Malawi cichlids as a group experience a very active and dynamic diversification process. Current effective population size estimates range form 2,000 to near 40,000, and to >120,000 for estimates of ancestral population sizes. It appears that very recent speciation and gene flow are among the reasons why it has been difficult to discern the phylogenetic history of Malawi cichlids.

Using nuclear haplotypes with microsatellites to study gene flow between recently separated Cichlid species.

When populations or species have recently separated they often share genetic variation. However, it can be difficult to determine whether shared polymorphisms are the result of gene flow, the result of the persistence of variation in both populations since the time of common ancestry, or both of these factors. We have developed an empirical protocol for using loci that include unique nuclear DNA sequence haplotypes together with linked microsatellites or short tandem repeats (STRs). These 'HapSTRs' offer the potentially high resolution associated with the high mutation rate of STRs, together with the advantages of low homoplasy of unique sequence DNA. We also describe a new procedure for estimating the likelihood of HapSTR data under an Isolation with Migration model. An example using Cichlid fishes from Lake Malawi is described. The analysis suggests that the species have been exchanging genes since the time they began to diverge.