Loss of genetic variation and ancestral sex determination system in North American northern pike characterized by whole-genome resequencing.

The northern pike Esox lucius is a freshwater fish with low genetic diversity but ecological success throughout the Northern Hemisphere. Here, we generate an annotated chromosome-level genome assembly of 941 Mbp in length with 25 chromosome-length scaffolds. We then genotype 47 northern pike from Alaska through New Jersey at a genome-wide scale and characterize a striking decrease in genetic diversity along the sampling range. Individuals west of the North American Continental Divide have substantially higher diversity than those to the east (e.g. Interior Alaska and St. Lawrence River have on average 181 and 64K heterozygous SNPs per individual, or a heterozygous SNP every 5.2 and 14.6 kbp, respectively). Individuals clustered within each population with strong support, with numerous private alleles observed within each population. Evidence for recent population expansion was observed for a Manitoba hatchery and the St. Lawrence population (Tajima's D = -1.07 and -1.30, respectively). Several chromosomes have large regions with elevated diversity, including LG24, which holds amhby, the ancestral sex determining gene. As expected amhby was largely male-specific in Alaska and the Yukon and absent southeast to these populations, but we document some amhby(-) males in Alaska and amhby(+) males in the Columbia River, providing evidence for a patchwork of presence of this system in the western region. These results support the theory that northern pike recolonized North America from refugia in Alaska and expanded following deglaciation from west to east, with probable founder effects resulting in loss of both neutral and functional diversity (e.g. amhby).

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function.

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Population-size history inferences from the coho salmon (Oncorhynchus kisutch) genome.

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.

Modulation of the Host Cell Transcriptome and Epigenome by Fusobacterium nucleatum.

Fusobacterium nucleatum is a ubiquitous opportunistic pathogen with an emerging role as an oncomicrobe in colorectal cancer and other cancer settings. F. nucleatum can adhere to and invade host cells in a manner that varies across F. nucleatum strains and host cell phenotypes. Here, we performed pairwise cocultures between three F. nucleatum strains and two immortalized primary host cell types (human colonic epithelial [HCE] cells and human carotid artery endothelial [HCAE] cells) followed by transcriptome sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) to investigate transcriptional and epigenetic host cell responses. We observed that F. nucleatum-induced host cell transcriptional modulation involves strong upregulation of genes related to immune migration and inflammatory processes, such as TNF, CXCL8, CXCL1, and CCL20. Furthermore, we identified genes strongly upregulated in a cell line-specific manner. In HCE cells, overexpressed genes included UBD and DUOX2/DUOXA2, associated with p53 degradation-mediated proliferation and intestinal reactive oxygen species (ROS) production, respectively. In HCAE cells, overexpressed genes included EFNA1 and LIF, two genes commonly upregulated in colorectal cancer and associated with poor patient outcomes, and PTGS2 (COX2), a gene associated with the protective effect of aspirin in the colorectal cancer setting. Interestingly, we also observed downregulation of numerous histone modification genes upon F. nucleatum exposure. We used the ChIP-seq data to annotate chromatin states genome wide and found significant chromatin remodeling following F. nucleatum exposure in HCAE cells, with increased frequencies of active enhancer and low-signal/quiescent states. Thus, our results highlight increased inflammation and chemokine gene expression as conserved host cell responses to F. nucleatum exposure and extensive host cell epigenomic changes specific to host cell type. IMPORTANCE Fusobacterium nucleatum is a bacterium normally found in the healthy oral cavity but also has an emerging role in colorectal cancer and other cancer settings. The host-microbe interactions of F. nucleatum and its involvement in tumor initiation, progression, and treatment resistance are not fully understood. We explored host cell changes that occur in response to F. nucleatum. We identified key genes differentially expressed in response to various conditions of F. nucleatum exposure and determined that the conserved host cell response to F. nucleatum was dominated by increased inflammation and chemokine gene expression. Additionally, we found extensive host cell epigenomic changes as a novel aspect of host modulation associated with F. nucleatum exposure. These results extend our understanding of F. nucleatum as an emerging pathogen and highlight the importance of considering strain heterogeneity and host cell phenotypic variation when exploring pathogenic mechanisms of F. nucleatum.

A 200K SNP chip reveals a novel Pacific salmon louse genotype linked to differential efficacy of emamectin benzoate.

Antiparasitic drugs such as emamectin benzoate (EMB) are relied upon to reduce the parasite load, particularly of the sea louse Lepeophtheirus salmonis, on farmed salmon. The decline in EMB treatment efficacy for this purpose is an important issue for salmon producers around the world, and particularly for those in the Atlantic Ocean where widespread EMB tolerance in sea lice is recognized as a significant problem. Salmon farms in the Northeast Pacific Ocean have not historically experienced the same issues with treatment efficacy, possibly due to the relatively large population of endemic salmonid hosts that serve to both redistribute surviving lice and dilute populations potentially under selection by introducing naïve lice to farms. Frequent migration of lice among farmed and wild hosts should limit the effect of farm-specific selection pressures on changes to the overall allele frequencies of sea lice in the Pacific Ocean. A previous study using microsatellites examined L. salmonis oncorhynchi from 10 Pacific locations from wild and farmed hosts and found no population structure. Recently however, a farm population of sea lice was detected where EMB bioassay exposure tolerance was abnormally elevated. In response, we have developed a Pacific louse draft genome that complements the previously-released Atlantic louse sequence. These genomes were combined with whole-genome re-sequencing data to design a highly sensitive 201,279 marker SNP array applicable for both subspecies (90,827 validated Pacific loci; 153,569 validated Atlantic loci). Notably, kmer spectrum analysis of the re-sequenced samples indicated that Pacific lice exhibit a large within-individual heterozygosity rate (average of 1 in every 72 bases) that is markedly higher than that of Atlantic individuals (1 in every 173 bases). The SNP chip was used to produce a high-density map for Atlantic sea louse linkage group 5 that was previously shown to be associated with EMB tolerance in Atlantic lice. Additionally, 478 Pacific louse samples from farmed and wild hosts obtained between 2005 and 2014 were also genotyped on the array. Clustering analysis allowed us to detect the apparent emergence of an otherwise rare genotype at a high frequency among the lice collected from two farms in 2013 that had reported elevated EMB tolerance. This genotype was not observed in louse samples collected from the same farm in 2010, nor in any lice sampled from other locations prior to 2013. However, this genotype was detected at low frequencies in louse samples from farms in two locations reporting elevated EMB tolerance in 2014. These results suggest that a rare genotype present in Pacific lice may be locally expanded in farms after EMB treatment. Supporting this hypothesis, 437 SNPs associated with this genotype were found to be in a region of linkage group 5 that overlaps the region associated with EMB resistance in Atlantic lice. Finally, five of the top diagnostic SNPs within this region were used to screen lice that had been subjected to an EMB survival assay, revealing a significant association between these SNPs and EMB treatment outcome. To our knowledge this work is the first report to identify a genetic link to altered EMB efficacy in L. salmonis in the Pacific Ocean.