Search

Precise and rapid identification ofSalmonellaserotypes from suspect food matrices is critical for successful source attribution of illness outbreaks (Scallan et al., 2011). Currently, close to 3% of U.S. foodborneSalmonellaoutbreaks have been attributed to multiple Salmonellaserotypes (2.85%, 2000 - 2020) (CDC, 2022). Recent foodborne outbreaks that have been attributed to multipleSalmonellaserotypes force us to question whether these are rare events or if previous methods simply did not have the throughput to provide an accurate picture of the complex ecology that is connected to outbreak etiologies. (Hassan et al., 2019; FDA, 2021; Whitney et al., 2021).Anin-silicobenchmark dataset, comprising 29 uniqueSalmonella, 46 non-Salmonellabacterial and 10 viral genomes, was generated with varying read depths. For outbreak samples, analysis was performed on previously sequenced pre-enrichments and selective enrichments of papayas and peaches (fruits and leaves) that led to the identification of multiple serovars. Data analyses was performed using a custom-builtk-mertool, SeqSero2, Kallisto and bettercallsal.Thein-silicodataset analyzed with bettercallsal had accuracy, recall and specificity of 95%, 96% and 98 % respectively. In the papaya outbreak samples, bettercallsal identified multiple serovar presence in concordance with Bioplex assay results and the genome hits assigned to the samples areSalmonellaisolates from the papaya outbreak as evident by NCBI SNP cluster information. In peach outbreak samples, bettercallsal identified both the serovars (Alachua and Gaminara) in concordance withk-mer analysis and the Luminex xMap assay. bettercallsal outperformedk-mer, Kallisto and Seqsero2 in identifying multiple serovars from enrichment cultures using shotgun metagenomics sequencing.MostSalmonellasubtyping work has relied upon WGS methods which focuses on the high-resolution analysis of single genomes, or multiple single genomes picked from colonies on agar. Here we introduce laboratory and bioinformatics innovations for a metagenomic outbreak response workflow that accurately identifies multipleSalmonellaserovars at the same time in a much higher throughput approach. bettercallsal is one of the first analysis tools that can potentially identify multipleSalmonellaspp. serotypes from a metagenomic or quasi-metagenomic datasets with accuracy and can provide early insights into the etiology of the sample.

Insulin-like growth factor I (IGF-1) has been implicated in breast cancer due to its mitogenic and anti-apoptotic effects. Despite substantial research on the role of IGF-1 in tumor progression, the relationship of IGF-1 to tissue stem cells, particularly in mammary tissue, and the resulting tumor susceptibility has not been elucidated. Previous studies with the BK5.IGF-1 transgenic (Tg) mouse model reveals that IGF-1 does not act as a classical, post-carcinogen tumor promoter in the mammary gland. Pre-pubertal Tg mammary glands display increased numbers and enlarged sizes of terminal end buds, a niche for mammary stem cells (MaSCs). Here we show that MaSCs from both wild-type (WT) and Tg mice expressed IGF-1R and that overexpression of Tg IGF-1 increased numbers of MaSCs by undergoing symmetric division, resulting in an expansion of the MaSC and luminal progenitor (LP) compartments in pre-pubertal female mice. This expansion was maintained post-pubertally and validated by mammosphere assays in vitro and transplantation assays in vivo. The addition of recombinant IGF-1 promoted, and IGF-1R downstream inhibitors decreased mammosphere formation. Single-cell transcriptomic profiles generated from 2 related platforms reveal that IGF-1 stimulated quiescent MaSCs to enter the cell cycle and increased their expression of genes involved in proliferation, plasticity, tumorigenesis, invasion, and metastasis. This study identifies a novel, pro-tumorigenic mechanism, where IGF-1 increases the number of transformation-susceptible carcinogen targets during the early stages of mammary tissue development, and “primes” their gene expression profiles for transformation.

Rapid and comprehensive analysis of metagenomic data from any sample is of importance in food safety. Also important is the need for automated analysis pipelines allowing rapid and effective construction of metagenomic assembled genomes (MAGs) to enable bacterial source-tracking from metagenomic data. We developed a precision metagenomics approach for detecting and classifying shiga toxin-producingEscherichia coli (STEC) in enrichments of agricultural water using Oxford Nanopore long read sequencing where the bioinformatics data analysis employed many sequential manual steps (Maguire et al, 2021). Here we present centriflaken, a suite of automated data analysis workflows enabled by Nextflow, which takes metagenomic data and generates MAGs and performs in silico-based analysis as described in Maguire et al, 2021. The final summary plots and tables can be downloaded from the provided MultiQC HTML report generated as part of the pipeline. The centriflaken pipeline was validated with data from our previously published method and was able to replicate the detection and classification of STECs for each sample. We tested the pipeline with nanopore data obtained from 21 additional enriched samples from irrigation water and was able to perform the entire precision metagenomics analysis in less than 5 hours. We have further expanded this precision approach to include any user supplied taxa of interest (Listeria monocytogenes or Salmonella) not only STECs.

Antecedent viral infection may contribute to increased susceptibility to several neurological diseases, such as multiple sclerosis and Parkinson’s disease. Variation in clinical presentations of these diseases is often associated with gender, genetic background, or a combination of these and other factors. The complicated etiologies of these virally influenced diseases are difficult to study in conventional laboratory mouse models, which display a very limited number of phenotypes. We have used the genetically and phenotypically diverse Collaborative Cross mouse panel to examine complex neurological phenotypes after viral infection. Female and male mice from 18 CC strains were evaluated using a multifaceted phenotyping pipeline to define their unique disease profiles following infection with Theiler’s Murine Encephalomyelitis Virus, a neurotropic virus. We identified 4 distinct disease progression profiles based on limb-specific paresis and paralysis, tremors and seizures, and other clinical signs, along with separate gait profiles. We found that mice of the same strain had more similar profiles compared to those of different strains, and also identified strains and phenotypic parameters in which sex played a significant role in profile differences. These results demonstrate the value of using CC mice for studying complex disease subtypes influenced by sex and genetic background. Our findings will be useful for developing novel mouse models of virally induced neurological diseases with heterogenous presentation, an important step for designing personalized, precise treatments.

Background It is not unusual for stallions to have fertility problems. For many, artificial insemination with more dense spermatozoa (isolated by density gradient centrifugation) results in greater pregnancy rates compared with the rates when using unfractionated spermatozoa. RNAs in spermatozoa delivered to the oocyte at conception are required for embryo development. Novel molecular assays of spermatozoa that reflect function are needed to predict the fertility of stallions. Objectives To describe and compare the RNA populations in more dense and less dense spermatozoa from stallions. Materials and methods Spermatozoa from five stallions were separated into more dense and less dense populations by density gradient centrifugation. Complementary DNA libraries were made from each of the ten total RNA samples after ribosomal RNA removal. Next-generation sequencing characterized the RNA populations in more and less dense spermatozoa. Quantitative reverse transcription-PCR was used to confirm differential expression of selected RNAs. Results Stallion spermatozoa contain 11 215 RNAs, with the most prevalent RNA being a 1492 base long non-coding RNA. The levels of 159 RNAs were greater in more dense spermatozoa, while levels of seven other RNAs were greater in less dense spermatozoa. Quantitative reverse transcription-PCR confirmed the threefold greater levels of solute carrier family 26 member 8 (SLC26A8) mRNA in less dense spermatozoa, and sixfold and threefold greater expression levels of the SCP2 sterol binding domain containing 1 (SCP2D1) and spermatogenesis-associated protein 31D1 (SPATA31D1) mRNAs in more dense spermatozoa, respectively. Discussion and conclusion We identified 11 215 RNAs in stallion spermatozoa and 166 with differential expression between more dense and less dense fractions. Many prevalent RNAs were also found in bull, boar, and human spermatozoa. Many differentially expressed RNAs are known to be testis- or spermatozoa-specific. Our results may lead to identification of an RNA population in spermatozoa that is optimal for establishing successful pregnancies.

Significance Myeloid lineage cells are suspected of having an integral role in pathophysiological processes of multiple sclerosis (MS), but the molecular mechanism(s) governing their effector function remains incompletely understood. We show that STAT3 is activated in myeloid cells near active MS lesions. Conditional deletion of Stat3 in myeloid cells abolished symptoms of experimental autoimmune encephalomyelitis (EAE) by suppressing the generation of pathogenic T cells. Functional and transcriptomic analyses of myeloid cells from wild-type and conditional knockout mice implicated antigen processing/presentation and inflammatory cytokine production as the mechanism underlying the effects of STAT3 deletion on impaired T cell activation. Our data indicate that targeting STAT3 in myeloid cells might be a viable treatment option for autoimmune demyelinating disease., Multiple sclerosis (MS) is an autoimmune inflammatory demyelinating disease of the central nervous system. Dysregulation of STAT3, a transcription factor pivotal to various cellular processes including Th17 cell differentiation, has been implicated in MS. Here, we report that STAT3 is activated in infiltrating monocytic cells near active MS lesions and that activation of STAT3 in myeloid cells is essential for leukocyte infiltration, neuroinflammation, and demyelination in experimental autoimmune encephalomyelitis (EAE). Genetic disruption of Stat3 in peripheral myeloid lineage cells abrogated EAE, which was associated with decreased antigen-specific T helper cell responses. Myeloid cells from immunized Stat3 mutant mice exhibited impaired antigen-presenting functions and were ineffective in driving encephalitogenic T cell differentiation. Single-cell transcriptome analyses of myeloid lineage cells from preclinical wild-type and mutant mice revealed that loss of myeloid STAT3 signaling disrupted antigen-dependent cross-activation of myeloid cells and T helper cells. This study identifies a previously unrecognized requisite for myeloid cell STAT3 in the activation of myelin-reactive T cells and suggests myeloid STAT3 as a potential therapeutic target for autoimmune demyelinating disease.

Sex chromosomes and sex determining genes can evolve fast, with the sex-linked chromosomes often differing between closely related species. A substantial body of 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-linked chromosomes differ between some species pairs yet are relatively conserved in other taxa. Addressing this question will require comparing closely related taxa with conserved and divergent sex chromosomes and sex determination systems to identify biological features that could explain these rate 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 highly 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 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 help to identify features of muscid genomes that promote sex chromosome divergence.

Angelman syndrome is a devastating neurogenetic disorder for which there is currently no effective treatment. It is caused by mutations or epimutations affecting the expression or function of the maternally inherited allele of the ubiquitin-protein ligase E3A (UBE3A) gene. The paternalUBE3Aallele is imprinted in neurons of the central nervous system (CNS) by theUBE3Aantisense (UBE3A-AS) transcript, which represents the distal end of theSNHG14transcription unit. Reactivating the expression of the paternalUBE3Aallele in the CNS has long been pursued as a therapeutic option for Angelman syndrome. Here, we designed and optimized antisense oligonucleotides (ASO) targeting an evolutionarily conserved region demarcating the start of the humanUBE3A-AStranscript and show that ASOs targeting this region can reverse imprinting ofUBE3Ain cultured Angelman syndrome neurons and throughout the CNS of a non-human primate model. Findings from this study advanced the first investigational molecular therapy for Angelman syndrome into clinical development (ClinicalTrials.gov,NCT04259281).SUMMARYHere, we describe the preclinical studies supporting the first investigational molecular therapy for Angelman syndrome to advance into clinical development (ClinicalTrials.gov,NCT04259281).

Binge alcohol exposure during pregnancy results in diminished vessel function and altered proteome in the maternal uterine artery. We aimed to utilize high throughput RNA-seq deep- sequencing to characterize specific effects of binge alcohol exposure during pregnancy on the uterine artery transcriptome, and gain insight into mechanisms under lying alcohol-mediated uterine artery dysfunction. Pregnant Sprague-Dawley rats assigned to Pair-Fed Control or Alcohol groups, received a once-daily orogastric gavage in a binge paradigm. RNA-sequencing using Illumina NextSeq 500, identified 13,941 genes; 40 significantly altered genes were altered by log2(fold change) > 2; 27 genes were upregulated and 13 were downregulated in the Alcohol group. Transcripts altered included those which encode for aldehyde dehydrogenases, matrix metalloproteases, and molecules vital for vasodilation and vascular remodeling. Biological pathways that were disproportionally altered by alcohol were proline and citrulline biosynthesis/metabolism. Disruption of these pathways suggests candidate mechanism(s) for alcohol-mediated impairments to the proteome and vascular function.

Human alcoholic hepatitis (AH) carries a high mortality rate. AH is an acute-on-chronic form of liver injury characterized by hepatic steatosis, ballooned hepatocytes, neutrophil infiltration, and pericellular fibrosis. We aimed to study the pathogenesis of AH in an animal model which combines chronic hepatic fibrosis with intragastric alcohol administration. Adult male C57BL6/J mice were treated with CCl4 (0.2 ml/kg, 2×weekly by intraperitoneal injections for 6 weeks) to induce chronic liver fibrosis. Then, ethyl alcohol (up to 25 g/kg/day for 3 weeks) was administered continuously to mice via a gastric feeding tube, with or without one-half dose of CCl4. Liver and serum markers and liver transcriptome were evaluated to characterize acute-on-chronic-alcoholic liver disease in our model. CCl4 or alcohol treatment alone induced liver fibrosis or steatohepatitis, respectively, findings that were consistent with expected pathology. Combined treatment resulted in a marked exacerbation of liver injury, as evident by the development of inflammation, steatosis, and pericellular fibrosis, pathological features of human AH. E. coli and Candida were also detected in livers of mice cotreated with CCl4 and alcohol, indicating pathogen translocation from gut to liver, similar to human AH. Importantly, liver transcriptomic changes specific to combined treatment group demonstrated close concordance with pathways perturbed in patients with severe AH. Overall, mice treated with CCl4 and alcohol displayed key molecular and pathological characteristics of human AH—pericellular fibrosis, increased hepatic bacterial load, and dysregulation of the same molecular pathways. This model may be useful for developing therapeutics for AH.

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.

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.

Although epidermal growth factor receptor (EGFR)-targeted therapies are approved for colorectal cancer (CRC) treatment, only 15% of CRC patients respond to EGFR inhibition. Here, we show that colorectal cancers (CRC) can initiate and grow faster through an EGFR-independent mechanism, irrespective of the presence of EGFR, in two different mouse models using tissue-specific ablation of Egfr. The growth benefit in the absence of EGFR is also independent of Kras status. An EGFR-independent gene expression signature, also observed in human CRCs, revealed that anergy-inducing genes are overexpressed in EGFR-independent polyps, suggesting increased infiltration of anergic lymphocytes promotes an accelerated growth rate that is partially caused by escape from cell-mediated immune responses. Many genes in the EGFR-independent gene expression signature are downstream targets of interleukin 10 receptor alpha (IL10RA). We further show that IL10 is detectable in serum from mice with EGFR-independent colon polyps. Using organoids in vitro and Src ablation in vivo, we show that IL10 contributes to growth of EGFR-independent CRCs, potentially mediated by the well-documented role of SRC in IL10 signaling. Based on these data, we show that the combination of an EGFR inhibitor with an anti-IL10 neutralizing antibody results in decreased cell proliferation in organoids and in decreased polyp size in pre-clinical models harboring EGFR-independent CRCs, providing a new therapeutic intervention for CRCs resistant to EGFR inhibitor therapies.

This is a draft genome of an orf virus (ORFV) vaccine strain assembled via long- and short-read hybrid assembly. ORFV is a zoonotic pathogen that affects sheep and goats. The genome of the virus contained in the vaccine was found to have high similarity (98%) to those of other published strains.

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., 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.

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.

Mus musculus is the only known species from which embryonic stem cells (ESC) can be isolated under conditions requiring only leukemia inhibitory factor (LIF). Other species are non-permissive in LIF media, and form developmentally primed epiblast stem cells (EpiSC) similar to cells derived from post-implantation, egg cylinders. To evaluate whether non-permissiveness extends to induced pluripotent stem cells (iPSC), we derived iPSC from the eight founder strains of the mouse Collaborative Cross. Two strains, NOD/ShiLtJ and the WSB/EiJ, were non-permissive, consistent with the previous classification of NOD/ShiLtJ as non-permissive to ESC derivation. We determined non-permissiveness is recessive, and that non-permissive genomes do not compliment. We overcame iPSC non-permissiveness by using GSK3B and MEK inhibitors with serum, a technique we termed 2iS reprogramming. Although used for ESC derivation, GSK3B and MEK inhibitors have not been used during iPSC reprogramming because they inhibit survival of progenitor differentiated cells. iPSC derived in 2iS are more transcriptionally similar to ESC than EpiSC, indicating that 2iS reprogramming acts to overcome genetic background constraints. Finally, of species tested for ESC or iPSC derivation, only some M. musculus strains are permissive under LIF culture conditions suggesting that this is an evolutionarily derived characteristic in the M. musculus lineage.

Strangles, caused by Streptococcus equi subspecies equi (S. equi) is an infectious disease of horses with worldwide distribution, but there are limited data available regarding strain variation using whole genome sequencing among and within outbreaks in the United States (US), and how US isolates compare with S. equi isolated globally. To address this knowledge-gap, we compared the whole genomes of 54 S. equi isolates from Texas and Kentucky and those of 230 publicly available sequences of S. equi isolates collected from other countries. Our results show that despite minimal variation among isolates within an outbreak some mutations do occur among individual outbreak isolates. Some S. equi strains from the US are closely related to S. equi isolates from other countries, likely reflecting international dissemination of isolates. Collectively, these data improve our understanding of phenotypic and genotypic variation of isolates within an outbreak, and the international distribution of S. equi. We also identify a novel variant of the S. equi M-protein, and observed cases of strangles that were caused by the modified-live vaccine but that were not recognized as vaccine-associated at the time of clinical sample submission.

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.

Recent progress in DNA sequencing technologies and advanced bioinformatic tools have enabled researchers to rapidly decipher the tick microbiome. To date, however, a number of microbiome studies performed on Dermacentor reticulatus ticks is still quite limited. Despite the importance of this ixodid tick for veterinary and human medicine, only two investigations have examined its microbiome. Moreover, these studies analyzed only a limited number of ticks/tick pools. Given the scarcity of microbiome data for D. reticulatus in general and the lack of microbiome studies on tick species from Eastern Europe in particular, the objective of the current investigation was to analyze the microbiome of D. reticulatus ticks collected from three geographical regions of Ukraine. A total of 88 individual tick microbiomes were analyzed by sequencing the V6 region of 16S rRNA. As a result, numerous significant differences in the bacterial relative abundance were detected between males and females of D. reticulatus for each region. The alpha diversity measures indicate that microbiomes were significantly different between females of D. reticulatus inter-regionally. In contrast, the collective results for male ticks are more suggestive of inter-regional microbiome homogeneity. The overall findings indicate that the composition and diversity of the D. reticulatus microbiome can be impacted by geographical and sex-related factors.

Nonalcoholic fatty liver disease (NAFLD) is the most prevalent pathological liver condition in developed countries. NAFLD results in severe alterations in liver function, including xenobiotic metabolism. Perchloroethylene (PERC) is a ubiquitous environmental pollutant, a known hepatotoxicant in rodents, and a probable human carcinogen. It is known that PERC disposition and metabolism are affected by NAFLD in mice; here, we examined how NAFLD changes PERC-associated liver effects. Male C57Bl6/J mice were fed a low-fat diet (LFD), high-fat diet (HFD), or methionine/folate/choline-deficient diet (MCD) to model a healthy liver, or mild and severe forms of NAFLD, respectively. After 8 weeks on diets, mice were orally administered PERC (300 mg/kg/day) or vehicle (5% aqueous Alkamuls-EL620) for 5 days. PERC-induced liver effects were exacerbated in both NAFLD groups. PERC exposure was associated with up-regulation of genes involved in xenobiotic, lipid, and glutathione metabolism, and down-regulation of the complement and coagulation cascades, regardless of the diet. Interestingly, HFD-fed mice, not MCD-fed mice, were generally more sensitive to PERC-induced liver effects. This was indicated by histopathology and transcriptional responses, where induction of genes associated with cell cycle and inflammation were prominent. Liver effects positively correlated with diet-specific differences in liver concentrations of PERC. We conclude that NAFLD alters the toxicodynamics of PERC and that NAFLD is a susceptibility factor that should be considered in future risk management decisions for PERC and other chlorinated solvents.

Background Trichloroethylene (TCE) is a known carcinogen in humans and rodents. Previous studies of inter-strain variability in TCE metabolism were conducted in multi-strain panels of classical inbred mice with limited genetic diversity to identify gene-environment interactions associated with chemical exposure. Objectives To evaluate inter-strain variability in TCE metabolism and identify genetic determinants that are associated with TCE metabolism and effects using Collaborative Cross (CC), a large panel of genetically diverse strains of mice. Methods We administered a single oral dose of 0, 24, 80, 240, or 800 mg/kg of TCE to mice from 50 CC strains, and collected organs 24 h post-dosing. Levels of trichloroacetic acid (TCA), a major oxidative metabolite of TCE were measured in multiple tissues. Protein expression and activity levels of TCE-metabolizing enzymes were evaluated in the liver. Liver transcript levels of known genes perturbed by TCE exposure were also quantified. Genetic association mapping was performed on the acquired phenotypes. Results TCA levels varied in a dose- and strain-dependent manner in liver, kidney, and serum. The variability in TCA levels among strains did not correlate with expression or activity of a number of enzymes known to be involved in TCE oxidation. Peroxisome proliferator-activated receptor alpha (PPARα)-responsive genes were found to be associated with strain-specific differences in TCE metabolism. Conclusions This study shows that CC mouse population is a valuable tool to quantitatively evaluate inter-individual variability in chemical metabolism and to identify genes and pathways that may underpin population differences.

BACKGROUND Ideal management of sexually transmitted infections (STI) may require risk markers for pathology or vaccine development. Previously, we identified common genetic variants associated with chlamydial pelvic inflammatory disease (PID) and reduced fecundity. As this explains only a proportion of the long-term morbidity risk, we used whole-exome sequencing to identify biological pathways that may be associated with STI-related infertility. METHODS We obtained stored DNA from 43 non-Hispanic black women with PID from the PID Evaluation and Clinical Health Study. Infertility was assessed at a mean of 84 months. Principal component analysis revealed no population stratification. Potential covariates did not significantly differ between groups. Sequencing kernel association test was used to examine associations between aggregates of variants on a single gene and infertility. The results from the sequencing kernel association test were used to choose "focus genes" (P < 0.01; n = 150) for subsequent Ingenuity Pathway Analysis to identify "gene sets" that are enriched in biologically relevant pathways. RESULTS Pathway analysis revealed that focus genes were enriched in canonical pathways including, IL-1 signaling, P2Y purinergic receptor signaling, and bone morphogenic protein signaling. CONCLUSIONS Focus genes were enriched in pathways that impact innate and adaptive immunity, protein kinase A activity, cellular growth, and DNA repair. These may alter host resistance or immunopathology after infection. Targeted sequencing of biological pathways identified in this study may provide insight into STI-related infertility.

Genetic background is known to play a role in the ability to derive pluripotent, embryonic stem cells (ESC), a trait referred to as permissiveness. Previously we have demonstrated that induced pluripotent stem cells (iPSC) can be readily derived from non-permissive mouse strains by addition of serum-based media supplemented with GSK3B and MEK inhibitors, termed 2iS media, three days into reprogramming. Here we describe the derivation of second type of iPSC colony from non-permissive mouse strains that can be stably maintained independently of 2iS media. The resulting cells display transcriptional heterogeneity similar to that observed in ESC from permissive genetic backgrounds derived in conventional serum containing media supplemented with leukemia inhibitor factor. However, unlike previous studies that report exclusive subpopulations, we observe both exclusive and simultaneous expression of naive and primed cell surface markers. Herein, we explore shifts in pluripotency in the presence of 2iS and characterize heterogenous subpopulations to determine their pluripotent state and role in heterogenous iPSCs derived from the non-permissive NOD/ShiLtJ strain. We conclude that heterogeneity is a naturally occurring, necessary quality of stem cells that allows for the maintenance of pluripotency. This study further demonstrates the efficacy of the 2iS reprogramming technique. It is also the first study to derive stable ESC-like stem cells from the non-permissive NOD/ShiLtJ and WSB/EiJ strains, enabling easier and broader research possibilities into pluripotency for these and similar non-permissive mouse strains and species.

Here, we report the complete and draft genome sequences of 8 Staphylococcus pseudintermedius isolates, 4 from human bacteremia infections and 4 from canine bacteremia infections. This species is recognized primarily as an important canine pathogen, but it is increasingly being identified in human infections.

Metabolism of 1,3-butadiene, a known human and rodent carcinogen, results in formation of reactive epoxides, a key event in its carcinogenicity. Although mice exposed to 1,3-butadiene present DNA adducts in all tested tissues, carcinogenicity is limited to liver, lung, and lymphoid tissues. Previous studies demonstrated that strain- and tissue-specific epigenetic effects in response to 1,3-butadiene exposure may influence susceptibly to DNA damage and serve as a potential mechanism of tissue-specific carcinogenicity. This study aimed to investigate interindividual variability in the effects of 1,3-butadiene using a population-based mouse model. Male mice from 20 Collaborative Cross strains were exposed to 0 or 635 ppm 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. We evaluated DNA damage and epigenetic effects in target (lung and liver) and nontarget (kidney) tissues of 1,3-butadiene-induced carcinogenesis. DNA damage was assessed by measuring N-7-(2,3,4-trihydroxybut-1-yl)-guanine (THB-Gua) adducts. To investigate global histone modification alterations, we evaluated the trimethylation and acetylation of histones H3 and H4 across tissues. Changes in global cytosine DNA methylation were evaluated from the levels of methylation of LINE-1 and SINE B1 retrotransposons. We quantified the degree of variation across strains, deriving a chemical-specific human variability factor to address population variability in carcinogenic risk, which is largely ignored in current cancer risk assessment practice. Quantitative trait locus mapping identified four candidate genes related to chromatin remodeling whose variation was associated with interstrain susceptibility. Overall, this study uses 1,3-butadiene to demonstrate how the Collaborative Cross mouse population can be used to identify the mechanisms for and quantify the degree of interindividual variability in tissue-specific effects that are relevant to chemically induced carcinogenesis.

PURPOSE: This study aims to determine if intra-ovarian injection of bone marrow–derived mesenchymal stem cells (MSCs) improves or restores ovarian function in aged females. METHODS: Prospective randomized study of eight aged mares and six young mares receiving intra-ovarian injection of MSCs or vehicle. Main outcome measures were antral follicle count and serum anti-Müllerian hormone (AMH) (aged and young mares), and for aged mares, oocyte meiotic and developmental competence; gross and histological ovarian assessment; evaluation of presence of chimerism in recovered granulosa cells and in ovarian tissue samples; and gene expression in ovarian tissue as assessed by RNA sequencing. RESULTS: Injection of MSCs was not associated with significant changes in follicle number, oocyte recovery rate on follicle aspiration, oocyte maturation rate, or blastocyst rate after ICSI in aged mares, or in changes in follicle number in young mares. There were no significant changes in peripheral AMH concentrations, indicating a lack of effect on growing follicles. MSC donor DNA was not recovered in granulosa cells or in ovarian tissue, indicating lack of persistence of injected MSC. RNA sequencing revealed significant differences in gene expression between MSC- and vehicle-injected ovaries. CONCLUSIONS: Intra-ovarian injection of bone marrow–derived MSCs altered gene expression but did not improve ovarian function in aged mares. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s10815-018-1371-6) contains supplementary material, which is available to authorized users.

Summary In seeds and other parts of cultivated, tetraploid cotton (Gossypium hirsutum L.), multicellular groups of cells lysigenously form dark glands containing toxic terpenoids such as gossypol that defend the plant against pests and pathogens. Using RNA‐seq analysis of embryos from near‐isogenic glanded (Gl 2 Gl 2 Gl 3 Gl 3) versus glandless (gl 2 gl 2 gl 3 gl 3) plants, we identified 33 genes that expressed exclusively or at higher levels in embryos just prior to gland formation in glanded plants. Virus‐induced gene silencing against three gene pairs led to significant reductions in the number of glands in the leaves, and significantly lower levels of gossypol and related terpenoids. These genes encode transcription factors and have been designated the ‘Cotton Gland Formation’ (CGF) genes. No sequence differences were found between glanded and glandless cotton for CGF1 and CGF2 gene pairs. The glandless cotton has a transposon insertion within the coding sequence of the GoPGF (synonym CGF3) gene of the A subgenome and extensive mutations in the promoter of D subgenome homeolog. Overexpression of GoPGF (synonym CGF3) led to a dramatic increase in gossypol and related terpenoids in cultured cells, whereas CRISPR/Cas9 knockout of GoPGF (synonym CGF3) genes resulted in glandless phenotype. Taken collectively, the results show that the GoPGF (synonym CGF3) gene plays a critical role in the formation of glands in the cotton plant. Seed‐specific silencing of CGF genes, either individually or in combination, could eliminate glands, thus gossypol, from the cottonseed to render it safe as food or feed for monogastrics.

Multi-parental recombinant inbred populations, such as the Collaborative Cross (CC) mouse genetic reference population, are increasingly being used for analysis of quantitative trait loci (QTL). However specialized analytic software for these complex populations is typically built in R that works only on command-line, which limits the utility of these powerful resources for many users. To overcome analytic limitations, we developed gQTL, a web accessible, simple graphical user interface application based on the DOQTL platform in R to perform QTL mapping using data from CC mice.

Trueperella pyogenes is an opportunistic pathogen that causes suppurative infections in animals including humans. Data on phenotypic and genotypic properties of T. pyogenes isolated from ruminants, particularly goats and sheep, are lacking. We characterized, by phenotypic and genotypic means, T. pyogenes of caprine and ovine origin, and established their phylogenetic relationship with isolates from other ruminants. T. pyogenes isolates ( n = 50) from diagnostic specimens of bovine ( n = 25), caprine ( n = 19), and ovine ( n = 6) origin were analyzed. Overall, variable biochemical activities were observed among the T. pyogenes isolates. The fimbriae-encoding gene, fimE, and neuraminidase-encoding gene, nanH, were, respectively, more frequently detected in the large ( p = 0.0006) and small ( p = 0.0001) ruminant isolates. Moreover, genotype V ( plo/ nanH/ nanP/ fimA/ fimC) was only detected in the caprine and ovine isolates, whereas genotype IX ( plo/ nanP/ fimA/ fimC/ fimE) was solely present in the isolates of bovine origin ( p = 0.0223). The 16S rRNA gene sequences of all T. pyogenes isolates were clustered with the reference T. pyogenes strain ATCC 19411 and displayed a high degree of identity to each other. Our results highlight phenotypic and genotypic diversity among ruminant isolates of T. pyogenes and reinforce the importance of characterization of more clinical isolates to better understand the pathogenesis of this bacterium in different animal species.

Studies of gene expression are common in toxicology and provide important clues to mechanistic understanding of adverse effects of chemicals. Most prior studies have been performed in a single strain or cell line; however, gene expression is heavily influenced by the genetic background, and these genotype-expression differences may be key drivers of inter-individual variation in response to chemical toxicity. In this study, we hypothesized that the genetically-diverse Collaborative Cross (CC) mouse population can be used to gain insight and suggest mechanistic hypotheses for the dose- and genetic background-dependent effects of chemical exposure. This hypothesis was tested using a model liver toxicant trichloroethylene (TCE). Liver transcriptional responses to TCE exposure were evaluated 24 hours after dosing. Transcriptomic dose-responses were examined for both TCE and its major oxidative metabolite trichloroacetic acid (TCA). As expected, peroxisome- and fatty acid metabolism-related pathways were among the most dose-responsive enriched pathways in all strains. However, nearly half of the TCE-induced liver transcriptional perturbation was strain-dependent, with abundant evidence of strain/dose interaction, including in the peroxisomal signaling-associated pathways. These effects were highly concordant between the administered TCE dose and liver levels of TCA. Dose-response analysis of gene expression at the pathway level yielded points of departure similar to those derived from the traditional toxicology studies for both non-cancer and cancer effects. Mapping of expression-genotype-dose relationships revealed some significant associations; however, the effects of TCE on gene expression in liver appear to be highly polygenic traits that are challenging to positionally map. This study highlights the usefulness of mouse population-based studies in assessing inter-individual variation in toxicological responses, but cautions that genetic mapping may be challenging because of the complexity in gene-exposure-dose relationships.

Supplemental material, DS1_JVDI_10.1177_1040638718762479 for Phenotypic and genotypic characteristics of Trueperella pyogenes isolated from ruminants by Artem S. Rogovskyy, Sara Lawhon, Kathryn Kuczmanski, David C. Gillis, Jing Wu, Helen Hurley, Yuliya V. Rogovska, Kranti Konganti, Ching-Yuan Yang, Kay Duncan in Journal of Veterinary Diagnostic Investigation