Analyses of Expressed Sequence Tags from Chironomus riparius Using Pyrosequencing : Molecular Ecotoxicology Perspective

Protein Fe-S Centers as a Molecular Target of Toxicity of a Complex Transition Metal Oxide Nanomaterial with Downstream Impacts on Metabolism and Growth

Oxidative stress is frequently identified as a mechanism of toxicity of nanomaterials. However, rarely have the specific underlying molecular targets responsible for these impacts been identified. We previously demonstrated significant negative impacts of transition metal oxide (TMO) lithium-ion battery cathode nanomaterial, lithium cobalt oxide (LCO), on the growth, development, hemoglobin, and heme synthesis gene expression in the larvae of a model sediment invertebrate Chironomus riparius. Here, we propose that alteration of the Fe-S protein function by LCO is a molecular initiating event leading to these changes. A 10 mg/L LCO exposure causes significant oxidation of the aconitase 4Fe-4S center after 7 d as determined from the electron paramagnetic resonance spectroscopy measurements of intact larvae and a significant reduction in the aconitase activity of larval protein after 48 h (p < 0.05). Next-generation RNA sequencing identified significant changes in the expression of genes involved in 4Fe-4S center binding, Fe-S center synthesis, iron ion binding, and metabolism for 10 mg/L LCO at 48 h (FDR-adjusted, p < 0.1). We propose an adverse outcome pathway, where the oxidation of metabolic and regulatory Fe-S centers of proteins by LCO disrupts metabolic homeostasis, which negatively impacts the growth and development, a mechanism that may apply for these conserved proteins across species and for other TMO nanomaterials.

A High-Quality Genome Assembly from Short and Long Reads for the Non-biting Midge Chironomus riparius (Diptera)

Chironomus riparius is of great importance as a study species in various fields like ecotoxicology, molecular genetics, developmental biology and ecology. However, only a fragmented draft genome exists to date, hindering the recent rush of population genomic studies in this species. Making use of 50 NGS datasets, we present a hybrid genome assembly from short and long sequence reads that make C. riparius’ genome one of the most contiguous Dipteran genomes published, the first complete mitochondrial genome of the species, and the respective recombination rate among the first insect recombination rates at all. The genome assembly and associated resources will be highly valuable to the broad community working with dipterans in general and chironomids in particular. The estimated recombination rate will help evolutionary biologists gaining a better understanding of commonalities and differences of genomic patterns in insects.

The differential expression of Chironomus spp genes as useful tools in the search for pollution biomarkers in freshwater ecosystems

Insects of the Chironomidae family are characterized by a wide ecological diversity in freshwater ecosystems. The larvae have the physiological potential to tolerate environmental stress even when there is a low concentration of oxygen, the presence of toxic substances or when there are changes in temperature and salinity. On the other hand, it is important to consider that at a cellular level, when individual insects are exposed to environmental changes, it induces responses of groups of genes that govern the molecular mechanisms related to such tolerance. In this review, using fourth instar larvae of Chironomus spp. in natural conditions and of Chironomus columbiensis under controlled conditions, we will discuss the genetic expression of a group of genes that respond to detoxification and also the biological functions involved and impacted on by mining stressors. The study of macroinvertebrate bioindicator species and their gene expression as a result of mining activity opens a window on the search for genetic biomarkers that could be used in environmental pollution assessments in freshwater ecosystems.

Molecular Characterization and Expression Analysis of P38 MAPK Gene and Protein in Aquatic Midge, Chironomus riparius (Diptera: Chironomidae), Exposed to Environmental Contaminants

P38 Mitogen-activated protein kinase (MAPK), an important signaling protein involved in various cellular processes, including stress responses, has been well characterized in model organisms. P38 has been identified in a number of insects, including the genus Drosophila; however, its homologue in Chironomus riparius has not yet been identified. In this study, we identified and characterized p38 MAPK (Crp38) gene in C. riparius using a transcriptome database that was previously generated 454 GS-FLX pyrosequencing. Comparative and phylogenetic analyses were performed using the p38 homologue of other species, such as Drosophila melanogaster, Aedes aegypti, Bombyx mori, Caenorhabditis elegans, Homo sapiens, etc. Furthermore, to test its potential as a biomarker of environmental contamination, Crp38 gene expression was analyzed upon exposure to nonylphenol (NP), silver nanoparticles (AgNPs), and cadmium (Cd). Crp38 gene expression was up- or down-regulated depending on the concentration and exposure duration of chemicals. These results show the role of Crp38 gene in defense against environmental stresses, as well as its potential use as a biomarker for various environmental pollutants. We further synthesized p38 antibody based on the predicted amino acid sequence deduced from Crp38 cDNA and, using this customized antibody, examined p38 protein expression in Cd exposed C. riparius. Although transcriptional alteration was not translated to the protein level, this result showed the possible application of a protein level functional study using cDNA sequence information from next-generation sequencing database in nonmodel organisms.

Characterization and expression of cytochrome p450 cDNA (CYP9AT2) in Chironomus riparius fourth instar larvae exposed to multiple xenobiotics

We identified and characterized a CYP9 family gene, CrCYP9AT2, from Chironomus riparius, an eco-toxicologically important model organism. The 1978 base pair (bp) length CrCYP9AT2 cDNA has an open reading frame of 1587 bp encoding a putative 528 amino acid protein. There was 267 bp 5' and 123 bp 3' untranslated region with a polyadenylation signal site (AATAAA). The putative heme-binding cysteine at position 471 and the typical p450 signature sequence of 463-FGIGPRNCIG-473 were also present. The CrCYP9AT2 transcript was present in all life stages with the highest expression in larvae. The modulation of CrCYP9AT2 was studied using real-time polymerase chain reaction after 24h exposure to cadmium chloride, benzo(a)pyrene; bisphenol A; nonylphenol; chlorpyrifos and ethinylestradiol. Significant up-regulation of CrCYP9AT2 gene was observed after exposure to Cd, B(a)P and CP. However, CrCYP9AT2 was significantly down-regulated after exposure to BPA, NP and EE.

Transcriptional regulation of glutathione biosynthesis genes, γ-glutamyl-cysteine ligase and glutathione synthetase in response to cadmium and nonylphenol in Chironomus riparius

We characterized Chironomus riparius glutathione (GSH) biosynthesis genes, γ-glutamyl-cysteine ligase catalytic subunit (cr-gcl) and glutathione synthetase (cr-gs) and studied their expression after cadmium (Cd) and nonylphenol (NP) exposure. The full length cDNA of the Cr-GCL catalytic subunit was 2185 base pair (bp) in length containing an open reading frame of 1905bp, a 13bp 5' and 267bp 3' untranslated regions. The theoretical molecular mass of the deduced amino acid sequence (633) was 72.65kDa with an estimated pI of 5.42. The partial cDNA of Cr-GS was 739bp in length consisting 221 amino acids. The deduced amino acid sequence of Cr-GCL and Cr-GS cDNAs showed high conservation with homologs from other species. In phylogenetic analysis Cr-GCL and Cr-GS were grouped with equivalent genes from insects belonging to the dipteran order. The expression of cr-gcl and cr-gs was measured using quantitative real-time PCR after exposure to sub lethal concentrations of Cd (2, 10 and 20mg/L) and NP (10, 50 and 100μg/L) for 12, 24, 48 and 72h using real-time PCR methods. The mRNA expression of Cr-GCL and Cr-GS was significantly modulated after exposure to different concentrations of Cd and NP for different time periods. Total GSH levels showed a non-significant decrease after exposure to Cd for 24h. However, no change in GSH levels was observed after exposure to NP for 24h. These results suggest that Cr-GS and Cr-GCL expression is modulated by Cd and NP stress and may play an important role in detoxification of xenobiotics and antioxidant defense. We conclude that Cr-GS and Cr-GCL could be used as biomarkers of Cd and NP stress in aquatic environment for the studied species.

Characterization and expression of superoxide dismutase genes in Chironomus riparius (Diptera, Chironomidae) larvae as a potential biomarker of ecotoxicity

Superoxide dismutase (SOD, EC 1.15.1.1) is an enzyme involved in the scavenging of reactive oxygen species (ROS) into molecular oxygen and hydrogen peroxide. In this study, a copper-zinc superoxide dismutase (Cu-ZnSOD) gene and a manganese superoxide dismutase (MnSOD) gene in aquatic midge, Chironomus riparius (CrSODs) was identified using an Expressed Sequence Tag (EST) database generated by 454 pyrosequencing. A multiple sequence alignment of C. riparius sequences revealed high homology with other insect sequences in terms of the amino acid level. Phylogenetic analysis of the CrSODs revealed that they were grouped with SODs of other organisms, such as Polypedilum vanderplanki, Drosophila melanogaster, Aedes aegypti, Anopheles gambiae, Culex quinquefasciatus and Bombyx mori. Expression of the corresponding CrSODs was analyzed during different developmental stages and following exposure to various environmental contaminants with different mode of actions i.e., paraquat, cadmium, benzo[a]pyrene, and chloropyrifos. CrSOD gene expression was significantly up or down regulated in response to exposure to the chemicals tested. The overall results suggested that SOD gene expression provided a platform for the understanding of oxidative stress responses caused by exposure to various environmental contaminants, and the SOD genes could be used as biomarkers for environmental disturbances such as oxidative stress initiated by xenobiotics.

Combining next-generation sequencing and microarray technology into a transcriptomics approach for the non-model organism Chironomus riparius

Whole-transcriptome gene-expression analyses are commonly performed in species that have a sequenced genome and for which microarrays are commercially available. To do such analyses in species with no or limited genome data, i.e. non-model organisms, necessary transcriptomics resources, i.e. an annotated transcriptome and a validated gene-expression microarray, must first be developed. The aim of the present study was to establish an advanced approach for developing transcriptomics resources for non-model organisms by combining next-generation sequencing (NGS) and microarray technology. We applied our approach to the non-biting midge Chironomus riparius, an ecologically relevant species that is widely used in sediment ecotoxicity testing. We sampled extensively covering all C. riparius developmental stages as well as toxicant exposed larvae and obtained from a normalized cDNA library 1.5 M NGS reads totalling 501 Mbp. Using the NGS data we developed transcriptomics resources in several steps. First, we designed 844 k probes directly on the NGS reads, as well as 76 k probes targeting expressed sequence tags of related species. These probes were tested for their affinity to C. riparius DNA and mRNA, by performing two biological experiments with a 1 M probe-selection microarray that contained the entire probe-library. Subsequently, the 1.5 M NGS reads were assembled into 23,709 isotigs and 135,082 singletons, which were associated to ∼55 k, respectively, ∼61 k gene ontology terms and which corresponded together to 22,593 unique protein accessions. An algorithm was developed that took the assembly and the probe affinities to DNA and mRNA into account, what resulted in 59 k highly-reliable probes that targeted uniquely 95% of the isotigs and 18% of the singletons. Concluding, our approach allowed the development of high-quality transcriptomics resources for C. riparius, and is applicable to any non-model organism. It is expected, that these resources will advance ecotoxicity testing with C. riparius as whole-transcriptome gene-expression analysis are now possible with this species.

Effects of cadmium chloride and nonylphenol on the expression of StAR-related lipid transfer domain containing protein (START1) gene in aquatic midge, Chironomus riparius

We identified and characterized a partial cDNA of StAR-related lipid transfer domain containing protein gene from Chironomus riparius (CrSTART1) having homology with human MLN64 and Drosophila melanogaster START1 (DmSTART1) and evaluated the effects of cadmium chloride (Cd) and nonylphenol (NP) on its expression. Pfam analysis identified the presence of two StAR-related lipid transfer (START) domains in CrSTART1 having several conserved amino acid residues, characteristic of the MLN64 and DmSTART1. The mRNA expression of CrSTART1 was observed in all developmental stages. The modulation in the mRNA expression of CrSTART1 was investigated after exposure to different concentrations Cd (0, 2, 10, and 20 mg/L) and NP (0, 10, 50, and 100 μg/L) for different time intervals in fourth instar larvae of C. riparius. Significant downregulation of CrSTART1 mRNA was observed after exposure to 2, 10 and 20 mg/L of Cd for 24, 48 and 72 h. Significant upregulation of CrSTART1 was observed after exposure to 10 and 50 μg/L of NP for 24, and 48 h period. At 100 μg/L of NP significant upregulation of CrSTART1 was observed after 12 h and downregulated after 24, 48 and 72 h.

Characterization and transcriptional regulation of thioredoxin reductase 1 on exposure to oxidative stress inducing environmental pollutants in Chironomus riparius

We characterized thioredoxin reductase 1 (TrxR1) from Chironomus riparius (CrTrxR1) and studied its expression under oxidative stress. The full-length cDNA is 1820bp long and contains an open reading frame (ORF) of 1488bp. The deduced CrTrxR1 protein has 495 amino acids and a calculated molecular mass of 54.41kDa and an isoelectric point of 6.15. There was a 71bp 5' and a 261bp 3' untranslated region with a polyadenylation signal site (AATAAA). Homologous alignments showed the presence of conserved catalytic domain Cys-Val-Asn-Val-Gly-Cys (CVNVGC), the C-terminal amino acids 'CCS' and conserved amino acids required in catalysis. The expression of CrTrxR1 is measured using quantitative real-time PCR after exposure to 50 and 100mg/L of paraquat (PQ) and 2, 10 and 20mg/L of cadmium chloride (Cd). CrTrxR1 mRNA was upregulated after PQ exposure at all conditions tested. The highest level of CrTrxR1 expression was observed after exposure to 10mg/L of Cd for 24h followed by 20mg/L for 48h. Significant downregulation of CrTrxR1 was observed after exposure to 10 and 20mg/L of Cd for 72h. This study shows that the CrTrxR1 could be potentially used as a biomarker of oxidative stress inducing environmental contaminants.