Baocheng Guo

Population genomic evidence for adaptive differentiation in Baltic Sea three-spined sticklebacks

BackgroundThe degree of genetic differentiation among populations experiencing high levels of gene flow is expected to be low for neutral genomic sites, but substantial divergence can occur in sites subject to directional selection. Studies of highly mobile marine fish populations provide an opportunity to investigate this kind of heterogeneous genomic differentiation, but most studies to this effect have focused on a relatively low number of genetic markers and/or few populations. Hence, the patterns and extent of genomic divergence in high-gene-flow marine fish populations remain poorly understood.ResultsWe here investigated genome-wide patterns of genetic variability and differentiation in ten marine populations of three-spined stickleback (Gasterosteus aculeatus) distributed across a steep salinity and temperature gradient in the Baltic Sea, by utilizing >30,000 single nucleotide polymorphisms obtained with a pooled RAD-seq approach. We found that genetic diversity and differentiation varied widely across the genome, and identified numerous fairly narrow genomic regions exhibiting signatures of both divergent and balancing selection. Evidence was uncovered for substantial genetic differentiation associated with both salinity and temperature gradients, and many candidate genes associated with local adaptation in the Baltic Sea were identified.ConclusionsThe patterns of genetic diversity and differentiation, as well as candidate genes associated with adaptation, in Baltic Sea sticklebacks were similar to those observed in earlier comparisons between marine and freshwater populations, suggesting that similar processes may be driving adaptation to brackish and freshwater environments. Taken together, our results provide strong evidence for heterogenic genomic divergence driven by local adaptation in the face of gene flow along an environmental gradient in the post-glacially formed Baltic Sea.

Construction of Ultradense Linkage Maps with Lep-MAP2: Stickleback F2 Recombinant Crosses as an Example

High-density linkage maps are important tools for genome biology and evolutionary genetics by quantifying the extent of recombination, linkage disequilibrium, and chromosomal rearrangements across chromosomes, sexes, and populations. They provide one of the best ways to validate and refine de novo genome assemblies, with the power to identify errors in assemblies increasing with marker density. However, assembly of high-density linkage maps is still challenging due to software limitations. We describe Lep-MAP2, a software for ultradense genome-wide linkage map construction. Lep-MAP2 can handle various family structures and can account for achiasmatic meiosis to gain linkage map accuracy. Simulations show that Lep-MAP2 outperforms other available mapping software both in computational efficiency and accuracy. When applied to two large F2-generation recombinant crosses between two nine-spined stickleback (Pungitius pungitius) populations, it produced two high-density (∼6 markers/cM) linkage maps containing 18,691 and 20,054 single nucleotide polymorphisms. The two maps showed a high degree of synteny, but female maps were 1.5–2 times longer than male maps in all linkage groups, suggesting genome-wide recombination suppression in males. Comparison with the genome sequence of the three-spined stickleback (Gasterosteus aculeatus) revealed a high degree of interspecific synteny with a low frequency (<5%) of interchromosomal rearrangements. However, a fairly large (ca. 10 Mb) translocation from autosome to sex chromosome was detected in both maps. These results illustrate the utility and novel features of Lep-MAP2 in assembling high-density linkage maps, and their usefulness in revealing evolutionarily interesting properties of genomes, such as strong genome-wide sex bias in recombination rates.

Population genomic evidence for adaptive differentiation in the Baltic Sea herring

Detecting and estimating the degree of genetic differentiation among populations of highly mobile marine fish having pelagic larval stages is challenging because their effective population sizes can be large, and thus, little genetic drift and differentiation is expected in neutral genomic sites. However, genomic sites subject to directional selection stemming from variation in local environmental conditions can still show substantial genetic differentiation, yet these signatures can be hard to detect with low‐throughput approaches. Using a pooled RAD‐seq approach, we investigated genomewide patterns of genetic variability and differentiation within and among 20 populations of Atlantic herring in the Baltic Sea (and adjacent Atlantic sites), where previous low‐throughput studies and/or studies based on few populations have found limited evidence for genetic differentiation. Stringent quality control was applied in the filtering of 1 791 254 SNPs, resulting in a final data set of 68 182 polymorphic loci. Clear differentiation was identified between Atlantic and Baltic populations in many genomic sites, while differentiation within the Baltic Sea area was weaker and geographically less structured. However, outlier analyses – whether including all populations or only those within the Baltic Sea – uncovered hundreds of directionally selected loci in which variability was associated with either salinity, temperature or both. Hence, our results support the view that although the degree of genetic differentiation among Baltic Sea herring populations is low, there are many genomic regions showing elevated divergence, apparently as a response to temperature‐ and salinity‐related natural selection. As such, the results add to the increasing evidence of local adaptation in highly mobile marine organisms, and those in the young Baltic Sea in particular.

Genomic divergence between nine- and three-spined sticklebacks.

BackgroundComparative genomics approaches help to shed light on evolutionary processes that shape differentiation between lineages. The nine-spined stickleback (Pungitius pungitius) is a closely related species of the ecological ‘supermodel’ three-spined stickleback (Gasterosteus aculeatus). It is an emerging model system for evolutionary biology research but has garnered less attention and lacks extensive genomic resources. To expand on these resources and aid the study of sticklebacks in a phylogenetic framework, we characterized nine-spined stickleback transcriptomes from brain and liver using deep sequencing.ResultsWe obtained nearly eight thousand assembled transcripts, of which 3,091 were assigned as putative one-to-one orthologs to genes found in the three-spined stickleback. These sequences were used for evaluating overall differentiation and substitution rates between nine- and three-spined sticklebacks, and to identify genes that are putatively evolving under positive selection. The synonymous substitution rate was estimated to be 7.1 × 10-9 per site per year between the two species, and a total of 165 genes showed patterns of adaptive evolution in one or both species. A few nine-spined stickleback contigs lacked an obvious ortholog in three-spined sticklebacks but were found to match genes in other fish species, suggesting several gene losses within 13 million years since the divergence of the two stickleback species. We identified 47 SNPs in 25 different genes that differentiate pond and marine ecotypes. We also identified 468 microsatellites that could be further developed as genetic markers in nine-spined sticklebacks.ConclusionWith deep sequencing of nine-spined stickleback cDNA libraries, our study provides a significant increase in the number of gene sequences and microsatellite markers for this species, and identifies a number of genes showing patterns of adaptive evolution between nine- and three-spined sticklebacks. We also report several candidate genes that might be involved in differential adaptation between marine and freshwater nine-spined sticklebacks. This study provides a valuable resource for future studies aiming to identify candidate genes underlying ecological adaptation in this and other stickleback species.

The Roles and Evolutionary Patterns of Intronless Genes in Deuterostomes

Genes without introns are a characteristic feature of prokaryotes, but there are still a number of intronless genes in eukaryotes. To study these eukaryotic genes that have prokaryotic architecture could help to understand the evolutionary patterns of related genes and genomes. Our analyses revealed a number of intronless genes that reside in 6 deuterostomes (sea urchin, sea squirt, zebrafish, chicken, platypus, and human). We also determined the conservation for each intronless gene in archaea, bacteria, fungi, plants, metazoans, and other eukaryotes. Proportions of intronless genes that are inherited from the common ancestor of archaea, bacteria, and eukaryotes in these species were consistent with their phylogenetic positions, with more proportions of ancient intronless genes residing in more primitive species. In these species, intronless genes belong to different cellular roles and gene ontology (GO) categories, and some of these functions are very basic. Part of intronless genes is derived from other intronless genes or multiexon genes in each species. In conclusion, we showed that a varying number and proportion of intronless genes reside in these 6 deuterostomes, and some of them function importantly. These genes are good candidates for subsequent functional and evolutionary analyses specifically.

Hox Genes of the Japanese Eel Anguilla japonica and Hox Cluster Evolution in Teleosts

Compared with other diploid teleosts (2n=48), anguilloid fish have a specialized karyotype (2n=38) and remarkable morphological variation, and represent one basal group species of teleosts. To investigate the Hox gene/cluster inventory in basal teleosts, a PCR-based survey of Hox genes in the Japanese eel (Anguilla japonica) was conducted with both gene-specific and homeobox-targeted degenerate primers. Our data provide evidence that at least 34 distinct Hox genes exist in the Japanese eel genome and that they represent eight Hox clusters. Duplication of Hox genes in the Japanese eel appears to be the result of the fish-specific genome duplication (FSGD) event. The Japanese eel shared the FSGD event with other teleosts such as zebrafish and pufferfish. A member of Hox paralog group one (HoxA1b) was preserved in the Japanese eel but was lost in other teleosts. Available Hox data revealed that the Hox cluster evolved distinctly in different teleost lineages. All duplicated Hox clusters were retained after the FSGD event in basal teleosts like in the Japanese eel, whereas crown teleosts lost one cluster (HoxCb or HoxDb). Based on current teleostean phylogeny, the HoxDb cluster was lost independently in the teleost lineages Otocephala and Euteleostei.

Integrating multi-origin expression data improves the resolution of deep phylogeny of ray-finned fish (Actinopterygii)

The actinopterygians comprise nearly one-half of all extant vertebrate species and are very important for human well-being. However, the phylogenetic relationships among certain groups within the actinopterygians are still uncertain, and debates about these relationships have continued for a long time. Along with the progress achieved in sequencing technologies, phylogenetic analyses based on multi-gene sequences, termed phylogenomic approaches, are becoming increasingly common and often result in well-resolved and highly supported phylogenetic hypotheses. Based on the transcriptome sequences generated in this study and the extensive expression data currently available from public databases, we obtained alignments of 274 orthologue groups for 26 scientifically and commercially important actinopterygians, representing 17 out of 44 orders within the class Actinopterygii. Using these alignments and probabilistic methods, we recovered relationships between basal actinopterygians and teleosts, among teleosts within protacanthopterygians and related lineages, and also within acanthomorphs. These relationships were recovered with high confidence.

Duplicated Gene Evolution Following Whole-Genome Duplication in Teleost Fish

Baocheng Guo1,2,3, Andreas Wagner1,2 and Shunping He3* 1 Institute of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich 2 The Swiss Institute of Bioinformatics, Quartier Sorge-Batiment Genopode, Lausanne 3Fish Phylogenetics and Biogeography Group, Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 1,2Switzerland 3PR China

Sox genes evolution in closely related young tetraploid cyprinid fishes and their diploid relative

Previous study and analysis of cytochrome b suggested that polyploidization event in the genus Tor occurred about 10 Mya ago. In order to understand evolutionary fates of Sox gene in the early stage of genome duplication at the nucleotide level, PCR surveys for Sox genes in three closely related cyprinid fishes T. douronensis (2n=100), T. qiaojiensis (2n=?), T. sinensis (2n=100) and their relative T. brevifilis (2n=50) were performed. Totally, 52 distinct Sox genes were obtained in these four species, representing SoxB, SoxC, and SoxE group. As expected, isoforms of some Sox genes correspond with the ploidy of species, such as two copies of Sox9a exist in tetraploid species. Analysis indicated that duplicated Sox gene pairs caused by polyploidization evolved independently of each other within polyploid species. Results of substitution rate showed nearly equal rate of nonsynonymous substitution of duplicated Sox orthologs among different polyploid species and their diploid relative orthologs, suggesting at the early stage of genome duplicated Sox orthologs are under similar selective constraints in different polyploidy species and their diploid relative at the amino acid level. All PCR fragments of Sox genes obtained in this study are not accompanied by obvious increase in mutations and pseudogene formation which means that they are under strong purifying selection, suggesting that they are functional at the DNA level. Genealogical analysis revealed that T. qiaojiensis was tetraploid, and T. douronensis, T. qiaojiensis as well as T. sinensis had an allotetraploid ancestor.

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F2 interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h2 ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.