Gill raker and pyloric caeca counts differ between Arctic char (Salvelinus alpinus) and Dolly Varden (S. malma) populations across their ranges.

Meristic characters are often used to differentiate between closely related forms, morphs, and species of fishes, and lend insight into ecology and post-glacial recolonization in taxa with complicated or contentious phylogenies, including the genus Salvelinus. Previous studies of meristics in Salvelinus have focused mostly on individual populations. We collated data from 456 populations/systems across the North American and Russian Arctic and sub-Arctic, and found that counts of pyloric caeca and gill rakers differed consistently between fish visually and/or genetically identified as Arctic char and Dolly Varden across their distributional ranges.

Characterization and distribution of a 14-Mb chromosomal inversion in native populations of rainbow trout (Oncorhynchus mykiss).

Multiple studies in a range of taxa have found links between structural variants and the development of ecologically important traits. Such variants are becoming easier to find due, in large part, to the increase in the amount of genome-wide sequence data in nonmodel organisms. The salmonids (salmon, trout, and charr) are a taxonomic group with abundant genome-wide datasets due to their importance in aquaculture, fisheries, and variation in multiple ecologically important life-history traits. Previous research on rainbow trout (Oncorhynchus mykiss) has documented a large pericentric (∼55 Mb) chromosomal inversion (CI) on chromosome 5 (Omy05) and a second smaller (∼14 Mb) chromosome inversion on Omy20. While the Omy05 inversion appears to be associated with multiple adaptive traits, the inversion on Omy20 has received far less attention. In this study, we re-analyze RAD-seq and amplicon data from several populations of rainbow trout (O. mykiss) to better document the structure and geographic distribution of variation in the Omy20 CI. Moreover, we utilize phylogenomic techniques to characterize both the age- and the protein-coding gene content of the Omy20 CI. We find that the age of the Omy20 inversion dates to the early stages of O. mykiss speciation and predates the Omy05 inversion by ∼450,000 years. The 2 CIs differ further in terms of the frequency of the homokaryotypes. While both forms of the Omy05 CI are found across the eastern Pacific, the ancestral version of the Omy20 CI is restricted to the southern portion of the species range in California. Furthermore, the Omy20 inverted haplotype is comparable in genetic diversity to the ancestral form, whereas derived CIs typically show substantially reduced genetic diversity. These data contribute to our understanding of the age and distribution of a large CI in rainbow trout and provide a framework for researchers looking to document CIs in other nonmodel species.

Genomic structural variation in Barramundi Perch Lates calcarifer and potential roles in speciation and adaptation.

Advances in genome sequencing and assembly techniques have increased the documentation of structural variants in wild organisms. Of these variants, chromosomal inversions are especially prominent due to their large size and active recombination suppression between alternative homokaryotypes. This suppression enables the two forms of the inversion to be maintained and allows preservation of locally adapted alleles. The Barramundi Perch (Lates calcarifer) is a widespread species complex with three main genetic lineages located in the biogeographic regions of Australia and New Guinea (AUS+NG), Southeast Asia (SEA), and the Indian Subcontinent (IND). Barramundi Perch are typically considered to be a protandrous sequential hermaphrodite species that exhibits catadromy. Freshwater occupancy and intraspecific variation in life history (e.g., partially migratory populations) exist and provide opportunities for strongly divergent selection associated with, for example, salinity tolerance, swimming ability, and marine dispersal. Herein, we utilize genomic data generated from all three genetic lineages to identify and describe three polymorphic candidate chromosomal inversions. These candidate chromosomal inversions appear to be fixed for ancestral variants in the IND lineage, fixed for inverted versions in the AUS+NG lineage, and exhibit variation in all three inversions in the SEA lineage. Barramundi Perch have a diverse portfolio of life history options that includes migratory strategy as well as sexual system (i.e., hermaphroditism and gonochorism). We propose that the some of the life history variability observed in Barramundi Perch may be linked to inversions and, in so doing, we present genetic data that might be useful in enhancing aquaculture production and population management.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha.

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Comparative genomics of rainbow trout (Oncorhynchus mykiss): Is the genetic architecture of migratory behavior conserved among populations?

Rainbow trout (Oncorhynchus mykiss) are a partially migratory species wherein some individuals undergo long-distance anadromous migrations, and others stay as residents in their native freshwater streams. The decision to migrate is known to be highly heritable, and yet, the underlying genes and alleles associated with migration are not fully characterized. Here we used a pooled approach of whole-genome sequence data from migratory and resident trout of two native populations-Sashin Creek, Alaska and Little Sheep Creek, Oregon-to obtain a genome-wide perspective of the genetic architecture of resident and migratory life history. We calculated estimates of genetic differentiation, genetic diversity, and selection between the two phenotypes to locate regions of interest and then compared these associations between populations. We identified numerous genes and alleles associated with life history development in the Sashin Creek population with a notable area on chromosome 8 that may play a critical role in the development of the migratory phenotype. However, very few alleles appeared to be associated with life history development in the Little Sheep Creek system, suggesting population-specific genetic effects are likely important in the development of anadromy. Our results indicate that a migratory life history is not controlled by a singular gene or region but supports the idea that there are many independent ways for a migratory phenotype to emerge in a population. Therefore, conserving and promoting genetic diversity in migratory individuals is paramount to conserving these populations. Ultimately, our data add to a growing body of literature that suggests that population-specific genetic effects, likely mediated through environmental variation, contribute to life history development in rainbow trout.

Genomics reveal the origins and current structure of a genetically depauperate freshwater species in its introduced Alaskan range.

Invasive species are a major threat to global biodiversity, yet also represent large-scale unplanned ecological and evolutionary experiments to address fundamental questions in nature. Here we analyzed both native and invasive populations of predatory northern pike (Esox lucius) to characterize landscape genetic variation, determine the most likely origins of introduced populations, and investigate a presumably postglacial population from Southeast Alaska of unclear provenance. Using a set of 4329 SNPs from 351 individual Alaskan northern pike representing the most widespread geographic sampling to date, our results confirm low levels of genetic diversity in native populations (average 𝝅 of 3.18 × 10-4) and even less in invasive populations (average 𝝅 of 2.68 × 10-4) consistent with bottleneck effects. Our analyses indicate that invasive northern pike likely came from multiple introductions from different native Alaskan populations and subsequently dispersed from original introduction sites. At the broadest scale, invasive populations appear to have been founded from two distinct regions of Alaska, indicative of two independent introduction events. Genetic admixture resulting from introductions from multiple source populations may have mitigated the negative effects associated with genetic bottlenecks in this species with naturally low levels of genetic diversity. Genomic signatures strongly suggest an excess of rare, population-specific alleles, pointing to a small number of founding individuals in both native and introduced populations consistent with a species' life history of limited dispersal and gene flow. Lastly, the results strongly suggest that a small isolated population of pike, located in Southeast Alaska, is native in origin rather than stemming from a contemporary introduction event. Although theory predicts that lack of genetic variation may limit colonization success of novel environments, we detected no evidence that a lack of standing variation limited the success of this genetically depauperate apex predator.

A candidate chromosome inversion in Arctic charr (Salvelinus alpinus) identified by population genetic analysis techniques.

The "genomics era" has allowed questions to be asked about genome organization and genome architecture of non-model species at a rate not previously seen. Analyses of these genome-wide datasets have documented many examples of novel structural variants (SVs) such as chromosomal inversions, copy number variants, and chromosomal translocations, many of which have been linked to adaptation. The salmonids are a taxonomic group with abundant genome-wide datasets due to their importance in aquaculture and fisheries. However, the number of documented SVs in salmonids is surprisingly low and is most likely due to removing loci in high linkage disequilibrium when analyzing structure and gene flow. Here we re-analyze RAD-seq data from several populations of Arctic charr (Salvelinus alpinus) and document a novel ∼1.2 MB SV at the distal end of LG12. This variant contains 15 protein-coding genes connected to a wide-range of functions including cell adhesion and signal transduction. Interestingly, we studied the frequency of this polymorphism in four disjointed populations of charr-one each from Nunavut, Newfoundland, Eastern Russia, and Scotland-and found evidence of the variant only in Nunavut, Canada, suggesting the polymorphism is novel and recently evolved.

Differential gene expression associated with behavioral variation in ecotypes of Lake Superior brook trout (Salvelinus fontinalis).

Associations between behaviors and the development of different life history tactics have been documented in several species of salmon, trout, and charr. While it is well known that such behaviors are heritable the genes and molecular pathways connected to these behaviors remain unknown. We used an RNA-seq approach to identify genes and molecular pathways differentially regulated in brain tissue between "shy" and "bold" brook trout (Salvelinus fontinalis). A small number of genes were differentially expressed between the behavioral types at several months after hatching and two years of age. Pathway analysis revealed that EIF2 signaling differed consistently between shy and bold individuals suggesting large-scale differences in protein synthesis between behavioral types in the brain. Additionally, the RNA-seq data were used to find polymorphisms within the brook trout genome and a GWAS approach was used to test for statistical associations between genetic variants and behavior type. One allele located in a transcription factor (TSHZ3) contained a protein-coding non-synonymous SNP suggesting that functional variation within TSHZ3 is connected to the development of different behaviors. These results suggest that the molecular basis of behavioral development is complex and due to the differential expression of many genes involved in a wide-range of different molecular pathways.

Comparative Genomic Analyses and a Novel Linkage Map for Cisco (Coregonus artedi) Provide Insights into Chromosomal Evolution and Rediploidization Across Salmonids.

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco (Coregonus artedi), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts.

A large-scale chromosomal inversion is not associated with life history development in rainbow trout from Southeast Alaska.

In studying the causative mechanisms behind migration and life history, the salmonids-salmon, trout, and charr-are an exemplary taxonomic group, as life history development is known to have a strong genetic component. A double inversion located on chromosome 5 in rainbow trout (Oncorhynchus mykiss) is associated with life history development in multiple populations, but the importance of this inversion has not been thoroughly tested in conjunction with other polymorphisms in the genome. To that end, we used a high-density SNP chip to genotype 192 F1 migratory and resident rainbow trout and focused our analyses to determine whether this inversion is important in life history development in a well-studied population of rainbow trout from Southeast Alaska. We identified 4,994 and 436 SNPs-predominantly outside of the inversion region-associated with life history development in the migrant and resident familial lines, respectively. Although F1 samples showed genomic patterns consistent with the double inversion on chromosome 5 (reduced observed and expected heterozygosity and an increase in linkage disequilibrium), we found no statistical association between the inversion and life history development. Progeny produced by crossing resident trout and progeny produced by crossing migrant trout both consisted of a mix of migrant and resident individuals, irrespective of the individuals' inversion haplotype on chromosome 5. This suggests that although the inversion is present at a low frequency, it is not strongly associated with migration as it is in populations of Oncorhynchus mykiss from lower latitudes.

Long-Term Conservation of Ohnologs Through Partial Tetrasomy Following Whole-Genome Duplication in Salmonidae.

Whole-genome duplications (WGDs) have occurred repeatedly and broadly throughout the evolutionary history of eukaryotes. However, the effects of WGD on genome function and evolution remain unclear. The salmonid WGD that occurred approximately 88 million years ago presents an excellent opportunity for studying the effects of WGD as ∼10-15% of each salmonid genome still exhibits tetrasomic inheritance. Herein, we utilized the rainbow trout (Oncorhynchus mykiss) genome assembly and brain transcriptome data to examine the fate of gene pairs (ohnologs) following the salmonid whole-genome duplication. We find higher sequence identity between ohnologs located within known tetrasomic regions than between ohnologs found in disomic regions, and that tetrasomically inherited ohnologs showed greater similarity in patterns of gene expression and per ohnolog were lower expressed, than disomically inherited ohnologs. Enrichment testing for Gene Ontology terms identified 49 over-represented terms in tetrasomically inherited ohnologs compared to disomic ohnologs. However, why these ohnologs are retained as tetrasomic is difficult to answer. It could be that we have identified salmonid specific "dangerous duplicates", that is, genes that cannot take on new roles following WGD. Alternatively, there may be adaptive advantages for retaining genes as functional duplicates in tetrasomic regions, as presumably, movement of these genes into disomic regions would affect both their sequence identity and their gene expression patterns.

Transcriptional profiling of the clpX mutant in Bacillus anthracis reveals regulatory connection with the lrgAB operon.

ClpX functions as either an independent chaperone or a component of the ClpXP protease, a conserved intracellular protease that acts as a global regulator in the bacterial cell by degrading regulatory proteins, stress response proteins and rate-limiting enzymes. Previously, we found that loss of clpX in Bacillus anthracis Sterne leads to increased susceptibility to antimicrobial agents that target the cell envelope. The aim of this study was to identify genes within the regulatory network of clpX that contribute to antimicrobial resistance. Using microarray analysis, we found 119 genes that are highly differentially expressed in the ∆clpX mutant, with the majority involved in metabolic, transport or regulatory functions. Several of these differentially expressed genes, including glpF, sigM, mrsA, lrgA and lrgB, are associated with cell wall-active antibiotics in other bacterial species. We focused on lrgA and lrgB, which form the lrgAB operon and are downregulated in ∆clpX, because loss of lrgAB increases autolytic activity and penicillin susceptibility in Staphylococcus aureus. While we observed no changes in autolytic activity in either ∆clpX or ∆lrgAB B. anthracis Sterne, we find that both mutants have increased susceptibility to the antimicrobial peptide LL-37 and daptomycin. However, phenotypes between ∆clpX and ∆lrgAB are not identical as ∆clpX also displays increased susceptibility to penicillin and nisin but ∆lrgAB does not. Therefore, while decreased expression of lrgAB may be partially responsible for the increased antimicrobial susceptibility seen in the ∆clpX mutant, disruption of other pathways must also contribute to this phenotype.

Evidence of sex-bias in gene expression in the brain transcriptome of two populations of rainbow trout (Oncorhynchus mykiss) with divergent life histories.

Sex-bias in gene expression is a mechanism that can generate phenotypic variance between the sexes, however, relatively little is known about how patterns of sex-bias vary during development, and how variable sex-bias is between different populations. To that end, we measured sex-bias in gene expression in the brain transcriptome of rainbow trout (Oncorhynchus mykiss) during the first two years of development. Our sampling included from the fry stage through to when O. mykiss either migrate to the ocean or remain resident and undergo sexual maturation. Samples came from two F1 lines: One from migratory steelhead trout and one from resident rainbow trout. All samples were reared in a common garden environment and RNA sequencing (RNA-seq) was used to estimate patterns of gene expression. A total of 1,716 (4.6% of total) genes showed evidence of sex-bias in gene expression in at least one time point. The majority (96.7%) of sex-biased genes were differentially expressed during the second year of development, indicating that patterns of sex-bias in expression are tied to key developmental events, such as migration and sexual maturation. Mapping of differentially expressed genes to the O. mykiss genome revealed that the X chromosome is enriched for female upregulated genes, and this may indicate a lack of dosage compensation in rainbow trout. There were many more sex-biased genes in the migratory line than the resident line suggesting differences in patterns of gene expression in the brain between populations subjected to different forces of selection. Overall, our results suggest that there is considerable variation in the extent and identity of genes exhibiting sex-bias during the first two years of life. These differentially expressed genes may be connected to developmental differences between the sexes, and/or between adopting a resident or migratory life history.

Using RNA-seq to determine patterns of sex-bias in gene expression in the brain of the sex-role reversed Gulf Pipefish (Syngnathus scovelli).

Sex-bias in gene expression is a widespread mechanism for controlling the development of phenotypes that differ between males and females. Most studies on sex-bias in gene expression have focused on species that exhibit traditional sex-roles (male-male competition and female parental care). By contrast the Syngnathid fishes (sea horses, pipefish, and sea dragons) are a group of organisms where many species exhibit male brooding and sex-role reversal (female-female competition for mates and paternal parental care), and little is known about how patterns of sex-bias in gene expression vary in species with sex-role reversal. Here we utilize RNA-seq technology to investigate patterns of sex-bias in gene expression in the brain tissue of the Gulf Pipefish (Syngnathus scovelli) a species that exhibits sex-role reversal. Gene expression analysis identified 73 sex-biased genes, 26 genes upregulated in females and 47 genes upregulated in males. Gene ontology analysis found 52 terms enriched for the sex-biased genes in a wide range of pathways suggesting that multiple functions and processes differ between the sexes. We focused on two areas of interest: sex steroids/hormones and circadian rhythms, both of which exhibited sex-bias in gene expression, and are known to influence sexual development in other species. Lastly, the work presented herein contributes to a growing body of genome data available for the Syngnathids, increasing our knowledge on patterns of gene expression in these unusual fishes.

Using Linkage Maps as a Tool To Determine Patterns of Chromosome Synteny in the Genus Salvelinus.

Next generation sequencing techniques have revolutionized the collection of genome and transcriptome data from nonmodel organisms. This manuscript details the application of restriction site-associated DNA sequencing (RADseq) to generate a marker-dense genetic map for Brook Trout (Salvelinus fontinalis). The consensus map was constructed from three full-sib families totaling 176 F1 individuals. The map consisted of 42 linkage groups with a total female map size of 2502.5 cM, and a total male map size of 1863.8 cM. Synteny was confirmed with Atlantic Salmon for 38 linkage groups, with Rainbow Trout for 37 linkage groups, Arctic Char for 36 linkage groups, and with a previously published Brook Trout linkage map for 39 linkage groups. Comparative mapping confirmed the presence of 8 metacentric and 34 acrocentric chromosomes in Brook Trout. Six metacentric chromosomes seem to be conserved with Arctic Char suggesting there have been at least two species-specific fusion and fission events within the genus Salvelinus In addition, the sex marker (sdY; sexually dimorphic on the Y chromosome) was mapped to Brook Trout BC35, which is homologous with Atlantic Salmon Ssa09qa, Rainbow Trout Omy25, and Arctic Char AC04q. Ultimately, this linkage map will be a useful resource for studies on the genome organization of Salvelinus, and facilitates comparisons of the Salvelinus genome with Salmo and Oncorhynchus.

RNA-seq reveals differential gene expression in the brains of juvenile resident and migratory smolt rainbow trout (Oncorhynchus mykiss).

Many migratory traits are heritable, but there is a paucity of evidence identifying the molecular mechanisms underlying differentiation in alternative migratory tactics, or in linking variation in gene expression to migratory behaviors. To that end, we examined differential gene expression in the brain transcriptome between young steelhead trout that had undergone the smoltification process, and resident rainbow trout (Oncorhynchus mykiss) from Sashin Creek, Alaska. Samples were sequenced from two time points: immediately before (at 20months of age) and during (2years of age) the presumed peak of smoltification. Smolt and resident individuals came from two genetic crosses, one where both parents were migratory, and another where both parents were residents. A total of 533 (1.9%) genes were differentially expressed between crosses, or between smolt and resident samples. These genes include some candidate migratory genes (such as POMC), as well as genes with no previous known involvement in the migratory process. Progeny from resident parents showed more upregulated genes than progeny from migrant parents at both time points. Pathway analysis showed enrichment in 227 biological pathways between cross type, and 171 biological pathways were enriched between residents and smolts. Enriched pathways had connections to many biofunctions, and most were only enriched in one contrast. However, pathways connected to phototransduction were enriched between both cross type and migratory tactics in 11 out of 12 contrasts, suggesting there are fundamental differences in how smolts and residents process light in the brain. The genes and pathways described herein constitute an a priori candidate list for future studies of migration in other populations of O. mykiss, and other migratory species.

Ontogenetic changes in embryonic and brain gene expression in progeny produced from migratory and resident Oncorhynchus mykiss.

Little information has been gathered regarding the ontogenetic changes that contribute to differentiation between resident and migrant individuals, particularly before the onset of gross morphological and physiological changes in migratory individuals. The aim of this study was to evaluate gene expression during early development in Oncorhynchus mykiss populations with different life histories, in a tissue known to integrate environmental cues to regulate complex developmental processes and behaviours. We sampled offspring produced from migrant and resident parents, collecting whole embryos prior to the beginning of first feeding, and brain tissue at three additional time points over the first year of development. RNA sequencing for 32 individuals generated a reference transcriptome of 30 177 genes that passed count thresholds. Differential gene expression between migrant and resident offspring was observed for 1982 genes. The greatest number of differentially expressed genes occurred at 8 months of age, in the spring a full year before the obvious physiological transformation from stream-dwelling parr to sea water-adaptable smolts begins for migrant individuals. Sex and age exhibited considerable effects on differential gene expression between migrants and resident offspring. Differential gene expression was observed in genes previously associated with migration, but also in genes previously unassociated with early life history divergence. Pathway analysis revealed coordinated differential expression in genes related to phototransduction, which could modulate photoperiod responsiveness and variation in circadian rhythms. The role for early differentiation in light sensitivity and biological rhythms is particularly intriguing in understanding early brain processes involved in differentiation of migratory and resident life history types.

Mapping and Expression of Candidate Genes for Development Rate in Rainbow Trout (Oncorhynchus mykiss).

Development rate has important implications for individual fitness and physiology. In salmonid fishes, development rate correlates with many traits later in life, including life-history diversity, growth, and age and size at sexual maturation. In rainbow trout (Oncorhynchus mykiss), a quantitative trait locus for embryonic development rate has been detected on chromosome 5 across populations. However, few candidate genes have been identified within this region. In this study, we use gene mapping, gene expression, and quantitative genetic methods to further identify the genetic basis of embryonic developmental rate in O. mykiss Among the genes located in the region of the major development rate quantitative trait locus (GHR1, Clock1a, Myd118-1, and their paralogs), all were expressed early in embryonic development (fertilization through hatch), but none were differentially expressed between individuals with the fast- or slow-developing alleles for a major embryonic development rate quantitative trait locus. In a follow-up study of migratory and resident rainbow trout from natural populations in Alaska, we found significant additive variation in development rate and, moreover, found associations between development rate and allelic variation in all 3 candidate genes within the quantitative trait locus for embryonic development. The mapping of these genes to this region and associations in multiple populations provide positional candidates for further study of their roles in growth, development, and life-history diversity in this model salmonid.

Microsatellite analyses across three diverse vertebrate transcriptomes (Acipenser fulvescens, Ambystoma tigrinum, and Dipodomys spectabilis).

Historically, many population genetics studies have utilized microsatellite markers sampled at random from the genome and presumed to be selectively neutral. Recent studies, however, have shown that microsatellites can occur in transcribed regions, where they are more likely to be under selection. In this study, we mined microsatellites from transcriptomes generated by 454-pyrosequencing for three vertebrate species: lake sturgeon (Acipenser fulvescens), tiger salamander (Ambystoma tigrinum), and kangaroo rat (Dipodomys spectabilis). We evaluated (i) the occurrence of microsatellites across species; (ii) whether particular gene ontology terms were over-represented in genes that contained microsatellites; (iii) whether repeat motifs were located in untranslated regions or coding sequences of genes; and (iv) in silico polymorphism. Microsatellites were less common in tiger salamanders than in either lake sturgeon or kangaroo rats. Across libraries, trinucleotides were found more frequently than any other motif type, presumably because they do not cause frameshift mutations. By evaluating variation across reads assembled to a given contig, we were able to identify repeat motifs likely to be polymorphic. Our study represents one of the first comparative data sets on the distribution of vertebrate microsatellites within expressed genes. Our results reinforce the idea that microsatellites do not always occur in noncoding DNA, but commonly occur in expressed genes.

Evaluating adaptive divergence between migratory and nonmigratory ecotypes of a salmonid fish, Oncorhynchus mykiss.

Next-generation sequencing and the application of population genomic and association approaches have made it possible to detect selection and unravel the genetic basis to variable phenotypic traits. The use of these two approaches in parallel is especially attractive in nonmodel organisms that lack a sequenced and annotated genome, but only works well when population structure is not confounded with the phenotype of interest. Herein, we use population genomics in a nonmodel fish species, rainbow trout (Oncorhynchus mykiss), to better understand adaptive divergence between migratory and nonmigratory ecotypes and to further our understanding about the genetic basis of migration. Restriction site-associated DNA (RAD) tag sequencing was used to identify single-nucleotide polymorphisms (SNPs) in migrant and resident O. mykiss from two systems, one in Alaska and the other in Oregon. A total of 7920 and 6755 SNPs met filtering criteria in the Alaska and Oregon data sets, respectively. Population genetic tests determined that 1423 SNPs were candidates for selection when loci were compared between resident and migrant samples. Previous linkage mapping studies that used RAD DNA tag SNPs were available to determine the position of 1990 markers. Several significant SNPs are located in genome regions that contain quantitative trait loci for migratory-related traits, reinforcing the importance of these regions in the genetic basis of migration/residency. Annotation of genome regions linked to significant SNPs revealed genes involved in processes known to be important in migration (such as osmoregulatory function). This study adds to our growing knowledge on adaptive divergence between migratory and nonmigratory ecotypes of this species; across studies, this complex trait appears to be controlled by many loci of small effect, with some in common, but many loci not shared between populations studied.