Finding Nemo’s Genes: A chromosome‐scale reference assembly of the genome of the orange clownfish Amphiprion percula

The iconic orange clownfish, Amphiprion percula, is a model organism for studying the ecology and evolution of reef fishes, including patterns of population connectivity, sex change, social organization, habitat selection and adaptation to climate change. Notably, the orange clownfish is the only reef fish for which a complete larval dispersal kernel has been established and was the first fish species for which it was demonstrated that antipredator responses of reef fishes could be impaired by ocean acidification. Despite its importance, molecular resources for this species remain scarce and until now it lacked a reference genome assembly. Here, we present a de novo chromosome‐scale assembly of the genome of the orange clownfish Amphiprion percula. We utilized single‐molecule real‐time sequencing technology from Pacific Biosciences to produce an initial polished assembly comprised of 1,414 contigs, with a contig N50 length of 1.86 Mb. Using Hi‐C‐based chromatin contact maps, 98% of the genome assembly were placed into 24 chromosomes, resulting in a final assembly of 908.8 Mb in length with contig and scaffold N50s of 3.12 and 38.4 Mb, respectively. This makes it one of the most contiguous and complete fish genome assemblies currently available. The genome was annotated with 26,597 protein‐coding genes and contains 96% of the core set of conserved actinopterygian orthologs. The availability of this reference genome assembly as a community resource will further strengthen the role of the orange clownfish as a model species for research on the ecology and evolution of reef fishes.     

Overview on the Diversity of Sounds Produced by Clownfishes (Pomacentridae): Importance of Acoustic Signals in Their Peculiar Way of Life

Background

Clownfishes (Pomacentridae) are brightly colored coral reef fishes well known for their mutualistic symbiosis with tropical sea anemones. These fishes live in social groups in which there is a size-based dominance hierarchy. In this structure where sex is socially controlled, agonistic interactions are numerous and serve to maintain size differences between individuals adjacent in rank. Clownfishes are also prolific callers whose sounds seem to play an important role in the social hierarchy. Here, we aim to review and to synthesize the diversity of sounds produced by clownfishes in order to emphasize the importance of acoustic signals in their way of life.

Methodology/Principal Findings

Recording the different acoustic behaviors indicated that sounds are divided into two main categories: aggressive sounds produced in conjunction with threat postures (charge and chase), and submissive sounds always emitted when fish exhibited head shaking movements (i.e. a submissive posture). Both types of sounds showed size-related intraspecific variation in dominant frequency and pulse duration: smaller individuals produce higher frequency and shorter duration pulses than larger ones, and inversely. Consequently, these sonic features might be useful cues for individual recognition within the group. This observation is of significant importance due to the size-based hierarchy in clownfish group. On the other hand, no acoustic signal was associated with the different reproductive activities.

Conclusions/Significance

Unlike other pomacentrids, sounds are not produced for mate attraction in clownfishes but to reach and to defend the competition for breeding status, which explains why constraints are not important enough for promoting call diversification in this group.     

Normal development of the tomato clownfish Amphiprion frenatus: live imaging and in situ hybridization analyses of mesodermal and neurectodermal development

The normal embryonic development of the tomato clownfish Amphiprion frenatus was analysed using live imaging and by in situ hybridization for detection of mesodermal and neurectodermal development. Both morphology of live embryos and tissue‐specific staining revealed significant differences in the gross developmental programme of A. frenatus compared with better‐known teleost fish models, in particular, initiation of somitogenesis before complete epiboly, initiation of narrowing of the neurectoderm (neurulation) before somitogenesis, relatively early pigmentation of melanophores at the 10–15 somite stage and a distinctive pattern of melanophore distribution. These results suggest evolutionary adaptability of the teleost developmental programme. The ease of obtaining eggs, in vitro culture of the embryo, in situ staining analyses and these reported characteristics make A. frenatus a potentially important model marine fish species for studying embryonic development, physiology, ecology and evolution.     

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.     

Genomic variation at the tips of the adaptive radiation of Darwin's finches

Adaptive radiation unfolds as selection acts on the genetic variation underlying functional traits. The nature of this variation can be revealed by studying the tips of an ongoing adaptive radiation. We studied genomic variation at the tips of the Darwin's finch radiation; specifically focusing on polymorphism within, and variation among, three sympatric species of the genus Geospiza. Using restriction site‐associated DNA (RAD‐seq), we characterized 32 569 single‐nucleotide polymorphisms (SNPs), from which 11 outlier SNPs for beak and body size were uncovered by a genomewide association study (GWAS). Principal component analysis revealed that these 11 SNPs formed four statistically linked groups. Stepwise regression then revealed that the first PC score, which included 6 of the 11 top SNPs, explained over 80% of the variation in beak size, suggesting that selection on these traits influences multiple correlated loci. The two SNPs most strongly associated with beak size were near genes associated with beak morphology across deeper branches of the radiation: delta‐like 1 homologue (DLK1) and high‐mobility group AT‐hook 2 (HMGA2). Our results suggest that (i) key adaptive traits are associated with a small fraction of the genome (11 of 32 569 SNPs), (ii) SNPs linked to the candidate genes are dispersed throughout the genome (on several chromosomes), and (iii) micro‐ and macro‐evolutionary variation (roots and tips of the radiation) involve some shared and some unique genomic regions.     

A de novo chromosome‐level genome assembly of Coregonus sp. “Balchen”: One representative of the Swiss Alpine whitefish radiation

Salmonids are of particular interest to evolutionary biologists due to their incredible diversity of life‐history strategies and the speed at which many salmonid species have diversified. In Switzerland alone, over 30 species of Alpine whitefish from the subfamily Coregoninae have evolved since the last glacial maximum, with species exhibiting a diverse range of morphological and behavioural phenotypes. This, combined with the whole genome duplication which occurred in the ancestor of all salmonids, makes the Alpine whitefish radiation a particularly interesting system in which to study the genetic basis of adaptation and speciation and the impacts of ploidy changes and subsequent rediploidization on genome evolution. Although well‐curated genome assemblies exist for many species within Salmonidae, genomic resources for the subfamily Coregoninae are lacking. To assemble a whitefish reference genome, we carried out PacBio sequencing from one wild‐caught Coregonus sp. “Balchen” from Lake Thun to ~90× coverage. PacBio reads were assembled independently using three different assemblers, falcon , canu and wtdbg2 and subsequently scaffolded with additional Hi‐C data. All three assemblies were highly contiguous, had strong synteny to a previously published Coregonus linkage map, and when mapping additional short‐read data to each of the assemblies, coverage was fairly even across most chromosome‐scale scaffolds. Here, we present the first de novo genome assembly for the Salmonid subfamily Coregoninae. The final 2.2‐Gb wtdbg2 assembly included 40 scaffolds, an N50 of 51.9 Mb and was 93.3% complete for BUSCOs. The assembly consisted of ~52% transposable elements and contained 44,525 genes.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

A transcriptomic scan for positively selected genes in two closely related marine fishes: Sebastes caurinus and S. rastrelliger

Comparative genomic analyses can provide valuable insight into functional evolutionary divergence among closely related species. Here we employ a comparative evolutionary analysis of expressed sequence tags (ESTs) from two closely related species of marine fishes (genus Sebastes--rockfish). Sebastes is a highly diverse group of marine fishes that inhabit a wide array of marine habitats and the study of this group can provide insights into speciation in the marine environment. ESTs were developed for S. caurinus (23,668 from brain, kidney, and spleen tissues) and S. rastrelliger (11,207 from brain and pituitary tissues). Following assembly we were able to identify, with high confidence, 257 orthologous sequence pairs between the two species through a reciprocal best hit blast search. An analysis of functional divergence between orthologs revealed that 19.46% had Ka/Ks values greater than 0.5 and 8.17% had Ka/Ks values greater than one, identifying a large pool of candidate genes to further study adaptive divergence in the group. Genes with elevated Ka/Ks values belonged to the following functional categories: immune function, metabolism, longevity, and reproductive behavior, indicating that adaptive divergence in these functional groups may be important in the diversification of this group of fishes. This study provides the ground work to better understand the molecular evolution of genes involved in a radiation of marine fishes.     

The adaptive genomic landscape of beak morphology in Darwin's finches

Beak shape in Darwin's ground finches (Geospiza) is emblematic of natural selection and adaptive radiation, yet our understanding of the genetic basis of beak shape variation, and thus the genetic target of natural selection, is still evolving. Here we reveal the genomic architecture of beak shape variation using genomewide comparisons of four closely related and hybridizing species across 13 islands subject to parallel natural selection. Pairwise contrasts among species were used to identify a large number of genomic loci that are consistently related to species differences across a complex landscape. These loci are associated with hundreds of genes that have enriched GO categories significantly associated with development. One genomic region of particular interest is a section of Chromosome 1A with many candidate genes and increased linkage. The distinct, pointed beak shape of the cactus finch is linked to an excess of intermediate frequency alleles and increased heterozygosity in significant SNPs, but not across the rest of the genome. Alleles associated with pointier beaks among species were associated with pointier‐beaked populations within each species, thus establishing a common basis for natural selection, species divergence and adaptive radiation. The adaptive genomic landscape for Darwin's finches mirrors theoretical expectations based on morphological variation. The implication that a large number of genes are actively maintained to facilitate beak variation across parallel populations with documented interspecies admixture challenges our understanding of evolutionary processes in the wild.     

Diversification and adaptive evolution of putative sweet taste receptors in threespine stickleback

The threespine stickleback Gasterosteus aculeatus is known to include several morphologically and ecologically divergent forms. Its phenotypic traits related to feeding vary among forms, and are considered to be a result of adaptations to various environments to find foods effectively. To examine whether the diversification of feeding modes in the stickleback involves genetic changes of the sense of taste, taste receptor family 1 (T1R) genes in stickleback were analyzed and compared with those in other model fishes. Ten T1R genes and 2 pseudogenes were identified from the stickleback genomic sequences. In particular, putative sweet taste receptors (T1R2s) highly increased in number in stickleback (8 genes and 2 pseudogenes) compared to other fishes (2-3 genes). Maximum likelihood estimations of nonsynonymous-synonymous nucleotide substitution rate have indicated that stickleback T1R2 are under positive selection. Expression analysis by RT-PCR revealed that most stickleback T1R genes were expressed in the taste organs; however, at least two T1R2 genes were not expressed in the taste organs, suggesting that the expression levels of these T1R2 genes may be fluctuated through the life history. In addition, sequencing analysis showed that several T1R2 genes in an anadromous form stickleback individual collected from the western Pacific (Japan) were substantially different from those in genomic data derived from a freshwater form individual collected in North America. This suggested that intra-specific variations of stickleback T1R2 genes were considerably large. Our results imply that, in stickleback, T1R2s have diversified through adaptation to various environments, probably to perceive substances important for its survival and reproduction.     

温度和体重对克氏双锯鱼仔鱼代谢率的影响

采用实验生态学的方法研究了不同温度和体重对克氏双锯鱼仔鱼呼吸和排泄的影响。结果表明,在23、26和29℃下,仔鱼个体日常呼吸率和排氨率与鱼体重呈幂函数(R=aWb)关系,即鱼体重越大其呼吸率和排氨率越大,b值随温度升高而增大,分别为0.8873,0.9033和0.9323(呼吸率),以及0.7625,0.8012和0.8278(排氨率)。温度和体重对克氏双锯鱼仔鱼个体呼吸率和排氨率的影响可用复合线性公式表示:RR=0.042(±0.007)W0.889(±0.026)×e0.122(±0.005)T,ER=0.002(±0.000)W0.797(±0.029)×e0.115(±0.007)T。比体重呼吸率和排氨率在相同温度条件下随个体增长而降低;在整个仔鱼期,比体重呼吸率和排氨率随温度升高而增加。克氏双锯鱼仔鱼呼吸和排泄Q10值在26—29℃较低,其可能是克氏双锯鱼仔鱼生长发育的最佳温度范围。克氏双锯鱼仔鱼在温度23—29℃时O/N范围为52—57,表明在温度23—29℃时克氏双锯鱼仔鱼代谢底物除了蛋白质外,脂肪和碳水化合物为能源也占了比较大的比例。     

Genomic analyses provide comprehensive insights into the domestication of bast fiber crop ramie (Boehmeria nivea)

Ramie (Boehmeria nivea) is an economically important natural fiber-producing crop that has been cultivated for thousands of years in China; however, the evolution of this crop remains largely unknown. Here, we report a ramie domestication analysis based on genome assembly and resequencing of cultivated and wild accessions. Two chromosome-level genomes representing wild and cultivated ramie were assembled de novo. Numerous structural variations between two assemblies, together with the genetic variations from population resequencing, constituted a comprehensive genomic variation map for ramie. Domestication analysis identified 71 high-confidence selective sweeps comprising 320 predicted genes, and 29 genes from sweeps were associated with fiber growth in the expression. In addition, we identified seven genetic loci associated with the fiber yield trait in the segregated population derived from the crossing of two assembled accessions, and two of which showed an overlap with the selective sweeps. These findings indicated that bast fiber traits were focused on during the domestication history of ramie. This study sheds light on the domestication of ramie and provides a valuable resource for biological and breeding studies of this important crop.     

River hydrological seasonality influences life history strategies of tropical riverine fishes

Under a particular set of selective forces, specific combinations of traits (strategies) will be favored in a given population, within the particular constraints of the considered species. For fishes, three demographic strategies have been suggested to result from adaptive responses to environmental predictability (i.e., seasonality): periodic, opportunistic and equilibrium [Winemiller KO, Rose KA (1992) Patterns of life-history diversification in North American fishes: implications for population regulation. Can J Fish Aquat Sci 49:2196–2218]. These strategies optimize fitness within predictable, unpredictable and stable systems, respectively. We tested these predictions of life history trait distribution along a gradient of hydrologic seasonality in West African tropical rivers at the drainage basin scale. We used logistic regression of species presence–absence data to test whether dominant life history traits of species caused community compositional change in response to a gradient of seasonality in hydrologic regime across basins. After accounting for taxonomic relatedness, species body size and statistical redundancy inherent to related traits, we found a higher proportion of species producing a great number of small oocytes, reproducing within a short period of time and presenting a low degree of parental care (the periodic strategy) in highly seasonal drainage basins (e.g., rivers with a short and predictable favorable season). Conversely, in more stable drainage basins (e.g., rivers with a wet season of several months), we observed a greater proportion of species producing small numbers of large oocytes, reproducing within a long period of time and providing parental care to their offspring (the equilibrium strategy). Our results suggest that distributions of tropical freshwater fishes at the drainage basin scale can be partly explained by the match between life history strategies and seasonality gradients in hydrological conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Complete Mitochondrial DNA Sequences of the Threadfin Cichlid (Petrochromis trewavasae) and the Blunthead Cichlid (Tropheus moorii) and Patterns of Mitochondrial Genome Evolution in Cichlid Fishes

The cichlid fishes of the East African Great Lakes represent a model especially suited to study adaptive radiation and speciation. With several African cichlid genome projects being in progress, a promising set of closely related genomes is emerging, which is expected to serve as a valuable data base to solve questions on genotype-phenotype relations. The mitochondrial (mt) genomes presented here are the first results of the assembly and annotation process for two closely related but eco-morphologically highly distinct Lake Tanganyika cichlids, Petrochromis trewavasae and Tropheus moorii. The genomic sequences comprise 16,588 bp (P. trewavasae) and 16,590 bp (T. moorii), and exhibit the typical mitochondrial structure, with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a non-coding control region. Analyses confirmed that the two species are very closely related with an overall sequence similarity of 96%. We analyzed the newly generated sequences in the phylogenetic context of 21 published labroid fish mitochondrial genomes. Consistent with other vertebrates, the D-loop region was found to evolve faster than protein-coding genes, which in turn are followed by the rRNAs; the tRNAs vary greatly in the rate of sequence evolution, but on average evolve the slowest. Within the group of coding genes, ND6 evolves most rapidly. Codon usage is similar among examined cichlid tribes and labroid families; although a slight shift in usage patterns down the gene tree could be observed. Despite having a clearly different nucleotide composition, ND6 showed a similar codon usage. C-terminal ends of Cox1 exhibit variations, where the varying number of amino acids is related to the structure of the obtained phylogenetic tree. This variation may be of functional relevance for Cox1 synthesis.     

Host specialist clownfishes are environmental niche generalists

Why generalist and specialist species coexist in nature is a question that has interested evolutionary biologists for a long time. While the coexistence of specialists and generalists exploiting resources on a single ecological dimension has been theoretically and empirically explored, biological systems with multiple resource dimensions (e.g. trophic, ecological) are less well understood. Yet, such systems may provide an alternative to the classical theory of stable evolutionary coexistence of generalist and specialist species on a single resource dimension. We explore such systems and the potential trade-offs between different resource dimensions in clownfishes. All species of this iconic clade are obligate mutualists with sea anemones yet show interspecific variation in anemone host specificity. Moreover, clownfishes developed variable environmental specialization across their distribution. In this study, we test for the existence of a relationship between host-specificity (number of anemones associated with a clownfish species) and environmental-specificity (expressed as the size of the ecological niche breadth across climatic gradients). We find a negative correlation between host range and environmental specificities in temperature, salinity and pH, probably indicating a trade-off between both types of specialization forcing species to specialize only in a single direction. Trade-offs in a multi-dimensional resource space could be a novel way of explaining the coexistence of generalist and specialists.