A chromosomal inversion may facilitate adaptation despite periodic gene flow in a freshwater fish

Differences in genomic architecture between populations, such as chromosomal inversions, may play an important role in facilitating adaptation despite opportunities for gene flow. One system where chromosomal inversions may be important for eco‐evolutionary dynamics is in freshwater fishes, which often live in heterogenous environments characterized by varying levels of connectivity and varying opportunities for gene flow. In the present study, reduced representation sequencing was used to study possible adaptation in n = 345 walleye (Sander vitreus) from three North American waterbodies: Cedar Bluff Reservoir (Kansas, USA), Lake Manitoba (Manitoba, Canada), and Lake Winnipeg (Manitoba, Canada). Haplotype and outlier‐based tests revealed a putative chromosomal inversion that contained three expressed genes and was nearly fixed in walleye assigned to Lake Winnipeg. These patterns exist despite the potential for high gene flow between these proximate Canadian lakes, suggesting that the inversion may be important for facilitating adaptive divergence between the two lakes despite gene flow. However, a specific adaptive role for the putative inversion could not be tested with the present data. Our study illuminates the importance of genomic architecture consistent with local adaptation in freshwater fishes. Furthermore, our results provide additional evidence that inversions may facilitate local adaptation in many organisms that inhabit connected but heterogenous environments.     

Phylogeography and adaptation genetics of stickleback from the Haida Gwaii archipelago revealed using genome‐wide single nucleotide polymorphism genotyping

Threespine stickleback populations are model systems for studying adaptive evolution and the underlying genetics. In lakes on the Haida Gwaii archipelago (off western Canada), stickleback have undergone a remarkable local radiation and show phenotypic diversity matching that seen throughout the species distribution. To provide a historical context for this radiation, we surveyed genetic variation at >1000 single nucleotide polymorphism (SNP) loci in stickleback from over 100 populations. SNPs included markers evenly distributed throughout genome and candidate SNPs tagging adaptive genomic regions. Based on evenly distributed SNPs, the phylogeographic pattern differs substantially from the disjunct pattern previously observed between two highly divergent mtDNA lineages. The SNP tree instead shows extensive within watershed population clustering and different watersheds separated by short branches deep in the tree. These data are consistent with separate colonizations of most watersheds, despite underlying genetic connections between some independent drainages. This supports previous suppositions that morphological diversity observed between watersheds has been shaped independently, with populations exhibiting complete loss of lateral plates and giant size each occurring in several distinct clades. Throughout the archipelago, we see repeated selection of SNPs tagging candidate freshwater adaptive variants at several genomic regions differentiated between marine–freshwater populations on a global scale (e.g. EDA, Na/K ATPase). In estuarine sites, both marine and freshwater allelic variants were commonly detected. We also found typically marine alleles present in a few freshwater lakes, especially those with completely plated morphology. These results provide a general model for postglacial colonization of freshwater habitat by sticklebacks and illustrate the tremendous potential of genome‐wide SNP data sets hold for resolving patterns and processes underlying recent adaptive divergences.     

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long‐term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene–temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene–environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.     

Less is more: extreme genome complexity reduction with ddRAD using Ion Torrent semiconductor technology

Massively parallel sequencing a small proportion of the whole genome at high coverage enables answering a wide range of questions from molecular evolution and evolutionary biology to animal and plant breeding and forensics. In this study, we describe the development of restriction‐site associated DNA (RAD) sequencing approach for Ion Torrent PGM platform. Our protocol results in extreme genome complexity reduction using two rare‐cutting restriction enzymes and strict size selection of the library allowing sequencing of a relatively small number of genomic fragments with high sequencing depth. We applied this approach to a common freshwater fish species, the Eurasian perch (Perca fluviatilis L.), and generated over 2.2 MB of novel sequence data consisting of ~17 000 contigs, identified 1259 single nucleotide polymorphisms (SNPs). We also estimated genetic differentiation between the DNA pools from freshwater (Lake Peipus) and brackish water (the Baltic Sea) populations and identified SNPs with the strongest signal of differentiation that could be used for robust individual assignment in the future. This work represents an important step towards developing genomic resources and genetic tools for the Eurasian perch. We expect that our ddRAD sequencing protocol for semiconductor sequencing technology will be useful alternative for currently available RAD protocols.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

The availability of diverse ecological niches can promote adaptation of trophic specializations and related traits, as has been repeatedly observed in evolutionary radiations of freshwater fish. The role of genetics, environment, and history in ecologically driven divergence and adaptation, can be studied on adaptive radiations or populations showing ecological polymorphism. Salmonids, especially the Salvelinus genus, are renowned for both phenotypic diversity and polymorphism. Arctic charr (Salvelinus alpinus) invaded Icelandic streams during the glacial retreat (about 10,000 years ago) and exhibits many instances of sympatric polymorphism. Particularly, well studied are the four morphs in Lake Þingvallavatn in Iceland. The small benthic (SB), large benthic (LB), planktivorous (PL), and piscivorous (PI) charr differ in many regards, including size, form, and life history traits. To investigate relatedness and genomic differentiation between morphs, we identified variable sites from RNA‐sequencing data from three of those morphs and verified 22 variants in population samples. The data reveal genetic differences between the morphs, with the two benthic morphs being more similar and the PL‐charr more genetically different. The markers with high differentiation map to all linkage groups, suggesting ancient and pervasive genetic separation of these three morphs. Furthermore, GO analyses suggest differences in collagen metabolism, odontogenesis, and sensory systems between PL‐charr and the benthic morphs. Genotyping in population samples from all four morphs confirms the genetic separation and indicates that the PI‐charr are less genetically distinct than the other three morphs. The genetic separation of the other three morphs indicates certain degree of reproductive isolation. The extent of gene flow between the morphs and the nature of reproductive barriers between them remain to be elucidated.     

Hierarchical population structure and genetic diversity of lake trout (Salvelinus namaycush) in a dendritic system in Northern Labrador

相似文献   

Double‐digest RAD sequencing using Ion Proton semiconductor platform (ddRADseq‐ion) with nonmodel organisms

Research in evolutionary biology involving nonmodel organisms is rapidly shifting from using traditional molecular markers such as mtDNA and microsatellites to higher throughput SNP genotyping methodologies to address questions in population genetics, phylogenetics and genetic mapping. Restriction site associated DNA sequencing (RAD sequencing or RADseq) has become an established method for SNP genotyping on Illumina sequencing platforms. Here, we developed a protocol and adapters for double‐digest RAD sequencing for Ion Torrent (Life Technologies; Ion Proton, Ion PGM) semiconductor sequencing. We sequenced thirteen genomic libraries of three different nonmodel vertebrate species on Ion Proton with PI chips: Arctic charr Salvelinus alpinus, European whitefish Coregonus lavaretus and common lizard Zootoca vivipara. This resulted in ~962 million single‐end reads overall and a mean of ~74 million reads per library. We filtered the genomic data using Stacks, a bioinformatic tool to process RAD sequencing data. On average, we obtained ~11 000 polymorphic loci per library of 6–30 individuals. We validate our new method by technical and biological replication, by reconstructing phylogenetic relationships, and using a hybrid genetic cross to track genomic variants. Finally, we discuss the differences between using the different sequencing platforms in the context of RAD sequencing, assessing possible advantages and disadvantages. We show that our protocol can be used for Ion semiconductor sequencing platforms for the rapid and cost‐effective generation of variable and reproducible genetic markers.     

Genetic variation within and between populations of hermaphroditic <Emphasis Type="Italic">Bulinus truncatus</Emphasis> tetraploid freshwater snails of the Albertine Rift,East Africa

Genetic variability within and among Bulinus truncatus of the Albertine Rift freshwater bodies were assessed to investigate the degree of inbreeding and gene flow in the snail populations. The effect of ploidy on the genetic structuring of B. truncatus is also described. We characterized the genetic structure of seven B. truncatus populations from Lake Albert, Lake Kivu, and Katosho swamp in Tanzania using five polymorphic microsatellite loci. Genetic differentiation was quantified using pairwise FST values and Nei’s standard genetic distances. Different alleles were observed across all loci and genetic diversity was low although it varied greatly across populations; observed heterozygosity was, however, higher than the expected heterozygosity in three of the populations studied. Significant heterozygote deficiencies were observed coupled with significant linkage disequilibria in five populations for all the five loci examined in this study. We found significant genetic differentiation among the seven freshwater bodies; private alleles were observed across all loci indicating restricted or absence of gene flow between populations. Limited snail dispersal and the reproductive biology of B. truncatus are the major forces shaping the genetic variation observed. Low genetic variation within B. truncatus populations exposes them to a high parasite infection risk as predicted in the Red Queen hypothesis.     

Migratory behaviour shapes spatial genetic structure of cyprinid fishes within the Lake Malawi catchment

• Genetic differences among freshwater fish populations are dependent on life‐history characteristics of the species, including the range of adult dispersal and the extent of homing to natal breeding grounds. However, the effects of variation in such characteristics on population genetic connectivity are rarely studied comparatively among closely related species.
• We studied population genetic structure within three congeneric cyprinid species from the Lake Malawi catchment that differ substantially in life‐history traits and conservation status, using a combination of microsatellite and mitochondrial DNA markers. Mpasa (Opsaridium microlepis) is a large (70 cm total length) migratory species that spawns in rivers, but as an adult is exclusively known from the main lake body. Sanjika (Opsaridium microcephalum), is a medium size (30 cm total length) species that exists in lake breeding, river‐lake migratory and apparently landlocked populations. Dwarf sanjika (Opsaridium tweddleorum) is a small non‐migratory species (15 cm total length) that persists in small tributaries surrounding the main lake and adjoining rivers.
• The results revealed striking differences among the three species in spatial genetic structuring. The river‐lake migratory mpasa showed only weak yet significant population genetic structure within the main Lake Malawi catchment, suggesting that there is no strong natal homing. The habitat‐generalist sanjika showed only weak spatial genetic differentiation at microsatellite loci within the Lake Malawi catchment, but moderate structure in mitochondrial DNA, potentially reflecting male‐biased dispersal. The river‐restricted dwarf sanjika showed strong genetic structure in both microsatellite and mitochondrial DNA, suggesting strictly limited dispersal at both adult and juvenile stages.
• We conclude that contrasting migration life histories have resulted in dramatically different patterns of population genetic structure among these congeneric species. The observed patterns demonstrate how divergent life‐history evolution may strongly influence broader patterns of population genetic connectivity in freshwater fish, with consequences for management and conservation. Specifically the results suggesting gene flow among Lake Malawi populations of mpasa, an IUCN red‐listed ‘Endangered’ species endemic to the lake catchment, imply that conservation initiatives operating at both local and catchment scales are needed to reverse local population decline.

     

Genomewide scan for adaptive differentiation along altitudinal gradient in the Andrew's toad Bufo andrewsi

Recent studies of humans, dogs and rodents have started to discover the genetic underpinnings of high altitude adaptations, yet amphibians have received little attention in this respect. To identify possible signatures of adaptation to altitude, we performed a genome scan of 15 557 single nucleotide polymorphisms (SNPs) obtained with restriction site‐associated DNA sequencing of pooled samples from 11 populations of Andrew's toad (Bufo andrewsi) from the edge of the Tibetan Plateau, spanning an altitudinal gradient from 1690 to 2768 m.a.s.l. We discovered significant geographic differentiation among all sites, with an average F_ST = 0.023 across all SNPs. Apart from clear patterns of isolation by distance, we discovered numerous outlier SNPs showing strong associations with variation in altitude (1394 SNPs), average annual temperature (1859 SNPs) or both (1051 SNPs). Levels and patterns of genetic differentiation in these SNPs were consistent with the hypothesis that they have been subject to directional selection and reflect adaptation to altitudinal variation among the study sites. Genes with footprints of selection were significantly enriched in binding and metabolic processes. Several genes potentially related to high altitude adaptation were identified, although the identity and functional significance of most genomic targets of selection remain unknown. In general, the results provide genomic support for results of earlier common garden and low coverage genetic studies that have uncovered substantial adaptive differentiation along altitudinal and latitudinal gradients in amphibians.     

Comparative analyses of actinobacterial genomic fragments from Lake Kinneret

The high genomic G+C group of Actinobacteria possesses a variety of physiological and metabolic properties, and exhibits diverse lifestyles and ecological distribution. In recent years, Actinobacteria have been found to frequently dominate samples obtained from freshwater samples. Furthermore, phylogenetic analyses have shown that 16S rRNA genes from uncultured actinobacterial freshwater samples cluster in four distinct lineages. While these lineages are abundant, little is known about them and currently no pure‐culture representatives or genomic fragments of them are available. In a screen of a genomic library from the moderately eutrophic freshwater Lake Kinneret, five fosmid clones containing actinobacterial genomic fragments were found. Three ～40 kb genomic fragments were chosen for sequencing. Fosmids K003 and K005 showed high similarity and were affiliated with the acIV actinobacterial freshwater lineage. Fosmid K004 was affiliated with the highly abundant acI lineage. A comparative genomic analysis revealed high synteny between the two freshwater clones K003 and K005 but a lower synteny between these two and the K004 fosmid. Fosmids K003 and K005 share an identical arrangement of arginine biosynthesis gene while K004 showed a slightly different arrangement by lacking the argF gene. Fosmid Ant4E12, an Antarctic actinobacterial clone, showed a higher synteny with K003/5 than K004 and a similar arginine operon, but in a different genomic context. The Clusters of Orthologous Groups categories assignment of the three fosmids yielded genes that were mostly involved in amino acid and nucleotide metabolism, as well as transport and ribosomal RNA translation, structure and biogenesis. These genomic fragments represent the first sequences to be published from these lineages, providing a cornerstone for future work on this environmentally dominant group.     

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.     

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the “eucalypts”) are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole‐genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome‐wide, interspecific differentiation (F_ST = 0.15) and intraspecific differentiation between localities (F_ST ≈ 0.01–0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.     

High levels of interspecific gene flow in an endemic cichlid fish adaptive radiation from an extreme lake environment

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (F_ST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.     

Ecological genomics in full colour

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning et al. ( 2014 ) detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson et al. ( 2014 ) identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson et al. ( 2014 ) go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics – efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.     

Structural genomic variation leads to genetic differentiation in Lake Tanganyika's sardines

Identifying patterns in genetic structure and the genetic basis of ecological adaptation is a core goal of evolutionary biology and can inform the management and conservation of species that are vulnerable to population declines exacerbated by climate change. We used reduced‐representation genomic sequencing methods to gain a better understanding of genetic structure among and within populations of Lake Tanganyika's two sardine species, Limnothrissa miodon and Stolothrissa tanganicae. Samples of these ecologically and economically important species were collected across the length of Lake Tanganyika, as well as from nearby Lake Kivu, where L. miodon was introduced in 1959. Our results reveal differentiation within both S. tanganicae and L. miodon that is not explained by geography. Instead, this genetic differentiation is due to the presence of large sex‐specific regions in the genomes of both species, but involving different polymorphic sites in each species. Our results therefore indicate rapidly evolving XY sex determination in the two species. Additionally, we found evidence of a large chromosomal rearrangement in L. miodon, creating two homokaryotypes and one heterokaryotype. We found all karyotypes throughout Lake Tanganyika, but the frequencies vary along a north–south gradient and differ substantially in the introduced Lake Kivu population. We do not find evidence for significant isolation by distance, even over the hundreds of kilometres covered by our sampling, but we do find shallow population structure.     

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.