Search

Because the liver plays a vital role in the clearance of exogenous chemical compounds, it is susceptible to chemical-induced toxicity. Animal-based testing is routinely used to assess the hepatotoxic potential of chemicals. While large-scale high-throughput sequencing data can indicate the genes affected by chemical exposures, we need system-level approaches to interpret these changes. To this end, we developed an updated rat genome-scale metabolic model to integrate large-scale transcriptomics data and utilized a chemical structure similarity-based ToxProfiler tool to identify chemicals that bind to specific toxicity targets to understand the mechanisms of toxicity. We used high-throughput transcriptomics data from a 5-day in vivo study where rats were exposed to different non-toxic and hepatotoxic chemicals at increasing concentrations and investigated how liver metabolism was differentially altered between the non-toxic and hepatotoxic chemical exposures. Our analysis indicated that the genes identified via toxicity target analysis and those mapped to the metabolic model showed a distinct gene expression pattern, with the majority showing upregulation for hepatotoxicants compared to non-toxic chemicals. Similarly, when we mapped the metabolic genes at the pathway level, we identified several pathways in carbohydrate, amino acid, and lipid metabolism that were significantly upregulated for hepatotoxic chemicals. Furthermore, using our system-level integration of gene expression data with the rat metabolic model, we could differentiate metabolites in these pathways that were systematically elevated or suppressed due to hepatotoxic versus non-toxic chemicals. Thus, using our combined approach, we were able to identify a set of potential gene signatures that clearly differentiated liver toxic responses from non-toxic chemicals, which helped us identify potential metabolic pathways and metabolites that are systematically associated with the toxicant exposure., (Published by Oxford University Press 2025.)

Toxicants with the potential to bioaccumulate in humans and animals have long been a cause for concern, particularly due to their association with multiple diseases and organ injuries. Per- and polyfluoro alkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAH) are two such classes of chemicals that bioaccumulate and have been associated with steatosis in the liver. Although PFAS and PAH are classified as chemicals of concern, their molecular mechanisms of toxicity remain to be explored in detail. In this study, we aimed to identify potential mechanisms by which an acute exposure to PFAS and PAH chemicals can induce lipid accumulation and whether the responses depend on chemical class, dose, and sex. To this end, we analyzed mechanisms beginning with the binding of the chemical to a molecular initiating event (MIE) and the consequent transcriptomic alterations. We collated potential MIEs using predictions from our previously developed ToxProfiler tool and from published steatosis adverse outcome pathways. Most of the MIEs are transcription factors, and we collected their target genes by mining the TRRUST database. To analyze the effects of PFAS and PAH on the steatosis mechanisms, we performed a computational MIE-target gene analysis on high-throughput transcriptomic measurements of liver tissue from male and female rats exposed to either a PFAS or PAH. The results showed peroxisome proliferator-activated receptor (PPAR)-α targets to be the most dysregulated, with most of the genes being upregulated. Furthermore, PFAS exposure disrupted several lipid metabolism genes, including upregulation of fatty acid oxidation genes ( Acadm , Acox1 , Cpt2 , Cyp4a1 - 3 ) and downregulation of lipid transport genes ( Apoa1 , Apoa5 , Pltp ). We also identified multiple genes with sex-specific behavior. Notably, the rate-limiting genes of gluconeogenesis ( Pck1 ) and bile acid synthesis ( Cyp7a1 ) were specifically downregulated in male rats compared to female rats, while the rate-limiting gene of lipid synthesis ( Scd ) showed a PFAS-specific upregulation. The results suggest that the PPAR signaling pathway plays a major role in PFAS-induced lipid accumulation in rats. Together, these results show that PFAS exposure induces a sex-specific multi-factorial mechanism involving rate-limiting genes of gluconeogenesis and bile acid synthesis that could lead to activation of an adverse outcome pathway for steatosis., Competing Interests: Authors AH, MA, and VP were employed by The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc. Authors MB-M, DM, DP, ES, and RS were employed by Sciome LLC. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Hari, AbdulHameed, Balik-Meisner, Mav, Phadke, Scholl, Shah, Casey, Auerbach, Wallqvist and Pannala.)

The exposome collectively refers to all exposures, beginning in utero and continuing throughout life, and comprises not only standard environmental exposures such as point source pollution and ozone levels but also exposures from diet, medication, lifestyle factors, stress, and occupation. The exposome interacts with individual genetic and epigenetic characteristics to affect human health and disease, but large-scale studies that characterize the exposome and its relationships with human disease are limited. To address this gap, we used extensive questionnaire data from the diverse North Carolina-based Personalized Environment and Genes Study (PEGS, n = 9, 429) to evaluate exposure associations in relation to common diseases. We performed an exposome-wide association study (ExWAS) to examine single exposure models and their associations with 11 common complex diseases, namely allergic rhinitis, asthma, bone loss, fibroids, high cholesterol, hypertension, iron-deficient anemia, ovarian cysts, lower GI polyps, migraines, and type 2 diabetes. Across diseases, we found associations with lifestyle factors and socioeconomic status as well as asbestos, various dust types, biohazardous material, and textile-related exposures. We also found disease-specific associations such as fishing with lead weights and migraines. To differentiate between a replicated result and a novel finding, we used an AI-based literature search and database tool that allowed us to examine the current literature. We found both replicated findings, especially for lifestyle factors such as sleep and smoking across diseases, and novel findings, especially for occupational exposures and multiple diseases., Competing Interests: None declared., (Published by Oxford University Press 2024.)

The correlations among individual exposures in the exposome, which refers to all exposures an individual encounters throughout life, are important for understanding the landscape of how exposures co-occur, and how this impacts health and disease. Exposome-wide association studies (ExWAS), which are analogous to genome-wide association studies (GWAS), are increasingly being used to elucidate links between the exposome and disease. Despite increased interest in the exposome, tools and publications that characterize exposure correlations and their relationships with human disease are limited, and there is a lack of data and results sharing in resources like the GWAS catalog. To address these gaps, we developed the PEGS Explorer web application to explore exposure correlations in data from the diverse North Carolina-based Personalized Environment and Genes Study (PEGS) that were rigorously calculated to account for differing data types and previously published results from ExWAS. Through globe visualizations, PEGS Explorer allows users to explore correlations between exposures found to be associated with complex diseases. The exposome data used for analysis includes not only standard environmental exposures such as point source pollution and ozone levels but also exposures from diet, medication, lifestyle factors, stress, and occupation. The web application addresses the lack of accessible data and results sharing, a major challenge in the field, and enables users to put results in context, generate hypotheses, and, importantly, replicate findings in other cohorts. PEGS Explorer will be updated with additional results as they become available, ensuring it is an up-to-date resource in exposome science., Competing Interests: None declared., (Published by Oxford University Press 2024.)

The purpose of this study was to determine the toxicological and pharmacokinetic properties of sucralose-6-acetate, a structural analog of the artificial sweetener sucralose. Sucralose-6-acetate is an intermediate and impurity in the manufacture of sucralose, and recent commercial sucralose samples were found to contain up to 0.67% sucralose-6-acetate. Studies in a rodent model found that sucralose-6-acetate is also present in fecal samples with levels up to 10% relative to sucralose which suggest that sucralose is also acetylated in the intestines. A MultiFlow® assay, a high-throughput genotoxicity screening tool, and a micronucleus (MN) test that detects cytogenetic damage both indicated that sucralose-6-acetate is genotoxic. The mechanism of action was classified as clastogenic (produces DNA strand breaks) using the MultiFlow® assay. The amount of sucralose-6-acetate in a single daily sucralose-sweetened drink might far exceed the threshold of toxicological concern for genotoxicity (TTC genotox ) of 0.15 µg/person/day. The RepliGut® System was employed to expose human intestinal epithelium to sucralose-6-acetate and sucralose, and an RNA-seq analysis was performed to determine gene expression induced by these exposures. Sucralose-6-acetate significantly increased the expression of genes associated with inflammation, oxidative stress, and cancer with greatest expression for the metallothionein 1 G gene (MT1G). Measurements of transepithelial electrical resistance (TEER) and permeability in human transverse colon epithelium indicated that sucralose-6-acetate and sucralose both impaired intestinal barrier integrity. Sucralose-6-acetate also inhibited two members of the cytochrome P450 family (CYP1A2 and CYP2C19). Overall, the toxicological and pharmacokinetic findings for sucralose-6-acetate raise significant health concerns regarding the safety and regulatory status of sucralose itself.

Streblospio benedicti is a common marine annelid that has become an important model for developmental evolution. It is the only known example of poecilogony (where two distinct developmental modes occur within a single species) that is due to a heritable difference in egg size. The dimorphic developmental programs and life-histories exhibited in this species depend on differences within the genome, making it an optimal model for understanding the genomic basis of developmental divergence. Studies using S. benedicti have begun to uncover the genetic and genomic principles that underlie developmental uncoupling, but until now they have been limited by the lack of availability of genomic tools. Here, we present an annotated chromosomal-level genome assembly of S. benedicti generated from a combination of Illumina reads, Nanopore long reads, Chicago and Hi-C chromatin interaction sequencing, and a genetic map from experimental crosses. At 701.4 Mb, the S. benedicti genome is the largest annelid genome to date that has been assembled to chromosomal scaffolds. The complete genome of S. benedicti is valuable for functional genomic analyses of development and evolution, as well as phylogenetic comparison within the annelida and the Lophotrochozoa. Despite having two developmental modes, there is no evidence of genome duplication or substantial gene number expansions. Instead, lineage-specific repeats account for much of the expansion of this genome compared with other annelids., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The common green bottle blow fly Lucilia sericata (family, Calliphoridae) is widely used for maggot debridement therapy, which involves the application of sterile maggots to wounds. The larval excretions and secretions are important for consuming necrotic tissue and inhibiting bacterial growth in wounds of patients. Lucilia sericata is also of importance as a pest of sheep and in forensic studies to estimate a postmortem interval. Here we report the assembly of a 565.3 Mb genome from long read PacBio DNA sequencing of genomic DNA. The genome contains 14,704 predicted protein coding genes and 1709 non-coding genes. Targeted annotation and transcriptional analyses identified genes that are highly expressed in the larval salivary glands (secretions) and Malpighian tubules (excretions) under normal growth conditions and following heat stress. The genomic resources will underpin future genetic studies and in development of engineered strains for genetic control of L. sericata and for biotechnology-enhanced maggot therapy., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Of the 12 full-length feline leukocyte antigen class I (FLAI) loci, 3 are presumed to be classical: FLAI-E, FLAI-H, and FLAI-K. As diversity is a class Ia hallmark, multi-allelism is an important surrogate supporting a classical designation, in the absence of direct demonstration of T-cell restriction. Conversely, limited polymorphism at an expressed locus suggests regulation of immune effectors with invariant receptors, and non-classical status. FLAI-A, FLAI-J, FLAI-L, and FLAI-O are putative class Ib genes in cats. For both classes, identifying prevalent variants across outbred populations can illuminate specific genotypes to be prioritized for immune studies, as shared alleles direct shared responses. Since variation is concentrated in exons 2 and 3, which encode the antigen-binding domains, partial-length cloning/sequencing can be used for allele discovery, but is laborious and occasionally ambiguous. Here we develop a targeted approach to FLAI genotyping, using the single-molecule real-time (SMRT) platform, which allows full-length (3.4-kb) reads without assembly. Consensus sequences matched full-length Sanger references. Thirty-one new class Ia genes were found in 17 cats. Alleles segregated strongly by loci, and the origins of formerly difficult-to-assign sequences were resolved. Although not targeted, FLAI-L and FLAI-J, and the pseudogene FLAI-F, were also returned. Eighteen class Ib alleles were identified. Diversity was restricted and outside hypervariable regions. Both class Ib genes were transcriptionally active. Novel alternative splicing of FLAI-L was observed. SMRT sequencing of FLAI amplicons is useful for full-length genotyping at feline class Ia loci. High-throughput sequencing could allow highly accurate allele surveys in large cat cohorts., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

We report the complete circular genome sequences of six Lactobacillus strains and their plasmids, if any, from the fecal material of quarter horses at different ages., (Copyright © 2020 Meinders et al.)

We report the complete circular genome sequences of Lactobacillus crispatus strain C25, its plasmid, and Lactobacillus animalis strain P38; both strains were isolated from the cecum of 4-week-old chickens. These isolates represent potential probiotic strains for poultry., (Copyright © 2020 Hassan et al.)

The New World Screwworm fly, Cochliomyia hominivorax, is a major pest of livestock in South America and Caribbean. However, few genomic resources have been available for this species. A genome of 534 Mb was assembled from long read PacBio DNA sequencing of DNA from a highly inbred strain. Analysis of molecular evolution identified 40 genes that are likely under positive selection. Developmental RNA-seq analysis identified specific genes associated with each stage. We identify and analyze the expression of genes that are likely important for host-seeking behavior (chemosensory), development of larvae in open wounds in warm-blooded animals (heat shock protein, immune response) and for building transgenic strains for genetic control programs including gene drive (sex determination, germline). This study will underpin future experiments aimed at understanding the parasitic lifestyle of the screwworm fly and greatly facilitate future development of strains for efficient systems for genetic control of screwworm.

Background: Cytauxzoonosis is a disease of felids in North America caused by the tick-transmitted apicomplexan parasite Cytauxzoon felis . Cytauxzoonosis is particularly virulent for domestic cats, but no vaccine currently exists. The parasite cannot be cultivated in vitro, presenting a significant limitation for vaccine development., Methods: Recent sequencing of the C. felis genome has identified over 4300 putative protein-encoding genes. From this pool we constructed a protein microarray containing 673 putative C. felis proteins. This microarray was probed with sera from C. felis -infected and naïve cats to identify differentially reactive antigens which were incorporated into two expression library vaccines, one polyvalent and one monovalent. We assessed the efficacy of these vaccines to prevent of infection and/or disease in a tick-challenge model., Results: Probing of the protein microarray resulted in identification of 30 differentially reactive C. felis antigens that were incorporated into the two expression library vaccines. However, expression library immunization failed to prevent infection or disease in cats challenged with C. felis ., Conclusions: Protein microarray facilitated high-throughput identification of novel antigens, substantially increasing the pool of characterized C. felis antigens. These antigens should be considered for development of C. felis vaccines, diagnostics, and therapeutics.

Four bacterial strains identified as members of the Acidovorax genus were isolated from two geographically distinct but similarly contaminated soils in North Carolina, USA, characterized, and their genomes sequenced. Their 16S rRNA genes were highly similar to those previously recovered during stable-isotope probing (SIP) of one of the soils with the polycyclic aromatic hydrocarbon (PAH) phenanthrene. Heterotrophic growth of all strains occurred with a number of organic acids, as well as phenanthrene, but no other tested PAHs. Optimal growth occurred aerobically under mesophilic temperature, neutral pH, and low salinity conditions. Predominant fatty acids were C 16:1 ω7c/C 16:1 ω6c, C 16:0 , and C 18:1 ω7c, and were consistent with the genus. Genomic G+C contents ranged from 63.6 to 64.2%. A combination of whole genome comparisons and physiological analyses indicated that these four strains likely represent a single species within the Acidovorax genus. Chromosomal genes for phenanthrene degradation to phthalate were nearly identical to highly conserved regions in phenanthrene-degrading Delftia, Burkholderia, Alcaligenes, and Massilia species in regions flanked by transposable or extrachromosomal elements. The lower degradation pathway for phenanthrene metabolism was inferred by comparisons to described genes and proteins. The novel species Acidovorax carolinensis sp. nov. is proposed, comprising the four strains described in this study with strain NA3 T as the type strain (=LMG 30136, =DSM 105008)., (Copyright © 2018 Elsevier GmbH. All rights reserved.)

Carboxylesterase 1 (CES1) metabolizes methylphenidate and other drugs. CES1 gene variation only partially explains pharmacokinetic (PK) variability. Biomarkers predicting the PKs of drugs metabolized by CES1 are needed. We identified lipids in plasma from 44 healthy subjects that correlated with CES1 activity as determined by PK parameters of methylphenidate including a ceramide (q value = 0.001) and a phosphatidylcholine (q value = 0.005). Carriers of the CES1 143E allele had decreased methylphenidate metabolism and altered concentration of this phosphatidylcholine (q value = 0.040) and several high polyunsaturated fatty acid lipids (PUFAs). The half-maximal inhibitory concentration (IC 50 ) values of chenodeoxycholate and taurocholate were 13.55 and 19.51 μM, respectively, consistent with a physiological significance. In silico analysis suggested that bile acid inhibition of CES1 involved both binding to the active and superficial sites of the enzyme. We initiated identification of metabolites predicting PKs of drugs metabolized by CES1 and suggest lipids to regulate or be regulated by this enzyme., (© 2018 The Authors CPT: Pharmacometrics & Systems Pharmacology published by Wiley Periodicals, Inc. on behalf of American Society for Clinical Pharmacology and Therapeutics.)

Background: Modern societies are exposed to vast numbers of potentially hazardous chemicals. Despite demonstrated linkages between chemical exposure and severe health effects, there are limited, often conflicting, data on how adverse health effects of exposure differ across individuals., Objectives: We tested the hypothesis that population variability in response to certain chemicals could elucidate a role for gene-environment interactions (GxE) in differential susceptibility., Methods: High-throughput screening (HTS) data on thousands of chemicals in genetically heterogeneous zebrafish were leveraged to identify a candidate chemical (Abamectin) with response patterns indicative of population susceptibility differences. We tested the prediction by generating genome-wide sequence data for 276 individual zebrafish displaying susceptible (Affected) vs. resistant (Unaffected) phenotypes following identical chemical exposure., Results: We found GxE associated with differential susceptibility in the sox7 promoter region and then confirmed gene expression differences between phenotypic response classes., Conclusions: The results for Abamectin in zebrafish demonstrate that GxE associated with naturally occurring, population genetic variation play a significant role in mediating individual response to chemical exposure. https://doi.org/10.1289/EHP2662.

It has been hypothesized that the Drosophila 4 th  chromosome is derived from an ancient X chromosome [1]. In the Australian sheep blowfly, Lucilia cuprina, the heterochromatic X chromosome contains few active genes and orthologs of Drosophila X-linked genes are autosomal. Of 8 X-linked genes identified previously in L. cuprina, 6 were orthologs of Drosophila 4 th -chromosome genes [2]. The X-linked genes were expressed equally in males and females. Here we identify an additional 51 X-linked genes and show that most are dosage compensated. Orthologs of 49 of the 59 X-linked genes are on the 4 th  chromosome in D. melanogaster. Because painting of fourth (Pof) is important for expression of Drosophila 4 th -chromosome genes [3], we used Cas9 to make a loss-of-function knockin mutation in an L. cuprina Pof ortholog we call no blokes (nbl). Homozygous nbl males derived from homozygous nbl mothers die at the late pupal stage. Homozygous nbl females are viable, fertile, and live longer than heterozygous nbl females. RNA expression of most X-linked genes was reduced in homozygous nbl male pupae and to a lesser extent in nbl females compared to heterozygous siblings. The results suggest that NBL could be important for X chromosome dosage compensation in L. cuprina. NBL may also facilitate gene expression in the heterochromatic environment of the X chromosome in both sexes. This study supports the hypothesis on the origin of the Drosophila 4 th chromosome and that a POF-like protein was required for normal gene expression on the ancient X chromosome., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Toxicological and pharmacological researchers have seized upon the many benefits of zebrafish, including the short generation time, well-characterized development, and early maturation as clear embryos. A major difference from many model organisms is that standard husbandry practices in zebrafish are designed to maintain population diversity. While this diversity is attractive for translational applications in human and ecological health, it raises critical questions on how interindividual genetic variation might contribute to chemical exposure or disease susceptibility differences. Findings from pooled samples of zebrafish support this supposition of diversity yet cannot directly measure allele frequencies for reference versus alternate alleles. Using the Tanguay lab Tropical 5D zebrafish line (T5D), we performed whole genome sequencing on a large group (n = 276) of individual zebrafish embryos. Paired-end reads were collected on an Illumina 3000HT, then aligned to the most recent zebrafish reference genome (GRCz10). These data were used to compare observed population genetic variation across species (humans, mice, zebrafish), then across lines within zebrafish. We found more single nucleotide polymorphisms (SNPs) in T5D than have been reported in SNP databases for any of the WIK, TU, TL, or AB lines. We theorize that some subset of the novel SNPs may be shared with other zebrafish lines but have not been identified in other studies due to the limitations of capturing population diversity in pooled sequencing strategies. We establish T5D as a model that is representative of diversity levels within laboratory zebrafish lines and demonstrate that experimental design and analysis can exert major effects when characterizing genetic diversity in heterogeneous populations.

Trichloroethylene (TCE) and tetrachloroethylene (PCE) are ubiquitous environmental contaminants and occupational health hazards. Recent health assessments of these agents identified several critical data gaps, including lack of comparative analysis of their effects. This study examined liver and kidney effects of TCE and PCE in a dose-response study design. Equimolar doses of TCE (24, 80, 240, and 800 mg/kg) or PCE (30, 100, 300, and 1000 mg/kg) were administered by gavage in aqueous vehicle to male B6C3F1/J mice. Tissues were collected 24 h after exposure. Trichloroacetic acid (TCA), a major oxidative metabolite of both compounds, was measured and RNA sequencing was performed. PCE had a stronger effect on liver and kidney transcriptomes, as well as greater concentrations of TCA. Most dose-responsive pathways were common among chemicals/tissues, with the strongest effect on peroxisomal β-oxidation. Effects on liver and kidney mitochondria-related pathways were notably unique to PCE. We performed dose-response modeling of the transcriptomic data and compared the resulting points of departure (PODs) to those for apical endpoints derived from long-term studies with these chemicals in rats, mice, and humans, converting to human equivalent doses using tissue-specific dosimetry models. Tissue-specific acute transcriptional effects of TCE and PCE occurred at human equivalent doses comparable to those for apical effects. These data are relevant for human health assessments of TCE and PCE as they provide data for dose-response analysis of the toxicity mechanisms. Additionally, they provide further evidence that transcriptomic data can be useful surrogates for in vivo PODs, especially when toxicokinetic differences are taken into account., (© The Author 2017. Published by Oxford University Press on behalf of the Society of Toxicology.All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

To replicate efficiently, viruses must create favorable cell conditions and overcome cell antiviral responses. We previously reported that the reovirus protein μ2 from strain T1L, but not strain T3D, represses one antiviral response: alpha/beta interferon signaling. We report here that T1L, but not T3D, μ2 localizes to nuclear speckles, where it forms a complex with the mRNA splicing factor SRSF2 and alters its subnuclear localization. Reovirus replicates in cytoplasmic viral factories, and there is no evidence that reovirus genomic or messenger RNAs are spliced, suggesting that T1L μ2 might target splicing of cell RNAs. Indeed, RNA sequencing revealed that reovirus T1L, but not T3D, infection alters the splicing of transcripts for host genes involved in mRNA posttranscriptional modifications. Moreover, depletion of SRSF2 enhanced reovirus replication and cytopathic effect, suggesting that T1L μ2 modulation of splicing benefits the virus. This provides the first report of viral antagonism of the splicing factor SRSF2 and identifies the viral protein that determines strain-specific differences in cell RNA splicing. IMPORTANCE Efficient viral replication requires that the virus create favorable cell conditions. Many viruses accomplish this by repressing specific antiviral responses. We demonstrate here that some mammalian reoviruses, RNA viruses that replicate strictly in the cytoplasm, express a protein variant that localizes to nuclear speckles, where it targets a cell mRNA splicing factor. Infection with a reovirus strain that targets this splicing factor alters splicing of cell mRNAs involved in the maturation of many other cell mRNAs. Depletion of this cell splicing factor enhances reovirus replication and cytopathic effect. Our results provide the first evidence of viral antagonism of this splicing factor and suggest that downstream consequences to the cell are global and benefit the virus., (Copyright © 2017 American Society for Microbiology.)

The taxonomy of the order Piroplasmida, which includes a number of clinically and economically relevant organisms, is a hotly debated topic amongst parasitologists. Three genera (Babesia, Theileria, and Cytauxzoon) are recognized based on parasite life cycle characteristics, but molecular phylogenetic analyses of 18S sequences have suggested the presence of five or more distinct Piroplasmida lineages. Despite these important advancements, a few studies have been unable to define the taxonomic relationships of some organisms (e.g. C. felis and T. equi) with respect to other Piroplasmida. Additional evidence from mitochondrial genome sequences and synteny should aid in the inference of Piroplasmida phylogeny and resolution of taxonomic uncertainties. In this study, we have amplified, sequenced, and annotated seven previously uncharacterized mitochondrial genomes (Babesia canis, Babesia vogeli, Babesia rossi, Babesia sp. Coco, Babesia conradae, Babesia microti-like sp., and Cytauxzoon felis) and identified additional ribosomal fragments in ten previously characterized mitochondrial genomes. Phylogenetic analysis of concatenated mitochondrial and 18S sequences as well as cox1 amino acid sequence identified five distinct Piroplasmida groups, each of which possesses a unique mitochondrial genome structure. Specifically, our results confirm the existence of four previously identified clades (B. microti group, Babesia sensu stricto, Theileria equi, and a Babesia sensu latu group that includes B. conradae) while supporting the integration of Theileria and Cytauxzoon species into a single fifth taxon. Although known biological characteristics of Piroplasmida corroborate the proposed phylogeny, more investigation into parasite life cycles is warranted to further understand the evolution of the Piroplasmida. Our results provide an evolutionary framework for comparative biology of these important animal and human pathogens and help focus renewed efforts toward understanding the phylogenetic relationships within the group., Competing Interests: The authors have declared that no competing interests exist.

Foodborne infections caused by Salmonella enterica serovars are a significant problem worldwide. Presented here is the genome sequence of the nontyphoidal S. enterica serovar Typhimurium mutant strain NC983, a potential vaccine candidate., (Copyright © 2016 Troxell et al.)

The bacterial strain TR3.2, representing a novel deeply branching lineage within the Gammaproteobacteria, was isolated and its genome sequenced. This isolate is the first cultivated representative of the previously described "Pyrene Group 2" (PG2) and represents a variety of environmental sequences primarily associated with petrochemical contamination and aromatic hydrocarbon degradation., (Copyright © 2016 Singleton et al.)

An early-stage, population-wide biomarker for ovarian cancer (OVC) is essential to reverse its high mortality rate. Aberrant glycosylation by OVC has been reported, but studies have yet to identify an N-glycan with sufficiently high specificity. We curated a human biorepository of 82 case-control plasma samples, with 27%, 12%, 46%, and 15% falling across stages I-IV, respectively. For relative quantitation, glycans were analyzed by the individuality normalization when labeling with glycan hydrazide tags (INLIGHT) strategy for enhanced electrospray ionization, MS/MS analysis. Sixty-three glycan cancer burden ratios (GBRs), defined as the log10 ratio of the case-control extracted ion chromatogram abundances, were calculated above the limit of detection. The final GBR models, built using stepwise forward regression, included three significant terms: OVC stage, normalized mean GBR, and tag chemical purity; glycan class, fucosylation, or sialylation were not significant variables. After Bonferroni correction, seven N-glycans were identified as significant (p < 0.05), and after false discovery rate correction, an additional four glycans were determined to be significant (p < 0.05), with one borderline (p = 0.05). For all N-glycans, the vectors of the effects from stages II-IV were sequentially reversed, suggesting potential biological changes in OVC morphology or in host response.

Enterococcus cecorum (EC) is the dominant enteric commensal of adult chickens and contributes to the gut consortia of many avian and mammalian species. While EC infection is an uncommon zoonosis, like other enterococcal species it can cause life-threating nosocomial infection in people. In contrast to other enterococci which are considered opportunistic pathogens, emerging pathogenic strains of EC cause outbreaks of musculoskeletal disease in broiler chickens. Typical morbidity and mortality is comparable to other important infectious diseases of poultry. In molecular epidemiologic studies, pathogenic EC strains were found to be genetically clonal. These findings suggested acquisition of specific virulence determinants by pathogenic EC. To identify divergent genomic features and acquired virulence determinants in pathogenic EC; comparative genomic analysis was performed on genomes of 3 pathogenic and 3 commensal strains of EC. Pathogenic isolates had smaller genomes with a higher GC content, and they demonstrated large regions of synteny compared to commensal isolates. A molecular phylogenetic analysis demonstrated sequence divergence in pathogenic EC genomes. At a threshold of 98% identity, 414 predicted proteins were identified that were highly conserved in pathogenic EC but not in commensal EC. Among these, divergent CRISPR-cas defense loci were observed. In commensal EC, the type IIA arrangement typical for enterococci was present; however, pathogenic EC had a type IC locus, which is novel in enterococci but commonly observed in streptococci. Potential mediators of virulence identified in this analysis included a polysaccharide capsular locus similar to that recently described for E. faecium, an epa-like locus, and cell wall associated proteins which may bind host extracellular matrix. This analysis identified specific genomic regions, coding sequences, and predicted proteins which may be related to the divergent evolution and increased virulence of emerging pathogenic strains of EC.

A polycyclic aromatic hydrocarbon-degrading bacterium designated strain Ca6, a member of the family Rhodocyclaceae and a representative of the uncharacterized pyrene group 1 (PG1), was isolated and its genome sequenced. The presence of several genes suspected to be associated with PG1 was confirmed, and additional genes for aromatic compound metabolism were detected., (Copyright © 2015 Singleton et al.)

Cytauxzoonosis is an emerging infectious disease of domestic cats (Felis catus) caused by the apicomplexan protozoan parasite Cytauxzoon felis. The growing epidemic, with its high morbidity and mortality points to the need for a protective vaccine against cytauxzoonosis. Unfortunately, the causative agent has yet to be cultured continuously in vitro, rendering traditional vaccine development approaches beyond reach. Here we report the use of comparative genomics to computationally and experimentally interpret the C. felis genome to identify a novel candidate vaccine antigen for cytauxzoonosis. As a starting point we sequenced, assembled, and annotated the C. felis genome and the proteins it encodes. Whole genome alignment revealed considerable conserved synteny with other apicomplexans. In particular, alignments with the bovine parasite Theileria parva revealed that a C. felis gene, cf76, is syntenic to p67 (the leading vaccine candidate for bovine theileriosis), despite a lack of significant sequence similarity. Recombinant subdomains of cf76 were challenged with survivor-cat antiserum and found to be highly seroreactive. Comparison of eleven geographically diverse samples from the south-central and southeastern USA demonstrated 91-100% amino acid sequence identity across cf76, including a high level of conservation in an immunogenic 226 amino acid (24 kDa) carboxyl terminal domain. Using in situ hybridization, transcription of cf76 was documented in the schizogenous stage of parasite replication, the life stage that is believed to be the most important for development of a protective immune response. Collectively, these data point to identification of the first potential vaccine candidate antigen for cytauxzoonosis. Further, our bioinformatic approach emphasizes the use of comparative genomics as an accelerated path to developing vaccines against experimentally intractable pathogens.

Root-knot nematodes (Meloidogyne spp.) cause major yield losses to many of the world's crops, but efforts to understand how these pests recognize and interact with their hosts have been hampered by a lack of genetic resources. Starting with progeny of a cross between inbred strains (VW8 and VW9) of Meloidogyne hapla that differed in host range and behavioral traits, we exploited the novel, facultative meiotic parthenogenic reproductive mode of this species to produce a genetic linkage map. Molecular markers were derived from SNPs identified between the sequenced and annotated VW9 genome and de novo sequence of VW8. Genotypes were assessed in 183 F2 lines. The colinearity of the genetic and physical maps supported the veracity of both. Analysis of local crossover intervals revealed that the average recombination rate is exceptionally high compared with that in other metazoans. In addition, F2 lines are largely homozygous for markers flanking crossover points, and thus resemble recombinant inbred lines. We suggest that the unusually high recombination rate may be an adaptation to generate within-population genetic diversity in this organism. This work presents the most comprehensive linkage map of a parasitic nematode to date and, together with genomic and transcript sequence resources, empowers M. hapla as a tractable model. Alongside the molecular map, these progeny lines can be used for analyses of genome organization and the inheritance of phenotypic traits that have key functions in modulating parasitism, behavior, and survival and for the eventual identification of the responsible genes.

Background: The striped bass and its relatives (genus Morone) are important fisheries and aquaculture species native to estuaries and rivers of the Atlantic coast and Gulf of Mexico in North America. To open avenues of gene expression research on reproduction and breeding of striped bass, we generated a collection of expressed sequence tags (ESTs) from a complementary DNA (cDNA) library representative of their ovarian transcriptome., Results: Sequences of a total of 230,151 ESTs (51,259,448 bp) were acquired by Roche 454 pyrosequencing of cDNA pooled from ovarian tissues obtained at all stages of oocyte growth, at ovulation (eggs), and during preovulatory atresia. Quality filtering of ESTs allowed assembly of 11,208 high-quality contigs ≥ 100 bp, including 2,984 contigs 500 bp or longer (average length 895 bp). Blastx comparisons revealed 5,482 gene orthologues (E-value < 10-3), of which 4,120 (36.7% of total contigs) were annotated with Gene Ontology terms (E-value < 10-6). There were 5,726 remaining unknown unique sequences (51.1% of total contigs). All of the high-quality EST sequences are available in the National Center for Biotechnology Information (NCBI) Short Read Archive (GenBank: SRX007394). Informative contigs were considered to be abundant if they were assembled from groups of ESTs comprising ≥ 0.15% of the total short read sequences (≥ 345 reads/contig). Approximately 52.5% of these abundant contigs were predicted to have predominant ovary expression through digital differential display in silico comparisons to zebrafish (Danio rerio) UniGene orthologues. Over 1,300 Gene Ontology terms from Biological Process classes of Reproduction, Reproductive process, and Developmental process were assigned to this collection of annotated contigs., Conclusions: This first large reference sequence database available for the ecologically and economically important temperate basses (genus Morone) provides a foundation for gene expression studies in these species. The predicted predominance of ovary gene expression and assignment of directly relevant Gene Ontology classes suggests a powerful utility of this dataset for analysis of ovarian gene expression related to fundamental questions of oogenesis. Additionally, a high definition Agilent 60-mer oligo ovary 'UniClone' microarray with 8 × 15,000 probe format has been designed based on this striped bass transcriptome (eArray Group: Striper Group, Design ID: 029004).

Sequencing of expressed genes has shown that nematodes, particularly the plant-parasitic nematodes, have genes purportedly acquired from other kingdoms by horizontal gene transfer. The prevailing orthodoxy is that such transfer has been a driving force in the evolution of niche specificity, and a recent paper in BMC Evolutionary Biology that presents a detailed phylogenetic analysis of cellulase genes in the free-living nematode Pristionchus pacificus at the species, genus and family levels substantiates this hypothesis.

On the basis of the complete genome sequence of the root-knot nematode Melodogyne hapla, we have deduced and annotated the entire proteome of this plant-parasite to create a database of 14,420 proteins. We have made this database, termed HapPep3, available from the Superfamily repository of model organism proteomes (http://supfam.mrc-lmb.cam.ac.uk/SUPERFAMILY). To experimentally confirm the HapPep3 assignments using proteomics, we applied a data-independent LC/MS(E) analysis to M. hapla protein extracts fractionated by SDS-PAGE. A total of 516 nonredundant proteins were identified with an average of 9 unique peptides detected per protein. Some proteins, including examples with complex gene organization, were defined by more than 20 unique peptide matches, thus, providing experimental confirmation of computational predictions of intron/exon structures. On the basis of comparisons of the broad physicochemical properties of the experimental and computational proteomes, we conclude that the identified proteins reflect a true and unbiased sampling of HapPep3. Conversely, HapPep3 appears to broadly cover the protein space able to be experimentally sampled. To estimate the false discovery rate, we queried human, plant, and bacterial databases for matches to the LC/MS(E)-derived peptides, revealing fewer than 1% of matches, most of which were to highly conserved proteins. To provide a functional comparison of the acquired and deduced proteomes, each was subjected to higher order annotation, including comparisons of Gene Ontology, protein domains, signaling, and localization predictions, further indicating concordance, although those proteins that did deviate seem to be highly significant. Approximately 20% of the experimentally sampled proteome was predicted to be secreted, and thus potentially play a role at the host-parasite interface. We examined reference pathways to determine the extent of proteome similarity of M. hapla to that of the free-living nematode, Caenorhabditis elegans, revealing significant similarities and differences. Collectively, the analyzed protein set provides an initial foundation to experimentally dissect the basis of plant parasitism by M. hapla.

Root-knot nematode (RKN; Meloidogyne spp.) is a major crop pathogen worldwide. Effective resistance exists for a few plant species, including that conditioned by Mi in tomato (Solanum lycopersicum). We interrogated the root transcriptome of the resistant (Mi+) and susceptible (Mi-) cultivars 'Motelle' and 'Moneymaker,' respectively, during a time-course infection by the Mi-susceptible RKN species Meloidogyne incognita and the Mi-resistant species Meloidogyne hapla. In the absence of RKN infection, only a single significantly regulated gene, encoding a glycosyltransferase, was detected. However, RKN infection influenced the expression of broad suites of genes; more than half of the probes on the array identified differential gene regulation between infected and uninfected root tissue at some stage of RKN infection. We discovered 217 genes regulated during the time of RKN infection corresponding to establishment of feeding sites, and 58 genes that exhibited differential regulation in resistant roots compared to uninfected roots, including the glycosyltransferase. Using virus-induced gene silencing to silence the expression of this gene restored susceptibility to M. incognita in 'Motelle,' indicating that this gene is necessary for resistance to RKN. Collectively, our data provide a picture of global gene expression changes in roots during compatible and incompatible associations with RKN, and point to candidates for further investigation.

Sequence-based phylogenetic analyses typically are based on a small number of character sets and report gene trees which may not reflect the true species tree. We employed an EST mining strategy to suppress such incongruencies, and recovered the most robust phylogeny for five species of plant-parasitic nematode (Meloidogyne arenaria, M. chitwoodi, M. hapla, M. incognita, and M. javanica), three closely related tylenchid taxa (Heterodera glycines, Globodera pallida, and G. rostochiensis) and a distant taxon, Caenorhabditis elegans. Our multiple-gene approach is based on sampling more than 80,000 publicly available tylenchid EST sequences to identify phylum-wide orthologues. Bayesian inference, minimum evolution, maximum likelihood and protein distance methods were employed for phylogenetic reconstruction and hypothesis tests were constructed to elucidate differential selective pressures across the phylogeny for each gene. Our results place M. incognita and M. javanica as sister taxa, with M. arenaria as the next closely related nematode. Significant differences in selective pressure were revealed for some genes under some hypotheses, though all but one gene are exclusively under purifying selection, indicating conservation across the orthologous groups. This EST-based multi-gene analysis is a first step towards accomplishing genome-wide coverage for tylenchid evolutionary analyses.

Background: Published accounts of horizontally acquired genes in plant-parasitic nematodes have not been the result of a specific search for gene transfer per se, but rather have emerged from characterization of individual genes. We present a method for a high-throughput genome screen for horizontally acquired genes, illustrated using expressed sequence tag (EST) data from three species of root-knot nematode, Meloidogyne species., Results: Our approach identified the previously postulated horizontally transferred genes and revealed six new candidates. Screening was partially dependent on sequence quality, with more candidates identified from clustered sequences than from raw EST data. Computational and experimental methods verified the horizontal gene transfer candidates as bona fide nematode genes. Phylogenetic analysis implicated rhizobial ancestors as donors of horizontally acquired genes in Meloidogyne., Conclusions: High-throughput genomic screening is an effective way to identify horizontal gene transfer candidates. Transferred genes that have undergone amelioration of nucleotide composition and codon bias have been identified using this approach. Analysis of these horizontally transferred gene candidates suggests a link between horizontally transferred genes in Meloidogyne and parasitism.