Search

The heritability and actionability of variants in homologous recombination-related genes in biliary tract cancers (BTCs) are uncertain. Although associations between BTC and BRCA germline variants have been reported, homologous recombination deficiency has not been investigated in BTCs.We sequenced germline variants in 27 cancer-predisposing genes in 1,292 BTC cases and 37,583 controls without a personal nor family history of cancer. We compared pathogenic germline variant frequencies between cases and controls and documented the demographic and clinical characteristics of carriers. In addition, whole-genome sequencing of 45 BTC tissues was performed to evaluate homologous recombination deficiency status.Targeted sequencing identified 5,018 germline variants, which were classified into 317 pathogenic, 3,611 variants of uncertain significance, and 1,090 benign variants. Seventy-one BTC cases (5.5%) had at least one pathogenic variant among 27 cancer-predisposing genes. Pathogenic germline variants enriched in BTCs were present in BRCA1, BRCA2, APC, and MSH6 (p 0.00185). PALB2 variants were marginally associated with BTC (p = 0.01). APC variants were predominantly found in ampulla of Vater carcinomas. Whole-genome sequencing demonstrated that three BTCs with pathogenic germline variants in BRCA2 and PALB2, accompanied by loss of heterozygosity, displayed homologous recombination deficiency. Conversely, pathogenic germline variants in without a second hit or other homologous recombination-realted genes such as ATM and BRIP1 showed homologous recombination-proficient phenotypes.In this study, we describe the heritability and actionability of variants in homologous recombination-related genes, which could be used to guide screening and therapeutic strategies for BTCs.We found that 5.5% of biliary tract cancers (BTCs) in a Japanese population possessed hereditary cancer-predisposing gene alterations, including in BRCA and genes associated with colorectal cancer. Two hits in homologous recombination-related genes were required to confer a homologous recombination-deficient phenotype. PARP inhibitors and DNA-damaging regimens may be effective strategies against BTCs exhibiting homologous recombination deficiency. Hence, in this study, genome-wide sequencing has revealed a potential new therapeutic strategy that could be applied to a subset of BTCs.

Von Hippel–Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80–90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5–22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study’s goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs’ driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We perform whole-genome sequencing and RNA sequencing analysis of 141 ESCC biopsy specimens before NAC treatment to generate a machine-learning-based diagnostic model to predict NAC reactivity in ESCC and analyzed the association between immunogenomic features and NAC response. Neutrophil infiltration may play an important role in ESCC response to NAC. We also demonstrate that specific copy-number alterations and copy-number signatures in the ESCC genome are significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a good indicator of therapeutic capability and a therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

Introduction: Treatment options for biliary tract cancer (BTC) are very limited. It is necessary to investigate actionable genes and candidate drugs using a sophisticated knowledgebase (KB) and characterize BTCs immunologically for evaluating the actionability of molecular and immune therapies. Materials and Methods: The genomic and transcriptome data of 219 patients with BTC who underwent surgery were analyzed. Actionable mutations and candidate drugs were annotated using the largest available KB of the Asian population (CancerSCAN®). Predictive biomarkers of immune checkpoint inhibitors were analyzed using DNA and RNA sequencing data. Results: Twenty-two actionable genes and 43 candidate drugs were annotated in 74 patients (33.8%). The most frequent actionable genes were PTEN (7.3%), CDKN2A (6.8%), KRAS (6.4%). BRCA2, CDKN2A, and FGFR2 mutations were most frequently identified in case of intrahepatic cholangiocarcinoma. PTEN and CDKN2A mutations were associated with significantly shorter overall survival. PD-L1 and PD-1 expression was significantly higher in case of extrahepatic cholangiocarcinoma and T-cell-high expression. In total, 49.7% of cases were evaluated as having actionability for molecular therapy or immune checkpoint inhibitors. Conclusions: Identifying actionable genes and candidate drugs using the KB contribute to the development of therapeutic drugs and personalized treatment for BTC.

Esophageal squamous cell carcinoma (ESCC) is the predominant type of esophageal cancer in the Asian region, including Japan. A previous study reported mutational landscape of Japanese ESCCs by using exome sequencing. However, somatic structural alterations were yet to be explored. To provide a comprehensive mutational landscape, we performed whole genome sequencing (WGS) analysis of biopsy specimens from 20 ESCC patients in a Japanese population. WGS analysis identified non-silent coding mutations of TP53, ZNF750 and FAT1 in ESCC. We detected six mutational signatures in ESCC, one of which showed significant association with smoking status. Recurrent structural variations, many of which were chromosomal deletions, affected genes such as LRP1B, TTC28, CSMD1, PDE4D, SDK1 and WWOX in 25%–30% of tumors. Somatic copy number amplifications at 11q13.3 (CCND1), 3q26.33 (TP63/SOX2), and 8p11.23 (FGFR1) and deletions at 9p21.3 (CDKN2A) were identified. Overall, these multi-dimensional view of genomic alterations improve the understanding of the ESCC development at molecular level and provides future prognosis and therapeutic implications for ESCC in Japan.

Simple Summary Gallbladder cancer (GBC) is a rare but lethal cancer. Molecular characterization of GBC is insufficient so far, and a comprehensive molecular portrait is warranted to uncover new targets and classify GBC. Clustering analysis of RNA expression revealed two subclasses of 36 GBCs, which reflects the status of the tumor microenvironment (TME) and poor prognosis of GBC, including epithelial-mesenchymal transition (EMT), immune suppression, and the TGF-beta signaling pathway. The knockout of miR125B1 in GBC cell lines decreased its invasion ability and altered the EMT pathway. Mutations of the genes related to the TGF-beta signaling pathway were enriched in the poor-prognosis/TME-rich cluster of GBCs. This comprehensive molecular analysis provides a new classification of GBCs based on the TME activity, which is involved with EMT and immune suppression for poor prognosis of GBC. This information may be useful for GBC prognosis and therapeutic decision-making. Gallbladder cancer (GBC), a rare but lethal disease, is often diagnosed at advanced stages. So far, molecular characterization of GBC is insufficient, and a comprehensive molecular portrait is warranted to uncover new targets and classify GBC. We performed a transcriptome analysis of both coding and non-coding RNAs from 36 GBC fresh-frozen samples. The results were integrated with those of comprehensive mutation profiling based on whole-genome or exome sequencing. The clustering analysis of RNA-seq data facilitated the classification of GBCs into two subclasses, characterized by high or low expression levels of TME (tumor microenvironment) genes. A correlation was observed between gene expression and pathological immunostaining. TME-rich tumors showed significantly poor prognosis and higher recurrence rate than TME-poor tumors. TME-rich tumors showed overexpression of genes involved in epithelial-to-mesenchymal transition (EMT) and inflammation or immune suppression, which was validated by immunostaining. One non-coding RNA, miR125B1, exhibited elevated expression in stroma-rich tumors, and miR125B1 knockout in GBC cell lines decreased its invasion ability and altered the EMT pathway. Mutation profiles revealed TP53 (47%) as the most commonly mutated gene, followed by ELF3 (13%) and ARID1A (11%). Mutations of ARID1A, ERBB3, and the genes related to the TGF-beta signaling pathway were enriched in TME-rich tumors. This comprehensive analysis demonstrated that TME, EMT, and TGF-beta pathway alterations are the main drivers of GBC and provides a new classification of GBCs that may be useful for therapeutic decision-making.

With the high cost and the long-term assessment of developmental toxicity testing in mammals, the vertebrate zebrafish has become a useful alternative model organism for high-throughput developmental toxicity testing. Zebrafish is also very favorable for the 3R perspective in toxicology; however, the methodologies used by research groups vary greatly, posing considerable challenges to integrative analysis. In this review, we discuss zebrafish developmental toxicity testing, focusing on the methods of chemical exposure, the assessment of morphological abnormalities, housing conditions and their effects on the production of healthy embryos, and future directions. Zebrafish as a systems toxicology model has the potential to elucidate developmental toxicity pathways, and to provide a sound basis for human health risk assessments.

Gastroenteropancreatic (GEP) neuroendocrine neoplasm (NEN) that consists of neuroendocrine tumor and neuroendocrine carcinoma (NEC) is a lethal but under-investigated disease owing to its rarity. To fill the scarcity of clinically relevant models of GEP-NEN, we here established 25 lines of NEN organoids and performed their comprehensive molecular characterization. GEP-NEN organoids recapitulated pathohistological and functional phenotypes of the original tumors. Whole-genome sequencing revealed frequent genetic alterations in TP53 and RB1 in GEP-NECs, and characteristic chromosome-wide loss of heterozygosity in GEP-NENs. Transcriptome analysis identified molecular subtypes that are distinguished by the expression of distinct transcription factors. GEP-NEN organoids gained independence from the stem cell niche irrespective of genetic mutations. Compound knockout of TP53 and RB1, together with overexpression of key transcription factors, conferred on the normal colonic epithelium phenotypes that are compatible with GEP-NEN biology. Altogether, our study not only provides genetic understanding of GEP-NEN, but also connects its genetics and biological phenotypes.

We performed genomic and transcriptomic sequencing of 133 combined hepatocellular and intrahepatic cholangiocarcinoma (cHCC-ICC) cases, including separate, combined, and mixed subtypes. Integrative comparison of cHCC-ICC with hepatocellular carcinoma and intrahepatic cholangiocarcinoma revealed that combined and mixed type cHCC-ICCs are distinct subtypes with different clinical and molecular features. Integrating laser microdissection, cancer cell fraction analysis, and single nucleus sequencing, we revealed both mono- and multiclonal origins in the separate type cHCC-ICCs, whereas combined and mixed type cHCC-ICCs were all monoclonal origin. Notably, cHCC-ICCs showed significantly higher expression of Nestin, suggesting Nestin may serve as a biomarker for diagnosing cHCC-ICC. Our results provide important biological and clinical insights into cHCC-ICC.

The developing brain is extremely sensitive to many chemicals. Exposure to neurotoxicants during development has been implicated in various neuropsychiatric and neurological disorders, including autism spectrum disorder, attention deficit hyperactive disorder, schizophrenia, Parkinson's disease, and Alzheimer's disease. Although rodents have been widely used for developmental neurotoxicity testing, experiments using large numbers of rodents are time-consuming, expensive, and raise ethical concerns. Using alternative non-mammalian animal models may relieve some of these pressures by allowing testing of large numbers of subjects while reducing expenses and minimizing the use of mammalian subjects. In this review, we discuss some of the advantages of using zebrafish in developmental neurotoxicity testing, focusing on central nervous system development, neurobehavior, toxicokinetics, and toxicodynamics in this species. We also describe some important examples of developmental neurotoxicity testing using zebrafish combined with gene expression profiling, neuroimaging, or neurobehavioral assessment. Zebrafish may be a systems toxicology model that has the potential to reveal the pathways of developmental neurotoxicity and to provide a sound basis for human risk assessments.

Age-related macular degeneration (AMD) is the most common cause of vision loss in elderly individuals throughout the developed world. Inhibitors of vascular endothelial growth factor have been successfully used to treat choroidal neovascularization in late-stage AMD. The pathogenesis of early-stage AMD, however, remains largely unknown, impairing efforts to develop effective therapies that prevent progression to late-stage AMD. To address this, we performed comparative transcriptomics of macular and extramacular retinal pigmented epithelium-choroid (RPE-choroid) tissue from early-stage AMD patients. We found that expression of fatty acid desaturase 1 (FADS1), FADS2, and acetyl-CoA acetyltransferase 2 (ACAT2) is increased in macular but not extramacular tissue, possibly through activation of sterol regulatory element binding factor 1 (SREBF1). Consistent with this, we also found that expression of Fads1 is increased in RPE-choroid in a mouse model of early-stage AMD. In zebrafish, deletion of fads2, which encodes a protein that functions as both Fads1 and Fads2 in other species, enhanced apoptosis in the retina upon exposure to intense light. Similarly, pharmacological inhibition of Srebf1 enhanced apoptosis and reduced fads2 expression in zebrafish exposed to intense light. These results suggest that the SREBF1–FADS1/2 axis may be activated in macular RPE-choroid as a protective response during early-stage AMD and could thus be a therapeutic target for early-stage AMD.

Oligodendrocytes are major myelin-producing cells and play essential roles in the function of a healthy nervous system. However, they are also one of the most vulnerable neural cell types in the central nervous system (CNS), and myelin abnormalities in the CNS are found in a wide variety of neurological disorders, including multiple sclerosis, adrenoleukodystrophy, and schizophrenia. There is an urgent need to identify small molecular weight compounds that can stimulate myelination. In this study, we performed comparative transcriptome analysis to identify pharmacodynamic effects common to miconazole and clobetasol, which have been shown to stimulate myelination by mouse oligodendrocyte progenitor cells (OPCs). Of the genes differentially expressed in both miconazole- and clobetasol-treated mouse OPCs compared with untreated cells, we identified differentially expressed genes (DEGs) common to both drug treatments. Gene ontology analysis revealed that these DEGs are significantly associated with the sterol biosynthetic pathway, and further bioinformatics analysis suggested that sterol regulatory element binding factors (SREBFs) might be key upstream regulators of the DEGs. In silico screening of a public database for chemicals associated with SREBF activation identified fenofibrate, a peroxisome proliferator-activated receptor α (PPARα) agonist, as a drug that increases the expression of known SREBF targets, raising the possibility that fenofibrate may also stimulate myelination. To test this, we performed in vivo imaging of zebrafish expressing a fluorescent reporter protein under the control of the myelin basic protein (mbp) promoter. Treatment of zebrafish with fenofibrate significantly increased expression of the fluorescent reporter compared with untreated zebrafish. This increase was attenuated by co-treatment with fatostatin, a specific inhibitor of SREBFs, confirming that the fenofibrate effect was mediated via SREBFs. Furthermore, incubation of zebrafish with another PPARα agonist, gemfibrozil, also increased expression of the mbp promoter-driven fluorescent reporter in an SREBF-dependent manner. These results suggest that activation of SREBFs by small molecular weight compounds may be a feasible therapeutic approach to stimulate myelination.

Pulmonary arterial hypertension (PAH) is a heterogeneous disorder associated with a progressive increase in pulmonary artery resistance and pressure. Although various therapies have been developed, the 5-year survival rate of PAH patients remains low. There is thus an important need to identify novel genes that are commonly dysregulated in PAH of various etiologies and could be used as biomarkers and/or therapeutic targets. In this study, we performed comparative transcriptome analysis of five mammalian PAH datasets downloaded from a public database. We identified 228 differentially expressed genes (DEGs) from a rat PAH model caused by inhibition of vascular endothelial growth factor receptor under hypoxic conditions, 379 DEGs from a mouse PAH model associated with systemic sclerosis, 850 DEGs from a mouse PAH model associated with schistosomiasis, 1598 DEGs from one cohort of human PAH patients, and 4260 DEGs from a second cohort of human PAH patients. Gene-by-gene comparison identified four genes that were differentially upregulated or downregulated in parallel in all five sets of DEGs. Expression of coiled-coil domain containing 80 (CCDC80) and anterior gradient 2 genes was significantly increased in the five datasets, whereas expression of SMAD family member 6 and granzyme A was significantly decreased. Weighted gene co-expression network analysis revealed a connection between CCDC80 and collagen type I alpha 1 (COL1A1) expression. To validate the function of CCDC80 in vivo, we knocked out ccdc80 in zebrafish using the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system. In vivo imaging of zebrafish expressing a fluorescent protein in endothelial cells showed that ccdc80 deletion significantly increased the diameter of the ventral artery, a vessel supplying blood to the gills. We also demonstrated that expression of col1a1 and endothelin-1 mRNA was significantly decreased in the ccdc80-knockout zebrafish. Finally, we examined Ccdc80 immunoreactivity in a rat PAH model and found increased expression in the hypertrophied media and adventitia of the pre-acinar pulmonary arteries and in the thickened intima, media, and adventitia of the obstructed intra-acinar pulmonary arteries. These results suggest that increased expression of CCDC80 may be involved in the pathogenesis of PAH, potentially by modulating the expression of endothelin-1 and COL1A1.

Mammals exhibit poor recovery after spinal cord injury (SCI), whereas non-mammalian vertebrates exhibit significant spontaneous recovery after SCI. The mechanisms underlying this difference have not been fully elucidated; therefore, the purpose of this study was to investigate these mechanisms. Using comparative transcriptome analysis, we demonstrated that genes related to cell cycle were significantly enriched in the genes specifically dysregulated in zebrafish SCI. Most of the cell cycle-related genes dysregulated in zebrafish SCI were down-regulated, possibly through activation of e2f4. Using a larval zebrafish model of SCI, we demonstrated that the recovery of locomotive function and neuronal regeneration after SCI were significantly inhibited in zebrafish treated with an E2F4 inhibitor. These results suggest that activation of e2f4 after SCI may be responsible, at least in part, for the significant recovery in zebrafish. This provides novel insight into the lack of recovery after SCI in mammals and informs potential therapeutic strategies.

Hypertrophic cardiomyopathy (HCM) is characterized by left ventricular hypertrophy and is associated with a number of potential outcomes, including impaired diastolic function, heart failure, and sudden cardiac death. Various etiologies have been described for HCM, including pressure overload and mutations in sarcomeric and non-sarcomeric genes. However, the molecular pathogenesis of HCM remains incompletely understood. In this study, we performed comparative transcriptome analysis to identify dysregulated genes common to five mouse HCM models of differing etiology: (i) mutation of myosin heavy chain 6, (ii) mutation of tropomyosin 1, (iii) expressing human phospholamban on a null background, (iv) knockout of frataxin, and (v) transverse aortic constriction. Gene-by-gene comparison identified five genes dysregulated in all five HCM models. Glutathione S-transferase kappa 1 (Gstk1) was significantly downregulated in the five models, whereas myosin heavy chain 7 (Myh7), connective tissue growth factor (Ctgf), periostin (Postn), and reticulon 4 (Rtn4) were significantly upregulated. Gene ontology comparison revealed that 51 cellular processes were significantly enriched in genes dysregulated in each transcriptome dataset. Among them, six processes (oxidative stress, aging, contraction, developmental process, cell differentiation, and cell proliferation) were related to four of the five genes dysregulated in all HCM models. GSTK1 was related to oxidative stress only, whereas the other four genes were related to all six cell processes except MYH7 for oxidative stress. Gene-gene functional interaction network analysis suggested correlative expression of GSTK1, MYH7, and actin alpha 2 (ACTA2). To investigate the implications of Gstk1 downregulation for cardiac function, we knocked out gstk1 in zebrafish using the clustered regularly interspaced short palindromic repeats/Cas9 system. We found that expression of the zebrafish homologs of MYH7, ACTA2, and actin alpha 1 were increased in the gstk1-knockout zebrafish. In vivo imaging of zebrafish expressing a fluorescent protein in cardiomyocytes showed that gstk1 deletion significantly decreased the end diastolic volume and, to a lesser extent, end systolic volume. These results suggest that downregulation of GSTK1 may be a common mechanism underlying HCM of various etiologies, possibly through increasing oxidative stress and the expression of sarcomere genes.

Intestinal helminths cause iron-deficiency anemia in pregnant women, associated with premature delivery, low birth weight, maternal ill health, and maternal death. Although benzimidazole compounds such as mebendazole (MBZ) are highly efficacious against helminths, there are limited data on its use during pregnancy. In this study, we performed in vivo imaging of the retinas of zebrafish larvae exposed to MBZ, and found that exposure to MBZ during 2 and 3 days post-fertilization caused malformation of the retinal layers. To identify the molecular mechanism underlying the developmental toxicity of MBZ, we performed transcriptome analysis of zebrafish eyes. The analysis revealed that the DNA damage response was involved in the developmental toxicity of MBZ. We were also able to demonstrate that inhibition of ATM significantly attenuated the apoptosis induced by MBZ in the zebrafish retina. These results suggest that MBZ causes developmental toxicity in the zebrafish retina at least partly by activating the DNA damage response, including ATM signaling, providing a potential adverse outcome pathway in the developmental toxicity of MBZ in mammals.

Mitochondrial dysfunction has been implicated in various drug-induced toxicities and genetic disorders. Recently, the zebrafish has emerged as a versatile animal model for both chemical and genetic screenings. Taking advantage of its transparency, various in vivo fluorescent imaging methods have been developed to identify novel functions of chemicals and genes in zebrafish. However, there have not been fluorescent probes that can detect mitochondrial membrane potential in living zebrafish. In this study, we identified a novel cyanine dye called ZMJ214 that detects mitochondrial membrane potential in living zebrafish from 4 to 8 days post fertilization and is administered by simple immersion. The fluorescence intensity of ZMJ214 in zebrafish was increased and decreased by oligomycin and FCCP, respectively, suggesting a positive correlation between ZMJ214 fluorescence and mitochondrial membrane potential. In vivo imaging of zebrafish stained with ZMJ214 allowed for the detection of altered mitochondrial membrane potential induced by the antidiabetic drug troglitazone and the antiepileptic drug tolcapone, both of which have been withdrawn from the market due to mitochondrial toxicity. In contrast, pioglitazone and entacapone, which are similar to troglitazone and tolcapone, respectively, and have been used commercially, did not cause a change in mitochondrial membrane potential in zebrafish stained with ZMJ214. Live imaging of zebrafish stained with ZMJ214 also revealed that knock-down of slc25a12, a mitochondrial carrier protein associated with autism, dysregulated the mitochondrial membrane potential. These results suggest that ZMJ214 can be a useful tool to identify chemicals and genes that cause mitochondrial dysfunction in vivo.

With the high cost and the long-term assessment of developmental toxicity testing in mammals, the vertebrate zebrafish has become a useful alternative model organism for high-throughput developmental toxicity testing. Zebrafish is also very favorable for the 3R perspective in toxicology; however, the methodologies used by research groups vary greatly, posing considerable challenges to integrative analysis. In this review, we discuss zebrafish developmental toxicity testing, focusing on the methods of chemical exposure, the assessment of morphological abnormalities, housing conditions and their effects on the production of healthy embryos, and future directions. Zebrafish as a systems toxicology model has the potential to elucidate developmental toxicity pathways, and to provide a sound basis for human health risk assessments.