Phylogeny of salmonine fishes based on growth hormone introns: Atlantic (Salmo) and Pacific (Oncorhynchus) salmon are not sister taxa

Molecular and functional characterization of two IgM subclasses in large yellow croaker (Larimichthys crocea)

As the predominant immunoglobulin (Ig) isotype, IgM plays a crucial role in the acquired immunity of vertebrates. There is only one Igμ gene in mammals, except cattle, while the number of Igμ gene varies among teleost fish. In the current study, we found two functional Igμ genes (Igμ1 and Igμ2) and a pseudo Cμ gene (ψIgμ) in large yellow croaker (Larimichthys crocea). Both Igμ1 and Igμ2 genes possessed two transcript variants, which encoded the heavy chains of secreted (sIgM1 and sIgM2) and membrane-bound IgM1 and IgM2 (mIgM1 and mIgM2), respectively. Both the heavy chains of sIgM1 and sIgM2 consisted of a variable Ig domain, four constant Ig domains (CH1, CH2, CH3 and CH4) and a secretory tail, while those of mIgM1 and mIgM2 consisted of a variable Ig domain, three constant Ig domains (CH1, CH2 and CH3), a transmembrane domain and a short cytoplasmic tail. Cysteine residues that are necessary for the formation of intrachain and interchain disulfide bonds and tryptophan residues that are important for the folding of the Ig superfamily domain were well conserved in large yellow croaker IgM1 and IgM2. Interestingly, large yellow croaker IgM2 had an extra cysteine (C94) in the CH1 domain compared with IgM1, which may cause the structural difference between IgM1 and IgM2. A liquid chromatography-tandem mass spectrometry analysis revealed that both IgM1 and IgM2 were present at the protein level in large yellow croaker serum. Both the Igμ1 and Igμ2 genes were mainly expressed in systemic immune tissues, such as head kidney and spleen, but the expression level of Igμ2 was much lower than that of Igμ1. After Pseudomonas plecoglossicida infection, the expression levels of Igμ1 and Igμ2 in both the spleen and head kidney were significantly upregulated, with a higher upregulation of Igμ2 than that of Igμ1. These results suggested that Igμ1 and Igμ2 may play a differential role in the immune response of large yellow croaker against bacterial infection.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Molecular Phylogeny and Adaptive Mitochondrial DNA Evolution of Salmonids (Pisces: Salmonidae)

Salmonids are composed of anadromous and freshwater fishes, which is an important model for studying adaptive evolution. Herein, 49 salmonid complete mitochondrial genomes and those of two outgroups were used to infer a robust phylogeny for the family Salmonidae. The BI and RAxML phylogenetic trees based on 13 concatenated mitochondrial protein-coding genes showed well-supported nodes, and topologies were highly congruent. The concatenated 13 mitochondrial protein-coding genes, ND2, ND3, and ND5 genes were shown to have significantly larger dN/dS ratios in anadromous species than in freshwater species of Salmonidae, but the CYTB gene had significantly smaller dN/dS in anadromous species. The FEL analysis identified positively selected sites and negatively selected sites in each mitochondrial protein-coding gene separately. The RELAX program revealed that the ATP8 and CYTB genes supported intensified selection of the anadromous lineages. Our results demonstrated the phylogeny of Salmonidae and explored the mitochondrial DNA evolution pattern between anadromous and freshwater salmonids.

Immunoglobulins in teleosts

Immunoglobulins are glycoproteins which are produced as membrane-bound receptors on B-cells or in a secreted form, known as antibodies. In teleosts, three immunoglobulin isotypes, IgM, IgT, and IgD, are present, each comprising two identical heavy and two identical light polypeptide chains. The basic mechanisms for generation of immunoglobulin diversity are similar in teleosts and higher vertebrates. The B-cell pre-immune repertoire is diversified by VDJ recombination, junctional flexibility, addition of nucleotides, and combinatorial association of light and heavy chains, while the post-immune repertoire undergoes somatic hypermutation during clonal expansion. Typically, the teleost immunoglobulin heavy chain gene complex has a modified translocon arrangement where the Dτ-Jτ-Cτ cluster of IgT is generally located between the variable heavy chain (VH) region and the Dμ/δ-Jμ/δ-Cμ-Cδ gene segments, or within the set of VH gene segments. However, multiple genome duplication and deletion events and loss of some individual genes through evolution has complicated the IgH gene organization. The IgH gene arrangement allows the expression of either IgT or IgM/IgD. Alternative splicing is responsible for the regulation of IgM/IgD expression and the secreted versus transmembrane forms of IgT, IgD, and IgM. The overall structure of IgM and IgT is usually conserved across species, whereas IgD has a large variety of structures. IgM is the main effector molecule in both systemic and mucosal immunity and shows a broad range of concentrations in different teleost species. Although IgM is usually present in higher concentrations under normal conditions, IgT is considered the main mucosal Ig.

Polymorphisms of the growth hormone gene and their association with growth traits and sex in Sarcocheilichthys sinensis

The growth hormone gene (gh) of Sarcocheilichthys sinensis was cloned and characterized in this study. The cDNA length of gh was 973 bp, containing a 5'-UTR of 15 bp, a 3'-UTR of 325 bp and an open reading frame of 633 bp. The genomic DNA of gh was 2135 bp in length containing five exons and four introns. The precursor peptide of gh contained 210 amino acids (aa), including a signal peptide of 22 aa (Met¹-Ala²²) and a mature region of 188 aa (Ser²³-Leu²¹⁰). The similarity and identity ranges of the gh precursor peptide with those of other cyprinids were 88.6%-99.0% and 84.8%-98.6%, respectively. The gh of S. sinensis expressed at the highest level in the pituitary, and its expression was also detected in muscle and brain. Six polymorphic sites were detected in intron 1 (g.51InDel, g.64InDel and g.242InDel), intron 2 (g.864T>C), intron 3 (g.1017InDel) and intron 4 (g.1541A>G). Among these sites, g.242InDel was significantly associated with condition factor, g.1541A>G was associated with all six growth traits, while g.864T>C was associated with sex. The data obtained herein provide useful information for further studies on the regulation mechanisms of growth and sexual growth differences in S. sinensis.

Addressing incomplete lineage sorting and paralogy in the inference of uncertain salmonid phylogenetic relationships

Recent and continued progress in the scale and sophistication of phylogenetic research has yielded substantial advances in knowledge of the tree of life; however, segments of that tree remain unresolved and continue to produce contradicting or unstable results. These poorly resolved relationships may be the product of methodological shortcomings or of an evolutionary history that did not generate the signal traits needed for its eventual reconstruction. Relationships within the euteleost fish family Salmonidae have proven challenging to resolve in molecular phylogenetics studies in part due to ancestral autopolyploidy contributing to conflicting gene trees. We examine a sequence capture dataset from salmonids and use alternative strategies to accommodate the effects of gene tree conflict based on aspects of salmonid genome history and the multispecies coalescent. We investigate in detail three uncertain relationships: (1) subfamily branching, (2) monophyly of Coregonus and (3) placement of Parahucho. Coregoninae and Thymallinae are resolved as sister taxa, although conflicting topologies are found across analytical strategies. We find inconsistent and generally low support for the monophyly of Coregonus, including in results of analyses with the most extensive dataset and complex model. The most consistent placement of Parahucho is as sister lineage of Salmo.

Choosing and Using Introns in Molecular Phylogenetics

Introns are now commonly used in molecular phylogenetics in an attempt to recover gene trees that are concordant with species trees, but there are a range of genomic, logistical and analytical considerations that are infrequently discussed in empirical studies that utilize intron data. This review outlines expedient approaches for locus selection, overcoming paralogy problems, recombination detection methods and the identification and incorporation of LVHs in molecular systematics. A range of parsimony and Bayesian analytical approaches are also described in order to highlight the methods that can currently be employed to align sequences and treat indels in subsequent analyses. By covering the main points associated with the generation and analysis of intron data, this review aims to provide a comprehensive introduction to using introns (or any non-coding nuclear data partition) in contemporary phylogenetics.

Patterns of variation and covariation in the shapes of mandibular bones of juvenile salmonids in the genus Oncorhynchus

What is the nature of evolutionary divergence of the jaw skeleton within the genus Oncorhynchus? How can two associated bones evolve new shapes and still maintain functional integration? Here, we introduce and test a "concordance" hypothesis, in which an extraordinary matching of the evolutionary shape changes of the dentary and angular articular serves to preserve their fitting together. To test this hypothesis, we examined morphologies of the dentary and angular articular at parr (juvenile) stage, and at three levels of biological organization—between salmon and trout, between sister species within both salmon and trout, and among three types differing in life histories within one species, Oncorhynchus mykiss. The comparisons show bone shape divergences among the groups at each level; morphological divergence between salmon and trout is marked even at this relatively early life history stage. We observed substantial matching between the two mandibular bones in both pattern and amount of shape variation, and in shape covariation across species. These findings strongly support the concordance hypothesis, and reflect functional and/or developmental constraint on morphological evolution. We present evidence for developmental modularity within both bones. The locations of module boundaries were predicted from the patterns of evolutionary divergences, and for the dentary, at least, would appear to facilitate its functional association with the angular articular. The modularity results suggest that development has biased the course of evolution.

Comparative Analysis of the Shared Sex-Determination Region (SDR) among Salmonid Fishes

Salmonids present an excellent model for studying evolution of young sex-chromosomes. Within the genus, Oncorhynchus, at least six independent sex-chromosome pairs have evolved, many unique to individual species. This variation results from the movement of the sex-determining gene, sdY, throughout the salmonid genome. While sdY is known to define sexual differentiation in salmonids, the mechanism of its movement throughout the genome has remained elusive due to high frequencies of repetitive elements, rDNA sequences, and transposons surrounding the sex-determining regions (SDR). Despite these difficulties, bacterial artificial chromosome (BAC) library clones from both rainbow trout and Atlantic salmon containing the sdY region have been reported. Here, we report the sequences for these BACs as well as the extended sequence for the known SDR in Chinook gained through genome walking methods. Comparative analysis allowed us to study the overlapping SDRs from three unique salmonid Y chromosomes to define the specific content, size, and variation present between the species. We found approximately 4.1 kb of orthologous sequence common to all three species, which contains the genetic content necessary for masculinization. The regions contain transposable elements that may be responsible for the translocations of the SDR throughout salmonid genomes and we examine potential mechanistic roles of each one.

Recurrent evolution of life history ecotypes in sockeye salmon: implications for conservation and future evolution

We examine the evolutionary history and speculate about the evolutionary future of three basic life history ecotypes that contribute to the biocomplexity of sockeye salmon (Oncorhynchus nerka). The ‘recurrent evolution’ (RE) hypothesis claims that the sea/river ecotype is ancestral, a ‘straying’ form with poorly differentiated (meta)population structure, and that highly structured populations of lake-type sockeye and kokanee have evolved repeatedly in parallel adaptive radiations between recurrent glaciations of the Pleistocene Epoch. Basic premises of this hypothesis are consistent with new, independent evidence from recent surveys of genetic variation in mitochondrial and microsatellite DNA: (1) sockeye salmon are most closely related to pink (O. gorbuscha) and chum (O. keta) salmon with sea-type life histories; (2) the sockeye life history ecotypes exist as polyphyletic lineages within large drainages and geographic regions; (3) the sea/river ecotype exhibits less genetic differentiation among populations than the lake or kokanee ecotypes both within and among drainages; and (4) genetic diversity is typically higher in the sea/river ecotype than in the lake and kokanee ecotypes. Anthropogenic modification of estuarine habitat and intensive coastal fisheries have likely reduced and fragmented historic metapopulations of the sea/river ecotype, particularly in southern areas. In contrast, the kokanee ecotype appears to be favoured by marine fisheries and predicted changes in climate.

Phylogeny of suckermouth catfishes (Mochokidae: Chiloglanis) from Kenya: the utility of Growth Hormone introns in species level phylogenies

African suckermouth catfishes (Mochokidae: Chiloglanis) occur in freshwater throughout tropical Africa. Specimens from all major drainages across Kenya were collected over three field seasons. Here we present a phylogeny inferred from both mitochondrial cytochrome b (cyt b) and introns of the nuclear Growth Hormone gene (GH). The phylogeny inferred from introns is largely congruent with the results from an analysis of cyt b. The length and variability of GH introns make them ideal species level nuclear markers without the problem of introgression commonly encountered with mitochondrial genes. This analysis confirmed the presence of two previously known undescribed Chiloglanis species and also suggests the presence of previously unknown diversity within the Athi River system. The resulting phylogeny also indicates the presence of two separate lineages within C. brevibarbis. The historical biogeography of Chiloglanis within Kenya is discussed. The utility of GH intron for species level phylogenies of Siluriformes is compared to that in other groups.

Complete mitochondrial genome of Sakhalin taimen Parahucho perryi (Salmoniformes, Salmonidae) without two frame-disrupting indels in the ND4 gene

The complete mitochondrial genome was sequenced in two individuals of Sakhalin taimen Parahucho perryi. The genome sizes were 16,652 bp and 16,653 bp in the two isolates. Gene arrangement, base composition, and size of the two sequenced haplotypes are very similar to the P. perryi genome previously published (KC897653), but do not contain two frame-disrupting indels in the ND4 gene.

Molecular phylogeny of Salmo of the western Balkans, based upon multiple nuclear loci

Background

Classification of species within the genus Salmo is still a matter of discussion due to their high level of diversity and to the low power of resolution of mitochondrial (mt)DNA-based phylogeny analyses that have been traditionally used in evolutionary studies of the genus. We apply a new marker system based on nuclear (n)DNA loci to present a novel view of the phylogeny of Salmo representatives and we compare it with the mtDNA-based phylogeny.

Methods

Twenty-two nDNA loci were sequenced for 76 individuals of the brown trout complex: Salmo trutta (Danubian, Atlantic, Adriatic, Mediterranean and Duero mtDNA lineages), Salmo marmoratus (marble trout), Salmo obtusirostris (softmouth trout), and Salmo ohridanus (Ohrid belvica or belushka). Sequences were phylogenetically analyzed using maximum-likelihood and Bayesian Inference methods. The divergence time of the major clades was estimated using the program BEAST.

Results

The existence of five genetic units i.e. S. salar, S. ohridanus, S. obtusirostris, S. marmoratus and the S. trutta complex, including its major phylogenetic lineages was confirmed. Contrary to previous observations, S. obtusirostris was found to be sister to the S. trutta complex and the S. marmoratus clade rather than to the S. ohridanus clade. Reticulate evolution of S. obtusirostris was confirmed and a time for its pre-glacial origin suggested. S. marmoratus was found to be a separate species as S. trutta and S. obtusirostris. Relationships among lineages within the S. trutta complex were weakly supported and remain largely unresolved.

Conclusions

Nuclear DNA-based results showed a fairly good match with the phylogeny of Salmo inferred from mtDNA analyses. The comparison of nDNA and mtDNA data revealed at least four cases of mitochondrial–nuclear DNA discordance observed that were all confined to the Adriatic basin of the Western Balkans. Together with the well-known extensive morphological and genetic variability of Balkan trouts, this observation highlights an interesting and variegated evolutionary history of Salmo in this area.

Genome duplication and multiple evolutionary origins of complex migratory behavior in Salmonidae

Multiple rounds of whole genome duplication have repeatedly marked the evolution of vertebrates, and correlate strongly with morphological innovation. However, less is known about the behavioral, physiological and ecological consequences of genome duplication, and whether these events coincide with major transitions in vertebrate complexity. The complex behavior of anadromy - where adult fishes migrate up rivers from the sea to their natal site to spawn - is well known in salmonid fishes. Some hypotheses suggest that migratory behavior evolved as a consequence of an ancestral genome duplication event, which permitted salinity tolerance and osmoregulatory plasticity. Here we test whether anadromy evolved multiple times within salmonids, and whether genome duplication coincided with the evolution of anadromy. We present a method that uses ancestral character simulation data to plot the frequency of character transitions over a time calibrated phylogenetic tree to provide estimates of the absolute timing of character state transitions. Furthermore, we incorporate extinct and extant taxa to improve on previous estimates of divergence times. We present the first phylogenetic evidence indicating that anadromy evolved at least twice from freshwater salmonid ancestors. Results suggest that genome duplication did not coincide in time with changes in migratory behavior, but preceded a transition to anadromy by 55-50 million years. Our study represents the first attempt to estimate the absolute timing of a complex behavioral trait in relation to a genome duplication event.

Framing the Salmonidae family phylogenetic portrait: a more complete picture from increased taxon sampling

Considerable research efforts have focused on elucidating the systematic relationships among salmonid fishes; an understanding of these patterns of relatedness will inform conservation- and fisheries-related issues, as well as provide a framework for investigating evolutionary mechanisms in the group. However, uncertainties persist in current Salmonidae phylogenies due to biological and methodological factors, and a comprehensive phylogeny including most representatives of the family could provide insight into the causes of these difficulties. Here we increase taxon sampling by including nearly all described salmonid species (n = 63) to present a time-calibrated and more complete portrait of Salmonidae using a combination of molecular markers and analytical techniques. This strategy improved resolution by increasing the signal-to-noise ratio and helped discriminate methodological and systematic errors from sources of difficulty associated with biological processes. Our results highlight novel aspects of salmonid evolution. First, we call into question the widely-accepted evolutionary relationships among sub-families and suggest that Thymallinae, rather than Coregoninae, is the sister group to the remainder of Salmonidae. Second, we find that some groups in Salmonidae are older than previously thought and that the mitochondrial rate of molecular divergence varies markedly among genes and clades. We estimate the age of the family to be 59.1 MY (CI: 63.2-58.1 MY) old, which likely corresponds to the timing of whole genome duplication in salmonids. The average, albeit highly variable, mitochondrial rate of molecular divergence was estimated as ~0.31%/MY (CI: 0.27-0.36%/MY). Finally, we suggest that some species require taxonomic revision, including two monotypic genera, Stenodus and Salvethymus. In addition, we resolve some relationships that have been notoriously difficult to discern and present a clearer picture of the evolution of the group. Our findings represent an important contribution to the systematics of Salmonidae, and provide a useful tool for addressing questions related to fundamental and applied evolutionary issues.

A different perspective on the phylogenetic relationships of the Moxostomatini (Cypriniformes: Catostomidae) based on cytochrome-b and Growth Hormone intron sequences

We have examined phylogenetic relationships of suckers of tribe Moxostomatini (Cypriniformes, Catostomidae) using cytochrome-b and Growth Hormone gene intron sequences. Phylogenies were significantly different from recent estimates of relationships based primarily on morphology (Smith, 1992) and cytochrome-b sequences (Harris et al., 2002). Overall, there was little support for many basal nodes in the phylogeny, however it was clear that Scartomyzon and Moxostoma were not monophyletic, despite morphological evidence provided Robins and Raney (1956, 1957), Jenkins (1970), and Smith (1992). Growth Hormone sequences provided good support for a monophyletic Western Scartomyzon lineage and thus suggested a single ancestral invasion of Scartomyzon-like fishes into drainages of Texas and Mexico. Phylogenetic relationships of Western Scartomyzon are structured geographically and do not conform well to current taxonomy.

Whole genome duplication: challenges and considerations associated with sequence orthology assignment in Salmoninae

To illustrate some of the challenges and considerations in assigning correct orthology necessary for any comparative genomic investigation among salmonids, sequence data from the non-coding regions of different chromosomes in three members of the subfamily Salmoninae, rainbow trout Oncorhynchus mykiss, Atlantic salmon Salmo salar and Arctic charr Salvelinus alpinus, were compared. By analysing c. 55 distinct loci, corresponding to c. 142 kbp sequence information per species, 18 duplicated patterns representative of the two sequential rounds of teleost-specific whole genome duplications (i.e. 3R and 4R WGD) were identified. Sequence similarities between the 4R paralogues were c. 90%, which was slightly lower than those of the 4R orthologues and c. 60% for the 3R products. Through careful examination of the sequence data, however, only 14 loci could reliably be assigned as true orthologues. Locus-specific trees were constructed through maximum parsimony, maximum likelihood and neighbour-joining methods and were rooted using the information from a close relative, lake whitefish Coregonus clupeaformis. All approaches generated congruent trees supporting the {Coregonus [Salmo (Oncorhynchus, Salvelinus)]} topology. The general phenotypic characteristics of sequences, however, were highly suggestive of the basal position of Oncorhynchus, raising the hypothesis of an accelerated rate of nucleotide evolution in this species.

Spectral and polarization sensitivity of juvenile Atlantic salmon (Salmo salar): phylogenetic considerations

We were interested in comparing the characteristics of polarization sensitivity in Atlantic salmon to those in Pacific salmon. Here we show that the common ancestor to the clade containing Salmo salar, Oncorhynchus mykiss, O. nerka, O. clarkii and Salvelinus fontinalis has the trait of ultraviolet polarization sensitivity. We examined spectral and polarization sensitivity of juvenile Atlantic salmon (Salmo salar) using both optic nerve compound action potential (CAP) and electroretinogram (ERG) recordings. Our experiments employed photic manipulation to adjust the sensitivity of the four cone mechanisms of Atlantic salmon. A spectrally broad background was used to ensure a contribution of all cone mechanisms to both spectral and polarization sensitivity. Chromatic adaptation was used to isolate the sensitivity of each of the four cone mechanisms for both spectral and polarization sensitivity. Under spectrally broad conditions, UV sensitive (UVS), mid wavelength sensitive (MWS) and long wavelength sensitive (LWS) cone mechanisms contributed to polarization sensitivity. CAP recordings produced the typical 'W' shaped polarization sensitivity curve reflecting two active polarization detectors with peaks at e-vector orientations of 0 deg, 90 deg and 180 deg, and troughs at 30 deg and 150 deg. ERG recordings produced a four-peaked polarization sensitivity curve reflecting two active polarization detectors and negative feedback activity, with peaks at e-vectors 0 deg, 45 deg, 90 deg, 135 deg and 180 deg, and troughs at 30 deg, 60 deg, 120 deg and 150 deg. Polarization-sensitivity measurements of isolated cone mechanisms revealed two orthogonal polarization detector mechanisms in Atlantic salmon, identical to that found in rainbow trout and other Pacific salmonid fishes. Moreover, under spectrally broad background conditions, CAP and ERG polarization sensitivity of Atlantic salmon did not differ significantly from that reported in Pacific salmonids.

Evolution of duplicated IgH loci in Atlantic salmon, Salmo salar

BACKGROUND

The Atlantic salmon (Salmo salar) immunoglobulin heavy chain (IgH) locus possesses two parallel IgH isoloci (IGH-A and IGH-B), that are related to the genomic duplication event in the family Salmonidae. These duplicated IgH loci in Atlantic salmon provide a unique opportunity to examine the mechanisms of genome diversity and genome evolution of the IgH loci in vertebrates. In this study, we defined the structure of these loci in Atlantic salmon, and sequenced 24 bacterial artificial chromosome (BAC) clones that were assembled into the IGH-A (1.1 Mb) and IGH-B (0.9 Mb) loci. In addition, over 7,000 cDNA clones from the IgH variable (VH) region have been sequenced and analyzed.

RESULTS

The present study shows that the genomic organization of the duplicated IgH loci in Atlantic salmon differs from that in other teleosts and other vertebrates. The loci possess multiple Cτ genes upstream of the Cμ region, with three of the Cτ genes being functional. Moreover, the duplicated loci possess over 300 VH segments which could be classified into 18 families. This is the largest number of VH families currently defined in any vertebrate. There were significant structural differences between the two loci, indicating that both IGH-A and -B loci have evolved independently in the short time after the recent genome duplication approximately 60 mya.

CONCLUSIONS

Our results indicate that the duplication of the IgH loci in Atlantic salmon significantly contributes to the increased diversity of the antibody repertoire, as compared with the single IgH locus in other vertebrates.

Multiplexed identification of different fish species by detection of parvalbumin, a common fish allergen gene: a DNA application of multi-analyte profiling (xMAP) technology

Fish are a common cause of allergic reactions associated with food consumption, with parvalbumin being the major allergenic protein. Some fish-hypersensitive patients tolerate some fish species while being allergic to others. Reliable detection methods for allergenic fish species in foods are necessary to ensure compliance with food allergen labeling guidelines to protect fish-allergic consumers. The objective of this project was to develop a multi-analyte detection method for the presence of fish in food. Therefore, conserved parvalbumin exon sequences were utilized for the design of universal PCR primers amplifying intron DNA and small regions of exons flanking the enclosed intron from even very distantly related fish species. An assay for the identification of eight fish species was developed using xMAP technology with probes targeting species-specific parvalbumin intron regions. Additionally, a universal fish probe was designed targeting a highly conserved exon region located between the intron and the reverse primer region. The universal fish assay showed no cross-reactivity with other species, such as beef, pork, lamb, chicken, turkey, and shrimp. Importantly, with the exception of one notable case with fish in the same subfamily, species-specific detection showed no cross-reactivity with other fish species. Limits of detection for these eight species were experimentally estimated to range from 0.01% to 0.04%, with potential to increase the detection sensitivity. This report introduces a newly developed method for the multiplex identification of at least eight allergenic fish species in food, which could conceivably be extended to detect up to 100 species simultaneously in one sample.

Evolutionary divergence of duplicate copies of the growth hormone gene in suckers (Actinopterygii: catostomidae)

Catostomid fishes (suckers) have duplicate copies of the growth hormone gene and other nuclear genes, due to a genome duplication event early in the group’s history. Yet, paralogs of GH in suckers are more than 90% conserved in nucleotide (nt) and amino acid (aa) sequence. Within paralogs across species, variation in nt and aa sequence averages 3.33% and 4.46% for GHI, and 3.22% and 2.43% for GHII, respectively. Selection tests suggest that the two GH paralogs are under strong purifying selection. Consensus trees from phylogenetic analysis of GH coding region data for 23 species of suckers, other cypriniform fishes and outgroups resolved cypriniform relationships and relationships among GHI sequences of suckers more or less consistently with analyses based on other molecular data. However, the analysis failed to resolve all sucker GHI and GHII sequences as monophyletic sister groups. This unexpected topology did not differ significantly from topologies constrained to make all GH sequences monophyletic. We attribute this result either to limitations in our GHII data set or convergent adaptive changes in GHII of tribe Catostomini.

Nibea coibor growth hormone gene: its phylogenetic significance, microsatellite variation and expression analysis

The growth hormone (GH) gene has been characterized for a number of fishes and used to establish phylogenetic relationships and as a candidate gene for studies of genetic variation in connection with growth traits. In this study, we report the genomic structure of Nibea coibor GH (designated as ncGH) including its 5'-flanking region, being cloned by homology-cloning and chromosome walking methods. The ncGH gene spans approximately 3.0 kb and consists of six exons and five introns, as found for all cloned teleost GH genes with the exception of carps and catfish. The 5'-flanking region contains consensus sequences for a TATA box, a CRE, a pit-1alpha, a TRE, two HNF-3, a ERE and a GRE. Five microsatellites are identified in the ncGH gene and three of them are polymorphic marker. The open reading frame (ORF) of ncGH is 615 bp in length encoding a polypeptide of 204 amino acids with an estimated molecular mass of 23.04 kDa and theoretical isoelectric point of 7.79. The precursor of ncGH consists of a 17 amino-acid signal peptide and a 187 amino-acid mature peptide. The four Cys residues are located at conserved positions (Cys(69), Cys(177), Cys(194) and Cys(202)), and One possible site for N-glycosylation (Asn-X-Ser/Thr motif) is present at Asn(201). The coding region sequence of ncGH is used to align with the sequences of 18 other species from Percoidei and one species from Anabantoidei using Clustal X. A matrix of 612 bp was used to construct the phylogenetic trees using neighbor-joining and maximum parsimony methods. The phylogenetic trees by two methods are identical in most of the clades with high bootstrap support. Every family all forms independent monophyly on the phylogenetic trees, in the family, the different species also forms the monophyly according to the different genera. The results are also identical to those from morphological data, and demonstrated that the GH gene is very suitable for phylogenetic relationship analysis of Percoidei. To validate the predicted exon/intron boundaries, ncGH cDNA is cloned using RT-PCR, and tissue distributions are investigated using semi-quantitative RT-PCR method. The results indicate that the predicted exon/intron is correct, the ncGH mRNA are mainly expressed in pituitary, and weakly expressed in ovary, brain, liver, gill, intestine, muscle and hear, but not expressed in spleen and kidney.

A salmonid EST genomic study: genes, duplications, phylogeny and microarrays

Background

Salmonids are of interest because of their relatively recent genome duplication, and their extensive use in wild fisheries and aquaculture. A comprehensive gene list and a comparison of genes in some of the different species provide valuable genomic information for one of the most widely studied groups of fish.

Results

298,304 expressed sequence tags (ESTs) from Atlantic salmon (69% of the total), 11,664 chinook, 10,813 sockeye, 10,051 brook trout, 10,975 grayling, 8,630 lake whitefish, and 3,624 northern pike ESTs were obtained in this study and have been deposited into the public databases. Contigs were built and putative full-length Atlantic salmon clones have been identified. A database containing ESTs, assemblies, consensus sequences, open reading frames, gene predictions and putative annotation is available. The overall similarity between Atlantic salmon ESTs and those of rainbow trout, chinook, sockeye, brook trout, grayling, lake whitefish, northern pike and rainbow smelt is 93.4, 94.2, 94.6, 94.4, 92.5, 91.7, 89.6, and 86.2% respectively. An analysis of 78 transcript sets show Salmo as a sister group to Oncorhynchus and Salvelinus within Salmoninae, and Thymallinae as a sister group to Salmoninae and Coregoninae within Salmonidae. Extensive gene duplication is consistent with a genome duplication in the common ancestor of salmonids. Using all of the available EST data, a new expanded salmonid cDNA microarray of 32,000 features was created. Cross-species hybridizations to this cDNA microarray indicate that this resource will be useful for studies of all 68 salmonid species.

Conclusion

An extensive collection and analysis of salmonid RNA putative transcripts indicate that Pacific salmon, Atlantic salmon and charr are 94–96% similar while the more distant whitefish, grayling, pike and smelt are 93, 92, 89 and 86% similar to salmon. The salmonid transcriptome reveals a complex history of gene duplication that is consistent with an ancestral salmonid genome duplication hypothesis. Genome resources, including a new 32 K microarray, provide valuable new tools to study salmonids.

The genetic basis of smoltification-related traits in Oncorhynchus mykiss

The timing and propensity for migration between fresh- and seawater is a key theme in the diversity of life histories within the salmonid fishes. Across salmonid species, life-history strategies range from wholly freshwater-resident populations, to migratory and nonmigratory variation within populations, to populations and species that are primarily migratory. Despite the central theme of migration to the evolution of these fishes, the genetic architecture of migration-related processes is poorly understood. Using a genetic cross of clonal lines derived from migratory and nonmigratory life-history types of Onchorhynchus mykiss (steelhead and rainbow trout, respectively), we have dissected the genetic architecture of the complex physiological and morphological transformation that occurs immediately prior to seaward migration (termed smoltification). Quantitative trait loci (QTL) analyses were used to identify the number, effects, and genomic location of loci associated with smoltification-related traits, including growth and condition factor, body coloration, morphology, and osmoregulatory enzymes during the smoltification period. Genetic analyses revealed numerous QTL, but one locus in particular is associated with multiple traits in single and joint analyses. Dissecting the genetic architecture of this highly complex trait has profound implications for understanding the genetic and evolutionary basis of life-history diversity within and among migratory fishes.

Retroposons of salmonoid fishes (Actinopterygii: Salmonoidei) and their evolution

Short and long retroposons, or non-LTR retrotransposons (SINEs and LINEs, respectively) are two groups of interspersed repetitive elements amplifying in the genome via RNA and cDNA-mediated reverse transcription. In this process, SINEs entirely depend on the enzymatic machinery of autonomous LINEs. The impact of retroposons on the host genome is difficult to overestimate: their sequences account for significant portion of the eukaryotic genome, while propagation of their active copies gradually reshapes it. In this way, the retropositional activity plays a role of important evolutionary factor. More than 100 LINE and nearly 100 SINE families have been described to date from the genomes of various eukaryotes, and it is salmonoid fishes (Actinopterygii: Salmonoidei) that are particularly noticeable for the diversity of transposons they host in their genomes, including two LINE and seven SINE families. Moreover, this group of ray-finned fish represents an excellent opportunity to study such a rare evolutionary phenomenon as lateral gene transfer, due to a great variety of transposons and other sequences salmons share with a blood fluke, Schistosoma japonicum (Trematoda: Strigeiformes)--a parasitic helminth infecting various vertebrates. The aim of the present review is to structure all knowledge accumulated about salmonoid retroposons by now, as well as to complement it with the new data pertaining to the distribution of some SINE families.

Patterns of Historical Balancing Selection on the Salmonid Major Histocompatibility Complex Class II β Gene

Allelic variation in the major histocompatibility class (MHC) IIB gene of salmonids is analyzed for patterns indicative of natural selection acting at the molecular level. Sequence data for the second exon of this MHC gene were generated for 11 species in three salmonid genera: Oncorhynchus, Salmo, and Salvelinus. Phylogenetic analysis of nucleotide sequences revealed: (1) monophyletic grouping of alleles from each genus, (2) transspecies evolution of alleles within Salmo and Salvelinus, and (3) differential patterns of transspecies evolution within the genus Oncorhynchus. Within Oncorhynchus, five of seven species had alleles that were species-specific or nearly so, while the remaining two, O. mykiss and O. clarkii, retained ancestral polymorphisms. The different patterns in Oncorhynchus and the other two genera could be due to historical demographic effects or functional differences in MHC molecules in the three genera, but the two hypotheses could not be distinguished with the current dataset. An analysis of recombination/gene conversion identified numerous recombinant alleles, which is consistent with what has been found in other vertebrate taxa. However, these gene conversion events could not account for the species-specific allelic lineages observed in five of the Oncorhynchus species. Analyses of the relative rates of nonsynonymous and synonymous substitutions revealed the signature of selection on the class IIB gene in all 11 of the salmonid species for both the ABS and the non-ABS codons. Codon-based analyses of selection identified seven codons that have experienced selection in the majority of the species. More than half of these sites were mammalian ABS codons, but several were not, suggesting subtle functional differences in the mammalian and teleost fish MHC molecules.

Novel SINE families from salmons validate Parahucho (Salmonidae) as a distinct genus and give evidence that SINEs can incorporate LINE-related 3'-tails of other SINEs

Short interspersed elements (SINEs) constitute a group of retroposons propagating in the genome via a mechanism of reverse transcription, in which they depend on the enzymatic machinery of long retroposons (LINEs). Over 70 SINE families have been described to date from the genomes of various eukaryotes. Here, we characterize two novel SINEs from salmons (Actinopterygii: Salmonoidei). The first family, termed SlmI, was shown to be widespread among all genera of the suborder. These SINEs have a tRNA(Leu)-related promoter region at their 5'-end, a unique central conserved domain with a subfamily-specific region, and an end with RSg-1-LINE-derived 3'-terminus preceding the A/T-rich tail. The same LINE-related segment is also shared by two other salmonid SINEs: HpaI and OS-SINE1. The structural peculiarities and overall sequence identity of the SlmI 3'-terminus suggest that it has been acquired from HpaI SINEs but not directly from the partner LINE. This region plays a crucial role in the process of retrotransposition of short interspersed elements, and the case of its SINE-to-SINE transmission is the first recorded to date. Possible scenarios and potential evolutionary implications of the observed interaction between short retroposons are discussed. Apart from the above, we found a copy of the SlmI SINE in the GenBank entry for the blood fluke, Schistosoma japonicum (Trematoda: Strigeiformes) -- a trematode causing one of the most important human helminth infections, with its genome known to host other groups of salmonoid retroposons. In the present article, we suggest our views with regard to possible ways in which such an intensive horizontal transfer of salmonoid retroposons to the schistosomal genome occurs. The second novel SINE family, termed SlmII, originates from one of the SlmI subfamilies, with which it shares the same tRNA-related region, central domain, and a part of RSg-1-derived segment, but has a different 3'-tail of unidentified origin. Its distribution among salmonids validates Parahucho (Japanese huchen) as a distinct monotypic genus.

Quantitative trait loci x maternal cytoplasmic environment interaction for development rate in Oncorhynchus mykiss

Effects of maternal cytoplasmic environment (MCE) on development rate in rainbow trout were evaluated within a quantitative trait loci (QTL) analysis framework. Previous research had identified QTL for development rate in doubled haploid (DH) progeny produced from a cross between the Oregon State University (OSU) and the Swanson (SW) River rainbow trout clonal lines. In this study, progeny for QTL mapping were produced from a cross between the OSU and Clearwater (CW) River clonal lines. Doubled haploids were produced from the OSU x CW F1 by androgenesis using eggs from different females (or MCEs); with androgenesis, the maternal nuclear genome was destroyed by irradiation and diploidy was restored by blocking the first embryonic cleavage by heat shock. All embryos were incubated at the same temperature and development rate quantified as time to hatch. Using a linkage map constructed primarily with AFLP markers, QTL mapping was performed, including MCE covariates and QTL x MCE effects in models for testing. The major QTL for development rate in the OSU x SW cross overlaps with the major QTL found in this OSU x CW cross; effects at this locus were the same across MCEs. Both MCE and QTL x MCE effects contribute to variability in development rate, but QTL x MCE were minor and detected only at small-effect QTL.

Isolation and characterization of salmonid telomeric and centromeric satellite DNA sequences

Satellite DNA clones with a 37 bp repeat unit were obtained from BglII-digested genomic DNA of Masu salmon (Oncorhynchus masou) and Chum salmon (O. keta). Fluorescence in situ hybridization (FISH) analysis with the isolated clones as a probe showed that these repetitive sequences were localized in the telomeric regions of chromosomes in both species. Southern and dot blot analyses suggested conservation of homologous sequences with similar repeat unit in other salmonids including the species of the genus Oncorhynchus and Salvelinus, but lack or scarcity of such sequences in the genus Hucho and Salmo. Similarly, polymerase chain reaction (PCR)-based cloning of satellite DNA referring to a reported Rainbow trout (O. mykiss) centromeric sequence was successful for the Oncorhynchus, Salvelinus and Hucho species. The obtained satellite DNA clones were localized with FISH in the centromeric regions of chromosomes of the species from these three genera. Although PCR cloning of the centromeric satellite DNA had failed in the Salmo species due to some base changes in the priming sites, dot blot hybridization analysis suggested conservation of homologous satellite DNA in the genus Salmo as in the other three genera. In the neighbor-joining tree of cloned centromeric satellite DNA sequences, the genus Oncorhynchus and Salvelinus formed adjacent clades, and the clade of the genus Hucho included the reported centromeric sequence of the genus Salmo. Conservation pattern and molecular phylogeny of the telomeric and centromeric satellite DNA sequences isolated herein support a close phylogenetic relationship between the genus Oncorhynchus and Salvelinus and between the Salmo and Hucho.

Single nucleotide polymorphism (SNP) discovery in duplicated genomes: intron-primed exon-crossing (IPEC) as a strategy for avoiding amplification of duplicated loci in Atlantic salmon (Salmo salar) and other salmonid fishes

Background

Single nucleotide polymorphisms (SNPs) represent the most abundant type of DNA variation in the vertebrate genome, and their applications as genetic markers in numerous studies of molecular ecology and conservation of natural populations are emerging. Recent large-scale sequencing projects in several fish species have provided a vast amount of data in public databases, which can be utilized in novel SNP discovery in salmonids. However, the suggested duplicated nature of the salmonid genome may hamper SNP characterization if the primers designed in conserved gene regions amplify multiple loci.

Results

Here we introduce a new intron-primed exon-crossing (IPEC) method in an attempt to overcome this duplication problem, and also evaluate different priming methods for SNP discovery in Atlantic salmon (Salmo salar) and other salmonids. A total of 69 loci with differing priming strategies were screened in S. salar, and 27 of these produced ~13 kb of high-quality sequence data consisting of 19 SNPs or indels (one per 680 bp). The SNP frequency and the overall nucleotide diversity (3.99 × 10^-4) in S. salar was lower than reported in a majority of other organisms, which may suggest a relative young population history for Atlantic salmon. A subset of primers used in cross-species analyses revealed considerable variation in the SNP frequencies and nucleotide diversities in other salmonids.

Conclusion

Sequencing success was significantly higher with the new IPEC primers; thus the total number of loci to screen in order to identify one potential polymorphic site was six times less with this new strategy. Given that duplication may hamper SNP discovery in some species, the IPEC method reported here is an alternative way of identifying novel polymorphisms in such cases.