Cell Simulation as Cell Segmentation.

Single-cell spatial transcriptomics promises a highly detailed view of a cell's transcriptional state and microenvironment, yet inaccurate cell segmentation can render this data murky by misattributing large numbers of transcripts to nearby cells or conjuring nonexistent cells. We adopt methods from ab initio cell simulation to rapidly infer morphologically plausible cell boundaries that preserve cell type heterogeneity. Benchmarking applied to datasets generated by three commercial platforms show superior performance and computational efficiency of this approach compared with existing methods. We show that improved accuracy in cell segmentation aids greatly in detection of difficult to accurately segment tumor infiltrating immune cells such as neutrophils and T cells. Lastly, through improvements in our ability to delineate subsets of tumor infiltrating T cells, we show that CXCL13-expressing CD8+ T cells tend to be more closely associated with tumor cells than their CXCL13-negative counterparts in data generated from renal cell carcinoma patient samples. Proseg is available under at open source license at https://github.com/dcjones/proseg.

Modulation of Influenza A virus NS1 expression reveals prioritization of host response antagonism at single-cell resolution.

Influenza A virus (IAV) is an important human respiratory pathogen that causes significant seasonal epidemics and potential devastating pandemics. As part of its life cycle, IAV encodes the multifunctional protein NS1, that, among many roles, prevents immune detection and limits interferon (IFN) production. As distinct host immune pathways exert different selective pressures against IAV, as replication progresses, we expect a prioritization in the host immune antagonism by NS1. In this work, we profiled bulk transcriptomic differences in a primary bronchial epithelial cell model facing IAV infections at distinct NS1 levels. We further demonstrated that, at single cell level, the intracellular amount of NS1 in-part shapes the heterogeneity of the host response. We found that modulation of NS1 levels reveal a ranking in its inhibitory roles: modest NS1 expression is sufficient to inhibit immune detection, and thus the expression of pro-inflammatory cytokines (including IFNs), but higher levels are required to inhibit IFN signaling and ISG expression. Lastly, inhibition of chaperones related to the unfolded protein response requires the highest amount of NS1, often associated with later stages of viral replication. This work demystifies some of the multiple functions ascribed to IAV NS1, highlighting the prioritization of NS1 in antagonizing the different pathways involved in the host response to IAV infection.

Invasion genetics of the silver carp Hypophthalmichthys molitrix across North America: Differentiation of fronts, introgression, and eDNA metabarcode detection.

In the 1970s, the introduced silver carp Hypophthalmichthys molitrix (which is indigenous to eastern Asia) escaped from southern U.S. aquaculture to spread throughout the Mississippi River basin, and since has steadily moved northward. This large, prolific filter-feeder reduces food availability for other fishes. It now has reached the threshold of the Laurentian Great Lakes, where it likely will significantly impact food chains and fisheries. Our study evaluates population genetic variability and differentiation of the silver carp using 10 nuclear DNA microsatellite loci, and sequences of two mitochondrial genes-cytochrome b and cytochrome c oxidase subunit 1, along with the nuclear ribosomal protein S7 gene intron 1. We analyze population samples from: two primary Great Lakes' invasion fronts (at the Illinois River outside of Chicago, IL in Lake Michigan and in the Wabash River, which leads into the Maumee River and western Lake Erie), the original establishment "core" in the Lower Mississippi River, and expansion areas in the Upper Mississippi and Missouri rivers. We analyze and compare our results with bighead and other invasive carps, and cyprinid relatives. Results reveal that the silver carp invasion possesses moderate levels of genetic diversity, with more mtDNA haplotypes and unique microsatellite alleles in the "core" Lower Mississippi River population, which also diverges the most. The two invasion fronts also significantly genetically differ. About 3% of individuals (including all populations except the Illinois River) contain a unique and very divergent mtDNA haplotype, which likely stems from historic introgression in Asia with female largescale silver carp H. harmandi. The nuclear microsatellites and S7 sequences of the introgressed individuals do not differ from silver carp and are very distant from bighead carp. These sequence variation data are employed to design and evaluate a targeted high-throughput metabarcoding sequence assay that identifies and distinguishes among species of invasive carps (i.e., silver, bighead, grass, black, and common carps, along with goldfish), as well as native cyprinids, using cytochrome b. Our assay further differentiates among selected silver carp haplotypes (including between H. molitrix and H. harmandi), for use in population genetics and future analyses of spread pathways. We test and evaluate this assay on environmental (e)DNA water samples from 48 bait shops in the Great Lakes' region (along the Lake Erie, Lake St. Clair, and Wabash River watersheds), using positive and negative controls and custom bioinformatic processing. Test results discern silver carp eDNA in four of the shops-three in Lake Erie and one in the Wabash River watershed-and bighead carp from one of the same Lake Erie venues, suggesting that retailers (who often source from established southerly populations) comprise another introduction vector. Our overall findings thus provide key population genetic and phylogenetic data for understanding and tracing introductions, vectors, and spread pathways for silver carp, their variants, and their relatives.

Human-mediated evolution in a threatened species? Juvenile life-history changes in Snake River salmon.

Evaluations of human impacts on Earth's ecosystems often ignore evolutionary changes in response to altered selective regimes. Freshwater habitats for Snake River fall Chinook salmon (SRFCS), a threatened species in the US, have been dramatically changed by hydropower development and other watershed modifications. Associated biological changes include a shift in juvenile life history: Historically essentially 100% of juveniles migrated to sea as subyearlings, but a substantial fraction have migrated as yearlings in recent years. In contemplating future management actions for this species should major Snake River dams ever be removed (as many have proposed), it will be important to understand whether evolution is at least partially responsible for this life-history change. We hypothesized that if this trait is genetically based, parents who migrated to sea as subyearlings should produce faster-growing offspring that would be more likely to reach a size threshold to migrate to sea in their first year. We tested this with phenotypic data for over 2,600 juvenile SRFCS that were genetically matched to parents of hatchery and natural origin. Three lines of evidence supported our hypothesis: (i) the animal model estimated substantial heritability for juvenile growth rate for three consecutive cohorts; (ii) linear modeling showed an association between juvenile life history of parents and offspring growth rate; and (iii) faster-growing juveniles migrated at greater speeds, as expected if they were more likely to be heading to sea. Surprisingly, we also found that parents reared a full year in a hatchery produced the fastest growing offspring of all-apparently an example of cross-generational plasticity associated with artificial propagation. We suggest that SRFCS is an example of a potentially large class of species that can be considered to be "anthro-evolutionary"-signifying those whose evolutionary trajectories have been profoundly shaped by altered selective regimes in human-dominated landscapes.

Large-Scale Genotyping-by-Sequencing Indicates High Levels of Gene Flow in the Deep-Sea Octocoral Swiftia simplex (Nutting 1909) on the West Coast of the United States.

Deep-sea corals are a critical component of habitat in the deep-sea, existing as regional hotspots for biodiversity, and are associated with increased assemblages of fish, including commercially important species. Because sampling these species is so difficult, little is known about the connectivity and life history of deep-sea octocoral populations. This study evaluates the genetic connectivity among 23 individuals of the deep-sea octocoral Swiftia simplex collected from Eastern Pacific waters along the west coast of the United States. We utilized high-throughput restriction-site associated DNA (RAD)-tag sequencing to develop the first molecular genetic resource for the deep-sea octocoral, Swiftia simplex. Using this technique we discovered thousands of putative genome-wide SNPs in this species, and after quality control, successfully genotyped 1,145 SNPs across individuals sampled from California to Washington. These SNPs were used to assess putative population structure across the region. A STRUCTURE analysis as well as a principal coordinates analysis both failed to detect any population differentiation across all geographic areas in these collections. Additionally, after assigning individuals to putative population groups geographically, no significant FST values could be detected (FST for the full data set 0.0056), and no significant isolation by distance could be detected (p = 0.999). Taken together, these results indicate a high degree of connectivity and potential panmixia in S. simplex along this portion of the continental shelf.

Comparative genome mapping between Chinook salmon (Oncorhynchus tshawytscha) and rainbow trout (O. mykiss) based on homologous microsatellite loci.

Comparative genome mapping can rapidly facilitate the transfer of DNA sequence information from a well-characterized species to one that is less described. Chromosome arm numbers are conserved between members of the teleost family Salmonidae, order Salmoniformes, permitting rapid alignment of large syntenic blocks of DNA between members of the group. However, extensive Robertsonian rearrangements after an ancestral whole-genome duplication event has resulted in different chromosome numbers across Salmonid taxa. In anticipation of the rapid application of genomic data across members of the Pacific salmon genus Oncorhynchus, we mapped the genome of Chinook salmon (O. tshawytscha) by using 361 microsatellite loci and compared linkage groups to those already derived for a well-characterized species rainbow trout (O. mykiss). The Chinook salmon female map length was 1526 cM, the male map 733 cM, and the consensus map between the two sexes was 2206 cM. The average female to male recombination ratio was 5.43 (range 1-42.8 across all pairwise marker comparisons). We detected 34 linkage groups that corresponded with all chromosome arms mapped with homologous loci in rainbow trout and inferred that 16 represented metacentric chromosomes and 18 represented acrocentric chromosomes. Up to 13 chromosomes were conserved between the two species, suggesting that their structure precedes the divergence between Chinook salmon and rainbow trout. However, marker order differed in one of these linkage groups. The remaining linkage group structures reflected independent Robertsonian chromosomal arrangements, possibly after divergence. The putative linkage group homologies presented here are expected to facilitate future DNA sequencing efforts in Chinook salmon.

Decreased mortality of Lake Michigan Chinook salmon after bacterial kidney disease challenge: evidence for pathogen-driven selection?

In the late 1960s, Chinook salmon Oncorhynchus tshawytscha from the Green River, Washington, were successfully introduced into Lake Michigan. During spring from 1988 to 1992, large fish die-offs affecting Chinook salmon occurred in the lake. Multiple ecological factors probably contributed to the severity of the fish kills, but the only disease agent found regularly was Renibacterium salmoninarum, the causative agent of bacterial kidney disease. In this study, survival after challenge by R. salmoninarum was compared between two Chinook salmon stocks: a Lake Michigan stock from Wisconsin (WI) and the progenitor stock from the Green River. We found that the WI stock had significantly greater survival than the Green River stock. Next, the WI and Green River stocks were exposed to the marine pathogen Listonella anguillarum (formerly Vibrio anguillarum), one of the causative agents of vibriosis; survival after this challenge was significantly poorer for the WI stock than for the Green River stock. A close genetic relationship between the Green River and WI stocks was confirmed by analyzing 13 microsatellite loci. These results collectively suggest that disease susceptibility of Lake Michigan Chinook salmon has diverged from that of the source population, possibly in response to pathogen-driven selection.