Developing molecular classifiers to detect environmental stressors, smolt stages and morbidity in coho salmon, Oncorhynchus kisutch.

Aquatic species are increasingly confronted with environmental stressors because of climate change. Although molecular technologies have advanced our understanding of how organisms respond to stressors in laboratory settings, the ability to detect physiological responses to specific stressors under complex field conditions remains underdeveloped. This research applied multi-stressor challenge trials on coho salmon, employing the "Salmon Fit-Chips" genomic tool and a random forest-based classification model to develop classifiers predictive for chronic thermal and hypoxic stress, as well as salinity acclimation, smolt stage and morbidity status. The study also examined how smolts and de-smolts (smolts not having entered SW during the smolt window) responded transcriptionally to exposure to saltwater. Using RF classifiers optimized with 4 to 12 biomarkers, we identified transcriptional signatures that accurately predicted the presence of each stressor and physiological state, achieving prediction accuracy rates between 86.8 % and 100 %, regardless of other background stressors present. Stressor recovery time was established by placing fish back into non-stressor conditions after stress exposure, providing important context to stressor detections in field applications. Recovery from thermal and hypoxic stress requires about 3 and 2 days, respectively, with >3 days needed for re-acclimation to freshwater for seawater acclimated fish. The study also found non-additive (synergistic) effects of multiple stressors on mortality risk. Importantly, osmotic stress associated with de-smolts was the most important predictor of mortality. In saltwater, de-smolts exposed to salinity, high temperature, and hypoxia experienced a 9-fold increase in mortality compared to those only exposed to saltwater, suggesting a synergistic response to multiple stressors. These findings suggest that delays in hatchery releases to support release of larger fish need to be carefully scrutinized to ensure fish are not being released as de-smolts, which are highly susceptible to additional climate-induced stressors like rising temperatures and reduced dissolved oxygen levels in the marine environment.

The core microbiome of cultured Pacific oyster spat is affected by age but not mortality.

The Pacific oyster is the most widely cultured shellfish worldwide, but production has been affected by mortality events, including in hatcheries that supply the seed for growers. Several pathogens cause disease in oysters, but in many cases, mortality events cannot be attributed to a single agent and appear to be multifactorial, involving environmental variables and microbial interactions. As an organism's microbiome can provide resilience against pathogens and environmental stressors, we investigated the microbiomes in cohorts of freshly settled oyster spat, some of which experienced notable mortality. Deep sequencing of 16S rRNA gene fragments did not show a significant difference among the microbiomes of cohorts experiencing different mortality levels, but revealed a characteristic core microbiome comprising 74 taxa. Irrespective of mortality, the relative abundance of taxa in the core microbiomes changed significantly as the spat aged, yet remained distinct from the microbial community in the surrounding water. The core microbiome was dominated by bacteria in the families Rhodobacteraceae, Nitrosomonadaceae, Flavobacteriaceae, Pirellulaeceae, and Saprospiraceae. Within these families, 14 taxa designated as the "Hard-Core Microbiome" were indicative of changes in the core microbiome as the spat aged. The variability in diversity and richness of the core taxa decreased with age, implying niche occupation. As well, there was exchange of microbes with surrounding water during development of the core microbiome. The shift in the core microbiome demonstrates the dynamic nature of the microbiome as oyster spat age.IMPORTANCEThe Pacific oyster (Magallana gigas, also known as Crassostrea gigas) is the most widely cultivated shellfish and is important to the economy of many coastal communities. However, high mortality of spat during the first few days following metamorphosis can affect the seed supply to oyster growers. Here, we show that the microbiome composition of recently settled oyster spat experiencing low or high mortality was not significantly different. Instead, development of the core microbiome was associated with spat aging and was partially driven by dispersal through the water. These findings imply the importance of early-stage rearing conditions for spat microbiome development in aquaculture facilities. Furthermore, shellfish growers could gain information about the developmental state of the oyster spat microbiome by assessing key taxa. Additionally, the study provides a baseline microbiome for future hypothesis testing and potential probiotic applications on developing spat.

Chinook and Coho salmon hybrids linked to habitat and climatic changes on Vancouver Island, British Columbia.

Between 2013 and 2019, 63 presumed Chinook salmon Oncorhynchus tshawytscha sampled primarily in the Strait of Georgia (0.63% of total sample) were identified as potential Chinook-Coho (Oncorhynchus kisutch) hybrids by the presence of anomalous microsatellite genotypes. Their hybrid origin was confirmed by single nucleotide polymorphism amplification of two species-specific amplicons. Mitochondrial DNA indicated that most of these fish resulted from the hybridization of Coho salmon females and Chinook salmon males. Although no diagnostic external features were identified, several individuals displayed an abnormal scale arrangement on the caudal peduncle. One hybrid juvenile examined for meristics exhibited a pyloric caeca count intermediate between published values for Chinook and Coho salmon. Most hybrids originated in the Cowichan River during the 2014 brood year. Their prevalence in the watershed is a naturally occurring event, likely exacerbated by prolonged low water levels which limit habitat and delay Chinook salmon spawning, in addition to the differential abundance of the parental species. This research is the first to document ongoing natural hybridization (Chinook-Coho salmon crosses) and link it to habitat and climatic changes, and includes the identification of eight F1 adults and two juvenile backcross or F2 hybrids. The potential negative impacts of hybridization, particularly in Coho salmon through potential introgression, warrant hybrid identification as an ecosystem monitoring tool within a survey program.

Aquaculture mediates global transmission of a viral pathogen to wild salmon.

Global expansion of aquaculture and agriculture facilitates disease emergence and catalyzes transmission to sympatric wildlife populations. The health of wild salmon stocks critically concerns Indigenous peoples, commercial and recreational fishers, and the general public. Despite potential impact of viral pathogens such as Piscine orthoreovirus-1 (PRV-1) on endangered wild salmon populations, their epidemiology in wild fish populations remains obscure, as does the role of aquaculture in global and local spread. Our phylogeographic analyses of PRV-1 suggest that development of Atlantic salmon aquaculture facilitated spread from Europe to the North and South East Pacific. Phylogenetic analysis and reverse transcription polymerase chain reaction surveillance further illuminate the circumstances of emergence of PRV-1 in the North East Pacific and provide strong evidence for Atlantic salmon aquaculture as a source of infection in wild Pacific salmon. PRV-1 is now an important infectious agent in critically endangered wild Pacific salmon populations, fueled by aquacultural transmission.

Descriptive multi-agent epidemiology via molecular screening on Atlantic salmon farms in the northeast Pacific Ocean.

Rapid expansion of salmon aquaculture has resulted in high-density populations that host diverse infectious agents, for which surveillance and monitoring are critical to disease management. Screening can reveal infection diversity from which disease arises, differential patterns of infection in live and dead fish that are difficult to collect in wild populations, and potential risks associated with agent transmission between wild and farmed hosts. We report results from a multi-year infectious-agent screening program of farmed salmon in British Columbia, Canada, using quantitative PCR to assess presence and load of 58 infective agents (viruses, bacteria, and eukaryotes) in 2931 Atlantic salmon (Salmo salar). Our analysis reveals temporal trends, agent correlations within hosts, and agent-associated mortality signatures. Multiple agents, most notably Tenacibaculum maritimum, were elevated in dead and dying salmon. We also report detections of agents only recently shown to infect farmed salmon in BC (Atlantic salmon calicivirus, Cutthroat trout virus-2), detection in freshwater hatcheries of two marine agents (Kudoa thyrsites and Tenacibaculum maritimum), and detection in the ocean of a freshwater agent (Flavobacterium psychrophilum). Our results provide information for farm managers, regulators, and conservationists, and enable further work to explore patterns of multi-agent infection and farm/wild transmission risk.

Discovery and surveillance of viruses from salmon in British Columbia using viral immune-response biomarkers, metatranscriptomics, and high-throughput RT-PCR.

The emergence of infectious agents poses a continual economic and environmental challenge to aquaculture production, yet the diversity, abundance, and epidemiology of aquatic viruses are poorly characterised. In this study, we applied salmon host transcriptional biomarkers to identify and select fish in a viral disease state, but only those that were negative for known viruses based on RT-PCR screening. These fish were selected for metatranscriptomic sequencing to discover potential viral pathogens of dead and dying farmed Atlantic (Salmo salar) and Chinook (Oncorhynchus tshawytscha) salmon in British Columbia (BC). We found that the application of the biomarker panel increased the probability of discovering viruses in aquaculture populations. We discovered two viruses that have not previously been characterised in Atlantic salmon farms in BC (Atlantic salmon calicivirus and Cutthroat trout virus-2), as well as partially sequenced three putative novel viruses. To determine the epidemiology of the newly discovered or emerging viruses, we conducted high-throughput reverse transcription polymerase chain reaction (RT-PCR) and screened over 9,000 farmed and wild salmon sampled over one decade. Atlantic salmon calicivirus and Cutthroat trout virus-2 were in more than half of the farmed Atlantic salmon we tested. Importantly we detected some of the viruses we first discovered in farmed Atlantic salmon in Chinook salmon, suggesting a broad host range. Finally, we applied in situ hybridisation to determine infection and found differing cell tropism for each virus tested. Our study demonstrates that continual discovery and surveillance of emerging viruses in these ecologically important salmon will be vital for management of both aquaculture and wild resources in the future.

Environmental DNA from multiple pathogens is elevated near active Atlantic salmon farms.

The spread of infection from reservoir host populations is a key mechanism for disease emergence and extinction risk and is a management concern for salmon aquaculture and fisheries. Using a quantitative environmental DNA methodology, we assessed pathogen environmental DNA in relation to salmon farms in coastal British Columbia, Canada, by testing for 39 species of salmon pathogens (viral, bacterial, and eukaryotic) in 134 marine environmental samples at 58 salmon farm sites (both active and inactive) over 3 years. Environmental DNA from 22 pathogen species was detected 496 times and species varied in their occurrence among years and sites, likely reflecting variation in environmental factors, other native host species, and strength of association with domesticated Atlantic salmon. Overall, we found that the probability of detecting pathogen environmental DNA (eDNA) was 2.72 (95% CI: 1.48, 5.02) times higher at active versus inactive salmon farm sites and 1.76 (95% CI: 1.28, 2.42) times higher per standard deviation increase in domesticated Atlantic salmon eDNA concentration at a site. If the distribution of pathogen eDNA accurately reflects the distribution of viable pathogens, our findings suggest that salmon farms serve as a potential reservoir for a number of infectious agents; thereby elevating the risk of exposure for wild salmon and other fish species that share the marine environment.

Endangered wild salmon infected by newly discovered viruses.

The collapse of iconic, keystone populations of sockeye (Oncorhynchus nerka) and Chinook (Oncorhynchus tshawytscha) salmon in the Northeast Pacific is of great concern. It is thought that infectious disease may contribute to declines, but little is known about viruses endemic to Pacific salmon. Metatranscriptomic sequencing and surveillance of dead and moribund cultured Chinook salmon revealed a novel arenavirus, reovirus and nidovirus. Sequencing revealed two different arenavirus variants which each infect wild Chinook and sockeye salmon. In situ hybridisation localised arenavirus mostly to blood cells. Population surveys of >6000 wild juvenile Chinook and sockeye salmon showed divergent distributions of viruses, implying different epidemiological processes. The discovery in dead and dying farmed salmon of previously unrecognised viruses that are also widely distributed in wild salmon, emphasizes the potential role that viral disease may play in the population dynamics of wild fish stocks, and the threat that these viruses may pose to aquaculture.

Distribution and Phylogeny of Erythrocytic Necrosis Virus (ENV) in Salmon Suggests Marine Origin.

Viral erythrocytic necrosis (VEN) affects over 20 species of marine and anadromous fishes in the North Atlantic and North Pacific Oceans. However, the distribution and strain variation of its viral causative agent, erythrocytic necrosis virus (ENV), has not been well characterized within Pacific salmon. Here, metatranscriptomic sequencing of Chinook salmon revealed that ENV infecting salmon was closely related to ENV from Pacific herring, with inferred amino-acid sequences from Chinook salmon being 99% identical to those reported for herring. Sequence analysis also revealed 89 protein-encoding sequences attributed to ENV, greatly expanding the amount of genetic information available for this virus. High-throughput PCR of over 19,000 fish showed that ENV is widely distributed in the NE Pacific Ocean and was detected in 12 of 16 tested species, including in 27% of herring, 38% of anchovy, 17% of pollock, and 13% of sand lance. Despite frequent detection in marine fish, ENV prevalence was significantly lower in fish from freshwater (0.03%), as assessed with a generalized linear mixed effects model (p = 5.5 × 10-8). Thus, marine fish are likely a reservoir for the virus. High genetic similarity between ENV obtained from salmon and herring also suggests that transmission between these hosts is likely.

Transcriptional shifts during juvenile Coho salmon (Oncorhynchus kisutch) life stage changes in freshwater and early marine environments.

There is a paucity of information on the physiological changes that occur over the course of salmon early marine migration. Here we aim to provide insight on juvenile Coho salmon (Oncorhynchus kisutch) physiology using the changes in gene expression (cGRASP 44K microarray) of four tissues (brain, gill, muscle, and liver) across the parr to smolt transition in freshwater and through the first eight months of ocean residence. We also examined transcriptome changes with body size as a covariate. The strongest shift in the transcriptome for brain, gill, and muscle occurred between summer and fall in the ocean, representing physiological changes that we speculate may be associated with migration preparation to feeding areas. Metabolic processes in the liver were positively associated with body length, generally consistent with enhanced feeding opportunities. However, a notable exception to this metabolic pattern was for spring post-smolts sampled soon after entry into the ocean, which showed a pattern of gene expression more likely associated with depressed feeding or recent fasting. Overall, this study has revealed life stages that may be the most critical developmentally (fall post-smolt) and for survival (spring post-smolt) in the early marine environment. These life stages may warrant further investigation.

The Genome of the Beluga Whale (Delphinapterus leucas).

The beluga whale is a cetacean that inhabits arctic and subarctic regions, and is the only living member of the genus Delphinapterus. The genome of the beluga whale was determined using DNA sequencing approaches that employed both microfluidic partitioning library and non-partitioned library construction. The former allowed for the construction of a highly contiguous assembly with a scaffold N50 length of over 19 Mbp and total reconstruction of 2.32 Gbp. To aid our understanding of the functional elements, transcriptome data was also derived from brain, duodenum, heart, lung, spleen, and liver tissue. Assembled sequence and all of the underlying sequence data are available at the National Center for Biotechnology Information (NCBI) under the Bioproject accession number PRJNA360851A.

Heart and skeletal muscle inflammation (HSMI) disease diagnosed on a British Columbia salmon farm through a longitudinal farm study.

Heart and skeletal muscle inflammation (HSMI) is an emerging disease of marine-farmed Atlantic Salmon (Salmo salar), first recognized in 1999 in Norway, and later also reported in Scotland and Chile. We undertook a longitudinal study involving health evaluation over an entire marine production cycle on one salmon farm in British Columbia (Canada). In previous production cycles at this farm site and others in the vicinity, cardiac lesions not linked to a specific infectious agent or disease were identified. Histologic assessments of both live and moribund fish samples collected at the farm during the longitudinal study documented at the population level the development, peak, and recovery phases of HSMI. The fish underwent histopathological evaluation of all tissues, Twort's Gram staining, immunohistochemistry, and molecular quantification in heart tissue of 44 agents known or suspected to cause disease in salmon. Our analysis showed evidence of HSMI histopathological lesions over an 11-month timespan, with the prevalence of lesions peaking at 80-100% in sampled fish, despite mild clinical signs with no associated elevation in mortalities reported at the farm level. Diffuse mononuclear inflammation and myodegeneration, consistent with HSMI, was the predominant histologic observation in affected heart and skeletal muscle. Infective agent monitoring identified three agents at high prevalence in salmon heart tissue, including Piscine orthoreovirus (PRV), and parasites Paranucleospora theridion and Kudoa thyrsites. However, PRV alone was statistically correlated with the occurrence and severity of histopathological lesions in the heart. Immunohistochemical staining further localized PRV throughout HSMI development, with the virus found mainly within red blood cells in early cases, moving into the cardiomyocytes within or, more often, on the periphery of the inflammatory reaction during the peak disease, and reducing to low or undetectable levels later in the production cycle. This study represents the first longitudinal assessment of HSMI in a salmon farm in British Columbia, providing new insights on the pathogenesis of the disease.

Infectious disease, shifting climates, and opportunistic predators: cumulative factors potentially impacting wild salmon declines.

Emerging diseases are impacting animals under high-density culture, yet few studies assess their importance to wild populations. Microparasites selected for enhanced virulence in culture settings should be less successful maintaining infectivity in wild populations, as once the host dies, there are limited opportunities to infect new individuals. Instead, moderately virulent microparasites persisting for long periods across multiple environments are of greatest concern. Evolved resistance to endemic microparasites may reduce susceptibilities, but as barriers to microparasite distributions are weakened, and environments become more stressful, unexposed populations may be impacted and pathogenicity enhanced. We provide an overview of the evolutionary and ecological impacts of infectious diseases in wild salmon and suggest ways in which modern technologies can elucidate the microparasites of greatest potential import. We present four case studies that resolve microparasite impacts on adult salmon migration success, impact of river warming on microparasite replication, and infection status on susceptibility to predation. Future health of wild salmon must be considered in a holistic context that includes the cumulative or synergistic impacts of multiple stressors. These approaches will identify populations at greatest risk, critically needed to manage and potentially ameliorate the shifts in current or future trajectories of wild populations.

Transcriptomics of environmental acclimatization and survival in wild adult Pacific sockeye salmon (Oncorhynchus nerka) during spawning migration.

Environmental shifts accompanying salmon spawning migrations from ocean feeding grounds to natal freshwater streams can be severe, with the underlying stress often cited as a cause of increased mortality. Here, a salmonid microarray was used to characterize changes in gene expression occurring between ocean and river habitats in gill and liver tissues of wild migrating sockeye salmon (Oncorhynchus nerka Walbaum) returning to spawn in the Fraser River, British Columbia, Canada. Expression profiles indicate that the transcriptome of migrating salmon is strongly affected by shifting abiotic and biotic conditions encountered along migration routes. Conspicuous shifts in gene expression associated with changing salinity, temperature, pathogen exposure and dissolved oxygen indicate that these environmental variables most strongly impact physiology during spawning migrations. Notably, transcriptional changes related to osmoregulation were largely preparatory and occurred well before salmon encountered freshwater. In the river environment, differential expression of genes linked with elevated temperatures indicated that thermal regimes within the Fraser River are approaching tolerance limits for adult salmon. To empirically correlate gene expression with survival, biopsy sampling of gill tissue and transcriptomic profiling were combined with telemetry. Many genes correlated with environmental variables were differentially expressed between premature mortalities and successful migrants. Parametric survival analyses demonstrated a broad-scale transcriptional regulator, cofactor required for Sp1 transcriptional activation (CRSP), to be significantly predictive of survival. As the environmental characteristics of salmon habitats continue to change, establishing how current environmental conditions influence salmon physiology under natural conditions is critical to conserving this ecologically and economically important fish species.

Development and application of DNA techniques for validating and improving pinniped diet estimates.

Polymerase chain reaction techniques were developed and applied to identify DNA from >40 species of prey contained in fecal (scat) soft-part matrix collected at terrestrial sites used by Steller sea lions (Eumetopias jubatus) in British Columbia and the eastern Aleutian Islands, Alaska. Sixty percent more fish and cephalopod prey were identified by morphological analyses of hard parts compared with DNA analysis of soft parts (hard parts identified higher relative proportions of Ammodytes sp., Cottidae, and certain Gadidae). DNA identified 213 prey occurrences, of which 75 (35%) were undetected by hard parts (mainly Salmonidae, Pleuronectidae, Elasmobranchii, and Cephalopoda), and thereby increased species occurrences by 22% overall and species richness in 44% of cases (when comparing 110 scats that amplified prey DNA). Prey composition was identical within only 20% of scats. Overall, diet composition derived from both identification techniques combined did not differ significantly from hard-part identification alone, suggesting that past scat-based diet studies have not missed major dietary components. However, significant differences in relative diet contributions across scats (as identified using the two techniques separately) reflect passage rate differences between hard and soft digesta material and highlight certain hypothesized limitations in conventional morphological-based methods (e.g., differences in resistance to digestion, hard part regurgitation, partial and secondary prey consumption), as well as potential technical issues (e.g., resolution of primer efficiency and sensitivity and scat subsampling protocols). DNA analysis of salmon occurrence (from scat soft-part matrix and 238 archived salmon hard parts) provided species-level taxonomic resolution that could not be obtained by morphological identification and showed that Steller sea lions were primarily consuming pink (Oncorhynchus gorbuscha) and chum (Oncorhynchus keta) salmon. Notably, DNA from Atlantic salmon (Salmo salar) that likely originated from a distant fish farm was also detected in two scats from one site in the eastern Aleutian Islands. Overall, molecular techniques are valuable for identifying prey in the fecal remains of marine predators. Combining DNA and hard-part identification will effectively alleviate certain predicted biases and will ultimately enhance measures of diet richness, fisheries interactions (especially salmon-related ones), and the ecological role of pinnipeds and other marine predators, to the benefit of marine wildlife conservationists and fisheries managers.

Salmon spawning migration: metabolic shifts and environmental triggers.

A large-scale functional genomics study revealed shifting metabolic processes in white muscle during the final 1300 km migration of wild sockeye salmon to their spawning grounds in the Fraser River, British Columbia. In 2006, Lower Adams stock sockeye salmon ceased feeding after passing the Queen Charlotte Islands, 850 km from the Fraser River. Enhanced protein turnover and reduced transcription of actin, muscle contractile and heme-related proteins were early starvation responses in saltwater. Arrival to the estuarine environment triggered massive protein turnover through induction of proteasomal and lysosomal proteolysis and protein biosynthesis, and a shift from anaerobic glycolysis to oxidative phosphorylation. Response to entry into freshwater was modest, with up-regulation of heat shock proteins and nitric oxide biosynthesis. High river temperatures resulted in a strong defense/immune response and high mortalities in 50% of fish. Arrival to the spawning grounds triggered further up-regulation of oxidative phosphorylation and proteolysis, down-regulation of protein biosynthesis and helicase activity, and continued down-regulation of muscle proteins and most glycolytic enzymes. However, sharp up-regulation of PFK-I indicated induction of glycolytic potential at the spawning grounds. The identification of potential environmental cues triggering genome-wide transcriptional shifts in white muscle associated with migration and the strong activation of proteasomal proteolysis were both novel findings.

The salmonid MHC class I: more ancient loci uncovered.

An unprecedented level of sequence diversity has been maintained in the salmonid major histocompatibility complex (MHC) class I UBA gene, with between lineage AA sequence identities as low as 34%. The derivation of deep allelic lineages may have occurred through interlocus exon shuffling or convergence of ancient loci with the UBA locus, but until recently, no such ancient loci were uncovered. Herein, we document the existence of eight additional MHC class I loci in salmon (UCA, UDA, UEA, UFA, UGA, UHA, ULA, and ZE), six of which share exon 2 and 3 lineages with UBA, and three of which have not been described elsewhere. Half of the UBA exon 2 lineages and all UBA exon 3 lineages are shared with other loci. Two loci, UGA and UEA, share only a single exon lineage with UBA, likely generated through exon shuffling. Based on sequence homologies, we hypothesize that most exchanges and duplications occurred before or during tetraploidization (50 to 100 Ma). Novel loci that share no relationship with other salmonid loci are also identified (UHA and ZE). Each locus is evaluated for its potential to function as a class Ia gene based on gene expression, conserved residues and polymorphism. UBA is the only locus that can indisputably be classified as a class Ia gene, although three of the eight loci (ZE, UCA, and ULA) conform in three out of four measures. We hypothesize that these additional loci are in varying states of degradation to class Ib genes.

Expansion and contraction of major histocompatibility complex genes: a teleostean example.

The MHC is a multigene family that has arisen through recurrent expansion and contraction of genes, and a continuum of the evolutionary process is observed in the teleost fishes. The number of duplicated genes observed in different phylogenetic groups of teleost fish varies from one to 42, with only a few genes observed in the primitive euteleost species, and greater numbers of genes observed in the more advanced neoteleost species. In this study, an attempt is made to isolate all of the Mhc class I genes of an early neoteleost species, Atlantic cod (Gadus morhua L.), in the superorder Paracanthopterygii. Eighty-three sequences were isolated from the cDNA of an individual G. morhua. The level of gene duplication observed within each of the lineages and sublineages was similar, and most contained an estimated two to four duplicated genes. Mhc class I gene duplication in G. morhua was independent of, and possibly more recent than, extensive duplication in the Acanthopterygian superorder. Only limited contraction of Mhc genes is observed in G. morhua. A low level of haplotype diversity is observed, with most individuals containing at least one copy of each of the lineages tested. Divergence of the conserved N- and C-terminal residues of the antigen recognition site is observed, indicative of the initial stage of degeneration from classical to non-classical genes. However, most or all of the lineages are still polymorphic, and degeneration is present both within and among lineages. Thus, the outcome (i.e., which genes will remain classical) is as yet undetermined.