High-volume hybridoma sequencing on the NeuroMabSeq platform enables efficient generation of recombinant monoclonal antibodies and scFvs for neuroscience research.

The Neuroscience Monoclonal Antibody Sequencing Initiative (NeuroMabSeq) is a concerted effort to determine and make publicly available hybridoma-derived sequences of monoclonal antibodies (mAbs) valuable to neuroscience research. Over 30 years of research and development efforts including those at the UC Davis/NIH NeuroMab Facility have resulted in the generation of a large collection of mouse mAbs validated for neuroscience research. To enhance dissemination and increase the utility of this valuable resource, we applied a high-throughput DNA sequencing approach to determine immunoglobulin heavy and light chain variable domain sequences from source hybridoma cells. The resultant set of sequences was made publicly available as a searchable DNA sequence database (neuromabseq.ucdavis.edu) for sharing, analysis and use in downstream applications. We enhanced the utility, transparency, and reproducibility of the existing mAb collection by using these sequences to develop recombinant mAbs. This enabled their subsequent engineering into alternate forms with distinct utility, including alternate modes of detection in multiplexed labeling, and as miniaturized single chain variable fragments or scFvs. The NeuroMabSeq website and database and the corresponding recombinant antibody collection together serve as a public DNA sequence repository of mouse mAb heavy and light chain variable domain sequences and as an open resource for enhancing dissemination and utility of this valuable collection of validated mAbs.

Long wavelength-sensing cones of zebrafish retina exhibit multiple layers of transcriptional heterogeneity.

Introduction: Understanding how photoreceptor genes are regulated is important for investigating retinal development and disease. While much is known about gene regulation in cones, the mechanism by which tandemly-replicated opsins, such as human long wavelength-sensitive and middle wavelength-sensitive opsins, are differentially regulated remains elusive. In this study, we aimed to further our understanding of transcriptional heterogeneity in cones that express tandemly-replicated opsins and the regulation of such differential expression using zebrafish, which express the tandemly-replicated opsins lws1 and lws2. Methods: We performed bulk and single cell RNA-Seq of LWS1 and LWS2 cones, evaluated expression patterns of selected genes of interest using multiplex fluorescence in situ hybridization, and used exogenous thyroid hormone (TH) treatments to test selected genes for potential control by thyroid hormone: a potent, endogenous regulator of lws1 and lws2 expression. Results: Our studies indicate that additional transcriptional differences beyond opsin expression exist between LWS1 and LWS2 cones. Bulk RNA-Seq results showed 95 transcripts enriched in LWS1 cones and 186 transcripts enriched in LWS2 cones (FC > 2, FDR < 0.05). In situ hybridization results also reveal underlying heterogeneity within the lws1- and lws2-expressing populations. This heterogeneity is evident in cones of mature zebrafish, and further heterogeneity is revealed in transcriptional responses to TH treatments. Discussion: We found some evidence of coordinate regulation of lws opsins and other genes by exogenous TH in LWS1 vs. LWS2 cones, as well as evidence of gene regulation not mediated by TH. The transcriptional differences between LWS1 and LWS2 cones are likely controlled by multiple signals, including TH.

NeuroMabSeq: high volume acquisition, processing, and curation of hybridoma sequences and their use in generating recombinant monoclonal antibodies and scFvs for neuroscience research.

The Neuroscience Monoclonal Antibody Sequencing Initiative (NeuroMabSeq) is a concerted effort to determine and make publicly available hybridoma-derived sequences of monoclonal antibodies (mAbs) valuable to neuroscience research. Over 30 years of research and development efforts including those at the UC Davis/NIH NeuroMab Facility have resulted in the generation of a large collection of mouse mAbs validated for neuroscience research. To enhance dissemination and increase the utility of this valuable resource, we applied a high-throughput DNA sequencing approach to determine immunoglobulin heavy and light chain variable domain sequences from source hybridoma cells. The resultant set of sequences was made publicly available as searchable DNA sequence database ( neuromabseq.ucdavis.edu ) for sharing, analysis and use in downstream applications. We enhanced the utility, transparency, and reproducibility of the existing mAb collection by using these sequences to develop recombinant mAbs. This enabled their subsequent engineering into alternate forms with distinct utility, including alternate modes of detection in multiplexed labeling, and as miniaturized single chain variable fragments or scFvs. The NeuroMabSeq website and database and the corresponding recombinant antibody collection together serve as a public DNA sequence repository of mouse mAb heavy and light chain variable domain sequences and as an open resource for enhancing dissemination and utility of this valuable collection of validated mAbs.

Sex disparate gut microbiome and metabolome perturbations precede disease progression in a mouse model of Rett syndrome.

Rett syndrome (RTT) is a regressive neurodevelopmental disorder in girls, characterized by multisystem complications including gut dysbiosis and altered metabolism. While RTT is known to be caused by mutations in the X-linked gene MECP2, the intermediate molecular pathways of progressive disease phenotypes are unknown. Mecp2 deficient rodents used to model RTT pathophysiology in most prior studies have been male. Thus, we utilized a patient-relevant mouse model of RTT to longitudinally profile the gut microbiome and metabolome across disease progression in both sexes. Fecal metabolites were altered in Mecp2e1 mutant females before onset of neuromotor phenotypes and correlated with lipid deficiencies in brain, results not observed in males. Females also displayed altered gut microbial communities and an inflammatory profile that were more consistent with RTT patients than males. These findings identify new molecular pathways of RTT disease progression and demonstrate the relevance of further study in female Mecp2 animal models.

Establishment of quantitative RNAi-based forward genetics in Entamoeba histolytica and identification of genes required for growth.

While Entamoeba histolytica remains a globally important pathogen, it is dramatically understudied. The tractability of E. histolytica has historically been limited, which is largely due to challenging features of its genome. To enable forward genetics, we constructed and validated the first genome-wide E. histolytica RNAi knockdown mutant library. This library allows for Illumina deep sequencing analysis for quantitative identification of mutants that are enriched or depleted after selection. We developed a novel analysis pipeline to precisely define and quantify gene fragments. We used the library to perform the first RNAi screen in E. histolytica and identified slow growth (SG) mutants. Among genes targeted in SG mutants, many had annotated functions consistent with roles in cellular growth or metabolic pathways. Some targeted genes were annotated as hypothetical or lacked annotated domains, supporting the power of forward genetics in uncovering functional information that cannot be gleaned from databases. While the localization of neither of the proteins targeted in SG1 nor SG2 mutants could be predicted by sequence analysis, we showed experimentally that SG1 localized to the cytoplasm and cell surface, while SG2 localized to the cytoplasm. Overexpression of SG1 led to increased growth, while expression of a truncation mutant did not lead to increased growth, and thus aided in defining functional domains in this protein. Finally, in addition to establishing forward genetics, we uncovered new details of the unusual E. histolytica RNAi pathway. These studies dramatically improve the tractability of E. histolytica and open up the possibility of applying genetics to improve understanding of this important pathogen.

Biosolids as a Source of Antibiotic Resistance Plasmids for Commensal and Pathogenic Bacteria.

Antibiotic resistance (AR) is a threat to modern medicine, and plasmids are driving the global spread of AR by horizontal gene transfer across microbiomes and environments. Determining the mobile resistome responsible for this spread of AR among environments is essential in our efforts to attenuate the current crisis. Biosolids are a wastewater treatment plant (WWTP) byproduct used globally as fertilizer in agriculture. Here, we investigated the mobile resistome of biosolids that are used as fertilizer. This was done by capturing resistance plasmids that can transfer to human pathogens and commensal bacteria. We used a higher-throughput version of the exogenous plasmid isolation approach by mixing several ESKAPE pathogens and a commensal Escherichia coli with biosolids and screening for newly acquired resistance to about 10 antibiotics in these strains. Six unique resistance plasmids transferred to Salmonella typhimurium, Klebsiella aerogenes, and E. coli. All the plasmids were self-transferable and carried 3-6 antibiotic resistance genes (ARG) conferring resistance to 2-4 antibiotic classes. These plasmids-borne resistance genes were further embedded in genetic elements promoting intracellular recombination (i.e., transposons or class 1 integrons). The plasmids belonged to the broad-host-range plasmid (BHR) groups IncP-1 or PromA. Several of them were persistent in their new hosts when grown in the absence of antibiotics, suggesting that the newly acquired drug resistance traits would be sustained over time. This study highlights the role of BHRs in the spread of ARG between environmental bacteria and human pathogens and commensals, where they may persist. The work further emphasizes biosolids as potential vehicles of highly mobile plasmid-borne antibiotic resistance genes.

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.

The Species-Specific Acquisition and Diversification of a K1-like Family of Killer Toxins in Budding Yeasts of the Saccharomycotina.

Killer toxins are extracellular antifungal proteins that are produced by a wide variety of fungi, including Saccharomyces yeasts. Although many Saccharomyces killer toxins have been previously identified, their evolutionary origins remain uncertain given that many of these genes have been mobilized by double-stranded RNA (dsRNA) viruses. A survey of yeasts from the Saccharomyces genus has identified a novel killer toxin with a unique spectrum of activity produced by Saccharomyces paradoxus. The expression of this killer toxin is associated with the presence of a dsRNA totivirus and a satellite dsRNA. Genetic sequencing of the satellite dsRNA confirmed that it encodes a killer toxin with homology to the canonical ionophoric K1 toxin from Saccharomyces cerevisiae and has been named K1-like (K1L). Genomic homologs of K1L were identified in six non-Saccharomyces yeast species of the Saccharomycotina subphylum, predominantly in subtelomeric regions of the genome. When ectopically expressed in S. cerevisiae from cloned cDNAs, both K1L and its homologs can inhibit the growth of competing yeast species, confirming the discovery of a family of biologically active K1-like killer toxins. The sporadic distribution of these genes supports their acquisition by horizontal gene transfer followed by diversification. The phylogenetic relationship between K1L and its genomic homologs suggests a common ancestry and gene flow via dsRNAs and DNAs across taxonomic divisions. This appears to enable the acquisition of a diverse arsenal of killer toxins by different yeast species for potential use in niche competition.

Diversification, Introgression, and Rampant Cytonuclear Discordance in Rocky Mountains Chipmunks (Sciuridae: Tamias).

Evidence from natural systems suggests that hybridization between animal species is more common than traditionally thought, but the overall contribution of introgression to standing genetic variation within species remains unclear for most animal systems. Here, we use targeted exon capture to sequence thousands of nuclear loci and complete mitochondrial genomes from closely related chipmunk species in the Tamias quadrivittatus group that are distributed across the Great Basin and the central and southern Rocky Mountains of North America. This recent radiation includes six overlapping, ecologically distinct species (Tamias canipes, Tamias cinereicollis, Tamias dorsalis, T. quadrivittatus, Tamias rufus, and Tamias umbrinus) that show evidence for widespread introgression across species boundaries. Such evidence has historically been derived from a handful of markers, typically focused on mitochondrial loci, to describe patterns of introgression; consequently, the extent of introgression of nuclear genes is less well characterized. We conducted a series of phylogenomic and species-tree analyses to resolve the phylogeny of six species in this group. In addition, we performed several population-genomic analyses to characterize nuclear genomes and infer coancestry among individuals. Furthermore, we used emerging quartets-based approaches to simultaneously infer the species tree (SVDquartets) and identify introgression (HyDe). We found that, in spite of rampant introgression of mitochondrial genomes between some species pairs (and sometimes involving up to three species), there appears to be little to no evidence for nuclear introgression. These findings mirror other genomic results where complete mitochondrial capture has occurred between chipmunk species in the absence of appreciable nuclear gene flow. The underlying causes of recurrent massive cytonuclear discordance remain unresolved in this group but mitochondrial DNA appears highly misleading of population histories as a whole. Collectively, it appears that chipmunk species boundaries are largely impermeable to nuclear gene flow and that hybridization, while pervasive with respect to mtDNA, has likely played a relatively minor role in the evolutionary history of this group. [Cytonuclear discordance; hyridization; introgression, phylogenomics; SVDquartets; Tamias.].

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.

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.

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.

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.

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias).

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

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.

Evolutionary Paths That Expand Plasmid Host-Range: Implications for Spread of Antibiotic Resistance.

The World Health Organization has declared the emergence of antibiotic resistance to be a global threat to human health. Broad-host-range plasmids have a key role in causing this health crisis because they transfer multiple resistance genes to a wide range of bacteria. To limit the spread of antibiotic resistance, we need to gain insight into the mechanisms by which the host range of plasmids evolves. Although initially unstable plasmids have been shown to improve their persistence through evolution of the plasmid, the host, or both, the means by which this occurs are poorly understood. Here, we sought to identify the underlying genetic basis of expanded plasmid host-range and increased persistence of an antibiotic resistance plasmid using a combined experimental-modeling approach that included whole-genome resequencing, molecular genetics and a plasmid population dynamics model. In nine of the ten previously evolved clones, changes in host and plasmid each slightly improved plasmid persistence, but their combination resulted in a much larger improvement, which indicated positive epistasis. The only genetic change in the plasmid was the acquisition of a transposable element from a plasmid native to the Pseudomonas host used in these studies. The analysis of genetic deletions showed that the critical genes on this transposon encode a putative toxin-antitoxin (TA) and a cointegrate resolution system. As evolved plasmids were able to persist longer in multiple naïve hosts, acquisition of this transposon also expanded the plasmid's host range, which has important implications for the spread of antibiotic resistance.

The population genomics of rapid adaptation: disentangling signatures of selection and demography in white sands lizards.

Understanding the process of adaptation during rapid environmental change remains one of the central focal points of evolutionary biology. The recently formed White Sands system of southern New Mexico offers an outstanding example of rapid adaptation, with a variety of species having rapidly evolved blanched forms on the dunes that contrast with their close relatives in the surrounding dark soil habitat. In this study, we focus on two of the White Sands lizard species, Sceloporus cowlesi and Aspidoscelis inornata, for which previous research has linked mutations in the melanocortin-1 receptor gene (Mc1r) to blanched coloration. We sampled populations both on and off the dunes and used a custom sequence capture assay based on probed fosmid libraries to obtain >50 kb of sequence around Mc1r and hundreds of other random genomic locations. We then used model-based statistical inference methods to describe the demographic and adaptive history characterizing the colonization of White Sands. We identified a number of similarities between the two focal species, including strong evidence of selection in the blanched populations in the Mc1r region. We also found important differences between the species, suggesting different colonization times, different genetic architecture underlying the blanched phenotype and different ages of the beneficial alleles. Finally, the beneficial allele is dominant in S. cowlesi and recessive in A. inornata, allowing for a rare empirical test of theoretically expected patterns of selective sweeps under these differing models.

Draft Genome Sequence of Pasteurella multocida Isolate P1062, Isolated from Bovine Respiratory Disease.

Here, we report the draft genome of Pasteurella multocida isolate P1062 recovered from pneumonic bovine lung in the United States in 1959.