Transcriptome of embryonic and neonatal mouse cortex by high-throughput RNA sequencing

Brain structure and function experience dramatic changes from embryonic to postnatal development. Microarray analyses have detected differential gene expression at different stages and in disease models, but gene expression information during early brain development is limited. We have generated >27 million reads to identify mRNAs from the mouse cortex for>16,000 genes at either embryonic day 18 (E18) or postnatal day 7 (P7), a period of significant synapto-genesis for neural circuit formation. In addition, we devised strategies to detect alternative splice forms and uncovered more splice variants. We observed differential expression of 3,758 genes between the 2 stages, many with known functions or predicted to be important for neural development. Neurogenesis-related genes, such as those encoding Sox4, Sox11, and zinc-finger proteins, were more highly expressed at E18 than at P7. In contrast, the genes encoding synaptic proteins such as synaptotagmin, complexin 2, and syntaxin were up-regulated from E18 to P7. We also found that several neurological disorder-related genes were highly expressed at E18. Our transcriptome analysis may serve as a blueprint for gene expression pattern and provide functional clues of previously unknown genes and disease-related genes during early brain development.

Characterization of the porcine differentially expressed PDK4 gene and association with meat quality

To investigate the differential expression of genes in the skeletal muscle between Yorkshire and Chinese indigenous breed Meishan pigs, suppression subtractive hybridization was carried out and many genes were proved to be expressed significantly different in the two breeds. One gene highly expressed in Meishan but lowly expressed in Yorkshire specific library, shared strong homology with human pyruvate dehydrogenase kinase 4 (PDK4). Using semi-quantity and quantity PCR, We confirmed its differential expression between the two breeds. Temporal and spatial expression analysis indicated that porcine PDK4 gene is highly expressed in skeletal muscle and the highest in neonatal pigs. Complete cDNA cloning and sequence analysis revealed that porcine PDK4 gene contains an open reading frame of 1,221 bp. The deduced amino acid sequence showed conservation in evolution. A G/A mutation in intron 9 was identified and association analysis showed that it was significantly associated with intramuscular fat, muscle water content.

Identification and characterization of hypoxia-induced genes in Carassius auratus blastulae embryonic cells using suppression subtractive hybridization

Organisms living in water are inevitably exposed to periods of hypoxia. Environmental hypoxia has been an important stressor having manifold effects on aquatic life. Many fish species have evolved behavioral, physiological, biochemical and molecular adaptations that enable them to cope with hypoxia. However, the molecular mechanisms of hypoxia tolerance in fish, remain unknown. in this study, we used suppression subtractive hybridization to examine the differential gene expression in CAB cells (Carassius auratus blastulae embryonic cells) exposed to hypoxia for 24 h. We isolated 2100 clones and identified 211 differentially expressed genes (e-value <= 5e-3; Identity > 45%). Among the genes whose expression is modified in cells, a vast majority involved in metabolism, signal transduction, cell defense, angiogenesis, cell growth and proliferation. Twelve genes encoding for ERO1-L, p53, CPO, HO-1, MKP2, PFK-2, cystatin B, GLUT1, BTG1, TGF beta 1, PGAM1, hypothetical protein F1508, were selected and identified to be hypoxia-induced using semi-quantitive RT-PCR and real-time PCR. Among the identified genes, two open reading frames (ORFs) encoding for CaBTG1 and Cacystatin B were obtained. The deduced amino acid sequence of CaBTG1 had 94.1%, 72.8%, 72.8%, 72.8%, 68.6% identity with that of DrBTG1, HsBTG1, BtBTG1, MmBTG1 and XIBTG1. Comparison of Cacystatin B with known cystatin B, the molecules exhibited 49.5 to 76.0% identity overall. These results may provide significant information for further understanding of the adaptive mechanism by which C. auratus responds to hypoxia. (c) 2008 Elsevier Inc. All rights reserved.

Characterization of development-related genes for the cestode Bothriocephalus acheilognathi

Differential gene expression of mature and immature Bothriocephalus acheilognathi cestodes was analyzed using the suppression subtractive hybridization technique. Five mature-associated cDNAs were isolated and characterized. Virtual Northern blot and RT-PCR analyses confirmed that four of the five genes were up-regulated in mature parasites. The sequence analysis revealed that one gene encoded the structural protein chorion precursor, and that three encoded functional proteins homologous to yolk ferritin, sodium/hydrogen exchanger and muscin-like protein. Another gene appeared to be specific to B. acheilognathi, encoding a putative metal-bound protein. Although results obtained in the present study are preliminary, the information about the five genes may provide clues for further investigation on the decline in parasite numbers during the maturation of B. acheilognathi.

Survivin expressions in human hepatoma HepG2 cells exposed to ionizing radiation of different LET

Survivin is a member of the inhibitors of apoptosis (IAP) protein family that interferes with post-mitochondrial events including activation of caspases. To examine the regulation of survivin expression in response to irradiation with different linear energy transfer (LET), human hepatoma HepG2 cells were irradiated in vitro with X-rays and carbon ions. Cellular sensitivities to low- and high-LET radiation were determined by colony formation. Survivin expression at mRNA and protein level were measured with RT-PCR and Western blot analyses, respectively. Radiation-induced cell cycle arrest and apoptosis were investigated with flow cytometry. We found that low-LET X-rays induced dose-dependent increases in survivin expression. After exposure to high-LET carbon ions, survivin expression gradually increased from 0 to 4 Gy, and then declined at 6 Gy. More pronounced survivin expression, stronger G(2)/M phase arrest was observed after exposure to carbon ions in comparison with X-rays at doses from 0 to 4 Gy. These observations indicate that there is a differential survivin expression in response to different LET radiations and the cycle arrest mechanism may be associated with it. In addition, our data on induction of apoptosis are compatible with the assumption that survivin expression induced by low-LET X-rays radiation may play a critical role in inhibiting apoptosis. However, after irradiation with ions, an anti-apoptotic function of survivin is not evident, possibly because of the serious damage produced by densely ionizing radiation.

Identifying novel molecular mechanisms of ghrelin receptor signalling underlying neural control of food intake: interaction with stress and impulsivity

The gut-hormone, ghrelin, activates the centrally expressed growth hormone secretagogue 1a (GHS-R1a) receptor, or ghrelin receptor. The ghrelin receptor is a G-protein coupled receptor (GPCR) expressed in several brain regions, including the arcuate nucleus (Arc), lateral hypothalamus (LH), ventral tegmental area (VTA), nucleus accumbens (NAcc) and amygdala. Activation of the GHS-R1a mediates a multitude of biological activities, including release of growth hormone and food intake. The ghrelin signalling system also plays a key role in the hedonic aspects of food intake and activates the dopaminergic mesolimbic circuit involved in reward signalling. Recently, ghrelin has been shown to be involved in mediating a stress response and to mediate stress-induced food reward behaviour via its interaction with the HPA-axis at the level of the anterior pituitary. Here, we focus on the role of the GHS-R1a receptor in reward behaviour, including the motivation to eat, its anxiogenic effects, and its role in impulsive behaviour. We investigate the functional selectivity and pharmacology of GHS-R1a receptor ligands as well as crosstalk of the GHS-R1a receptor with the serotonin 2C (5-HT2C) receptor, which represent another major target in the regulation of eating behaviour, stress-sensitivity and impulse control disorders. We demonstrate, to our knowledge for the first time, the direct impact of GHS-R1a signalling on impulsive responding in a 2-choice serial reaction time task (2CSRTT) and show a role for the 5-HT2C receptor in modulating amphetamine-associated impulsive action. Finally, we investigate differential gene expression patterns in the mesocorticolimbic pathway, specifically in the NAcc and PFC, between innate low- and high-impulsive rats. Together, these findings are poised to have important implications in the development of novel treatment strategies to combat eating disorders, including obesity and binge eating disorders as well as impulse control disorders, including, substance abuse and addiction, attention deficit hyperactivity disorder (ADHD) and mood disorders.

CAMP factor homologues in Propionibacterium acnes: a new protein family differentially expressed by types I and II

Analysis of the draft genome sequence of the opportunistic pathogen Propionibacterium acnes type strain NCTC 737 (=ATCC 6919) revealed five genes with sequence identity to the co-haemolytic Christie-Atkins-Munch-Peterson (CAMP) factor of Streptococcus agalactiae. The predicted molecular masses for the expressed proteins ranged from 28 to 30 kDa. The genes were present in each of the three recently identified recA-based phylogenetic groupings of P. acnes (IA, IB and 11), as assessed by PCR amplification. Conserved differences in CAMP factor gene sequences between these three groups were also consistent with their previous phylogenetic designations. All type IA, IB and 11 isolates were positive for the co-haemolytic; reaction on sheep blood agar. Immunoblotting and silver staining of SIDS-PAGE gels, however, revealed differential protein expression of CAMP factors amongst the different groups. Type IB and 11 isolates produced an abundance of CAMP factor 1, detectable by specific antibody labelling and silver staining of SDS-PAGE gels. In contrast, abundant CAMP factor production was lacking in type A isolates, although larger amounts of CAMP factor 2 were detectable by immunoblotting compared with type 11 isolates. While the potential role of the abundant CAMP factor 1 in host colonization or virulence remains to be determined, it should be noted that the type strain of P. acnes used in much of the published literature is a type A isolate and is, therefore, lacking in this attribute.

Hypoxia selects for androgen independent LNCaP cells with a more malignant geno- and phenotype.

Hypoxia confers resistance to common cancer therapies, however, it has also has been shown to result in genetic alterations which may allow a survival advantage and increase the tumorigenic properties of cancer cells. Additionally, it may exert a selection pressure, allowing expansion of tumor cells with a more aggressive phenotype. To further assess the role of hypoxia in malignant progression in prostate cancer we exposed human androgen dependent prostate cancer cells (LNCaP) to cycles of chronic hypoxia and isolated a subline, LNCaP-H1. This article describes the partial characterization of this cell line. The LNCaP-H1 subline showed altered growth characteristics and exhibited androgen independent growth both in vitro and in vivo. Furthermore, these cells were resistant to mitochondrial-mediated apoptosis, probably since the endogenous levels of Bax was lower and Bcl-2 higher than in the parental LNCaP cells. Microarray analysis revealed that a complex array of pathways had differential gene expression between the 2 cell lines, with LNCaP-H1 cells exhibiting a genetic profile which suggests that they may be more likely metastasize to distant organs, especially bone. This was supported by an in vitro invasion assay, and an in vivo metastasis study. This study shows that hypoxia can select for androgen independent prostate cancer cells which have a survival advantage and are more likely to invade and metastasize.

Androgen deprivation results in time-dependent hypoxia in LNCaP prostate tumours; informed scheduling of the bioreductive drug AQ4N improves treatment response

Androgen withdrawal induces hypoxia in androgen-sensitive tissue; this is important as in the tumour microenvironment hypoxia is known to drive malignant progression. This study examined the time-dependent effect of androgen deprivation therapy (ADT) on tumour oxygenation and investigated the role of ADT-induced hypoxia on malignant progression in prostate tumours. LNCaP xenografted tumours were treated with anti-androgens and tumour oxygenation measured. Dorsal skin fold chambers (DSF) were used to image tumour vasculature in vivo. Quantitative PCR (QPCR) identified differential gene expression following treatment with bicalutamide. Bicalutamide and vehicle-only treated tumours were re-established in vitro and invasion and sensitivity to docetaxel were measured. Tumour growth delay was calculated following treatment with bicalutamide combined with the bioreductive drug AQ4N. Tumour oxygenation measurements showed a precipitate decrease following initiation of ADT. A clinically relevant dose of bicalutamide (2mg/kg/day) decreased tumour oxygenation by 45% within 24h, reaching a nadir of 0.09% oxygen (0.67±0.06 mmHg) by day 7; this persisted until day 14 when it increased up to day 28. Using DSF chambers, LNCaP tumours treated with bicalutamide showed loss of small vessels at days 7 and 14 with revascularization occurring by day 21. QPCR showed changes in gene expression consistent with the vascular changes and malignant progression. Cells from bicalutamide-treated tumours were more malignant than vehicle-treated controls. Combining bicalutamide with AQ4N (50mg/kg; single dose) caused greater tumour growth delay than bicalutamide alone. This study shows that bicalutamide-induced hypoxia selects for cells that show malignant progression; targeting hypoxic cells may provide greater clinical benefit.

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

Novel Hypoxia-Associated Markers of Chemoresistance in Ovarian Cancer

Background:
Ovarian cancer is the fifth leading cause of cancer in women and has poor
long-term survival, in part, due to chemoresistance. Tumour hypoxia is associated with
chemoresistance in ovarian cancer. However, relatively little is known about the genes
activated in ovarian cancer which cause chemoresistance due to hypoxia. This study
aimed to firstly identify genes whose expression is associated with hypoxia-induced
chemoresistance, and secondly select hypoxia-associated biomarkers and evaluate their
expression in ovarian tumours.
Design:
Cisplatin-sensitive (A2780) and cisplatin-resistant (A2780cis) ovarian cancer
cell lines were exposed to combinations of hypoxia and/or cisplatin as part of a matrix
designed to reflect clinically relevant scenarios. RNA was extracted and interrogated
on Affymetrix Human Gene arrays. Differential gene expression was analysed for cells
exposed to hypoxia and/or treated with cisplatin. Potential markers of chemoresistance
were selected for evaluation in a cohort of ovarian tumour samples by R
T-PCR.
Results:
A wide range of genes associated with chemoresistance were differentially
expressed in cells exposed to hypoxia and/or cisplatin. Selected genes [ANGPTL4,
HER3 and HIF-1
α
] were chosen for further validation in a cohort of ovarian tumour
samples. High expression of ANGPTL4 trended towards reduced progression-free and
overall survival. High expression of HER3 trended to increased progression-free but
reduced overall survival, while high expression of HIF-1
α
trended towards reduced
progression-free and increased overall survival.
Conclusions:
In conclusion, this study has further characterized the relationship between
hypoxia and chemoresistance in an ovarian cancer model. We have also identified many
potential biomarkers of hypoxia and platinum resistance and provided initial validation
of a subset of these markers in ovarian cancer tissues.

Androgen receptor driven transcription in molecular apocrine breast cancer is mediated by FoxA1

Breast cancer is a heterogeneous disease and several distinct subtypes exist based on differential gene expression patterns. Molecular apocrine tumours were recently identified as an additional subgroup, characterised as oestrogen receptor negative and androgen receptor positive (ER- AR+), but with an expression profile resembling ER+ luminal breast cancer. One possible explanation for the apparent incongruity is that ER gene expression programmes could be recapitulated by AR. Using a cell line model of ER- AR+ molecular apocrine tumours (termed MDA-MB-453 cells), we map global AR binding events and find a binding profile that is similar to ER binding in breast cancer cells. We find that AR binding is a near-perfect subset of FoxA1 binding regions, a level of concordance never previously seen with a nuclear receptor. AR functionality is dependent on FoxA1, since silencing of FoxA1 inhibits AR binding, expression of the majority of the molecular apocrine gene signature and growth cell growth. These findings show that AR binds and regulates ER cis-regulatory elements in molecular apocrine tumours, resulting in a transcriptional programme reminiscent of ER-mediated transcription in luminal breast cancers.

Natural Killer Cell-like signature observed in locally advanced rectal cancer after neoadjuvant chemoradiotherapy in a tumour regression grade dependent manner

Background: Around 10-15% of patients with locally advanced rectal cancer (LARC) undergo a pathologically complete response (TRG4) to neoadjuvant chemoradiotherapy; the rest of patients exhibit a spectrum of tumour regression (TRG1-3). Understanding therapy-related genomic alterations may help us to identify underlying biology or novel targets associated with response that could increase the efficacy of therapy in patients that do not benefit from the current standard of care.
Methods: 48 FFPE rectal cancer biopsies and matched resections were analysed using the WG-DASL HumanHT-12_v4 Beadchip array on the illumina iScan. Bioinformatic analysis was conducted in Partek genomics suite and R studio. Limma and glmnet packages were used to identify genes differentially expressed between tumour regression grades. Validation of microarray results will be carried out using IHC, RNAscope and RT-PCR.
Results: Immune response genes were observed from supervised analysis of the biopsies which may have predictive value. Differential gene expression from the resections as well as pre and post therapy analysis revealed induction of genes in a tumour regression dependent manner. Pathway mapping and Gene Ontology analysis of these genes suggested antigen processing and natural killer mediated cytotoxicity respectively. The natural killer-like gene signature was switched off in non-responders and on in the responders. IHC has confirmed the presence of Natural killer cells through CD56+ staining.
Conclusion: Identification of NK cell genes and CD56+ cells in patients responding to neoadjuvant chemoradiotherapy warrants further investigation into their association with tumour regression grade in LARC. NK cells are known to lyse malignant cells and determining whether their presence is a cause or consequence of response is crucial. Interrogation of the cytokines upregulated in our NK-like signature will help guide future in vitro models.

Nucella lapillus: imposex transcriptome analysis and phenotypic plasticity

O conhecimento de mecanismos de genómica funcional tem sido maioritariamente adquirido pela utilização de organismos modelo que são mantidos em condições laboratoriais. Contudo, estes organismos não reflectem as respostas a alterações ambientais. Por outro lado, várias espécies, ecologicamente bem estudadas, reflectem bem as interacções entre genes e ambiente mas que, das quais não existem recursos genéticos disponíveis. O imposex, caracterizado pela superimposição de caracteres sexuais masculinos em fêmeas, é induzido pelo tributilestanho (TBT) e trifenilestanho (TPT) e representa um dos melhores exemplos de disrupção endócrina com causas antropogénicas no ambiente aquático. Com o intuito de elucidar as bases moleculares deste fenómeno, procedeu-se à combinação das metodologias de pirosequenciação (sequenciação 454 da Roche) e microarrays (Agilent 4*180K) de forma a contribuir para um melhor conhecimento desta interacção gene-ambiente no gastrópode Nucella lapillus, uma espécie sentinela para imposex. O trancriptoma de N. lapillus foi sequenciado, reconstruído e anotado e posteriormente utilizado para a produção de um “array” de nucleótidos. Este array foi então utilizado para explorar níveis de expressão génica em resposta à contaminação por TBT. Os resultados obtidos confirmaram as hipóteses anteriormente propostas (esteróidica, neuroendócrina, retinóica) e adicionalmente revelou a existência de potenciais novos mecanismos envolvidos no fenómeno imposex. Evidência para alvos moleculares de disrupção endócrina não relacionados com funções reprodutoras, tais como, sistema imunitário, apoptose e supressores de tumores, foram identificados. Apesar disso, tendo em conta a forte componente reprodutiva do imposex, esta componente funcional foi a mais explorada. Assim, factores de transcrição e receptores nucleares lipofílicos, funções mitocondriais e actividade de transporte celular envolvidos na diferenciação de géneros estão na base de potenciais novos mecanismos associados ao imposex em N. lapillus. Em particular, foi identificado como estando sobre-expresso, um possível homólogo do receptor nuclear “peroxisome proliferator-activated receptor gamma” (PPARγ), cuja função na indução de imposex foi confirmada experimentalmente in vivo após injecção dos animais com Rosiglitazone, um conhecido ligando de PPARγ em vertebrados. De uma forma geral, os resultados obtidos mostram que o fenómeno imposex é um mecanismo complexo, que possivelmente envolve a cascata de sinalização envolvendo o receptor retinoid X (RXR):PPARγ “heterodimer” que, até à data não foi descrito em invertebrados. Adicionalmente, os resultados obtidos apontam para alguma conservação de mecanismos de acção envolvidos na disrupção endócrina em invertebrados e vertebrados. Finalmente, a informação molecular produzida e as ferramentas moleculares desenvolvidas contribuem de forma significativa para um melhor conhecimento do fenómeno imposex e constituem importantes recursos para a continuação da investigação deste fenómeno e, adicionalmente, poderão vir a ser aplicadas no estudo de outras respostas a alterações ambientais usando N. lapillus como organismo modelo. Neste sentido, N. lapillus foi também utilizada para explorar a adaptação na morfologia da concha em resposta a alterações naturais induzidas por acção das ondas e pelo risco de predação por caranguejos. O contributo da componente genética, plástica e da sua interacção para a expressão fenotípica é crucial para compreender a evolução de caracteres adaptativos a ambientes heterogéneos. A contribuição destes factores na morfologia da concha de N. lapillus foi explorada recorrendo a transplantes recíprocos e experiências laboratoriais em ambiente comum (com e sem influência de predação) e complementada com análises genéticas, utilizando juvenis provenientes de locais representativos de costas expostas e abrigadas da acção das ondas. As populações estudadas são diferentes geneticamente mas possuem o mesmo cariótipo. Adicionalmente, análises morfométricas revelaram plasticidade da morfologia da concha em ambas as direcções dos transplantes recíprocos e também a retenção parcial, em ambiente comum, da forma da concha nos indivíduos da F2, indicando uma correlação positiva (co-gradiente) entre heritabilidade e plasticidade. A presença de estímulos de predação por caranguejos estimulou a produção de conchas com labros mais grossos, de forma mais evidente em animais recolhidos de costas expostas e também provocou alterações na forma da concha em animais desta proveniência. Estes dados sugerem contra-gradiente em alterações provocadas por predação na morfologia da concha, na produção de labros mais grossos e em níveis de crescimento. O estudo das interacções gene-ambiente descritas acima demonstram a actual possibilidade de produzir recursos e conhecimento genómico numa espécie bem caracterizada ecologicamente mas com limitada informação genómica. Estes recursos permitem um maior conhecimento biológico desta espécie e abrirão novas oportunidades de investigação, que até aqui seriam impossíveis de abordar.

Toxicogenomics of natural and anthropogenic effectors in eleutheroembryos

The introduction of chemicals into the environment by human activities may represent a serious risk to environmental and human health. Environmental risk assessment requires the use of efficient and sensitive tools to determine the impact of contaminants on the ecosystems. The use of zebrafish for the toxicity assessment of pharmaceuticals, drugs, and pollutants, is becoming well accepted due to zebrafish unique advantages for the screening of compounds for hazard identification. The aim of the present work is to apply toxicogenomic approaches to identify novel biomarkers and uncovered potential modes of action of classic and emergent contaminants able to disrupt endocrine systems, such as the Retinoic Acid Receptor, Retinoid X Receptor and the Aryl Hydrocarbon Receptor. This study relies on different nuclear and cytosolic protein receptors and other conditional (ligand- or stress- activated) transcriptional factors that are intimately involved in the regulation of defensome genes and in mechanisms of chemical toxicity. The transcriptomic effects of organic compounds, endogenous compounds, and nanoparticles were analysed during the early stages of zebrafish development. Studying the gene expression profiles of exposed and unexposed organisms to pollutants using microarrays allowed the identification of specific gene markers and to establish a "genetic code" for the tested compounds. Changes in gene expression were observed at toxicant concentrations that did not cause morphological effects. Even at low toxicant concentrations, the observed changes in transcript levels were robust for some target genes. Microarray responses of selected genes were further complemented by the real time quantitative polymerase chain reaction (qRT-PCR) methodology. The combination of bio-informatic, toxicological analyses of differential gene expression profiles, and biochemical and phenotypic responses across the treatments allowed the identification of uncovered potential mechanisms of action. In addition, this work provides an integrated set of tools that can be used to aid management-decision making by improving the predictive capability to measure environmental stress of contaminants in freshwater ecosystems. This study also illustrates the potential of zebrafish embryos for the systematic, large-scale analysis of chemical effects on developing vertebrates.