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.

Contagious Antibiotic Resistance: Plasmid Transfer among Bacterial Residents of the Zebrafish Gut.

By characterizing the trajectories of antibiotic resistance gene transfer in bacterial communities such as the gut microbiome, we will better understand the factors that influence this spread of resistance. Our aim was to investigate the host network of a multidrug resistance broad-host-range plasmid in the culturable gut microbiome of zebrafish. This was done through in vitro and in vivo conjugation experiments with Escherichia coli as the donor of the plasmid pB10::gfp When this donor was mixed with the extracted gut microbiome, only transconjugants of Aeromonas veronii were detected. In separate matings between the same donor and four prominent isolates from the gut microbiome, the plasmid transferred to two of these four isolates, A. veronii and Plesiomonas shigelloides, but not to Shewanella putrefaciens and Vibrio mimicus When these A. veronii and P. shigelloides transconjugants were the donors in matings with the same four isolates, the plasmid now also transferred from A. veronii to S. putrefaciensP. shigelloides was unable to donate the plasmid, and V. mimicus was unable to acquire it. Finally, when the E. coli donor was added in vivo to zebrafish through their food, plasmid transfer was observed in the gut, but only to Achromobacter, a rare member of the gut microbiome. This work shows that the success of plasmid-mediated antibiotic resistance spread in a gut microbiome depends on the donor-recipient species combinations and therefore their spatial arrangement. It also suggests that rare gut microbiome members should not be ignored as potential reservoirs of multidrug resistance plasmids from food.IMPORTANCE To understand how antibiotic resistance plasmids end up in human pathogens, it is crucial to learn how, where, and when they are transferred and maintained in members of bacterial communities such as the gut microbiome. To gain insight into the network of plasmid-mediated antibiotic resistance sharing in the gut microbiome, we investigated the transferability and maintenance of a multidrug resistance plasmid among the culturable bacteria of the zebrafish gut. We show that the success of plasmid-mediated antibiotic resistance spread in a gut microbiome can depend on which species are involved, as some are important nodes in the plasmid-host network and others are dead ends. Our findings also suggest that rare gut microbiome members should not be ignored as potential reservoirs of multidrug resistance plasmids from food.

Characterization and Evolution of the Spotted Gar Retina.

In this study, we characterize the retina of the spotted gar, Lepisosteus oculatus, a ray-finned fish. Gar did not undergo the whole genome duplication event that occurred at the base of the teleost fish lineage, which includes the model species zebrafish and medaka. The divergence of gars from the teleost lineage and the availability of a high-quality genome sequence make it a uniquely useful species to understand how genome duplication sculpted features of the teleost visual system, including photoreceptor diversity. We developed reagents to characterize the cellular organization of the spotted gar retina, including representative markers for all major classes of retinal neurons and Müller glia. We report that the gar has a preponderance of predicted short-wavelength shifted (SWS) opsin genes, including a duplicated set of SWS1 (ultraviolet) sensitive opsin encoding genes, a SWS2 (blue) opsin encoding gene, and two rod opsin encoding genes, all of which were expressed in retinal photoreceptors. We also report that gar SWS1 cones lack the geometric organization of photoreceptors observed in teleost fish species, consistent with the crystalline photoreceptor mosaic being a teleost innovation. Of note the spotted gar expresses both exo-rhodopsin (RH1-1) and rhodopsin (RH1-2) in rods. Exo-rhodopsin is an opsin that is not expressed in the retina of zebrafish and other teleosts, but rather is expressed in regions of the brain. This study suggests that exo-rhodopsin is an ancestral actinopterygian (ray finned fish) retinal opsin, and in teleosts its expression has possibly been subfunctionalized to the pineal gland.

Anxiolytic effects of fluoxetine and nicotine exposure on exploratory behavior in zebrafish.

Zebrafish (Danio rerio) have emerged as a popular model for studying the pharmacology and behavior of anxiety. While there have been numerous studies documenting the anxiolytic and anxiogenic effects of common drugs in zebrafish, many do not report or test for behavioral differences between the sexes. Previous studies have indicated that males and females differ in their baseline level of anxiety. In this study, we test for a sex interaction with fluoxetine and nicotine. We exposed fish to system water (control), 10 mg/L fluoxetine, or 1 mg/L nicotine for three minutes prior to being subjected to four minutes in an open-field drop test. Video recordings were tracked using ProAnalyst. Fish from both drug treatments reduced swimming speed, increased vertical position, and increased use of the top half of the open field when compared with the control, though fluoxetine had a larger effect on depth related behaviors while nicotine mostly affected swimming speed. A significant sex effect was observed where females swam at a slower and more constant speed than males, however neither drug produced a sex-dependent response.

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.

Oncorhynchus mykiss pax7 sequence variations with comparative analyses against other teleost species.

The paired box-7 (pax7) transcription factor expressed in satellite cells (SCs) is an essential regulator of skeletal muscle growth and regeneration in vertebrates including fish. Characterization of rainbow trout (Oncorhynchus mykiss) pax7 gene/s may offer novel insights into skeletal myogenesis by SCs in this indeterminate growth species. Further, evaluation of promoters for cis-regulatory regions may shed light on the evolutionary fate of the duplicated genes. Employing standard PCR, cloning and computational approach, we identified and report complete coding sequences of two pax7 paralogs of rainbow trout (rt); rtpax7α and rtpax7ß. Both genes show significant identity in the nucleotide (97%) and the predicted amino acid (98%) sequences, and bear the characteristic paired domain (PD), octapeptide (OP) and homeodomain (HD) motifs. We further report several splice variants of each gene and nucleotide differences in coding sequence that predicts six putative amino acid changes between the two genes. Additionally, we noted a trinucleotide deletion in rtpax7ß that results in putative serine elimination at the N-terminus and a single nucleotide polymorphism (SNP) in majority of the rtpax7ß variants (6/10) that predicts an arginine substitution for a lysine. We also deciphered the genomic organization up to the first three exons and the upstream putative promoter regions of both genes. Comparative in silico analysis of both the trout pax7 promoters with that of zebrafish pax7 duplicates; zfpax7a and zfpax7b; predicts several important cis-elements/transcription factor binding sites (TFBS) in these teleost pax7 promoter regions.

Maternal bisphenol a exposure impacts the fetal heart transcriptome.

Conditions during fetal development influence health and disease in adulthood, especially during critical windows of organogenesis. Fetal exposure to the endocrine disrupting chemical, bisphenol A (BPA) affects the development of multiple organ systems in rodents and monkeys. However, effects of BPA exposure on cardiac development have not been assessed. With evidence that maternal BPA is transplacentally delivered to the developing fetus, it becomes imperative to examine the physiological consequences of gestational exposure during primate development. Herein, we evaluate the effects of daily, oral BPA exposure of pregnant rhesus monkeys (Macaca mulatta) on the fetal heart transcriptome. Pregnant monkeys were given daily oral doses (400 µg/kg body weight) of BPA during early (50-100 ± 2 days post conception, dpc) or late (100 ± 2 dpc--term), gestation. At the end of treatment, fetal heart tissues were collected and chamber specific transcriptome expression was assessed using genome-wide microarray. Quantitative real-time PCR was conducted on select genes and ventricular tissue glycogen content was quantified. Our results show that BPA exposure alters transcription of genes that are recognized for their role in cardiac pathophysiologies. Importantly, myosin heavy chain, cardiac isoform alpha (Myh6) was down-regulated in the left ventricle, and 'A Disintegrin and Metalloprotease 12', long isoform (Adam12-l) was up-regulated in both ventricles, and the right atrium of the heart in BPA exposed fetuses. BPA induced alteration of these genes supports the hypothesis that exposure to BPA during fetal development may impact cardiovascular fitness. Our results intensify concerns about the role of BPA in the genesis of human metabolic and cardiovascular diseases.

Retinal regeneration is facilitated by the presence of surviving neurons.

Teleost fish regenerate their retinas after damage, in contrast to mammals. In zebrafish subjected to an extensive ouabain-induced lesion that destroys all neurons and spares Müller glia, functional recovery and restoration of normal optic nerve head (ONH) diameter take place at 100 days postinjury. Subsequently, regenerated retinas overproduce cells in the retinal ganglion cell (RGC) layer, and the ONH becomes enlarged. Here, we test the hypothesis that a selective injury, which spares photoreceptors and Müller glia, results in faster functional recovery and fewer long-term histological abnormalities. Following this selective retinal damage, recovery of visual function required 60 days, consistent with this hypothesis. In contrast to extensively damaged retinas, selectively damaged retinas showed fewer histological errors and did not overproduce neurons. Extensively damaged retinas had RGC axons that were delayed in pathfinding to the ONH, and showed misrouted axons within the ONH, suggesting that delayed functional recovery following an extensive lesion is related to defects in RGC axons exiting the eye and/or reaching their central targets. The atoh7, fgf8a, Sonic hedgehog (shha), and netrin-1 genes were differentially expressed, and the distribution of hedgehog protein was disrupted after extensive damage as compared with selective damage. Confirming a role for Shh signaling in supporting rapid regeneration, shha(t4) +/- zebrafish showed delayed functional recovery after selective damage. We suggest that surviving retinal neurons provide structural/molecular information to regenerating neurons, and that this patterning mechanism regulates factors such as Shh. These factors in turn control neuronal number, retinal lamination, and RGC axon pathfinding during retinal regeneration.

The quantitative genetic architecture of the bold-shy continuum in zebrafish, Danio rerio.

In studies of consistent individual differences (personality) along the bold-shy continuum, a pattern of behavioral correlations frequently emerges: individuals towards the bold end of the continuum are more likely to utilize risky habitat, approach potential predators, and feed under risky conditions. Here, we address the hypothesis that observed phenotypic correlations among component behaviors of the bold-shy continuum are a result of underlying genetic correlations (quantitative genetic architecture). We used a replicated three-generation pedigree of zebrafish (Danio rerio) to study three putative components of the bold-shy continuum: horizontal position, swim level, and feeding latency. We detected significant narrow-sense heritabilities as well as significant genetic and phenotypic correlations among all three behaviors, such that fish selected for swimming at the front of the tank swam closer to the observer, swam higher in the water column, and fed more quickly than fish selected for swimming at the back of the tank. Further, the lines varied in their initial open field behavior (swim level and activity level). The quantitative genetic architecture of the bold-shy continuum indicates that the multivariate behavioral phenotype characteristic of a "bold" personality type may be a result of correlated evolution via underlying genetic correlations.

Sex-specific transcriptional responses of the zebrafish (Danio rerio) brain selenoproteome to acute sodium selenite supplementation.

The potential benefits of selenium (Se) supplementation are currently under investigation for prevention of certain cancers and treatment of neurological disorders. However, little is known concerning the response of the brain to increased dietary Se under conditions of Se sufficiency, despite the majority of Se supplementation trials occurring in healthy, Se sufficient subjects. We evaluated the transcriptional response of Se-dependent genes, selenoproteins and the genes necessary for their synthesis (the selenoproteome), in the zebrafish (Danio rerio) brain to supplementation with nutritionally relevant levels of dietary Se (sodium selenite) during conditions of assumed Se sufficiency. We first used a microarray approach to analyze the response of the brain selenoproteome to dietary Se supplementation for 14 days and then assessed the immediacy and time-scale transcriptional response of the brain selenoproteome to 1, 7, and 14 days of Se supplementation by quantitative real-time PCR (qRT-PCR). The microarray approach did not indicate large-scale influences of Se on the brain transcriptome as a whole or the selenoproteome specifically; only one nonselenoproteome gene (si:ch73-44m9.2) was significantly differentially expressed. Our qRT-PCR results, however, indicate that increases of dietary Se cause small, but significant transcriptional changes within the brain selenoproteome, even after only 1 day of supplementation. These responses were dynamic over a short period of supplementation in a manner highly dependent on sex and the duration of Se supplementation. In nutritional intervention studies, it may be necessary to utilize methods such as qRT-PCR, which allow larger sample sizes, for detecting subtle transcriptional changes in the brain.

Is behavioral variation along the bold-shy continuum associated with variation in the stress axis in zebrafish?

We tested whether boldness is associated with attenuation of the physiological stress response in behaviorally selected lines of zebrafish Danio rerio. We measured three component behaviors of boldness: cortisol levels under control and stressed conditions, growth rate, and expression of key genes linked to the hypothalamic-pituitary-interrenal axis in the brain. Surprisingly, bold animals did not differ from shy animals with respect to cortisol levels. However, significant differences between these animals in the expression of glucocorticoid receptors and genes that regulate production of stress hormones indicate that there may still be a relationship between bold behavior and the stress axis. Perhaps the most surprising result of this study was the degree of sexual dimorphism: female zebrafish were bolder than male zebrafish, had significantly lower levels of cortisol, and differed significantly in the expression of several genes in the brain. Our data indicate that a bold behavioral type is associated with transcriptional attenuation of stress axis genes, but we do not yet know whether evolution along the bold-shy continuum is attributable to genetic changes in the stress axis. The bold and shy zebrafish lines will be valuable tools for additional research into the relationship between stress and behavior and the mechanisms regulating sexual dimorphism in these traits.

An improved algorithm for the detection of genomic variation using short oligonucleotide expression microarrays.

High-throughput microarray experiments often generate far more biological information than is required to test the experimental hypotheses. Many microarray analyses are considered finished after differential expression and additional analyses are typically not performed, leaving untapped biological information left undiscovered. This is especially true if the microarray experiment is from an ecological study of multiple populations. Comparisons across populations may also contain important genomic polymorphisms, and a subset of these polymorphisms may be identified with microarrays using techniques for the detection of single feature polymorphisms (SFP). SFPs are differences in microarray probe level intensities caused by genetic polymorphisms such as single-nucleotide polymorphisms and small insertions/deletions and not expression differences. In this study, we provide a new algorithm for the detection of SFPs, evaluate the algorithm using existing data from two publicly available Affymetrix Barley (Hordeum vulgare) microarray data sets and compare them to two previously published SFP detection algorithms. Results show that our algorithm provides more consistent and sensitive calling of SFPs with a lower false discovery rate. Simultaneous analysis of SFPs and differential expression is a low-cost method for the enhanced analysis of microarray data, enabling additional biological inferences to be made.

Brain transcriptome variation among behaviorally distinct strains of zebrafish (Danio rerio).

BACKGROUND: Domesticated animal populations often show profound reductions in predator avoidance and fear-related behavior compared to wild populations. These reductions are remarkably consistent and have been observed in a diverse array of taxa including fish, birds, and mammals. Experiments conducted in common environments indicate that these behavioral differences have a genetic basis. In this study, we quantified differences in fear-related behavior between wild and domesticated zebrafish strains and used microarray analysis to identify genes that may be associated with this variation. RESULTS: Compared to wild zebrafish, domesticated zebrafish spent more time near the water surface and were more likely to occupy the front of the aquarium nearest a human observer. Microarray analysis of the brain transcriptome identified high levels of population variation in gene expression, with 1,749 genes significantly differentially expressed among populations. Genes that varied among populations belonged to functional categories that included DNA repair, DNA photolyase activity, response to light stimulus, neuron development and axon guidance, cell death, iron-binding, chromatin reorganization, and homeobox genes. Comparatively fewer genes (112) differed between domesticated and wild strains with notable genes including gpr177 (wntless), selenoprotein P1a, synaptophysin and synaptoporin, and acyl-CoA binding domain containing proteins (acbd3 and acbd4). CONCLUSIONS: Microarray analysis identified a large number of genes that differed among zebrafish populations and may underlie behavioral domestication. Comparisons with similar microarray studies of domestication in rainbow trout and canids identified sixteen evolutionarily or functionally related genes that may represent components of shared molecular mechanisms underlying convergent behavioral evolution during vertebrate domestication. However, this conclusion must be tempered by limitations associated with comparisons among microarray studies and the low level of population-level replication inherent to these studies.

Retinal proliferation response in the buphthalmic zebrafish, bugeye.

The zebrafish retina regenerates in response to acute retinal lesions, replacing damaged neurons with new neurons. In this study we test the hypothesis that chronic stress to inner retinal neurons also triggers a retinal regeneration response in the bugeye zebrafish. Mutations in the lrp2 gene in zebrafish are associated with a progressive eye phenotype (bugeye) that models several risk factors for human glaucoma including buphthalmos (enlarged eyes), elevated intraocular pressure (IOP), and upregulation of genes related to retinal ganglion cell pathology. The retinas of adult bugeye zebrafish showed high rates of ongoing proliferation which resulted in the production of a small number of new retinal neurons, particularly photoreceptors. A marker of mechanical cell stress, Hsp27, was strongly expressed in inner retinal neurons and glia of bugeye retinas. The more enlarged eyes of individual bugeye zebrafish showed disrupted retinal lamination, and a persistent reduced density of neurons in the ganglion cell layer (GCL), although total numbers of GCL neurons were higher than in control eyes. Despite the presence of a proliferative response to damage, the adult bugeye zebrafish remained behaviorally blind. These findings suggest the existence of an unsuccessful regenerative response to a persistent pathological condition in the bugeye zebrafish.

Effect of aging on male reproduction in zebrafish (Danio rerio).

The study was designed to test the hypothesis that male aging is associated with a change in reproductive function in the zebrafish. Young (290 ± 37 d) and older (911 ± 48 d) males were combined with females (604 ± 24 d) to test the effect of male age on the number and fertility of eggs laid by their mates. 48% of breeding trials with young males and 25% of the trails with older males resulted in egg deposition. Although young males were associated with significantly more successful breeding attempts than older males, number of eggs laid per clutch, number and percent of fertilized eggs and the number and percent living embryos were not statistically different between young and older males. These data suggest that male aging is associated with altered reproductive behavior and/or female response but not in sperm quality per se. Consistent with this interpretation were the findings that percent motility and sperm motility characteristics did not differ between sperm from young and older males as assessed by computer-assisted sperm analysis. However, older males contained higher quantities of extractable sperm than did young males, perhaps associated with fewer successful breeding attempts. Age-related effects on male reproductive in the zebrafish may therefore be a consequence of behavioral or morphological features that play a role in female mate choice and/ or male sexual response.

Zebrafish (Danio rerio) vary by strain and sex in their behavioral and transcriptional responses to selenium supplementation.

We used the Nadia, Gaighatta, Scientific Hatcheries, and TM1 zebrafish (Danio rerio) strains to test the hypothesis that variation among populations influences the behavioral and transcriptional responses to selenium supplementation. When fed a diet with control levels of selenium, zebrafish strains differed significantly in behavior, characterized as their mean horizontal and vertical swimming positions within the tank. The four strains also differed in brain expression of selenoprotein P1a (sepp1a), glutathione peroxidase 3 (gpx3), thioredoxin reductase 1 (txnrd1), and tRNA selenocysteine associated protein 1 (secp43). Iodothyronine deiodinase 2 (dio2) did not differ among strains but showed a sex-specific expression pattern. When supplemented with selenium, all strains spent a greater proportion of time near the front of the tank, but the response of vertical swimming depth varied by strain. Selenium supplementation also caused changes in selenoprotein expression in the brain that varied by strain for sepp1a, secp43, and dio2, and varied by strain and sex for txnrd1. Expression of gpx3 was unaffected by selenium. Our data indicate that selenium homeostasis in the brain may be a regulator of behavior in zebrafish, and the strain-specific effects of selenium supplementation suggest that genetic heterogeneity among populations can influence the results of selenium supplementation studies.

Dietary carbohydrate level affects transcription factor expression that regulates skeletal muscle myogenesis in rainbow trout.

Understanding the effects of dietary carbohydrates on transcription factors that regulate myogenesis provides insight into the role of nutrient sensing by satellite cells towards myocyte differentiation. We evaluated the influence of dietary carbohydrate level (0, 15, 25 or 35%) on the temporal mRNA expression patterns (4, 8 or 12 weeks) of transcription factors that regulate satellite cell myocyte addition (MA) in rainbow trout (Oncorhynchus mykiss), a vertebrate with indeterminate growth. Relative to the 0% carbohydrate (NC) diet, 15 (IC-15) and 25% (IC-25) carbohydrate containing diets significantly up-regulate MyoD and Myf5, but not Pax7, after 12 weeks of feeding. Simultaneously, the Pax7/MyoD mRNA expression ratio declined significantly with both the IC diets. Myogenin mRNA expression also increased in rainbow trout (RBT) fed the IC-15 diet. The high carbohydrate (HC) diet (35%) attenuated the increased mRNA expression of these transcription factors. It is of note that the 4 and 8 week samples lacked the promyogenic expression patterns. The myogenic gene expression in fish fed the IC-15 diet for 12 weeks indicate a transcriptional signature that reflects increased satellite cell myogenesis. Our results suggest a potential role for satellite cells in the nutrient sensing ability of a vertebrate with indeterminate skeletal muscle growth.

Effect of starvation on transcriptomes of brain and liver in adult female zebrafish (Danio rerio).

We used microarray and quantitative real-time PCR (qRT-PCR) analyses in adult female zebrafish (Danio rerio) to identify metabolic pathways regulated by starvation in the liver and brain. The transcriptome of whole zebrafish brain showed little response to 21 days of starvation. Only agouti-related protein 1 (agrp1) significantly responded, with increased expression in brains of starved fish. In contrast, a 21-day period of starvation significantly downregulated 466 and upregulated 108 transcripts in the liver, indicating an overall decrease in metabolic activity, reduced lipid metabolism, protein biosynthesis, proteolysis, and cellular respiration, and increased gluconeogenesis. Starvation also regulated expression of many components of the unfolded protein response, the first such report in a species other than yeast (Saccharomyces cerevisiae) and mice (Mus musculus). The response of the zebrafish hepatic transcriptome to starvation was strikingly similar to that of rainbow trout (Oncorhynchus mykiss) and less similar to mouse, while the response of common carp (Cyprinus carpio) differed considerably from the other three species.

Sexual dimorphism in hepatic gene expression and the response to dietary carbohydrate manipulation in the zebrafish (Danio rerio).

In this study, we tested for the presence of sexual dimorphism in the hepatic transcriptome of the adult zebrafish and examined the effect of long term manipulation of dietary carbohydrate on gene expression in both sexes. Zebrafish were fed diets comprised of 0%, 15%, 25%, or 35% carbohydrate from the larval stage through sexual maturity, then sampled for hepatic tissue, growth, proximate body composition, and retention efficiencies. Using Affymetrix microarrays and qRT-PCR, we observed substantial sexual dimorphism in the hepatic transcriptome. Males up-regulated genes associated with oxidative metabolism, carbohydrate metabolism, energy production, and amelioration of oxidative stress, while females had higher expression levels of genes associated with translation. Restriction of dietary carbohydrate (0% diet) significantly affected hepatic gene expression, growth performance, retention efficiencies of protein and energy, and percentages of moisture, lipid, and ash. The response of some genes to dietary manipulation varied by sex; with increased dietary carbohydrate, males up-regulated genes associated with oxidative metabolism (e.g. hadhbeta) while females up-regulated genes associated with glucose phosphorylation (e.g. glucokinase). Our data support the use of the zebrafish model for the study of fish nutritional genomics, but highlight the importance of accounting for sexual dimorphism in these studies.