Dynamic transcriptomic profiles of zebrafish gills in response to zinc depletion

Abstract

<jats:title>Abstract</jats:title> <jats:sec> <jats:title>Background</jats:title> <jats:p>Zinc deficiency is detrimental to organisms, highlighting its role as an essential micronutrient contributing to numerous biological processes. To investigate the underlying molecular events invoked by zinc depletion we performed a temporal analysis of transcriptome changes observed within the zebrafish gill. This tissue represents a model system for studying ion absorption across polarised epithelial cells as it provides a major pathway for fish to acquire zinc directly from water whilst sharing a conserved zinc transporting system with mammals.</jats:p> </jats:sec> <jats:sec> <jats:title>Results</jats:title> <jats:p>Zebrafish were treated with either zinc-depleted (water = 2.61 μg L<jats:sup>-1</jats:sup>; diet = 26 mg kg<jats:sup>-1</jats:sup>) or zinc-adequate (water = 16.3 μg L<jats:sup>-1</jats:sup>; diet = 233 mg kg<jats:sup>-1</jats:sup>) conditions for two weeks. Gill samples were collected at five time points and transcriptome changes analysed in quintuplicate using a 16K oligonucleotide array. Of the genes represented the expression of a total of 333 transcripts showed differential regulation by zinc depletion (having a fold-change greater than 1.8 and an adjusted P-value less than 0.1, controlling for a 10% False Discovery Rate). Down-regulation was dominant at most time points and distinct sets of genes were regulated at different stages. Annotation enrichment analysis revealed that 'Developmental Process' was the most significantly overrepresented Biological Process GO term (<jats:italic>P</jats:italic> = 0.0006), involving 26% of all regulated genes. There was also significant bias for annotations relating to development, cell cycle, cell differentiation, gene regulation, butanoate metabolism, lysine degradation, protein tyrosin phosphatases, nucleobase, nucleoside and nucleotide metabolism, and cellular metabolic processes. Within these groupings genes associated with diabetes, bone/cartilage development, and ionocyte proliferation were especially notable. Network analysis of the temporal expression profile indicated that transcription factors <jats:italic>foxl1, wt1, nr5a1, nr6a1</jats:italic>, and especially, <jats:italic>hnf4a</jats:italic> may be key coordinators of the homeostatic response to zinc depletion.</jats:p> </jats:sec> <jats:sec> <jats:title>Conclusions</jats:title> <jats:p>The study revealed the complex regulatory pathways that allow the organism to subtly respond to the low-zinc condition. Many of the processes affected reflected a fundamental restructuring of the gill epithelium through reactivation of developmental programs leading to stem cell differentiation. The specific regulation of genes known to be involved in development of diabetes provides new molecular links between zinc deficiency and this disease. The present study demonstrates the importance of including the time-dimension in microarray studies.</jats:p> </jats:sec>