Linkage maps of the dwarf and Normal lake whitefish (Coregonus clupeaformis) species complex and their hybrids reveal the genetic architecture of population divergence

Genome assembly, structural variants, and genetic differentiation between lake whitefish young species pairs (Coregonus sp.) with long and short reads

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics.

Current genetic structure of European vendace Coregonus albula (L.) populations in Latvia after multiple past translocations

The European vendace Coregonus albula (L.), also known as the European cisco, is a widespread fish species in northern Europe, often regarded as an example of a glacial relict. It is an economically valuable fish and has been artificially propagated in Latvia since 1900. Despite past translocations of larvae and fry and its current protection status, it can be found in only 15 Latvian lakes. We used nine microsatellite markers to study vendace populations from nine Latvian lakes. A higher mean allelic richness and private allelic richness in Lake Riču suggest that this population may be indigenous. Three complementary clustering methods revealed similar grouping into three distinct genetic groups. According to the results, European vendace populations in the Latvian lakes studied may currently be a mixture of several other populations after multiple translocations.

Rapid buildup of sympatric species diversity in Alpine whitefish

Adaptive radiations in postglacial fish offer excellent settings to study the evolutionary mechanisms involved in the rapid buildup of sympatric species diversity from a single lineage. Here, we address this by exploring the genetic and ecological structure of the largest Alpine whitefish radiation known, that of Lakes Brienz and Thun, using microsatellite data of more than 2000 whitefish caught during extensive species-targeted and habitat-randomized fishing campaigns. We find six strongly genetically and ecologically differentiated species, four of which occur in both lakes, and one of which was previously unknown. These four exhibit clines of genetic differentiation that are paralleled in clines of eco-morphological and reproductive niche differentiation, consistent with models of sympatric ecological speciation along environmental gradients. In Lake Thun, we find two additional species, a profundal specialist and a species introduced in the 1930s from another Alpine whitefish radiation. Strong genetic differentiation between this introduced species and all native species of Lake Thun suggests that reproductive isolation can evolve among allopatric whitefish species within 15,000 years and persist in secondary sympatry. Consistent with speciation theory, we find stronger correlations between genetic and ecological differentiation for sympatrically than for allopatrically evolved species.

Thermal experience during embryogenesis contributes to the induction of dwarfism in whitefish Coregonus lavaretus

Ecotype pairs provide well-suited model systems for study of intraspecific phenotypical diversification of animals. However, little is still known about the processes that account for the development of different forms and sizes within a species, particularly in teleosts. Here, embryos of a normal-growing ‘large’ form and a dwarf form of whitefish Coregonus lavaretus were incubated at two temperatures that are usually experienced at their own spawning sites (2°C for the normal and 6°C for the dwarf form). All fish were subjected to similar thermal treatment after hatching. The present data demonstrate for the first time that different thermal experience in embryonic life has lasting effects on body and muscle growth of this ecotype pair and contributes to the development of the dwarf form. Thus, juvenile fish of the regular form are much smaller and have less muscle mass when pre-hatching thermal conditions were similar to those typical for the spawning sites of the dwarf form (6°C) than when subjected to conditions of their own spawning sites (2°C). Surprisingly, fish of the dwarf form exhibit a similar pattern of response to thermal history (2°-fish much larger than 6°-fish), indicating that in their case, normal spawning site temperature (6°C) is indeed likely to act as a growth limiting factor. Results also demonstrate that the hypertrophic and hyperplastic muscle growth modes are similarly affected by thermal history. Immunolabelling experiments for Pax7, H3P and Mef2 provide evidence that the cellular mechanisms behind the increased growth rates after cold incubation in both ecotypes are increased proliferation and reduced differentiation rates of muscle precursor cells. This is of major significance to aspects of ecological and developmental biology and from the evolutionary perspective.

Comparative High-Density Linkage Mapping Reveals Conserved Genome Structure but Variation in Levels of Heterochiasmy and Location of Recombination Cold Spots in the Common Frog

By combining 7077 SNPs and 61 microsatellites, we present the first linkage map for some of the early diverged lineages of the common frog, Rana temporaria, and the densest linkage map to date for this species. We found high homology with the published linkage maps of the Eastern and Western lineages but with differences in the order of some markers. Homology was also strong with the genome of the Tibetan frog Nanorana parkeri and we found high synteny with the clawed frog Xenopus tropicalis. We confirmed marked heterochiasmy between sexes and detected nonrecombining regions in several groups of the male linkage map. Contrary to the expectations set by the male heterogamety of the common frog, we did not find male heterozygosity excess in the chromosome previously shown to be linked to sex determination. Finally, we found blocks of loci showing strong transmission ratio distortion. These distorted genomic regions might be related to genetic incompatibilities between the parental populations, and are promising candidates for further investigation into the genetic basis of speciation and adaptation in the common frog.

High-density interspecific genetic linkage mapping provides insights into genomic incompatibility between channel catfish and blue catfish

Catfish is the leading aquaculture species in the United States. The interspecific hybrid catfish produced by mating female channel catfish with male blue catfish outperform both of their parent species in a number of traits. However, mass production of the hybrids has been difficult because of reproductive isolation. Investigations of genome structure and organization of the hybrids provide insights into the genetic basis for maintenance of species divergence in the face of gene flow, thereby helping develop strategies for introgression and efficient production of the hybrids for aquaculture. In this study, we constructed a high-density genetic linkage map using the hybrid catfish system with the catfish 250K SNP array. A total of 26,238 SNPs were mapped to 29 linkage groups, with 12,776 unique marker positions. The linkage map spans approximately 3240 cM with an average intermarker distance of 0.25 cM. A fraction of markers (986 of 12,776) exhibited significant deviation from the expected Mendelian ratio of segregation, and they were clustered in major genomic blocks across 15 LGs, most notably LG9 and LG15. The distorted markers exhibited significant bias for maternal alleles among the backcross progenies, suggesting strong selection against the blue catfish alleles. The clustering of distorted markers within genomic blocks should lend insights into speciation as marked by incompatibilities between the two species. Such findings should also have profound implications for understanding the genomic evolution of closely related species as well as the introgression of hybrid production programs in aquaculture.

Anthropogenic hybridization between endangered migratory and commercially harvested stationary whitefish taxa (Coregonus spp.)

Natural hybridization plays a key role in the process of speciation. However, anthropogenic (human induced) hybridization of historically isolated taxa raises conservation issues. Due to weak barriers to gene flow and the presence of endangered taxa, the whitefish species complex is an excellent study system to investigate the consequences of hybridization in conservation. We focused on three naturally reproductively isolated whitefish taxa in Germany: the endangered, anadromous North Sea houting (NSH) and Baltic houting (BH), which were reintroduced after local extinction, and the commercially stocked European whitefish (EW). To evaluate the genetic integrity of each taxon, source and reintroduced populations of NSH and BH, and EW populations were characterized based on two mitochondrial and 17 microsatellite loci. Additionally, we investigated gill raker counts as an adaptive phenotypic trait. Even though clear genetic and phenotypic differentiation confirmed the houtings as separate evolutionarily significant units, admixture analyses revealed an extensive hybrid zone. Hybridizations were introgressive, positively correlated with genetic diversity, and were reflected in the gill raker counts. The BH distribution range showed higher heterogeneity and stronger admixture than the NSH range. Erroneous stocking with non-native genotypes best explained these patterns, which pose challenges for the conservation of the endangered NSH and BH.

Hidden genetic variation evolves with ecological specialization: the genetic basis of phenotypic plasticity in Arctic charr ecomorphs

The genetic variance that determines phenotypic variation can change across environments through developmental plasticity and in turn play a strong role in evolution. Induced changes in genotype-phenotype relationships should strongly influence adaptation by exposing different sets of heritable variation to selection under some conditions, while also hiding variation. Therefore, the heritable variation exposed or hidden from selection is likely to differ among habitats. We used ecomorphs from two divergent populations of Arctic charr (Salvelinus alpinus) to test the prediction that genotype-phenotype relationships would change in relation to environment. If present over several generations this should lead to divergence in genotype-phenotype relationships under common conditions, and to changes in the amount and type of hidden genetic variance that can evolve. We performed a common garden experiment whereby two ecomorphs from each of two Icelandic lakes were reared under conditions that mimicked benthic and limnetic prey to induce responses in craniofacial traits. Using microsatellite based genetic maps, we subsequently detected QTL related to these craniofacial traits. We found substantial changes in the number and type of QTL between diet treatments and evidence that novel diet treatments can in some cases provide a higher number of QTL. These findings suggest that selection on phenotypic variation, which is both genetically and environmentally determined, has shaped the genetic architecture of adaptive divergence in Arctic charr. However, while adaptive changes are occurring in the genome there also appears to be an accumulation of hidden genetic variation for loci not expressed in the contemporary environment.

Insight into genomic changes accompanying divergence: genetic linkage maps and synteny of Lucania goodei and L. parva reveal a Robertsonian fusion

Linkage maps are important tools in evolutionary genetics and in studies of speciation. We performed a karyotyping study and constructed high-density linkage maps for two closely related killifish species, Lucania parva and L. goodei, that differ in salinity tolerance and still hybridize in their contact zone in Florida. Using SNPs from orthologous EST contigs, we compared synteny between the two species to determine how genomic architecture has shifted with divergence. Karyotyping revealed that L. goodei possesses 24 acrocentric chromosomes (1N) whereas L. parva possesses 23 chromosomes (1N), one of which is a large metacentric chromosome. Likewise, high-density single-nucleotide polymorphism-based linkage maps indicated 24 linkage groups for L. goodei and 23 linkage groups for L. parva. Synteny mapping revealed two linkage groups in L. goodei that were highly syntenic with the largest linkage group in L. parva. Together, this evidence points to the largest linkage group in L. parva being the result of a chromosomal fusion. We further compared synteny between Lucania with the genome of a more distant teleost relative medaka (Oryzias latipes) and found good conservation of synteny at the chromosomal level. Each Lucania LG had a single best match with each medaka chromosome. These results provide the groundwork for future studies on the genetic architecture of reproductive isolation and salinity tolerance in Lucania and other Fundulidae.

Evidence for genetic incompatibilities associated with post-zygotic reproductive isolation in the human fungal pathogen Cryptococcus neoformans

Hybridization is a potent mechanism for generating unique strains with broad host ranges and increased virulence in fungal pathogens. In the opportunistic basidiomycete pathogen Cryptococcus neoformans, intervarietal hybrids are commonly found infecting patients. The two parental varieties C. neoformans var. grubii and C. neoformans var. neoformans mate readily under laboratory conditions, but the hybrid basidiospores have germination rates about four times lower than those from intravarietal crosses. Here, we used microdissection to collect basidiospores from a hybrid cross and analysed the genotypes of germinated basidiospores to identify potentially antagonistic allelic combinations between loci that impact basidiospore germination. Our analyses showed clear evidence for Bateson–Dobzhansky–Muller (BDM) incompatibility affecting basidiospore viability. Antagonistic combinations of alleles from both two loci and three loci were found. Interestingly, most of the hybrid progeny showed segregation distortion in favour of the alleles from var. neoformans, consistent with large-scale epistatic interactions among loci affecting basidiospore viability. Our study presents the first evidence of BDM incompatibility between nuclear genes affecting post-zygotic reproductive isolation in this model basidiomycete yeast.

Chromosome rearrangements, recombination suppression, and limited segregation distortion in hybrids between Yellowstone cutthroat trout (Oncorhynchus clarkii bouvieri) and rainbow trout (O. mykiss)

Background

Introgressive hybridization is an important evolutionary process that can lead to the creation of novel genome structures and thus potentially new genetic variation for selection to act upon. On the other hand, hybridization with introduced species can threaten native species, such as cutthroat trout (Oncorhynchus clarkii) following the introduction of rainbow trout (O. mykiss). Neither the evolutionary consequences nor conservation implications of rainbow trout introgression in cutthroat trout is well understood. Therefore, we generated a genetic linkage map for rainbow-Yellowstone cutthroat trout (O. clarkii bouvieri) hybrids to evaluate genome processes that may help explain how introgression affects hybrid genome evolution.

Results

The hybrid map closely aligned with the rainbow trout map (a cutthroat trout map does not exist), sharing all but one linkage group. This linkage group (RYHyb20) represented a fusion between an acrocentric (Omy28) and a metacentric chromosome (Omy20) in rainbow trout. Additional mapping in Yellowstone cutthroat trout indicated the two rainbow trout homologues were fused in the Yellowstone genome. Variation in the number of hybrid linkage groups (28 or 29) likely depended on a Robertsonian rearrangement polymorphism within the rainbow trout stock. Comparison between the female-merged F₁ map and a female consensus rainbow trout map revealed that introgression suppressed recombination across large genomic regions in 5 hybrid linkage groups. Two of these linkage groups (RYHyb20 and RYHyb25_29) contained confirmed chromosome rearrangements between rainbow and Yellowstone cutthroat trout indicating that rearrangements may suppress recombination. The frequency of allelic and genotypic segregation distortion varied among parents and families, suggesting few incompatibilities exist between rainbow and Yellowstone cutthroat trout genomes.

Conclusions

Chromosome rearrangements suppressed recombination in the hybrids. This result supports several previous findings demonstrating that recombination suppression restricts gene flow between chromosomes that differ by arrangement. Conservation of synteny and map order between the hybrid and rainbow trout maps and minimal segregation distortion in the hybrids suggest rainbow and Yellowstone cutthroat trout genomes freely introgress across chromosomes with similar arrangement. Taken together, these results suggest that rearrangements impede introgression. Recombination suppression across rearrangements could enable large portions of non-recombined chromosomes to persist within admixed populations.

Natural selection and neutral evolution jointly drive population divergence between alpine and lowland ecotypes of the allopolyploid plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074-0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (FST  = 0.041-0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.

A hybrid genetic linkage map of two ecologically and morphologically divergent Midas cichlid fishes (Amphilophus spp.) obtained by massively parallel DNA sequencing (ddRADSeq)

Cichlid fishes are an excellent model system for studying speciation and the formation of adaptive radiations because of their tremendous species richness and astonishing phenotypic diversity. Most research has focused on African rift lake fishes, although Neotropical cichlid species display much variability as well. Almost one dozen species of the Midas cichlid species complex (Amphilophus spp.) have been described so far and have formed repeated adaptive radiations in several Nicaraguan crater lakes. Here we apply double-digest restriction-site associated DNA sequencing to obtain a high-density linkage map of an interspecific cross between the benthic Amphilophus astorquii and the limnetic Amphilophus zaliosus, which are sympatric species endemic to Crater Lake Apoyo, Nicaragua. A total of 755 RAD markers were genotyped in 343 F₂ hybrids. The map resolved 25 linkage groups and spans a total distance of 1427 cM with an average marker spacing distance of 1.95 cM, almost matching the total number of chromosomes (n = 24) in these species. Regions of segregation distortion were identified in five linkage groups. Based on the pedigree of parents to F₂ offspring, we calculated a genome-wide mutation rate of 6.6 × 10⁻⁸ mutations per nucleotide per generation. This genetic map will facilitate the mapping of ecomorphologically relevant adaptive traits in the repeated phenotypes that evolved within the Midas cichlid lineage and, as the first linkage map of a Neotropical cichlid, facilitate comparative genomic analyses between African cichlids, Neotropical cichlids and other teleost fishes.

Analysis of the meiotic segregation in intergeneric hybrids of tilapias

Tilapia species exhibit a large ecological diversity and an important propensity to interspecific hybridisation. This has been shown in the wild and used in aquaculture. However, despite its important evolutionary implications, few studies have focused on the analysis of hybrid genomes and their meiotic segregation. Intergeneric hybrids between Oreochromis niloticus and Sarotherodon melanotheron, two species highly differentiated genetically, ecologically, and behaviourally, were produced experimentally. The meiotic segregation of these hybrids was analysed in reciprocal second generation hybrid (F2) and backcross families and compared to the meiosis of both parental species, using a panel of 30 microsatellite markers. Hybrid meioses showed segregation in accordance to Mendelian expectations, independent from sex and the direction of crosses. In addition, we observed a conservation of linkage associations between markers, which suggests a relatively similar genome structure between the two parental species and the apparent lack of postzygotic incompatibility, despite their important divergence. These results provide genomics insights into the relative ease of hybridisation within cichlid species when prezygotic barriers are disrupted. Overall our results support the hypothesis that hybridisation may have played an important role in the evolution and diversification of cichlids.

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.

Coding Gene Single Nucleotide Polymorphism Mapping and Quantitative Trait Loci Detection for Physiological Reproductive Traits in Brook Charr, Salvelinus fontinalis

A linkage map of 40 linkage groups (LGs) was developed for brook charr, Salvelinus fontinalis, using an F(2) interstrain hybrid progeny (n = 171) and 256 coding gene SNP developed specifically for brook charr and validated from a large (>1000) subset of putative SNP, as well as 81 microsatellite markers. To identify quantitative trait loci (QTL) related to reproduction functions, these fish were also phenotyped at six physiological traits, including spermatozoid head diameter, sperm concentration, plasma testosterone, plasma 11-keto-testosterone, egg diameter, and plasma 17β-estradiol. Five significant QTL were detected over four LGs for egg diameter and plasma 17β-estradiol concentration in females, and sperm concentration as well as spermatozoid head diameter in males. In females, two different QTLs located on LG 11 and LG 34 were associated with the egg number, whereas one QTL was associated with plasma 17β-estradiol concentration (LG 8). Their total percent variance explained (PVE) was 26.7% and 27.6%, respectively. In males, two QTL were also detected for the sperm concentration, and their PVE were estimated at 18.58% and 14.95%, respectively. The low QTL number, associated with the high PVE, suggests that the variance in these reproductive physiological traits was either under the control of one major gene or a small number of genes. The QTL associated with sperm concentration, plasma 17β-estradiol, and egg diameter appeared to be under a dominance effect, whereas the two others were under a negative additive effect. These results show that genes underlying the phenotypic variance of these traits are under different modes of action (additive vs. dominance) and may be used to predict an increase or a decrease in their phenotypic values in subsequent generations of selective breeding. Moreover, this newly developed panel of mapped SNP located in coding gene regions will be useful for screening wild populations, especially in the context of investigating the genetic impact of massive stocking of domestic brook charr to support the angling industry throughout eastern North America.

Genomic scan as a tool for assessing the genetic component of phenotypic variance in wild populations

Methods for estimating quantitative trait heritability in wild populations have been developed in recent years which take advantage of the increased availability of genetic markers to reconstruct pedigrees or estimate relatedness between individuals, but their application to real-world data is not exempt from difficulties. This chapter describes a recent marker-based technique which, by adopting a genomic scan approach and focusing on the relationship between phenotypes and genotypes at the individual level, avoids the problems inherent to marker-based estimators of relatedness. This method allows the quantification of the genetic component of phenotypic variance ("degree of genetic determination" or "heritability in the broad sense") in wild populations and is applicable whenever phenotypic trait values and multilocus data for a large number of genetic markers (e.g., amplified fragment length polymorphisms, AFLPs) are simultaneously available for a sample of individuals from the same population. The method proceeds by first identifying those markers whose variation across individuals is significantly correlated with individual phenotypic differences ("adaptive loci"). The proportion of phenotypic variance in the sample that is statistically accounted for by individual differences in adaptive loci is then estimated by fitting a linear model to the data, with trait value as the dependent variable and scores of adaptive loci as independent ones. The method can be easily extended to accommodate quantitative or qualitative information on biologically relevant features of the environment experienced by each sampled individual, in which case estimates of the environmental and genotype × environment components of phenotypic variance can also be obtained.

The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing

Gene expression divergence is one of the mechanisms thought to be involved in the emergence of incipient species. Next-generation sequencing has become an extremely valuable tool for the study of this process by allowing whole transcriptome sequencing, or RNA-Seq. We have conducted a 454 GS-FLX pyrosequencing experiment to refine our understanding of adaptive divergence between dwarf and normal lake whitefish species (Coregonus clupeaformis spp.). The objectives were to: (i) investigate transcriptomic divergence as measured by liver RNA-Seq; (ii) test the correlation between divergence in expression and sequence polymorphism; and (iii) investigate the extent of allelic imbalance. We also compared the results of RNA-seq with those of a previous microarray study performed on the same fish. Following de novo assembly, results showed that normal whitefish overexpressed more contigs associated with protein synthesis while dwarf fish overexpressed more contigs related to energy metabolism, immunity and DNA replication and repair. Moreover, 63 SNPs showed significant allelic imbalance, and this phenomenon prevailed in the recently diverged dwarf whitefish. Results also showed an absence of correlation between gene expression divergence as measured by RNA-Seq and either polymorphism rate or sequence divergence between normal and dwarf whitefish. This study reiterates an important role for gene expression divergence, and provides evidence for allele-specific expression divergence as well as evolutionary decoupling of regulatory and coding sequences in the adaptive divergence of normal and dwarf whitefish. It also demonstrates how next-generation sequencing can lead to a more comprehensive understanding of transcriptomic divergence in a young species pair.

On the origin of species: insights from the ecological genomics of lake whitefish

In contrast to the large amount of ecological information supporting the role of natural selection as a main cause of population divergence and speciation, an understanding of the genomic basis underlying those processes is in its infancy. In this paper, we review the main findings of a long-term research programme that we have been conducting on the ecological genomics of sympatric forms of whitefish (Coregonus spp.) engaged in the process of speciation. We present this system as an example of how applying a combination of approaches under the conceptual framework of the theory of adaptive radiation has yielded substantial insight into evolutionary processes in a non-model species. We also discuss how the joint use of recent biotechnological developments will provide a powerful means to address issues raised by observations made to date. Namely, we present data illustrating the potential offered by combining next generation sequencing technologies with other genomic approaches to reveal the genomic bases of adaptive divergence and reproductive isolation. Given increasing access to these new genomic tools, we argue that non-model species studied in their ecological context such as whitefish will play an increasingly important role in generalizing knowledge of speciation.

Speciation genetics: current status and evolving approaches

The view of species as entities subjected to natural selection and amenable to change put forth by Charles Darwin and Alfred Wallace laid the conceptual foundation for understanding speciation. Initially marred by a rudimental understanding of hereditary principles, evolutionists gained appreciation of the mechanistic underpinnings of speciation following the merger of Mendelian genetic principles with Darwinian evolution. Only recently have we entered an era where deciphering the molecular basis of speciation is within reach. Much focus has been devoted to the genetic basis of intrinsic postzygotic isolation in model organisms and several hybrid incompatibility genes have been successfully identified. However, concomitant with the recent technological advancements in genome analysis and a newfound interest in the role of ecology in the differentiation process, speciation genetic research is becoming increasingly open to non-model organisms. This development will expand speciation research beyond the traditional boundaries and unveil the genetic basis of speciation from manifold perspectives and at various stages of the splitting process. This review aims at providing an extensive overview of speciation genetics. Starting from key historical developments and core concepts of speciation genetics, we focus much of our attention on evolving approaches and introduce promising methodological approaches for future research venues.

Understanding the onset of hybrid speciation

Natural hybridization between closely related taxa is a common phenomenon in both plants and animals. Hybridization has often been viewed as a destructive force that could erode established gene pools, but it is increasingly being recognized as a potentially creative force in evolution because it can lead to a mixture of novel genotypes, some of which have the potential for rapid adaptation to new environmental conditions. However, the evolutionary dynamics leading to the emergence of newly adapted gene pools after hybridization are largely unexplored. Here, we argue that the identification and analysis of the dynamic processes that occur after the first contact deserve specific attention, because this is the phase where hybrid speciation is most different from other forms of speciation.

Genomic linkage map of the human blood fluke Schistosoma mansoni

BACKGROUND

Schistosoma mansoni is a blood fluke that infects approximately 90 million people. The complete life cycle of this parasite can be maintained in the laboratory, making this one of the few experimentally tractable human helminth infections, and a rich literature reveals heritable variation in important biomedical traits such as virulence, host-specificity, transmission and drug resistance. However, there is a current lack of tools needed to study S. mansoni's molecular, quantitative, and population genetics. Our goal was to construct a genetic linkage map for S. mansoni, and thus provide a new resource that will help stimulate research on this neglected pathogen.

RESULTS

We genotyped grandparents, parents and 88 progeny to construct a 5.6 cM linkage map containing 243 microsatellites positioned on 203 of the largest scaffolds in the genome sequence. The map allows 70% of the estimated 300 Mb genome to be ordered on chromosomes, and highlights where scaffolds have been incorrectly assembled. The markers fall into eight main linkage groups, consistent with seven pairs of autosomes and one pair of sex chromosomes, and we were able to anchor linkage groups to chromosomes using fluorescent in situ hybridization. The genome measures 1,228.6 cM. Marker segregation reveals higher female recombination, confirms ZW inheritance patterns, and identifies recombination hotspots and regions of segregation distortion.

CONCLUSIONS

The genetic linkage map presented here is the first for S. mansoni and the first for a species in the phylum Platyhelminthes. The map provides the critical tool necessary for quantitative genetic analysis, aids genome assembly, and furnishes a framework for comparative flatworm genomics and field-based molecular epidemiological studies.

Gene expression divergence and hybrid misexpression between lake whitefish species pairs (Coregonus spp. Salmonidae)

Genomewide analyses of the transcriptome have confirmed that gene misexpression may underlie reproductive isolation mechanisms in interspecific hybrids. Here, using a 16,006 features cDNA microarray, we compared and contrasted gene expression divergence at two ontogenetic stages in incipient species of normal and dwarf whitefish (Coregonus clupeaformis) with that of first generation (normal x dwarf) and second-generation hybrid crosses (backcross: [normal x dwarf] x normal]. Our goal was to identify the main mode of action responsible for gene transcription and to discover key genes misexpressed in hybrids. Very few transcripts (five of 4,950 expressed) differed in mean expression level between parentals and hybrids at the embryonic stage, in contrast to 16-week-old juvenile fish for which 617 out of 5,359 transcripts differed significantly. We also found evidence for more misexpression in backcross hybrids whereby nonadditivity explained a larger fraction of hybrid inheritance patterns in backcross (54%) compared with F1-hybrids (9%). Gene expression in hybrids was more variable than in pure crosses and transgressive patterns of expression were ubiquitous in hybrids. In backcross embryos in particular, the expression of three key developmental genes involved in protein folding and mRNA translation was severely disrupted. Accordingly, gene misexpression in hybrids adds to other factors previously identified as contributing to the reproductive isolation of incipient species of lake whitefish.

Divergence along a steep ecological gradient in lake whitefish (Coregonus sp.)

To understand mechanisms structuring diversity in young adaptive radiations, quantitative and unbiased information about genetic and phenotypic diversity is much needed. Here, we present the first in-depth investigation of whitefish diversity in a Swiss lake, with continuous spawning habitat sampling in both time and space. Our results show a clear cline like pattern in genetics and morphology of populations sampled along an ecological depth gradient in Lake Neuchâtel. Divergent natural selection appears to be involved in shaping this cline given that trait specific P(ST)-values are significantly higher than F(ST)-values when comparing populations caught at different depths. These differences also tend to increase with increasing differences in depth, indicating adaptive divergence along a depth gradient, which persists despite considerable gene flow between adjacent demes. It however remains unclear, whether the observed pattern is a result of currently stable selection-gene flow balance, incipient speciation, or reverse speciation due to anthropogenic habitat alteration causing two formerly divergent species to collapse into a single gene pool.

The phenomics and expression quantitative trait locus mapping of brain transcriptomes regulating adaptive divergence in lake whitefish species pairs (Coregonus sp.)

We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between the sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development on the basis of colocalization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. Eighty-eight percent of eQTL-phenotypic QTL colocalization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots colocalized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc-finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of colocalization of brain eQTL with behavioral QTL, which provides insight into the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological, differences appear to play a key role in this divergence.

Pervasive sex-linked effects on transcription regulation as revealed by expression quantitative trait loci mapping in lake whitefish species pairs (Coregonus sp., Salmonidae)

Mapping of expression quantitative trait loci (eQTL) is a powerful means for elucidating the genetic architecture of gene regulation. Yet, eQTL mapping has not been applied toward investigating the regulation architecture of genes involved in the process of population divergence, ultimately leading to speciation events. Here, we conducted an eQTL mapping experiment to compare the genetic architecture of transcript regulation in adaptive traits, differentiating the recently evolved limnetic (dwarf) and benthic (normal) species pairs of lake whitefish. The eQTL were mapped in three data sets derived from an F(1) hybrid-dwarf backcrossed family: the entire set of 66 genotyped individuals and the two sexes treated separately. We identified strikingly more eQTL in the female data set (174), compared to both male (54) and combined (33) data sets. The majority of these genes were not differentially expressed between male and female progeny of the backcross family, thus providing evidence for a strong pleiotropic sex-linked effect in transcriptomic regulation. The subtelomeric region of a linkage group segregating in females encompassed >50% of all eQTL, which exhibited the most pronounced additive effects. We also conducted a direct comparison of transcriptomic profiles between pure dwarf and normal progeny reared in controlled conditions. We detected 34 differentially expressed transcripts associated with eQTL segregating only in sex-specific data sets and mostly belonging to functional groups that differentiate dwarf and normal whitefish in natural populations. Therefore, these eQTL are not related to interindividual variation, but instead to the adaptive and historical genetic divergence between dwarf and normal whitefish. This study exemplifies how the integration of genetic and transcriptomic data offers a strong means for dissecting the functional genomic response to selection by separating mapping family-specific effects from genetic factors under selection, potentially involved in the phenotypic divergence of natural populations.

Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionists

Recently, the amplified fragment length polymorphism (AFLP) technique has gained a lot of popularity, and is now frequently applied to a wide variety of organisms. Technical specificities of the AFLP procedure have been well documented over the years, but there is on the contrary little or scattered information about the statistical analysis of AFLPs. In this review, we describe the various methods available to handle AFLP data, focusing on four research topics at the population or individual level of analysis: (i) assessment of genetic diversity; (ii) identification of population structure; (iii) identification of hybrid individuals; and (iv) detection of markers associated with phenotypes. Two kinds of analysis methods can be distinguished, depending on whether they are based on the direct study of band presences or absences in AFLP profiles ('band-based' methods), or on allelic frequencies estimated at each locus from these profiles ('allele frequency-based' methods). We investigate the characteristics and limitations of these statistical tools; finally, we appeal for a wider adoption of methodologies borrowed from other research fields, like for example those especially designed to deal with binary data.

Development of silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) genetic maps using microsatellite and AFLP markers and a pseudo-testcross strategy

Silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) are two of the four most important pond-cultured fish species inhabiting the major river basins of China. In the present study, genetic maps of silver carp and bighead carp were constructed using microsatellite and AFLP markers and a two-way pseudo-testcross strategy. To create the maps, 60 individuals were obtained from a cross of a single bighead carp (female) and a single silver carp (male). The silver carp map consisted of 271 markers (48 microsatellites and 223 AFLPs) that were assembled into 27 linkage groups, of which 22 contained at least four markers. The total length of the silver carp map was 952.2 cM, covering 82.8% of the estimated genome size. The bighead carp map consisted of 153 markers (27 microsatellites and 126 AFLPs) which were organized into 30 linkage groups, of which 19 contained at least four markers. The total length of the bighead carp map was 852.0 cM, covering 70.5% of the estimated genome size. Eighteen microsatellite markers were common to both maps. These maps will contribute to discovery of genes and genetic regions controlling traits in the two species of carp.

Parallel evolution of lake whitefish dwarf ecotypes in association with limnological features of their adaptive landscape

Sympatric fish populations observed in many north temperate lakes are among the best models to study the processes of population divergence and adaptive radiation. Despite considerable research on such systems, little is known about the associations between ecological conditions and the extent of ecotypic divergence. In this study, we examined the biotic and abiotic properties of postglacial lakes in which lake whitefish, Coregonus clupeaformis, occur as a derived dwarf ecotype in sympatry with an ancestral normal ecotype. We compared 19 limnological variables between two groups of lakes known from previous studies to harbour sympatric dwarf and normal ecotypes with high and low levels of phenotypic and genetic differentiation respectively. We found clear environmental differences between the two lake groups. Namely, oxygen was the most discriminant variable, where lakes harbouring the most divergent populations were characterized by the greatest hypolimnetic oxygen depletion. These lakes also had lower zooplankton densities and a narrower distribution of zooplantonic prey length. These results suggest that the highest differentiation between sympatric ecotypes occurs in lakes with reduced habitat and prey availability that could increase competition for resources. This in turns supports the hypothesis that parallelism in the extent of phenotypic divergence among sympatric whitefish ecotypes is associated with parallelism in adaptive landscape in terms of differences in limnological characteristics, as well as availability and structure of the zooplanktonic community.