Characterisation and expression of the paired box protein 7 (Pax7) gene in polymorphic Arctic charr (Salvelinus alpinus)

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Differential gene expression during early development in recently evolved and sympatric Arctic charr morphs

Phenotypic differences between closely related taxa or populations can arise through genetic variation or be environmentally induced, leading to altered transcription of genes during development. Comparative developmental studies of closely related species or variable populations within species can help to elucidate the molecular mechanisms related to evolutionary divergence and speciation. Studies of Arctic charr (Salvelinus alpinus) and related salmonids have revealed considerable phenotypic variation among populations and in Arctic charr many cases of extensive variation within lakes (resource polymorphism) have been recorded. One example is the four Arctic charr morphs in the ∼10,000 year old Lake Thingvallavatn, which differ in numerous morphological and life history traits. We set out to investigate the molecular and developmental roots of this polymorphism by studying gene expression in embryos of three of the morphs reared in a common garden set-up. We performed RNA-sequencing, de-novo transcriptome assembly and compared gene expression among morphs during an important timeframe in early development, i.e., preceding the formation of key trophic structures. Expectedly, developmental time was the predominant explanatory variable. As the data were affected by some form of RNA-degradation even though all samples passed quality control testing, an estimate of 3'-bias was the second most common explanatory variable. Importantly, morph, both as an independent variable and as interaction with developmental time, affected the expression of numerous transcripts. Transcripts with morph effect, separated the three morphs at the expression level, with the two benthic morphs being more similar. However, Gene Ontology analyses did not reveal clear functional enrichment of transcripts between groups. Verification via qPCR confirmed differential expression of several genes between the morphs, including regulatory genes such as AT-Rich Interaction Domain 4A (arid4a) and translin (tsn). The data are consistent with a scenario where genetic divergence has contributed to differential expression of multiple genes and systems during early development of these sympatric Arctic charr morphs.

The duplicated paired box protein 7 (pax7) genes differentially transcribed during Japanese flounder (Paralichthys olivaceus) embryogenesis

PAX are important regulators of developmental processes. PAX7 plays crucial roles in patterning of the dorsal central nervous system (CNS), neural crest (NC), and skeletal muscle. Here, we identified six spliced isoforms of pax7a and one pax7b and characterized their expression patterns. All of flounder Pax7a-1, Pax7a-2, Pax7a-3, and Pax7b contain a conserved paired domain (PD), an octapeptide motif (OP), and a paired type homeodomain (HD). However, the PD of Pax7a-4 and the HD of Pax7a-5 are not intact, and there is no HD in Pax7a-4 and Pax7a-6. pax7a and pax7b show distinct spatiotemporal expression patterns during embryogenesis. Whole-mount in situ hybridization demonstrates that the expression patterns of pax7a and pax7b are overlapping but distinguishable in the dorsal central nervous system. pax7a is expressed in most part of the brain and the neural tube, while pax7b is expressed exclusively in the diencephalon and the midbrain. In addition, pax7a is also expressed in the cranial NC and the trunk NC. RT-PCR results show that there were different expression patterns between the different isoforms. These results indicate subfunction partitioning of the duplicated pax7 genes. The duplicated pax7 may provide additional flexibility in fine-tuning neurogenesis and somitogenesis.

Oncorhynchus mykiss pax7 sequence variations with comparative analyses against other teleost species

The paired box-7 (pax7) transcription factor expressed in satellite cells (SCs) is an essential regulator of skeletal muscle growth and regeneration in vertebrates including fish. Characterization of rainbow trout (Oncorhynchus mykiss) pax7 gene/s may offer novel insights into skeletal myogenesis by SCs in this indeterminate growth species. Further, evaluation of promoters for cis-regulatory regions may shed light on the evolutionary fate of the duplicated genes. Employing standard PCR, cloning and computational approach, we identified and report complete coding sequences of two pax7 paralogs of rainbow trout (rt); rtpax7α and rtpax7β. Both genes show significant identity in the nucleotide (97%) and the predicted amino acid (98%) sequences, and bear the characteristic paired domain (PD), octapeptide (OP) and homeodomain (HD) motifs. We further report several splice variants of each gene and nucleotide differences in coding sequence that predicts six putative amino acid changes between the two genes. Additionally, we noted a trinucleotide deletion in rtpax7β that results in putative serine elimination at the N-terminus and a single nucleotide polymorphism (SNP) in majority of the rtpax7β variants (6/10) that predicts an arginine substitution for a lysine. We also deciphered the genomic organization up to the first three exons and the upstream putative promoter regions of both genes. Comparative in silico analysis of both the trout pax7 promoters with that of zebrafish pax7 duplicates; zfpax7a and zfpax7b; predicts several important cis-elements/transcription factor binding sites (TFBS) in these teleost pax7 promoter regions.

Transcriptional dynamics of a conserved gene expression network associated with craniofacial divergence in Arctic charr

Background

Understanding the molecular basis of craniofacial variation can provide insights into key developmental mechanisms of adaptive changes and their role in trophic divergence and speciation. Arctic charr (Salvelinus alpinus) is a polymorphic fish species, and, in Lake Thingvallavatn in Iceland, four sympatric morphs have evolved distinct craniofacial structures. We conducted a gene expression study on candidates from a conserved gene coexpression network, focusing on the development of craniofacial elements in embryos of two contrasting Arctic charr morphotypes (benthic and limnetic).

Results

Four Arctic charr morphs were studied: one limnetic and two benthic morphs from Lake Thingvallavatn and a limnetic reference aquaculture morph. The presence of morphological differences at developmental stages before the onset of feeding was verified by morphometric analysis. Following up on our previous findings that Mmp2 and Sparc were differentially expressed between morphotypes, we identified a network of genes with conserved coexpression across diverse vertebrate species. A comparative expression study of candidates from this network in developing heads of the four Arctic charr morphs verified the coexpression relationship of these genes and revealed distinct transcriptional dynamics strongly correlated with contrasting craniofacial morphologies (benthic versus limnetic). A literature review and Gene Ontology analysis indicated that a significant proportion of the network genes play a role in extracellular matrix organization and skeletogenesis, and motif enrichment analysis of conserved noncoding regions of network candidates predicted a handful of transcription factors, including Ap1 and Ets2, as potential regulators of the gene network. The expression of Ets2 itself was also found to associate with network gene expression. Genes linked to glucocorticoid signalling were also studied, as both Mmp2 and Sparc are responsive to this pathway. Among those, several transcriptional targets and upstream regulators showed differential expression between the contrasting morphotypes. Interestingly, although selected network genes showed overlapping expression patterns in situ and no morph differences, Timp2 expression patterns differed between morphs.

Conclusion

Our comparative study of transcriptional dynamics in divergent craniofacial morphologies of Arctic charr revealed a conserved network of coexpressed genes sharing functional roles in structural morphogenesis. We also implicate transcriptional regulators of the network as targets for future functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-5-40) contains supplementary material, which is available to authorized users.

Regional expression of Pax7 in the brain of Xenopus laevis during embryonic and larval development

Pax7 is a member of the highly conserved Pax gene family that is expressed in restricted zones of the central nervous system (CNS) during development, being involved in early brain regionalization and the maintenance of the regional identity. Using sensitive immunohistochemical techniques we have analyzed the spatiotemporal pattern of Pax7 expression in the brain of the anuran amphibian Xenopus laevis, during development. Pax7 expression was first detected in early embryos in the basal plate of prosomere 3, roof and alar plates of prosomere 1 and mesencephalon, and the alar plate of rhombomere 1. As development proceeded, Pax7 cells were observed in the hypothalamus close to the catecholaminergic population of the mammillary region. In the diencephalon, Pax7 was intensely expressed in a portion of the basal plate of prosomere 3, in the roof plate and in scattered cells of the thalamus in prosomere 2, throughout the roof of prosomere 1, and in the commissural and juxtacommissural domains of the pretectum. In the mesencephalon, Pax7 cells were localized in the optic tectum and, to a lesser extent, in the torus semicircularis. The rostral portion of the alar part of rhombomere 1, including the ventricular layer of the cerebellum, expressed Pax7 and, gradually, some of these dorsal cells were observed to populate ventrally the interpeduncular nucleus and the isthmus (rhombomere 0). Additionally, Pax7 positive cells were found in the ventricular zone of the ventral part of the alar plate along the rhombencephalon and the spinal cord. The findings show that the strongly conserved features of Pax7 expression through development shared by amniote vertebrates are also present in the anamniote amphibians as a common characteristic of the brain organization of tetrapods.

Differentiation at the MHCIIα and Cath2 loci in sympatric Salvelinus alpinus resource morphs in Lake Thingvallavatn

Northern freshwater fish may be suitable for the genetic dissection of ecological traits because they invaded new habitats after the last ice age (∼10.000 years ago). Arctic charr (Salvelinus alpinus) colonizing streams and lakes in Iceland gave rise to multiple populations of small benthic morphotypes, often in sympatry with a pelagic morphotype. Earlier studies have revealed significant, but subtle, genetic differentiation between the three most common morphs in Lake Thingvallavatn. We conducted a population genetic screen on four immunological candidate genes Cathelicidin 2 (Cath2), Hepcidin (Hamp), Liver expressed antimicrobial peptide 2a (Leap-2a), and Major Histocompatibility Complex IIα (MHCIIα) and a mitochondrial marker (D-loop) among the three most common Lake Thingvallavatn charr morphs. Significant differences in allele frequencies were found between morphs at the Cath2 and MHCIIα loci. No such signal was detected in the D-loop nor in the other two immunological genes. In Cath2 the small benthic morph deviated from the other two (FST  = 0.13), one of the substitutions detected constituting an amino acid replacement polymorphism in the antimicrobial peptide. A more striking difference was found in the MHCIIα. Two haplotypes were very common in the lake, and their frequency differed greatly between the morphotypes (from 22% to 93.5%, FST  = 0.67). We then expanded our study by surveying the variation in Cath2 and MHCIIα in 9 Arctic charr populations from around Iceland. The populations varied greatly in terms of allele frequencies at Cath2, but the variation did not correlate with morphotype. At the MHCIIα locus, the variation was nearly identical to the variation in the two benthic morphs of Lake Thingvallavatn. The results are consistent with a scenario where parts of the immune systems have diverged substantially among Arctic charr populations in Iceland, after colonizing the island ∼10.000 years ago.

Validation of reference genes for expression studies during craniofacial development in arctic charr

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes.

Temporal gene expression variation associated with eyespot size plasticity in Bicyclus anynana

Seasonal polyphenism demonstrates an organism's ability to respond to predictable environmental variation with alternative phenotypes, each presumably better suited to its respective environment. However, the molecular mechanisms linking environmental variation to alternative phenotypes via shifts in development remain relatively unknown. Here we investigate temporal gene expression variation in the seasonally polyphenic butterfly Bicyclus anynana. This species shows drastic changes in eyespot size depending on the temperature experienced during larval development. The wet season form (larvae reared over 24°C) has large ventral wing eyespots while the dry season form (larvae reared under 19°C) has much smaller eyespots. We compared the expression of three proteins, Notch, Engrailed, and Distal-less, in the future eyespot centers of the two forms to determine if eyespot size variation is associated with heterochronic shifts in the onset of their expression. For two of these proteins, Notch and Engrailed, expression in eyespot centers occurred earlier in dry season than in wet season larvae, while Distal-less showed no temporal difference between the two forms. These results suggest that differences between dry and wet season adult wings could be due to a delay in the onset of expression of these eyespot-associated genes. Early in eyespot development, Notch and Engrailed may be functioning as repressors rather than activators of the eyespot gene network. Alternatively, temporal variation in the onset of early expressed genes between forms may have no functional consequences to eyespot size regulation and may indicate the presence of an 'hourglass' model of development in butterfly eyespots.

Activation of Pax7-positive cells in a non-contractile tissue contributes to regeneration of myogenic tissues in the electric fish S. macrurus

The ability to regenerate tissues is shared across many metazoan taxa, yet the type and extent to which multiple cellular mechanisms come into play can differ across species. For example, urodele amphibians can completely regenerate all lost tissues, including skeletal muscles after limb amputation. This remarkable ability of urodeles to restore entire limbs has been largely linked to a dedifferentiation-dependent mechanism of regeneration. However, whether cell dedifferentiation is the fundamental factor that triggers a robust regeneration capacity, and whether the loss or inhibition of this process explains the limited regeneration potential in other vertebrates is not known. Here, we studied the cellular mechanisms underlying the repetitive regeneration of myogenic tissues in the electric fish S. macrurus. Our in vivo microinjection studies of high molecular weight cell lineage tracers into single identified adult myogenic cells (muscle or noncontractile muscle-derived electrocytes) revealed no fragmentation or cellularization proximal to the amputation plane. In contrast, ultrastructural and immunolabeling studies verified the presence of myogenic stem cells that express the satellite cell marker Pax7 in mature muscle fibers and electrocytes of S. macrurus. These data provide the first example of Pax-7 positive muscle stem cells localized within a non-contractile electrogenic tissue. Moreover, upon amputation, Pax-7 positive cells underwent a robust replication and were detected exclusively in regions that give rise to myogenic cells and dorsal spinal cord components revealing a regeneration process in S. macrurus that is dependent on the activation of myogenic stem cells for the renewal of both skeletal muscle and the muscle-derived electric organ. These data are consistent with the emergent concept in vertebrate regeneration that different tissues provide a distinct progenitor cell population to the regeneration blastema, and these progenitor cells subsequently restore the original tissue.

The Pax3 and Pax7 paralogs cooperate in neural and neural crest patterning using distinct molecular mechanisms, in Xenopus laevis embryos

Pax3 and Pax7 paralogous genes have functionally diverged in vertebrate evolution, creating opportunity for a new distribution of roles between the two genes and the evolution of novel functions. Here we focus on the regulation and function of Pax7 in the brain and neural crest of amphibian embryos, which display a different pax7 expression pattern, compared to the other vertebrates already described. Pax7 expression is restricted to the midbrain, hindbrain and anterior spinal cord, and Pax7 activity is important for maintaining the fates of these regions, by restricting otx2 expression anteriorly. In contrast, pax3 displays broader expression along the entire neuraxis and Pax3 function is important for posterior brain patterning without acting on otx2 expression. Moreover, while both genes are essential for neural crest patterning, we show that they do so using two distinct mechanisms: Pax3 acts within the ectoderm which will be induced into neural crest, while Pax7 is essential for the inducing activity of the paraxial mesoderm towards the prospective neural crest.