Genomic surveillance identifies potential risk factors for SARS-CoV-2 transmission at a mid-sized university in a small rural town.

Understanding transmission dynamics of SARS-CoV-2 in institutions of higher education (IHEs) is important because these settings have potential for rapid viral spread. Here, we used genomic surveillance to retrospectively investigate transmission dynamics throughout the 2020-2021 academic year for the University of Idaho ("University"), a mid-sized IHE in a small rural town. We generated genome assemblies for 1168 SARS-CoV-2 samples collected during the academic year, representing 46.8% of positive samples collected from the University population and 49.8% of positive samples collected from the surrounding community ("Community") at the local hospital during this time. Transmission dynamics differed for the University when compared to the Community, with more infection waves that lasted shorter lengths of time, potentially resulting from high-transmission congregate settings along with mitigation efforts implemented by the University to combat outbreaks. We found evidence for low transmission rates between the University and Community, with approximately 8% of transmissions into the Community originating from the University, and approximately 6% of transmissions into the University originating from the Community. Potential transmission risk factors identified for the University included congregate settings such as sorority and fraternity events and residences, holiday travel, and high caseloads in the surrounding community. Knowledge of these risk factors can help the University and other IHEs develop effective mitigation measures for SARS-CoV-2 and similar pathogens.

High-Resolution Taxonomic Characterization Reveals Novel Human Microbial Strains with Potential as Risk Factors and Probiotics for Prediabetes and Type 2 Diabetes.

Alterations in the composition of the gut microbiota is thought to play a key role in causing type 2 diabetes, yet is not fully understood, especially at the strain level. Here, we used long-read DNA sequencing technology of 16S-ITS-23S rRNA genes for high-resolution characterization of gut microbiota in the development of type 2 diabetes. Gut microbiota composition was characterized from fecal DNA from 47 participants divided into 4 cohorts based on glycemic control: normal glycemic control (healthy; n = 21), reversed prediabetes (prediabetes/healthy; n = 8), prediabetes (n = 8), or type 2 diabetes (n = 10). A total of 46 taxa were found to be possibly related to progression from healthy state to type 2 diabetes. Bacteroides coprophilus DSM 18228, Bifidobacterium pseudocatenulatum DSM 20438, and Bifidobacterium adolescentis ATCC 15703 could confer resistance to glucose intolerance. On the other hand, Odoribacter laneus YIT 12061 may be pathogenic as it was found to be more abundant in type 2 diabetes participants than other cohorts. This research increases our understanding of the structural modulation of gut microbiota in the pathogenesis of type 2 diabetes and highlights gut microbiota strains, with the potential for targeted opportunistic pathogen control or consideration for probiotic prophylaxis and treatment.

Whole genome resequencing identifies local adaptation associated with environmental variation for redband trout.

Aquatic ectotherms are predicted to harbour genomic signals of local adaptation resulting from selective pressures driven by the strong influence of climate conditions on body temperature. We investigated local adaptation in redband trout (Oncorhynchus mykiss gairdneri) using genome scans for 547 samples from 11 populations across a wide range of habitats and thermal gradients in the interior Columbia River. We estimated allele frequencies for millions of single nucleotide polymorphism loci (SNPs) across populations using low-coverage whole genome resequencing, and used population structure outlier analyses to identify genomic regions under divergent selection between populations. Twelve genomic regions showed signatures of local adaptation, including two regions associated with genes known to influence migration and developmental timing in salmonids (GREB1L, ROCK1, SIX6). Genotype-environment association analyses indicated that diurnal temperature variation was a strong driver of local adaptation, with signatures of selection driven primarily by divergence of two populations in the northern extreme of the subspecies range. We also found evidence for adaptive differences between high-elevation desert vs. montane habitats at a smaller geographical scale. Finally, we estimated vulnerability of redband trout to future climate change using ecological niche modelling and genetic offset analyses under two climate change scenarios. These analyses predicted substantial habitat loss and strong genetic shifts necessary for adaptation to future habitats, with the greatest vulnerability predicted for high-elevation desert populations. Our results provide new insight into the complexity of local adaptation in salmonids, and important predictions regarding future responses of redband trout to climate change.

A new mouse SNP genotyping assay for speed congenics: combining flexibility, affordability, and power.

BACKGROUND: Speed congenics is an important tool for creating congenic mice to investigate gene functions, but current SNP genotyping methods for speed congenics are expensive. These methods usually rely on chip or array technologies, and a different assay must be developed for each backcross strain combination. "Next generation" high throughput DNA sequencing technologies have the potential to decrease cost and increase flexibility and power of speed congenics, but thus far have not been utilized for this purpose. RESULTS: We took advantage of the power of high throughput sequencing technologies to develop a cost-effective, high-density SNP genotyping assay that can be used across many combinations of backcross strains. The assay surveys 1640 genome-wide SNPs known to be polymorphic across > 100 mouse strains, with an expected average of 549 ± 136 SD diagnostic SNPs between each pair of strains. We demonstrated that the assay has a high density of diagnostic SNPs for backcrossing the BALB/c strain into the C57BL/6J strain (807-819 SNPs), and a sufficient density of diagnostic SNPs for backcrossing the closely related substrains C57BL/6N and C57BL/6J (123-139 SNPs). Furthermore, the assay can easily be modified to include additional diagnostic SNPs for backcrossing other closely related substrains. We also developed a bioinformatic pipeline for SNP genotyping and calculating the percentage of alleles that match the backcross recipient strain for each sample; this information can be used to guide the selection of individuals for the next backcross, and to assess whether individuals have become congenic. We demonstrated the effectiveness of the assay and bioinformatic pipeline with a backcross experiment of BALB/c-IL4/IL13 into C57BL/6J; after six generations of backcrosses, offspring were up to 99.8% congenic. CONCLUSIONS: The SNP genotyping assay and bioinformatic pipeline developed here present a valuable tool for increasing the power and decreasing the cost of many studies that depend on speed congenics. The assay is highly flexible and can be used for combinations of strains that are commonly used for speed congenics. The assay could also be used for other techniques including QTL mapping, standard F2 crosses, ancestry analysis, and forensics.

Wireworm (Coleoptera: Elateridae) genomic analysis reveals putative cryptic species, population structure, and adaptation to pest control.

The larvae of click beetles (Coleoptera: Elateridae), known as "wireworms," are agricultural pests that pose a substantial economic threat worldwide. We produced one of the first wireworm genome assemblies (Limonius californicus), and investigated population structure and phylogenetic relationships of three species (L. californicus, L. infuscatus, L. canus) across the northwest US and southwest Canada using genome-wide markers (RADseq) and genome skimming. We found two species (L. californicus and L. infuscatus) are comprised of multiple genetically distinct groups that diverged in the Pleistocene but have no known distinguishing morphological characters, and therefore could be considered cryptic species complexes. We also found within-species population structure across relatively short geographic distances. Genome scans for selection provided preliminary evidence for signatures of adaptation associated with different pesticide treatments in an agricultural field trial for L. canus. We demonstrate that genomic tools can be a strong asset in developing effective wireworm control strategies.

Regeneration associated transcriptional signature of retinal microglia and macrophages.

Zebrafish have the remarkable capacity to regenerate retinal neurons following a variety of damage paradigms. Following initial tissue insult and a period of cell death, a proliferative phase ensues that generates neuronal progenitors, which ultimately regenerate damaged neurons. Recent work has revealed that Müller glia are the source of regenerated neurons in zebrafish. However, the roles of another important class of glia present in the retina, microglia, during this regenerative phase remain elusive. Here, we examine retinal tissue and perform QuantSeq. 3'mRNA sequencing/transcriptome analysis to reveal localization and putative functions, respectively, of mpeg1 expressing cells (microglia/macrophages) during Müller glia-mediated regeneration, corresponding to a time of progenitor proliferation and production of new neurons. Our results indicate that in this regenerative state, mpeg1-expressing cells are located in regions containing regenerative Müller glia and are likely engaged in active vesicle trafficking. Further, mpeg1+ cells congregate at and around the optic nerve head. Our transcriptome analysis reveals several novel genes not previously described in microglia. This dataset represents the first report, to our knowledge, to use RNA sequencing to probe the microglial transcriptome in such context, and therefore provides a resource towards understanding microglia/macrophage function during successful retinal (and central nervous tissue) regeneration.

A Rapid Method for Sequencing Double-Stranded RNAs Purified from Yeasts and the Identification of a Potent K1 Killer Toxin Isolated from Saccharomyces cerevisiae.

Mycoviruses infect a large number of diverse fungal species, but considering their prevalence, relatively few high-quality genome sequences have been determined. Many mycoviruses have linear double-stranded RNA genomes, which makes it technically challenging to ascertain their nucleotide sequence using conventional sequencing methods. Different specialist methodologies have been developed for the extraction of double-stranded RNAs from fungi and the subsequent synthesis of cDNAs for cloning and sequencing. However, these methods are often labor-intensive, time-consuming, and can require several days to produce cDNAs from double-stranded RNAs. Here, we describe a comprehensive method for the rapid extraction and sequencing of dsRNAs derived from yeasts, using short-read next generation sequencing. This method optimizes the extraction of high-quality double-stranded RNAs from yeasts and 3' polyadenylation for the initiation of cDNA synthesis for next-generation sequencing. We have used this method to determine the sequence of two mycoviruses and a double-stranded RNA satellite present within a single strain of the model yeast Saccharomyces cerevisiae. The quality and depth of coverage was sufficient to detect fixed and polymorphic mutations within viral populations extracted from a clonal yeast population. This method was also able to identify two fixed mutations within the alpha-domain of a variant K1 killer toxin encoded on a satellite double-stranded RNA. Relative to the canonical K1 toxin, these newly reported mutations increased the cytotoxicity of the K1 toxin against a specific species of yeast.

Bacterial Diversity in Replicated Hydrogen Sulfide-Rich Streams.

Extreme environments typically require costly adaptations for survival, an attribute that often translates to an elevated influence of habitat conditions on biotic communities. Microbes, primarily bacteria, are successful colonizers of extreme environments worldwide, yet in many instances, the interplay between harsh conditions, dispersal, and microbial biogeography remains unclear. This lack of clarity is particularly true for habitats where extreme temperature is not the overarching stressor, highlighting a need for studies that focus on the role other primary stressors (e.g., toxicants) play in shaping biogeographic patterns. In this study, we leveraged a naturally paired stream system in southern Mexico to explore how elevated hydrogen sulfide (H2S) influences microbial diversity. We sequenced a portion of the 16S rRNA gene using bacterial primers for water sampled from three geographically proximate pairings of streams with high (> 20 µM) or low (~ 0 µM) H2S concentrations. After exploring bacterial diversity within and among sites, we compared our results to a previous study of macroinvertebrates and fish for the same sites. By spanning multiple organismal groups, we were able to illuminate how H2S may differentially affect biodiversity. The presence of elevated H2S had no effect on overall bacterial diversity (p = 0.21), a large effect on community composition (25.8% of variation explained, p < 0.0001), and variable influence depending upon the group-whether fish, macroinvertebrates, or bacteria-being considered. For bacterial diversity, we recovered nine abundant operational taxonomic units (OTUs) that comprised a core H2S-rich stream microbiome in the region. Many H2S-associated OTUs were members of the Epsilonproteobacteria and Gammaproteobacteria, which both have been implicated in endosymbiotic relationships between sulfur-oxidizing bacteria and eukaryotes, suggesting the potential for symbioses that remain to be discovered in these habitats.

Complete Genome Sequence of Broad-Host-Range Shiga Toxin-Converting Bacteriophage SH2026Stx1, Isolated from Escherichia coli O157:H7.

The Shiga toxin-encoding phage SH2026Stx1 was isolated from Escherichia coli O157:H7 strain 2026. SH2026Stx1 and its detoxified derivative can infect a broad range of E. coli strains, including commensal, enteropathogenic, and enteroaggregative strains. We report here the complete genome sequence of phage SH2026Stx1 and its important features.

Complete Mitochondrial Genome Sequence of Bighorn Sheep.

We report here the complete mitochondrial genome sequence of a Rocky Mountain bighorn sheep (Ovis canadensis) in the United States. The circular genome has a size of 16,466 bp and contains 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes.

High-Quality Complete Genome Sequences of Three Bovine Shiga Toxin-Producing Escherichia coli O177:H- (fliCH25) Isolates Harboring Virulent stx2 and Multiple Plasmids.

Shiga toxin-producing Escherichia coli (STEC) bacteria are zoonotic pathogens. We report here the high-quality complete genome sequences of three STEC O177:H- (fliCH25) strains, SMN152SH1, SMN013SH2, and SMN197SH3. The assembled genomes consisted of one optical map-verified circular chromosome for each strain, plus two plasmids for SMN013SH2 and three plasmids for SMN152SH1 and SMN197SH3, respectively.

Expression of bovine non-classical major histocompatibility complex class I proteins in mouse P815 and human K562 cells.

Major histocompatibility complex class I (MHC-I) proteins can be expressed as cell surface or secreted proteins. To investigate whether bovine non-classical MHC-I proteins are expressed as cell surface or secreted proteins, and to assess the reactivity pattern of monoclonal antibodies with non-classical MHC-I isoforms, we expressed the MHC proteins in murine P815 and human K562 (MHC-I deficient) cells. Following antibiotic selection, stably transfected cell lines were stained with H1A or W6/32 antibodies to detect expression of the MHC-I proteins by flow cytometry. Two non-classical proteins (BoLA-NC1*00501 and BoLA-NC3*00101) were expressed on the cell surface in both cell lines. Surprisingly, the BoLA-NC4*00201 protein was expressed on the cell membrane of human K562 but not mouse P815 cells. Two non-classical proteins (BoLA-NC1*00401, which lacks a transmembrane domain, and BoLA-NC2*00102) did not exhibit cell surface expression. Nevertheless, Western blot analyses demonstrated expression of the MHC-I heavy chain in all transfected cell lines. Ammonium-sulfate precipitation of proteins from culture supernatants showed that BoLA-NC1*00401 was secreted and that all surface expressed proteins where shed from the cell membrane by the transfected cells. Interestingly, the surface expressed MHC-I proteins were present in culture supernatants at a much higher concentration than BoLA-NC1*00401. This comprehensive study shows that bovine non-classical MHC-I proteins BoLA-NC1*00501, BoLA-NC3*00101, and BoLA-NC4*00201 are expressed as surface isoforms with the latter reaching the cell membrane only in K562 cells. Furthermore, it demonstrated that BoLA-NC1*00401 is a secreted isoform and that significant quantities of membrane associated MHC-I proteins can be shed from the cell membrane.

Mitochondrial Genome Sequence of the Galápagos Endemic Land Snail Naesiotus nux.

We report herein the draft mitochondrial genome sequence of Naesiotus nux, a Galápagos endemic land snail species of the genus Naesiotus. The circular genome is 15 kb and encodes 13 protein-coding genes, 2 rRNA genes, and 21 tRNA genes.

Full Mitochondrial Genome Sequence of the Sugar Beet Wireworm Limonius californicus (Coleoptera: Elateridae), a Common Agricultural Pest.

We report here the full mitochondrial genome sequence of Limonius californicus, a species of click beetle that is an agricultural pest in its larval form. The circular genome is 16.5 kb and contains 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes.