The Ahr2-Dependent wfikkn1 Gene Influences Zebrafish Transcriptome, Proteome, and Behavior.

The aryl hydrocarbon receptor (AHR) is required for vertebrate development and is also activated by exogenous chemicals, including polycyclic aromatic hydrocarbons (PAHs) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). AHR activation is well-understood, but roles of downstream molecular signaling events are largely unknown. From previous transcriptomics in 48 h postfertilization (hpf) zebrafish exposed to several PAHs and TCDD, we found wfikkn1 was highly coexpressed with cyp1a (marker for AHR activation). Thus, we hypothesized wfikkn1's role in AHR signaling, and showed that wfikkn1 expression was Ahr2 (zebrafish ortholog of human AHR)-dependent in developing zebrafish exposed to TCDD. To functionally characterize wfikkn1, we made a CRISPR-Cas9 mutant line with a 16-bp deletion in wfikkn1's exon, and exposed wildtype and mutants to dimethyl sulfoxide or TCDD. 48-hpf mRNA sequencing revealed over 700 genes that were differentially expressed (p < .05, log2FC > 1) between each pair of treatment combinations, suggesting an important role for wfikkn1 in altering both the 48-hpf transcriptome and TCDD-induced expression changes. Mass spectrometry-based proteomics of 48-hpf wildtype and mutants revealed 325 significant differentially expressed proteins. Functional enrichment demonstrated wfikkn1 was involved in skeletal muscle development and played a role in neurological pathways after TCDD exposure. Mutant zebrafish appeared morphologically normal but had significant behavior deficiencies at all life stages, and absence of Wfikkn1 did not significantly alter TCDD-induced behavior effects at all life stages. In conclusion, wfikkn1 did not appear to be significantly involved in TCDD's overt toxicity but is likely a necessary functional member of the AHR signaling cascade.

Single-cell RNA-seq Analysis Reveals That Prenatal Arsenic Exposure Results in Long-term, Adverse Effects on Immune Gene Expression in Response to Influenza A Infection.

Arsenic exposure via drinking water is a serious environmental health concern. Epidemiological studies suggest a strong association between prenatal arsenic exposure and subsequent childhood respiratory infections, as well as morbidity from respiratory diseases in adulthood, long after systemic clearance of arsenic. We investigated the impact of exclusive prenatal arsenic exposure on the inflammatory immune response and respiratory health after an adult influenza A virus (IAV) lung infection. C57BL/6J mice were exposed to 100 ppb sodium arsenite in utero, and subsequently infected with IAV (H1N1) after maturation to adulthood. Assessment of lung tissue and bronchoalveolar lavage fluid at various time points post-IAV infection reveals greater lung damage and inflammation in arsenic-exposed mice versus control mice. Single-cell RNA sequencing analysis of immune cells harvested from IAV-infected lungs suggests that the enhanced inflammatory response is mediated by dysregulation of innate immune function of monocyte-derived macrophages, neutrophils, natural killer cells, and alveolar macrophages. Our results suggest that prenatal arsenic exposure results in lasting effects on the adult host innate immune response to IAV infection, long after exposure to arsenic, leading to greater immunopathology. This study provides the first direct evidence that exclusive prenatal exposure to arsenic in drinking water causes predisposition to a hyperinflammatory response to IAV infection in adult mice, which is associated with significant lung damage.

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress.

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

Opportunities and Challenges for Dietary Arsenic Intervention.

The diet is emerging as the dominant source of arsenic exposure for most of the U.S. population. Despite this, limited regulatory efforts have been aimed at mitigating exposure, and the role of diet in arsenic exposure and disease processes remains understudied. In this brief, we discuss the evidence linking dietary arsenic intake to human disease and discuss challenges associated with exposure characterization and efforts to quantify risks. In light of these challenges, and in recognition of the potential longer-term process of establishing regulation, we introduce a framework for shorter-term interventions that employs a field-to-plate food supply chain model to identify monitoring, intervention, and communication opportunities as part of a multisector, multiagency, science-informed, public health systems approach to mitigation of dietary arsenic exposure. Such an approach is dependent on coordination across commodity producers, the food industry, nongovernmental organizations, health professionals, researchers, and the regulatory community. https://doi.org/10.1289/EHP3997.

Arsenic alters transcriptional responses to Pseudomonas aeruginosa infection and decreases antimicrobial defense of human airway epithelial cells.

Arsenic contamination of drinking water and food threatens the health of hundreds of millions of people worldwide by increasing the risk of numerous diseases. Arsenic exposure has been associated with infectious lung disease in epidemiological studies, but it is not yet understood how ingestion of low levels of arsenic increases susceptibility to bacterial infection. Accordingly, the goal of this study was to examine the effect of arsenic on gene expression in primary human bronchial epithelial (HBE) cells and to determine if arsenic altered epithelial cell responses to Pseudomonas aeruginosa, an opportunistic pathogen. Bronchial epithelial cells line the airway surface, providing a physical barrier and serving critical roles in antimicrobial defense and signaling to professional immune cells. We used RNA-seq to define the transcriptional response of HBE cells to Pseudomonas aeruginosa, and investigated how arsenic affected HBE gene networks in the presence and absence of the bacterial challenge. Environmentally relevant levels of arsenic significantly changed the expression of genes involved in cellular redox homeostasis and host defense to bacterial infection, and decreased genes that code for secreted antimicrobial factors such as lysozyme. Using pathway analysis, we identified Sox4 and Nrf2-regulated gene networks that are predicted to mediate the arsenic-induced decrease in lysozyme secretion. In addition, we demonstrated that arsenic decreased lysozyme in the airway surface liquid, resulting in reduced lysis of Microccocus luteus. Thus, arsenic alters the expression of genes and proteins in innate host defense pathways, thereby decreasing the ability of the lung epithelium to fight bacterial infection.

Acidosis increases the susceptibility of respiratory epithelial cells to Pseudomonas aeruginosa-induced cytotoxicity.

Bacterial infection can lead to acidosis of the local microenvironment, which is believed to exacerbate disease pathogenesis; however, the mechanisms by which changes in pH alter disease progression are poorly understood. We test the hypothesis that acidosis enhances respiratory epithelial cell death in response to infection with Pseudomonas aeruginosa Our findings support the idea that acidosis in the context of P. aeruginosa infection results in increased epithelial cell cytotoxicity due to ExoU intoxication. Importantly, enforced maintenance of neutral pH during P. aeruginosa infection demonstrates that cytotoxicity is dependent on the acidosis. Investigation of the underlying mechanisms revealed that host cell cytotoxicity correlated with increased bacterial survival during an acidic infection that was due to reduced bactericidal activity of host-derived antimicrobial peptides. These findings extend previous reports that the activities of antimicrobial peptides are pH-dependent and provide novel insights into the consequences of acidosis on infection-derived pathology. Therefore, this report provides the first evidence that physiological levels of acidosis increase the susceptibility of epithelial cells to acute Pseudomonas infection and demonstrates the benefit of maintaining pH homeostasis during a bacterial infection.

In Vivo Characterization of an AHR-Dependent Long Noncoding RNA Required for Proper Sox9b Expression.

Xenobiotic activation of the aryl hydrocarbon receptor (AHR) by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) prevents the proper formation of craniofacial cartilage and the heart in developing zebrafish. Downstream molecular targets responsible for AHR-dependent adverse effects remain largely unknown; however, in zebrafish sox9b has been identified as one of the most-reduced transcripts in several target organs and is hypothesized to have a causal role in TCDD-induced toxicity. The reduction of sox9b expression in TCDD-exposed zebrafish embryos has been shown to contribute to heart and jaw malformation phenotypes. The mechanisms by which AHR2 (functional ortholog of mammalian AHR) activation leads to reduced sox9b expression levels and subsequent target organ toxicity are unknown. We have identified a novel long noncoding RNA (slincR) that is upregulated by strong AHR ligands and is located adjacent to the sox9b gene. We hypothesize that slincR is regulated by AHR2 and transcriptionally represses sox9b. The slincR transcript functions as an RNA macromolecule, and slincR expression is AHR2 dependent. Antisense knockdown of slincR results in an increase in sox9b expression during both normal development and AHR2 activation, which suggests relief in repression. During development, slincR was expressed in tissues with sox9 essential functions, including the jaw/snout region, otic vesicle, eye, and brain. Reducing the levels of slincR resulted in altered neurologic and/or locomotor behavioral responses. Our results place slincR as an intermediate between AHR2 activation and the reduction of sox9b mRNA in the AHR2 signaling pathway.

Recent Advances and Emerging Applications in Text and Data Mining for Biomedical Discovery.

Precision medicine will revolutionize the way we treat and prevent disease. A major barrier to the implementation of precision medicine that clinicians and translational scientists face is understanding the underlying mechanisms of disease. We are starting to address this challenge through automatic approaches for information extraction, representation and analysis. Recent advances in text and data mining have been applied to a broad spectrum of key biomedical questions in genomics, pharmacogenomics and other fields. We present an overview of the fundamental methods for text and data mining, as well as recent advances and emerging applications toward precision medicine.

Monomethylarsonous Acid (MMAIII) Has an Adverse Effect on the Innate Immune Response of Human Bronchial Epithelial Cells to Pseudomonas aeruginosa.

UNLABELLED: Arsenic is the number one contaminant of concern with regard to human health according to the World Health Organization. Epidemiological studies on Asian and South American populations have linked arsenic exposure with an increased incidence of lung disease, including pneumonia, and chronic obstructive pulmonary disease, both of which are associated with bacterial infection. However, little is known about the effects of low dose arsenic exposure, or the contributions of organic arsenic to the innate immune response to bacterial infection. This study examined the effects on Pseudomonas aeruginosa (P. aeruginosa) induced cytokine secretion by human bronchial epithelial cells (HBEC) by inorganic sodium arsenite (iAsIII) and two major metabolites, monomethylarsonous acid (MMAIII) and dimethylarsenic acid (DMAV), at concentrations relevant to the U.S. POPULATION: Neither iAsIII nor DMAV altered P. aeruginosa induced cytokine secretion. By contrast, MMAIII increased P. aeruginosa induced secretion of IL-8, IL-6 and CXCL2. A combination of iAsIII, MMAIII and DMAV (10 pbb total) reduced IL-8 and CXCL1 secretion. These data demonstrate for the first time that exposure to MMAIII alone, and a combination of iAsIII, MMAIII and DMAV at levels relevant to the U.S. may have negative effects on the innate immune response of human bronchial epithelial cells to P. aeruginosa.

Structurally distinct polycyclic aromatic hydrocarbons induce differential transcriptional responses in developing zebrafish.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the environment as components of fossil fuels and by-products of combustion. These multi-ring chemicals differentially activate the aryl hydrocarbon receptor (AHR) in a structurally dependent manner, and induce toxicity via both AHR-dependent and -independent mechanisms. PAH exposure is known to induce developmental malformations in zebrafish embryos, and recent studies have shown cardiac toxicity induced by compounds with low AHR affinity. Unraveling the potentially diverse molecular mechanisms of PAH toxicity is essential for understanding the hazard posed by complex PAH mixtures present in the environment. We analyzed transcriptional responses to PAH exposure in zebrafish embryos exposed to benz(a)anthracene (BAA), dibenzothiophene (DBT) and pyrene (PYR) at concentrations that induced developmental malformations by 120 h post-fertilization (hpf). Whole genome microarray analysis of mRNA expression at 24 and 48 hpf identified genes that were differentially regulated over time and in response to the three PAH structures. PAH body burdens were analyzed at both time points using GC-MS, and demonstrated differences in PAH uptake into the embryos. This was important for discerning dose-related differences from those that represented unique molecular mechanisms. While BAA misregulated the least number of transcripts, it caused strong induction of cyp1a and other genes known to be downstream of the AHR, which were not induced by the other two PAHs. Analysis of functional roles of misregulated genes and their predicted regulatory transcription factors also distinguished the BAA response from regulatory networks disrupted by DBT and PYR exposure. These results indicate that systems approaches can be used to classify the toxicity of PAHs based on the networks perturbed following exposure, and may provide a path for unraveling the toxicity of complex PAH mixtures.

Comparative developmental toxicity of environmentally relevant oxygenated PAHs.

Oxygenated polycyclic aromatic hydrocarbons (OPAHs) are byproducts of combustion and photo-oxidation of parent PAHs. OPAHs are widely present in the environment and pose an unknown hazard to human health. The developing zebrafish was used to evaluate a structurally diverse set of 38 OPAHs for malformation induction, gene expression changes and mitochondrial function. Zebrafish embryos were exposed from 6 to 120h post fertilization (hpf) to a dilution series of 38 different OPAHs and evaluated for 22 developmental endpoints. AHR activation was determined via CYP1A immunohistochemistry. Phenanthrenequinone (9,10-PHEQ), 1,9-benz-10-anthrone (BEZO), xanthone (XAN), benz(a)anthracene-7,12-dione (7,12-B[a]AQ), and 9,10-anthraquinone (9,10-ANTQ) were evaluated for transcriptional responses at 48hpf, prior to the onset of malformations. qRT-PCR was conducted for a number of oxidative stress genes, including the glutathione transferase(gst), glutathione peroxidase(gpx), and superoxide dismutase(sod) families. Bioenergetics was assayed to measure in vivo oxidative stress and mitochondrial function in 26hpf embryos exposed to OPAHs. Hierarchical clustering of the structure-activity outcomes indicated that the most toxic of the OPAHs contained adjacent diones on 6-carbon moieties or terminal, para-diones on multi-ring structures. 5-carbon moieties with adjacent diones were among the least toxic OPAHs while the toxicity of multi-ring structures with more centralized para-diones varied considerably. 9,10-PHEQ, BEZO, 7,12-B[a]AQ, and XAN exposures increased expression of several oxidative stress related genes and decreased oxygen consumption rate (OCR), a measurement of mitochondrial respiration. Comprehensive in vivo characterization of 38 structurally diverse OPAHs indicated differential AHR dependency and a prominent role for oxidative stress in the toxicity mechanisms.

AHR2 mutant reveals functional diversity of aryl hydrocarbon receptors in zebrafish.

The aryl hydrocarbon receptor (AHR) is well known for mediating the toxic effects of TCDD and has been a subject of intense research for over 30 years. Current investigations continue to uncover its endogenous and regulatory roles in a wide variety of cellular and molecular signaling processes. A zebrafish line with a mutation in ahr2 (ahr2(hu3335)), encoding the AHR paralogue responsible for mediating TCDD toxicity in zebrafish, was developed via Targeting Induced Local Lesions IN Genomes (TILLING) and predicted to express a non-functional AHR2 protein. We characterized AHR activity in the mutant line using TCDD and leflunomide as toxicological probes to investigate function, ligand binding and CYP1A induction patterns of paralogues AHR2, AHR1A and AHR1B. By evaluating TCDD-induced developmental toxicity, mRNA expression changes and CYP1A protein in the AHR2 mutant line, we determined that ahr2(hu3335) zebrafish are functionally null. In silico modeling predicted differential binding of TCDD and leflunomide to the AHR paralogues. AHR1A is considered a non-functional pseudogene as it does not bind TCCD or mediate in vivo TCDD toxicity. Homology modeling, however, predicted a ligand binding conformation of AHR1A with leflunomide. AHR1A-dependent CYP1A immunohistochemical expression in the liver provided in vivo confirmation of the in silico docking studies. The ahr2(hu3335) functional knockout line expands the experimental power of zebrafish to unravel the role of the AHR during development, as well as highlights potential activity of the other AHR paralogues in ligand-specific toxicological responses.

Silver nanospheres are cytotoxic and genotoxic to fish cells.

Nanoparticles are being widely investigated for a range of applications due to their unique physical properties. For example, silver nanoparticles are used in commercial products for their antibacterial and antifungal properties. Some of these products are likely to result in silver nanoparticles reaching the aquatic environment. As such, nanoparticles pose a health concern for humans and aquatic species. We used a medaka (Oryzias latipes) cell line to investigate the cytotoxicity and genotoxicity of 30nm diameter silver nanospheres. Treatments of 0.05, 0.3, 0.5, 3 and 5microg/cm(2) induced 80, 45.7, 24.3, 1 and 0.1% survival, respectively, in a colony forming assay. Silver nanoparticles also induced chromosomal aberrations and aneuploidy. Treatments of 0, 0.05, 0.1 and 0.3microg/cm(2) induced damage in 8, 10.8, 16 and 15.8% of metaphases and 10.8, 15.6, 24 and 24 total aberrations in 100 metaphases, respectively. These data show that silver nanoparticles are cytotoxic and genotoxic to fish cells.

Medaka (Oryzias latipes) as a sentinel species for aquatic animals: Medaka cells exhibit a similar genotoxic response as North Atlantic right whale cells.

Hexavalent chromium (Cr(VI)) is emerging as a major concern for aquatic environments, particularly marine environments. Medaka (Oryzias latipes) has been used as a model species for human and aquatic health, including the marine environment, though few studies have directly compared toxicological responses in medaka to humans or other aquatic species. We used a medaka fin cell line to compare the genotoxic response of medaka to Cr(VI) to the response observed in North Atlantic right whale cells to see if responses in medaka were similar to those of other aquatic species, particularly aquatic mammals. We used the production of chromosomal aberrations as a measure of genotoxicity. We found that in medaka cells, concentrations of 1, 5 and 10 microM sodium chromate damaged 17, 32 and 43% of metaphases, respectively and these same concentrations 1, 2.5, 5 and 10 microM sodium chromate damaged 14, 24 and 49% of metaphases, respectively, in North Atlantic right whale lung cells and 11, 32 and 41% of metaphases, respectively, in North Atlantic right whale testes cells. These data show that genotoxic responses in medaka are comparable to those seen in North Atlantic right whale cells, consistent with the hypothesis that medaka are a useful model for other aquatic species.

The cytotoxicity and genotoxicity of hexavalent chromium in medaka (Oryzias latipes) cells.

Chromium is an increasing health concern for aquatic environments, however, the mechanism of chromium toxicity in aquatic species is yet unknown. We used a medaka (Oryzias latipes) fin cell line to investigate the cytotoxicity and genotoxicity of sodium chromate, a soluble form of hexavalent chromium. We used a clonogenic cytotoxicity assay to measure sodium chromate cytotoxicity, gamma-H2A.X immunofluoresence to measure DNA double-strand breaks, and chromosome damage to measure clastogenicity. We found that sodium chromate is cytotoxic to medaka fin cells, with toxicity increasing in a concentration-dependent manner. Treatments of 0.5, 1, 5, 10, 25, 50 and 100 microM sodium chromate caused 100, 103.5, 87.8, 77.5, 40.9, 15 and 2.7% survival, respectively, relative to the control. We visualized DNA double-strand breaks in medaka cells through the formation of gamma-H2A.X foci. Breaks could be detected at concentrations as low as 1 microM. We also found that sodium chromate induces chromosomal aberrations, causing chromatid lesions and exchanges that increase with concentration. Treatments of 0, 1, 5, 10 and 25 microM sodium chromate damaged 10.3, 17, 32.3, 43 and 51.6% of metaphases and induced 13, 23, 44, 69 and 118 total aberrations in 100 metaphases, respectively. These data show that hexavalent chromium is both cytotoxic and genotoxic to fish cells. Our results set the context for future work in the medaka cell culture model and provide important tools for investigating mechanisms of toxicity in aquatic organisms.

Homologous recombination repair protects against particulate chromate-induced chromosome instability in Chinese hamster cells.

Particulate hexavalent chromium [Cr(VI)] compounds are well-established human carcinogens. Cr(VI)-induced tumors are characterized by chromosomal instability (CIN); however, the mechanisms of this effect are unknown. We investigated the hypothesis that homologous recombination (HR) repair of DNA double-strand breaks protect cells from Cr(VI)-induced CIN by focusing on the XRCC3 and RAD51C genes, which play an important role in cellular resistance to DNA double-strand breaks. We used Chinese hamster cells defective in each HR gene (irs3 for RAD51C and irs1SF for XRCC3) and compared with their wildtype parental and cDNA-complemented controls. We found that the intracellular Cr ion levels varied among the cell lines after particulate chromate treatment. Importantly, accounting for differences in Cr ion levels, we discovered that XRCC3 and RAD51C cells treated with lead chromate had increased cytotoxicity and chromosomal aberrations, relative to wildtype and cDNA-complimented cells. We also observed the emergence of high levels of chromatid exchanges in the two mutant cell lines. For example, 1microg/cm(2) lead chromate induced 20 and 32 exchanges in XRCC3- and RAD51C-deficient cells, respectively, whereas no exchanges were detected in the wildtype and cDNA-complemented cells. These observations suggest that HR protects cells from Cr(VI)-induced CIN, consistent with the ability of particulate Cr(VI) to induce double-strand breaks.

Modelling the impacts of two exotic invasive species on a native butterfly: top-down vs. bottom-up effects.

1. Exotic invasive species can influence population dynamics of native species through top-down or bottom-up forces. The present study examined separate and interactive effects of multiple exotic species invasions on the native mustard white butterfly, Pieris napi oleracea Harris (Lepidoptera: Pieridae), using a stochastic simulation model. 2. P. n. oleracea populations in North America have decreased regionally since the 1860s. Competition with an exotic congener (P. rapae L.), loss of native host plants and parasitism by the introduced broconid wasp (Cotesia glomerata L.), have been suggested to be independently responsible for its decline. The present study examined these hypotheses, as well as an alternative, invasion by an exotic crucifer, garlic mustard (Alliaria petiolata[Bieb.] Cavara & Grande). 3. A stochastic simulation model of P. n. oleracea population dynamics revealed that decreasing the number of host plants available for oviposition and larval development (i.e. habitat loss), sharply reduced the probability of populations persistence and decreased population size for those that persisted. 4. Simulated invasion by garlic mustard also substantially decreased both probability of persistence (= 0 at approximately 50% cover) and mean population size. Persistence probability never reached zero under any C. glomerata scenarios, even when larval mortality in the second generation due to parasitism was 100%. The impact of garlic mustard was intensified by the addition of C. glomerata parasitism. 5. Results suggest that bottom-up forces, loss of host plants through forest understorey loss and/or garlic mustard invasion are the most important forces driving P. n. oleracea population decline. Parasitism by C. glomerata may interact to reduce P. n. oleracea populations more rapidly, but appears insufficient alone to cause local extinction.