Comparison of recombinant Culicoides allergens produced in different expression systems for IgE serology of insect bite hypersensitivity in horses of different origins.

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses, an IgE-mediated reaction to Culicoides midges. Causative Culicoides spp. are not indigenous in Iceland resulting in high prevalence of IBH in horses born in Iceland and exported as compared to Icelandic horses born in a Culicoides rich environment. The aims were (i) to compare IgE levels in sera of IBH-affected horses born in Iceland (n = 47) with horses of the Icelandic breed (n = 23) and of other breeds (n = 27) born in Culicoides infested area; (ii) to investigate if barley could be a useful production system of allergens for IBH immunoassays. IgE binding in sera was tested by ELISA on two recombinant Culicoides allergens, rCul n 3 and rCul n 4, each produced in E. coli, insect cells and barley. Significantly more IgE was detected against all allergens in sera from IBH-affected compared to healthy horses. Icelandic-born Icelandic horses stand out with higher IgE levels against the allergens and higher area under the curve (AUC) on rCul n 4 as compared to the European-born horses. The barley and E.coli produced allergens had very similar performance in distinguishing between IBH-affected and healthy horses.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs.

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

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems.

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

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 developmental transcriptome of contrasting Arctic charr ( Salvelinus alpinus) morphs.

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.   To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr. The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr. The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

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.

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.