Search

Most current Salmonella subtyping analyses rely on whole genome sequencing (WGS), which focuses on the high-resolution analysis of single genomes or multiple single genomes from the isolated colonies on microbiological agar plates. In this study, we introduce bioinformatics innovations for a metagenomic outbreak response workflow that accurately identifies multiple Salmonella serovars at the same time. bettercallsal is one of the first analysis tools to identify multiple Salmonella enterica serotypes from metagenomic or quasi-metagenomic datasets with high accuracy, allowing these isolate-independent methods to be incorporated into surveillance and root cause investigations. It was tested on an in silico benchmark dataset comprising 29 unique Salmonella serovars, 46 non-Salmonella bacterial genomes, and 10 viral genomes at varying read depths and on previously well-characterized and sequenced non-selective primary and selective enrichments of papaya and peach samples from separate outbreak investigations that resulted in the identification of multiple Salmonella serovars using traditional isolate culturing and WGS as well as nucleic acid assays. Analyses were also conducted on these datasets using a custom-built k-mer tool, SeqSero2, and Kallisto to compare serotype calling to bettercallsal. The in silico dataset analyzed with bettercallsal achieved the maximum precision, recall, and accuracy of 100, 83, and 94%, respectively. In the papaya outbreak samples, bettercallsal identified the presence of multiple serovars in agreement with the Luminex® xMAP assay results and also identified more serovars per sample, as evidenced by NCBI SNP clustering. In peach outbreak samples, bettercallsal identified two serovars in concordance with k-mer analysis and the Luminex xMAP assay. The genome hit reported by bettercallsal clustered with the chicken isolate genome, as reported by the FDA peach outbreak investigation from sequenced isolates (WGS). Overall, bettercallsal outperformed k-mer, Seqsero2, and Kallisto in identifying multiple serovars from enrichment cultures using shotgun metagenomic sequencing.

Multiple sclerosis (MS) is the most prevalent demyelinating disease of the central nervous system, characterized by myelin destruction, axonal degeneration, and progressive loss of neurological functions. Remyelination is considered an axonal protection strategy and may enable functional recovery, but the mechanisms of myelin repair, especially after chronic demyelination, remain poorly understood. Here, we used the cuprizone demyelination mouse model to investigate spatiotemporal characteristics of acute and chronic de- and remyelination and motor functional recovery following chronic demyelination. Extensive remyelination occurred after both the acute and chronic insults, but with less robust glial responses and slower myelin recovery in the chronic phase. Axonal damage was found at the ultrastructural level in the chronically demyelinated corpus callosum and in remyelinated axons in the somatosensory cortex. Unexpectedly, we observed the development of functional motor deficits after chronic remyelination. RNA sequencing of isolated brain regions revealed significantly altered transcripts across the corpus callosum, cortex and hippocampus. Pathway analysis identified selective upregulation of extracellular matrix/collagen pathways and synaptic signaling in the chronically de/remyelinating white matter. Our study demonstrates regional differences of intrinsic reparative mechanisms after a chronic demyelinating insult and suggests a potential link between long-term motor function alterations and continued axonal damage during chronic remyelination. Moreover, the transcriptome dataset of three brain regions and over an extended de/remyelination period provides a valuable platform for a better understanding of the mechanisms of myelin repair as well as the identification of potential targets for effective remyelination and neuroprotection for progressive MS.

ABSTRACT Chicks are ideal to follow the development of the intestinal microbiota and to understand how a pathogen perturbs this developing population. Taxonomic/metagenomic analyses captured the development of the chick microbiota in unperturbed chicks and in chicks infected with Salmonella enterica serotype Typhimurium (STm) during development. Taxonomic analysis suggests that colonization by the chicken microbiota takes place in several waves. The cecal microbiota stabilizes at day 12 posthatch with prominent Gammaproteobacteria and Clostridiales. Introduction of S. Typhimurium at day 4 posthatch disrupted the expected waves of intestinal colonization. Taxonomic and metagenomic shotgun sequencing analyses allowed us to identify species present in uninfected chicks. Untargeted metabolomics suggested different metabolic activities in infected chick microbiota. This analysis and gas chromatography-mass spectrometry on ingesta confirmed that lactic acid in cecal content coincides with the stable presence of enterococci in STm-infected chicks. Unique metabolites, including 2-isopropylmalic acid, an intermediate in the biosynthesis of leucine, were present only in the cecal content of STm-infected chicks. The metagenomic data suggested that the microbiota in STm-infected chicks contained a higher abundance of genes, from STm itself, involved in branched-chain amino acid synthesis. We generated an ilvC deletion mutant (STM3909) encoding ketol-acid-reductoisomerase, a gene required for the production of l-isoleucine and l-valine. ΔilvC mutants are disadvantaged for growth during competitive infection with the wild type. Providing the ilvC gene in trans restored the growth of the ΔilvC mutant. Our integrative approach identified biochemical pathways used by STm to establish a colonization niche in the chick intestine during development. IMPORTANCE Chicks are an ideal model to follow the development of the intestinal microbiota and to understand how a pathogen perturbs this developing population. Using taxonomic and metagenomic analyses, we captured the development of chick microbiota to 19 days posthatch in unperturbed chicks and in chicks infected with Salmonella enterica serotype Typhimurium (STm). We show that normal development of the microbiota takes place in waves and is altered in the presence of a pathogen. Metagenomics and metabolomics suggested that branched-chain amino acid biosynthesis is especially important for Salmonella growth in the infected chick intestine. Salmonella mutants unable to make l-isoleucine and l-valine colonize the chick intestine poorly. Restoration of the pathway for biosynthesis of these amino acids restored the colonizing ability of Salmonella. Integration of multiple analyses allowed us to correctly identify biochemical pathways used by Salmonella to establish a niche for colonization in the chick intestine during development.

ABSTRACT Coxiella burnetii requires a type IVB secretion system (T4SS) to promote intracellular replication and virulence. We hypothesized that Coxiella employs its T4SS to secrete effectors that enable stealthy colonization of immune cells. To address this, we used RNA sequencing to compare the transcriptional response of murine bone marrow-derived macrophages (BMDM) infected with those of wild-type Coxiella and a T4SS-null mutant at 8 and 24 h postinfection. We found a T4SS-independent upregulation of proinflammatory transcripts which was consistent with a proinflammatory polarization phenotype. Despite this, infected BMDM failed to completely polarize, as evidenced by modest surface expression of CD38 and CD11c, nitrate production, and reduced proinflammatory cytokine and chemokine secretion compared to positive controls. As these BMDM permitted replication of C. burnetii, we employed them to identify T4SS effectors that are essential in the specific cellular context of a primary macrophage. We found five Himar1 transposon mutants in T4SS effectors that had a replication defect in BMDM but not J774A.1 cells. The mutants were also attenuated in a SCID mouse model of infection. Among these candidate virulence factors, we found that CBU1639 contributed to the inhibition of macrophage proinflammatory responses to Coxiella infection. These data demonstrate that while T4SS is dispensable for the stealthy invasion of primary macrophages, Coxiella has evolved multiple T4SS effectors that specifically target macrophage function to proliferate within that specific cellular context. IMPORTANCE Coxiella burnetii, the causative agent of Q fever, preferentially infects macrophages of the respiratory tract when causing human disease. This work describes how primary macrophages respond to C. burnetii at the earliest stages of infection, before bacterial replication. We found that while infected macrophages increase expression of proinflammatory genes after bacterial entry, they fail to activate the accompanying antibacterial functions that might ultimately control the infection. This disconnect between initial response and downstream function was not mediated by the bacterium’s type IVB secretion system, suggesting that Coxiella has other virulence factors that dampen host responses early in the infection process. Nevertheless, we were able to identify several type IVB secreted effectors that were specifically required for survival in macrophages and mice. This work is the first to identify type IVB secretion effectors that are specifically required for infection and replication within primary macrophages.

Abstract Checkpoint blockade-mediated immunotherapy is emerging as an effective treatment modality for multiple cancer types. However, cancer cells frequently evade the immune system, compromising the effectiveness of immunotherapy. It is crucial to develop screening methods to identify the patients who would most benefit from these therapies because of the risk of the side effects and the high cost of treatment. Here we show that expression of the MHC class I transactivator (CITA), NLRC5, is important for efficient responses to anti-CTLA-4 and anti-PD1 checkpoint blockade therapies. Melanoma tumors derived from patients responding to immunotherapy exhibited significantly higher expression of NLRC5 and MHC class I-related genes compared to non-responding patients. In addition, multivariate analysis that included the number of tumor-associated non-synonymous mutations, predicted neo-antigen load and PD-L2 expression was capable of further stratifying responders and non-responders to anti-CTLA4 therapy. Moreover, expression or methylation of NLRC5 together with total somatic mutation number were significantly correlated with increased patient survival. These results suggest that NLRC5 tumor expression, alone or together with tumor mutation load constitutes a valuable predictive biomarker for both prognosis and response to anti-CTLA-4 and potentially anti-PD1 blockade immunotherapy in melanoma patients.

Sex chromosomes and sex determining genes can evolve fast, with the sex-linked chromosomes often differing between closely related species. Population genetics theory has been developed and tested to explain the rapid evolution of sex chromosomes and sex determination. However, we do not know why the sex chromosomes are divergent in some taxa and conserved in others. Addressing this question requires comparing closely related taxa with conserved and divergent sex chromosomes to identify biological features that could explain these differences. Cytological karyotypes suggest that muscid flies (e.g., house fly) and blow flies are such a taxonomic pair. The sex chromosomes appear to differ across muscid species, whereas they are conserved across blow flies. Despite the cytological evidence, we do not know the extent to which muscid sex chromosomes are independently derived along different evolutionary lineages. To address that question, we used genomic and transcriptomic sequence data to identify young sex chromosomes in two closely related muscid species, horn fly (Haematobia irritans) and stable fly (Stomoxys calcitrans). We provide evidence that the nascent sex chromosomes of horn fly and stable fly were derived independently from each other and from the young sex chromosomes of the closely related house fly (Musca domestica). We present three different scenarios that could have given rise to the sex chromosomes of horn fly and stable fly, and we describe how the scenarios could be distinguished. Distinguishing between these scenarios in future work could identify features of muscid genomes that promote sex chromosome divergence.

Summary: Prenatal alcohol exposure (PAE) results in cerebral cortical dysgenesis. Single-cell RNA sequencing was performed on murine fetal cerebral cortical cells from six timed pregnancies, to decipher persistent cell- and sex-specific effects of an episode of PAE during early neurogenesis. We found, in an analysis of 38 distinct neural subpopulations across 8 lineage subtypes, that PAE altered neural maturation and cell cycle and disrupted gene co-expression networks. Whereas most differentially regulated genes were inhibited, particularly in females, PAE also induced sex-independent neural expression of fetal hemoglobin, a presumptive epigenetic stress adaptation. PAE inhibited Bcl11a, Htt, Ctnnb1, and other upstream regulators of differentially expressed genes and inhibited several autism-linked genes, suggesting that neurodevelopmental disorders share underlying mechanisms. PAE females exhibited neural loss of X-inactivation, with correlated activation of autosomal genes and evidence for spliceosome dysfunction. Thus, episodic PAE persistently alters the developing neural transcriptome, contributing to sex- and cell-type-specific teratology.

Honey bee (Apis mellifera) queens have a remarkable organ, the spermatheca, which successfully stores sperm for years after a virgin queen mates. This study uniquely characterized and quantified the transcriptomes of the spermathecae from mated and virgin honey bee queens via RNA sequencing to identify differences in mRNA levels based on a queen's mating status. The transcriptome of drone semen was analyzed for comparison. Samples from three individual bees were independently analyzed for mated queen spermathecae and virgin queen spermathecae, and three pools of semen from ten drones each were collected from three separate colonies. In total, the expression of 11,233 genes was identified in mated queen spermathecae, 10,521 in virgin queen spermathecae, and 10,407 in drone semen. Using a cutoff log2 fold-change value of 2.0, we identified 212 differentially expressed genes between mated and virgin spermathecal queen tissues: 129 (1.4% of total) were up-regulated and 83 (0.9% of total) were down-regulated in mated queen spermathecae. Three genes in mated queen spermathecae, three genes in virgin queen spermathecae and four genes in drone semen that were more highly expressed in those tissues from the RNA sequencing data were further validated by real time quantitative PCR. Among others, expression of Kielin/chordin-like and Trehalase mRNAs was highest in the spermathecae of mated queens compared to virgin queen spermathecae and drone semen. Expression of the mRNA encoding Alpha glucosidase 2 was higher in the spermathecae of virgin queens. Finally, expression of Facilitated trehalose transporter 1 mRNA was greatest in drone semen. This is the first characterization of gene expression in the spermathecae of honey bee queens revealing the alterations in mRNA levels within them after mating. Future studies will extend to other reproductive tissues with the purpose of relating levels of specific mRNAs to the functional competence of honey bee queens and the colonies they head.

ABSTRACT Staphylococcus epidermidis is a leading cause of nosocomial infections in patients with a compromised immune system and/or an implanted medical device. Seventy to 90% of S. epidermidis clinical isolates are methicillin resistant and carry the mecA gene, present in a mobile genetic element (MGE) called the staphylococcal cassette chromosome mec (SCCmec) element. Along with the presence of antibiotic and heavy metal resistance genes, MGEs can also contain genes encoding secreted or cell wall-anchored virulence factors. In our earlier studies of S. epidermidis clinical isolates, we discovered S. epidermidis surface protein J (SesJ), a prototype of a recently discovered subfamily of the microbial surface component recognizing adhesive matrix molecule (MSCRAMM) group. MSCRAMMs are major virulence factors of pathogenic Gram-positive bacteria. Here, we report that the sesJ gene is always accompanied by two glycosyltransferase genes, gtfA and gtfB, and is present in two MGEs, called the arginine catabolic mobile element (ACME) and the staphylococcal cassette chromosome (SCC) element. The presence of the sesJ gene was associated with the left-hand direct repeat DR_B or DR_E. When inserted via DR_E, the sesJ gene was encoded in the SCC element. When inserted via DR_B, the sesJ gene was accompanied by the genes for the type 1 restriction modification system and was encoded in the ACME. Additionally, the SCC element and ACME carry different isoforms of the SesJ protein. To date, the genes encoding MSCRAMMs have been seen to be located in the bacterial core genome. Here, we report the presence of an MSCRAMM in an MGE in S. epidermidis clinical isolates. IMPORTANCE S. epidermidis is an opportunistic bacterium that has established itself as a successful nosocomial pathogen. The modern era of novel therapeutics and medical devices has extended the longevity of human life, but at the same time, we also witness the evolution of pathogens to adapt to newly available niches in the host. Increasing antibiotic resistance among pathogens provides an example of such pathogen adaptation. With limited opportunities to modify the core genome, most of the adaptation occurs by acquiring new genes, such as virulence factors and antibiotic resistance determinants present in MGEs. In this study, we describe that the sesJ gene, encoding a recently discovered cell wall-anchored protein in S. epidermidis, is present in both ACME and the SCC element. The presence of virulence factors in MGEs can influence the virulence potential of a specific strain. Therefore, it is critical to study the virulence factors found in MGEs in emerging pathogenic bacteria or strains to understand the mechanisms used by these bacteria to cause infections.

Haematobia irritans, commonly known as the horn fly, is a globally distributed blood-feeding pest of cattle that is responsible for significant economic losses to cattle producers. Chemical insecticides are the primary means for controlling this pest but problems with insecticide resistance have become common in the horn fly. To provide a foundation for identification of genomic loci for insecticide resistance and for discovery of new control technology, we report the sequencing, assembly, and annotation of the horn fly genome. The assembled genome is 1.14 Gb, comprising 76,616 scaffolds with N50 scaffold length of 23 Kb. Using RNA-Seq data, we have predicted 34,413 gene models of which 19,185 have been assigned functional annotations. Comparative genomics analysis with the Dipteran flies Musca domestica L., Drosophila melanogaster, and Lucilia cuprina, show that the horn fly is most closely related to M. domestica, sharing 8,748 orthologous clusters followed by D. melanogaster and L. cuprina, sharing 7,582 and 7,490 orthologous clusters respectively. We also identified a gene locus for the sodium channel protein in which mutations have been previously reported that confers target site resistance to the most common class of pesticides used in fly control. Additionally, we identified 276 genomic loci encoding members of metabolic enzyme gene families such as cytochrome P450s, esterases and glutathione S-transferases, and several genes orthologous to sex determination pathway genes in other Dipteran species.

Virus-induced neurological sequelae resulting from infection by Theiler’s murine encephalomyelitis virus (TMEV) are used for studying human conditions ranging from epileptic seizures to demyelinating disease. Mouse strains are typically considered susceptible or resistant to TMEV infection based on viral persistence and extreme phenotypes, such as demyelination. We have identified a broader spectrum of phenotypic outcomes by infecting strains of the genetically diverse Collaborative Cross (CC) mouse resource. We evaluated the chronic-infection gene expression profiles of hippocampi and thoracic spinal cords for 19 CC strains in relation to phenotypic severity and TMEV persistence. Strains were clustered based on similar phenotypic profiles and TMEV levels at 90 days post-infection, and we categorized distinct TMEV response profiles. The three most common profiles included “resistant” and “susceptible,” as before, as well as a “resilient” TMEV response group which experienced both TMEV persistence and mild neurological phenotypes even at 90 days post-infection. Each profile had a distinct gene expression signature, allowing the identification of pathways and networks specific to each TMEV response group. CC founder haplotypes for genes involved in these pathways/networks revealed candidate response-specific alleles. These alleles demonstrated pleiotropy and epigenetic (miRNA) regulation in long-term TMEV infection, with particular relevance for resilient mouse strains.

Abstract Background The ubiquitin protein E3A ligase gene (UBE3A) gene is imprinted with maternal-specific expression in neurons and biallelically expressed in all other cell types. Both loss-of-function and gain-of-function mutations affecting the dosage of UBE3A are associated with several neurodevelopmental syndromes and psychological conditions, suggesting that UBE3A is dosage-sensitive in the brain. The observation that loss of imprinting increases the dosage of UBE3A in brain further suggests that inactivation of the paternal UBE3A allele evolved as a dosage-regulating mechanism. To test this hypothesis, we examined UBE3A transcript and protein levels among cells, tissues, and species with different imprinting states of UBE3A. Results Overall, we found no correlation between the imprinting status and dosage of UBE3A. Importantly, we found that maternal Ube3a protein levels increase in step with decreasing paternal Ube3a protein levels during neurogenesis in mouse, fully compensating for loss of expression of the paternal Ube3a allele in neurons. Conclusions Based on our findings, we propose that imprinting of UBE3A does not function to reduce the dosage of UBE3A in neurons but rather to regulate some other, as yet unknown, aspect of gene expression or protein function.

Abstract Introduction Prenatal alcohol exposure can contribute to fetal alcohol spectrum disorders (FASD), characterized by a myriad of developmental impairments affecting behavior and cognition. Studies show that many of these functional impairments are associated with the hippocampus, a structure exhibiting exquisite vulnerability to developmental alcohol exposure and critically implicated in learning and memory; however, mechanisms underlying alcohol‐induced hippocampal deficits remain poorly understood. By utilizing a high‐throughput RNA‐sequencing (RNA‐seq) approach to address the neurobiological and molecular basis of prenatal alcohol‐induced hippocampal functional deficits, we hypothesized that chronic binge prenatal alcohol exposure alters gene expression and global molecular pathways in the fetal hippocampus. Methods Timed‐pregnant Sprague–Dawley rats were randomly assigned to a pair‐fed control (PF) or binge alcohol (ALC) treatment group on gestational day (GD) 4. ALC dams acclimatized from GDs 5–10 with a daily treatment of 4.5 g/kg alcohol and subsequently received 6 g/kg on GDs 11–20. PF dams received a once daily maltose dextrin gavage on GDs 5–20, isocalorically matching ALC counterparts. On GD 21, bilateral hippocampi were dissected, flash frozen, and stored at −80°C. Total RNA was then isolated from homogenized tissues. Samples were normalized to ~4nM and pooled equally. Sequencing was performed by Illumina NextSeq 500 on a 75 cycle, single‐end sequencing run. Results RNA‐seq identified 13,388 genes, of these, 76 genes showed a significant difference (p

In 2013, the U.S. Department of Agriculture Food Safety and Inspection Service (USDA-FSIS) began transitioning to whole genome sequencing (WGS) for foodborne disease outbreak- and recall-associated isolate identification of select bacterial species. While WGS offers greater precision, certain hurdles must be overcome before widespread application within the food industry is plausible. Challenges include diversity of sequencing platform outputs and lack of standardized bioinformatics workflows for data analyses. We sequenced DNA from USDA-FSIS approved, non-pathogenic E. coli surrogates and a derivative group of rifampicin-resistant mutants (rifR) via both Oxford Nanopore MinION and Illumina MiSeq platforms to generate and annotate complete genomes. Genome sequences from each clone were assembled separately so long-read, short-read, and combined sequence assemblies could be directly compared. The combined sequence data approach provides more accurate completed genomes. The genomes from these isolates were verified to lack functional key E. coli elements commonly associated with pathogenesis. Genetic alterations known to confer rifR were also identified. As the food industry adopts WGS within its food safety programs, these data provide completed genomes for commonly used surrogate strains, with a direct comparison of sequence platforms and assembly strategies relevant to research/testing workflows applicable for both processors and regulators.

We present SBEToolbox (Systems Biology and Evolution Toolbox), an open-source Matlab toolbox for biological network analysis. It takes a network file as input, calculates a variety of centralities and topological metrics, clusters nodes into modules, and displays the network using different graph layout algorithms. Straightforward implementation and the inclusion of high-level functions allow the functionality to be easily extended or tailored through developing custom plugins. SBEGUI, a menu-driven graphical user interface (GUI) of SBEToolbox, enables easy access to various network and graph algorithms for programmers and non-programmers alike. All source code and sample data are freely available at https://github.com/biocoder/SBEToolbox/releases .