Interspecific effects of 4A-DNT (4-amino-2,6-dinitrotoluene) and RDX (1,3,5-trinitro-1,3,5-triazine) in Japanese quail, Northern bobwhite, and Zebra finch

The seizure-inducing plastic explosive RDX inhibits the α1β2γ2 GABAA receptor

Objective

Royal demolition explosive (RDX) can induce seizures in wildlife and humans following release into the environment or after voluntary consumption. During the Vietnam War, RDX intoxication was the most common cause of generalized seizures in US service personnel, and in some sections of the armed forces, eating of RDX has continued as “a dare” to this day. After its mechanism of action was long unknown, RDX was recently shown to be a GABA_A receptor antagonist. We here determined the GABA_A receptor subtype‐selectivity of RDX and mapped its functional binding site.

Methods

We used whole‐cell patch‐clamp to determine the potency of RDX on 10 recombinantly expressed GABA_A receptors and mapped the RDX binding site using a combination of Rosetta molecular modeling and site‐directed mutagenesis.

Results

RDX was found to reversibly inhibit the α1β2γ2 GABA_A receptor with an IC₅₀ of 23 μmol/L (95% CI 15.1–33.3 μmol/L), whereas α4 and α6 containing GABA_A receptor combinations were 4–10‐fold less sensitive. RDX is binding to the noncompetitive antagonist (NCA) site in the pore. In a molecular model based on the cryo‐EM structure of the resting state of the α1β2γ2 receptor, RDX forms two hydrogen bonds with the threonines at the T6’ ring and makes hydrophobic interactions with the valine and alanine in 2′ position of the α1 or β2 subunits.

Interpretation

Our findings characterize the mechanism of action of RDX at the atomistic level and suggest that RDX‐induced seizures should be susceptible to treatment with GABA_A modulating drugs such as benzodiazepines, barbiturates, propofol, or neurosteroids.

Example of Adverse Outcome Pathway Concept Enabling Genome-to-Phenome Discovery in Toxicology

The following article represents a mini-review of an intensive 10-year progression of genome-to-phenome (G2P) discovery guided by the adverse outcome pathway (AOP) concept. This example is presented as a means to stimulate crossover of this toxicological concept to enhance G2P discovery within the broader biological sciences community. The case study demonstrates the benefits of the AOP approach for establishing causal linkages across multiple levels of biological organization ultimately linking molecular initiation (often at the genomic scale) to organism-level phenotypes of interest. The case study summarizes a US military effort to identify the mechanism(s) underlying toxicological phenotypes of lethargy and weight loss in response to nitroaromatic munitions exposures, such as 2,4,6-trinitrotoluene. Initial key discoveries are described including the toxicogenomic results that nitrotoluene exposures inhibited expression within the peroxisome proliferator activated receptor α (PPARα) pathway. We channeled the AOP concept to test the hypothesis that inhibition of PPARα signaling in nitrotoluene exposures impacted lipid metabolic processes, thus affecting systemic energy budgets, ultimately resulting in body weight loss. Results from a series of transcriptomic, proteomic, lipidomic, in vitro PPARα nuclear signaling, and PPARα knock-out investigations ultimately supported various facets of this hypothesis. Given these results, we next proceeded to develop a formalized AOP description of PPARα antagonism leading to body weight loss. This AOP was refined through intensive literature review and polished through multiple rounds of peer-review leading to final international acceptance as an Organisation for Economic Cooperation and Development-approved AOP. Briefly, that AOP identifies PPARα antagonist binding as the molecular initiating event (MIE) leading to a series of key events including inhibition of nuclear transactivation for genes controlling lipid metabolism and ketogenesis, inhibition of fatty acid beta-oxidation and ketogenesis dynamics, negative energy budget, and ultimately the adverse outcome (AO) of body-weight loss. Given that the PPARα antagonism MIE represented a reliable indicator of AO progression within the pathway, a phylogenetic analysis was conducted which indicated that PPARα amino acid relatedness generally tracked species relatedness. Additionally, PPARα amino acid relatedness analysis using the Sequence Alignment to Predict Across Species Susceptibility predicted susceptibility to the MIE across vertebrates providing context for AOP extrapolation across species. Overall, we hope this illustrative example of how the AOP concept has benefited toxicology sows a seed within the broader biological sciences community to repurpose the concept to facilitate enhanced G2P discovery in biology.

Systems toxicology identifies mechanistic impacts of 2-amino-4,6-dinitrotoluene (2A-DNT) exposure in Northern Bobwhite

Background

A systems toxicology investigation comparing and integrating transcriptomic and proteomic results was conducted to develop holistic effects characterizations for the wildlife bird model, Northern bobwhite (Colinus virginianus) dosed with the explosives degradation product 2-amino-4,6-dinitrotoluene (2A-DNT). A subchronic 60d toxicology bioassay was leveraged where both sexes were dosed via daily gavage with 0, 3, 14, or 30 mg/kg-d 2A-DNT. Effects on global transcript expression were investigated in liver and kidney tissue using custom microarrays for C. virginianus in both sexes at all doses, while effects on proteome expression were investigated in liver for both sexes and kidney in males, at 30 mg/kg-d.

Results

As expected, transcript expression was not directly indicative of protein expression in response to 2A-DNT. However, a high degree of correspondence was observed among gene and protein expression when investigating higher-order functional responses including statistically enriched gene networks and canonical pathways, especially when connected to toxicological outcomes of 2A-DNT exposure. Analysis of networks statistically enriched for both transcripts and proteins demonstrated common responses including inhibition of programmed cell death and arrest of cell cycle in liver tissues at 2A-DNT doses that caused liver necrosis and death in females. Additionally, both transcript and protein expression in liver tissue was indicative of induced phase I and II xenobiotic metabolism potentially as a mechanism to detoxify and excrete 2A-DNT. Nuclear signaling assays, transcript expression and protein expression each implicated peroxisome proliferator-activated receptor (PPAR) nuclear signaling as a primary molecular target in the 2A-DNT exposure with significant downstream enrichment of PPAR-regulated pathways including lipid metabolic pathways and gluconeogenesis suggesting impaired bioenergetic potential.

Conclusion

Although the differential expression of transcripts and proteins was largely unique, the consensus of functional pathways and gene networks enriched among transcriptomic and proteomic datasets provided the identification of many critical metabolic functions underlying 2A-DNT toxicity as well as impaired PPAR signaling, a key molecular initiating event known to be affected in di- and trinitrotoluene exposures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1798-4) contains supplementary material, which is available to authorized users.

A draft de novo genome assembly for the northern bobwhite (Colinus virginianus) reveals evidence for a rapid decline in effective population size beginning in the Late Pleistocene

Wild populations of northern bobwhites (Colinus virginianus; hereafter bobwhite) have declined across nearly all of their U.S. range, and despite their importance as an experimental wildlife model for ecotoxicology studies, no bobwhite draft genome assembly currently exists. Herein, we present a bobwhite draft de novo genome assembly with annotation, comparative analyses including genome-wide analyses of divergence with the chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genomes, and coalescent modeling to reconstruct the demographic history of the bobwhite for comparison to other birds currently in decline (i.e., scarlet macaw; Ara macao). More than 90% of the assembled bobwhite genome was captured within <40,000 final scaffolds (N50 = 45.4 Kb) despite evidence for approximately 3.22 heterozygous polymorphisms per Kb, and three annotation analyses produced evidence for >14,000 unique genes and proteins. Bobwhite analyses of divergence with the chicken and zebra finch genomes revealed many extremely conserved gene sequences, and evidence for lineage-specific divergence of noncoding regions. Coalescent models for reconstructing the demographic history of the bobwhite and the scarlet macaw provided evidence for population bottlenecks which were temporally coincident with human colonization of the New World, the late Pleistocene collapse of the megafauna, and the last glacial maximum. Demographic trends predicted for the bobwhite and the scarlet macaw also were concordant with how opposing natural selection strategies (i.e., skewness in the r-/K-selection continuum) would be expected to shape genome diversity and the effective population sizes in these species, which is directly relevant to future conservation efforts.