microPublication

Get Your Data Out, Be Cited

  • About
    • Editorial Policies
      • Editorial Staff
      • Editorial Board
      • Criteria For Publication
      • Publishing Information
      • Data Sharing Policy
    • For Authors
      • Preparation And Submission Of A Manuscript
      • Peer Review Process
      • Following Acceptance
      • Appeals
    • For Reviewers
    • Why micropublish?
  • Submit a microPublication
  • Journals
    • microPublication Biology
      • Editorial Board
  • microPublications
    • Biology
      • Species
        • Arabidopsis
        • C. elegans
        • D. discoideum
        • Drosophila
        • Human
        • Mouse
        • S. cerevisiae
        • S. pombe
        • Xenopus
        • Zebrafish
      • Categories
        • Phenotype Data
        • Methods
        • Expression Data
        • Genotype Data
        • Integrations
        • Genetic Screens
        • Models of Human Disease
        • Software
        • Interaction data
        • Database Updates
        • Electrophysiology Data
        • Phylogenetic Data
        • Science and Society
        • Biochemistry
  • Contact
  • More
    • Archives
    • FAQs
    • Newsletter
microPublication / Biology / The cadmium-responsive gene, cdr-6, does...
The cadmium-responsive gene, cdr-6, does not influence Caenorhabditis elegans responses to exogenous zinc
Kathryn S Evans1,2 and Erik C. Andersen1,3
1Molecular Biosciences, Northwestern University, Evanston, IL 60208
2Interdisciplinary Biological Sciences Program, Northwestern University, Evanston, IL 60208
3Robert H. Lurie Comprehensive Cancer Center, Northwestern University, Chicago, IL 60611
Correspondence to: Erik C. Andersen (erik.andersen@northwestern.edu)
Figure 1: Testing the role of cdr-6 in the C. elegans zinc response. A) Results from linkage mapping analysis using expression of cdr-6 as a quantitative trait. Genomic position (x-axis) is plotted against the logarithm of the odds (LOD) score (y-axis) for 13,003 genomic markers. Each significant QTL is indicated by a red triangle at the peak marker, and a blue rectangle shows the 95% confidence interval around the peak marker. The percentage of the total variance in the RIAIL population that can be explained by each QTL is shown above the QTL. B) For each QTL, the expression of cdr-6 (y-axis) in RIAILs split by genotype at the marker with the maximum LOD score (x-axis) are plotted as Tukey box plots. Each point corresponds to the relative expression of a unique recombinant strain. Strains with the N2 allele are colored orange, and strains with the CB4856 allele are colored blue. C) Gene model for cdr-6 is shown. Exons are represented by light purple rectangles and introns are represented by connecting lines. Location of the CRISPR-Cas9 deletion is shown with the grey box below the gene model. D) Strain genotypes are shown as colored rectangles (N2: orange, CB4856: blue) in detail for chromosome V (left) and in general for the rest of the chromosomes (right). The solid vertical line represents the peak marker of the zinc-response QTL, and the dashed vertical lines represent the confidence interval. Grey triangles represent cdr-6 deletions. E) Relative animal optical density in zinc (median.EXT, x-axis) is plotted as Tukey box plots against strain (see C, y-axis). Statistical significance of each strain compared to its parental strain (ECA1330 and ECA1332 to N2 and ECA1333 and ECA1334 to CB4856) is shown above each strain and colored by the parent strain it was tested against (ns = non-significant, p-value > 0.05; * = significant, p-value < 0.05).

Description

Regulation of the essential trace element zinc is necessary to avoid the toxic consequences caused by too little or too much of this metal (Vallee and Falchuk 1993; Rosen 2006). The zinc-response pathway has been extensively studied in the nematode roundworm Caenorhabditis elegans and several genes have been discovered that function to modulate sensitivity to both high and low zinc concentrations (Dietrich et al. 2016). Recently, we identified a quantitative trait locus (QTL) on the center of chromosome V, indicating that natural genetic variation between the laboratory strain, N2, and a genetically divergent wild isolate from Hawaii, CB4856, contributes to differential responses to excess zinc (Evans et al. 2020).

The confidence interval for this QTL is 1.6 Mb and contains 629 genes (WS263). Of these genes, 113 have one or more genetic variants predicted to modify the amino-acid sequence of the protein. However, protein-coding variation is just one of the ways that genetic variation can cause phenotypic variation. Another is variation in gene expression, which is hypothesized to be important in the majority of complex traits (Hindorff et al. 2009). We identified 83 expression QTL (eQTL) within this 1.6 Mb region using the eQTL dataset that mapped expression differences among a panel of recombinant inbred advanced intercross lines (RIAILs) also derived from N2 and CB4856 (Rockman et al. 2010; Evans and Andersen 2020). The most significant eQTL in this region caused a change in expression of the gene cdr-6 (Figure 1A). This gene, a homolog of the cadmium-response gene cdr-1, is highly expressed in intestinal and pharynx muscle cells during larval stages and is downregulated after treatments with arsenic, cadmium, or zinc (Dong et al. 2008). Furthermore, cdr-1 was previously shown to mitigate cadmium toxicity in C. elegans (Dong et al. 2005; Hall et al. 2012). Together, these data suggest that expression of cdr-6 might be toxic to C. elegans in the presence of heavy metals. Additionally, RIAILs with the CB4856 allele on chromosome V naturally express higher levels of cdr-6 (Figure 1B) and are also more sensitive to zinc than RIAILs with the N2 allele (Evans et al. 2020). This result indicates that strains with naturally low expression of cdr-6 are more resistant to excess zinc than strains with naturally high expression of this gene.

To test this hypothesis, we used CRISPR-Cas9 genome editing to create strains with large deletions of cdr-6 in both the N2 and CB4856 genetic backgrounds (Figure 1C,D). Because expression of cdr-6 was higher in RIAILs with the CB4856 allele (associated with zinc sensitivity) than in RIAILs with the N2 allele (associated with zinc resistance) (Figure 1B), we expected that a knockout of cdr-6 in the CB4856 genetic background would cause increased resistance to excess zinc. Alternatively, if variation in expression of cdr-6 underlies the zinc-response QTL on chromosome V, a knockout of cdr-6 in the N2 genetic background should not cause an increase in zinc resistance. We exposed N2, CB4856, and two strains with independently derived cdr-6 deletion alleles in each genetic background to elevated zinc and measured their optical densities using a high-throughput assay with the COPAS BIOSORT (Andersen et al. 2015; Evans and Andersen 2020; Evans et al. 2020). We found that strains with a deletion of cdr-6 phenocopied the strain with the same genetic background (Figure 1E), suggesting that differences in expression of cdr-6 do not underlie zinc responses.

This study not only provides evidence against cdr-6 as the causal gene underlying differences in zinc resistance between the N2 and CB4856 strains, but also indicates that cdr-6 does not influence zinc resistance in C. elegans. Although a previous study showed that cdr-6 was downregulated in response to zinc (Dong et al. 2008), our results do not contradict theirs. In fact, the authors also showed that the accumulation of fluid-filled droplets in the pseudocoelom, a phenotype observed in strains with inhibited CDR-6 expression using RNAi, is not increased in response to zinc or cadmium exposure (Dong et al. 2008). Taken together, we conclude that expression of cdr-6 decreases in response to zinc but animal development in the presence of zinc is not affected by cdr-6 function. It is likely that cdr-6 does not function in the nematode zinc response but rather is downregulated as an indirect effect of zinc exposure.

Methods

Request a detailed protocol

Strains

Animals were grown at 20°C on modified nematode growth media (NGMA) containing 1% agar and 0.7% agarose to prevent burrowing and fed OP50 (Andersen et al. 2014). The cdr-6 deletion strains are described below. All strains are available upon request.

Strain Genotype Deletion
ECA1330 cdr-6(ean181)

[N2 background]

730 bp deletion from V:12,413,681-12,412,951
ECA1332 cdr-6(ean183)

[N2 background]

811 bp deletion from V:12,413,723-12,412,912
ECA1333 cdr-6(ean184)

[CB4856 background]

632 bp deletion from V:12,413,034-12,413,665
ECA1334 cdr-6(ean185)

[CB4856 background]

788 bp deletion from V:12,412,864-12,413,651

Expression QTL mapping

Microarray expression for 14,107 probes were previously collected from 208 RIAILs (Rockman and Kruglyak 2009; Rockman et al. 2010), filtered (Andersen et al. 2014), and mapped using linkage mapping with 13,003 SNPs (Brady et al. 2019; Evans and Andersen 2020). The probe that measures expression for cdr-6 is A_12_P101644. All expression and eQTL data can be found in Evans and Anderson, 2020.

Generation of cdr-6 deletion strains

Deletion alleles for cdr-6 were generated using CRISPR-Cas9 genome editing as previously described (Hahnel et al. 2018; Brady et al. 2019). Briefly, 5’ and 3’ guide RNAs (Synthego, Redwood City, CA) were designed with the highest possible on-target and off-target scores (Doench et al. 2016). The CRISPR injection mix (1 µM dpy-10 sgRNA, 5 µM of each sgRNA for cdr-6, 0.5 µM of a single-stranded oligodeoxynucleotide template for homology-directed repair of dpy-10 (IDT, Skokie, IL), and 5 µM Cas9 protein (Q3B Berkeley, Berkeley, CA) in water) was assembled and incubated for an hour before injection. Young adults were mounted onto agar injection pads, injected in either the anterior or posterior arm of the gonad, and allowed to recover on 6 cm plates. Survivors were transferred to individual 6 cm plates 12 hours later and allowed to lay embryos. After two to three days, the F1 progeny were screened and individuals with Rol or Dpy phenotypes were selected and transferred to new 6 cm plates. After 48 hours, the F1 individuals were genotyped by PCR flanking the desired deletions. Individuals with deletions were propagated and genotyped for at least two additional generations to ensure homozygosity and to cross out the Rol mutation. Deletion amplicons were Sanger sequenced to identify the exact location of the deletion. Reagents used to generate and validate the deletion are found below:

crECA36 dpy-10 guide RNA:

GCUACCAUAGGCACCACGAG

crECA37 dpy-10 repair construct: CACTTGAACTTCAATACGGCAAGATGAGAATGACTGGAAACCGTACCGCATGCGGTGCCTA GGTAGCGGAGCTTCACATGGCTTCAGACCAACAGCCTAT

cdr-6 guide RNA:

crECA153: TCCAGTGACAACAACCCTAG

crECA115: AGTTTATAGTCGTTGTGTTG

External primers for deletion validation:

oECA1635: ACACACGTTCACTGGCTAGACT

oECA1636: AGCACGTCGTTGATATGCGAAC

Internal primers for deletion validation:

oECA1365: ACACACGTTCACTGGCTAGACT

oECA1637: TGTGTTCCCCATTGAGCTCGAT

High-throughput zinc-response assay

Zinc response was measured as a function of animal optical density (representing nematode development) as described previously (Brady et al. 2019; Evans and Andersen 2020). Briefly, strains were propagated for two generations, bleach-synchronized, and titered in K medium (Boyd et al. 2012) at a concentration of 25-50 embryos per well of a 96-well microtiter plate. The following day, arrested L1s were fed HB101 bacterial lysate (Pennsylvania State University Shared Fermentation Facility, State College, PA; (García-González et al. 2017)) at a final concentration of 5 mg/mL with either 1% water or 250 µM zinc in 1% water. Nematodes were allowed to grow at 20°C with constant shaking for 48 hours, treated with sodium azide (5 mM in water), and analyzed using the COPAS BIOSORT (Union Biometrica; Holliston, MA). Animal optical densities collected by the BIOSORT were processed and analyzed using the R package easysorter (Shimko and Andersen 2014) as described previously (Brady et al. 2019). Each strain has a minimum of 20 replicates and each replicate is defined by the median animal optical densities of six to 32 animals. Pairwise tests for statistically significant differences in the zinc response between strains were performed using the TukeyHSD function (R Core Team 2017) on an ANOVA model with the formula (phenotype ~ strain). For plotting purposes, these residual values were normalized from zero to one with zero being the well with the smallest value and one the well with the largest value.

Acknowledgments

We would like to thank members of the Andersen laboratory for their helpful comments on the manuscript. We would also like to thank WormBase and the C. elegans Natural Diversity Resource for data critical for our analysis.

References

Andersen, E. C., J. S. Bloom, J. P. Gerke, and L. Kruglyak, 2014 A variant in the neuropeptide receptor npr-1 is a major determinant of Caenorhabditis elegans growth and physiology. PLoS Genet. 10: e1004156.
PubMed
Andersen, E. C., T. C. Shimko, J. R. Crissman, R. Ghosh, J. S. Bloom et al., 2015 A powerful new quantitative genetics platform, combining Caenorhabditis elegans high-throughput fitness assays with a large collection of recombinant strains. G3 5: 911–920.
PubMed
Boyd, W. A., M. V. Smith, and J. H. Freedman, 2012 Caenorhabditis elegans as a model in developmental toxicology. Methods Mol. Biol. 889: 15–24.
PubMed
Brady, S. C., S. Zdraljevic, K. W. Bisaga, R. E. Tanny, D. E. Cook et al., 2019 A novel gene underlies bleomycin-response variation in Caenorhabditis elegans. Genetics 212: 1453–1468.
PubMed
Dietrich, N., C.-H. Tan, C. Cubillas, B. J. Earley, and K. Kornfeld, 2016 Insights into zinc and cadmium biology in the nematode Caenorhabditis elegans. Arch. Biochem. Biophys. 611: 120–133.
PubMed
Doench, J. G., N. Fusi, M. Sullender, M. Hegde, E. W. Vaimberg et al., 2016 Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9. Nat. Biotechnol. 34: 184–191.
PubMed
Dong, J., W. A. Boyd, and J. H. Freedman, 2008 Molecular characterization of two homologs of the Caenorhabditis elegans cadmium-responsive gene cdr-1: cdr-4 and cdr-6. J. Mol. Biol. 376: 621–633.
PubMed
Dong, J., M. O. Song, and J. H. Freedman, 2005 Identification and characterization of a family of Caenorhabditis elegans genes that is homologous to the cadmium-responsive gene cdr-1. Biochim. Biophys. Acta 1727: 16–26.
PubMed
Evans, K. S., and E. C. Andersen, 2020 The gene scb-1 underlies variation in Caenorhabditis elegans chemotherapeutic responses. G3 (Bethesda). 2020 Jul 7;10(7):2353-2364.
PubMed
Evans, K. S., S. Zdraljevic, L. Stevens, K. Collins, R. E. Tanny et al., 2020 Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans. bioRxiv.
2020.07.10.196857
García-González, A. P., A. D. Ritter, S. Shrestha, E. C. Andersen, L. S. Yilmaz et al., 2017 Bacterial metabolism affects the C. elegans response to cancer chemotherapeutics. Cell 169: 431–441.e8.
PubMed
Hahnel, S. R., S. Zdraljevic, B. C. Rodriguez, Y. Zhao, P. T. McGrath et al., 2018 Extreme allelic heterogeneity at a Caenorhabditis elegans beta-tubulin locus explains natural resistance to benzimidazoles. PLoS Pathog. 14: e1007226.
PubMed
Hall, J., K. L. Haas, and J. H. Freedman, 2012 Role of MTL-1, MTL-2, and CDR-1 in mediating cadmium sensitivity in Caenorhabditis elegans. Toxicol. Sci. 128: 418–426.
PubMed
Hindorff, L. A., P. Sethupathy, H. A. Junkins, E. M. Ramos, J. P. Mehta et al., 2009 Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc. Natl. Acad. Sci. U. S. A. 106: 9362–9367.
PubMed
R Core Team, 2017 R: A Language and Environment for Statistical Computing.
Rockman, M. V., and L. Kruglyak, 2009 Recombinational landscape and population genomics of Caenorhabditis elegans. PLoS Genet. 5: e1000419.
PubMed
Rockman, M. V., S. S. Skrovanek, and L. Kruglyak, 2010 Selection at linked sites shapes heritable phenotypic variation in C. elegans. Science 330: 372–376.
PubMed
Rosen, B. P., 2006 Metals in biology: past, present, and future, pp. 485–505 in Molecular Biology of Metal Homeostasis and Detoxification: From Microbes to Man, edited by M. J. Tamas and E. Martinoia. Springer Berlin Heidelberg, Berlin, Heidelberg.
Shimko, T. C., and E. C. Andersen, 2014 COPASutils: an R package for reading, processing, and visualizing data from COPAS large-particle flow cytometers. PLoS One 9: e111090.
PubMed
Vallee, B. L., and K. H. Falchuk, 1993 The biochemical basis of zinc physiology. Physiol. Rev. 73: 79–118.
PubMed

Funding

This work was supported by R01 ES029930 from the National Institute of Environmental Health Sciences to E.C.A. Additionally, K.S.E. received support from the NSF-Simons Center for Quantitative Biology at Northwestern University (awards Simons Foundation/SFARI 597491-RWC and the National Science Foundation 1764421) and the Cell and Molecular Basis of Disease Training grant (T32-GM008061).

Author Contributions

Kathryn S Evans: Conceptualization, Investigation, Methodology, Visualization, Writing - original draft, Writing - review and editing, Formal analysis
Erik C. Andersen: Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Writing - review and editing.

Reviewed By

Joel Meyer

History

Received: August 2, 2020
Revision received: September 4, 2020
Accepted: September 5, 2020
Published: September 14, 2020

Copyright

© 2020 by the authors. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Citation

Evans, KS; Andersen, EC (2020). The cadmium-responsive gene, cdr-6, does not influence Caenorhabditis elegans responses to exogenous zinc. microPublication Biology. 10.17912/micropub.biology.000305.
Download: RIS BibTeX
microPublication Biology is published by
1200 E. California Blvd. MC 1-43 Pasadena, CA 91125
The microPublication project is supported by
The National Institute of Health -- Grant #: 1U01LM012672-01
microPublication Biology:ISSN: 2578-9430