Zebrafish as a Model for Drug Screening in Genetic Kidney Diseases.

Genetic disorders account for a wide range of renal diseases emerging during childhood and adolescence. Due to the utilization of modern biochemical and biomedical techniques, the number of identified disease-associated genes is increasing rapidly. Modeling of congenital human disease in animals is key to our understanding of the biological mechanism underlying pathological processes and thus developing novel potential treatment options. The zebrafish (Danio rerio) has been established as a versatile small vertebrate organism that is widely used for studying human inherited diseases. Genetic accessibility in combination with elegant experimental methods in zebrafish permit modeling of human genetic diseases and dissecting the perturbation of underlying cellular networks and physiological processes. Beyond its utility for genetic analysis and pathophysiological and mechanistic studies, zebrafish embryos, and larvae are amenable for phenotypic screening approaches employing high-content and high-throughput experiments using automated microscopy. This includes large-scale chemical screening experiments using genetic models for searching for disease-modulating compounds. Phenotype-based approaches of drug discovery have been successfully performed in diverse zebrafish-based screening applications with various phenotypic readouts. As a result, these can lead to the identification of candidate substances that are further examined in preclinical and clinical trials. In this review, we discuss zebrafish models for inherited kidney disease as well as requirements and considerations for the technical realization of drug screening experiments in zebrafish.

Introduction

Modern genetic diagnostics allow the rapid discovery of human disease-associated mutations. Moreover, human genetic disorders can often be mimicked in animal models that can be exploited in large-scale chemical investigations for the search of modifiers of disease-associated phenotypes and potentially therapeutic compounds. The zebrafish (Danio rerio) has become an increasingly accepted vertebrate model organism for biomedical research (1, 2).

Despite being a member of the teleost class of fish species, there is great homology in development as well as cell- and organ-specific structural and physiological properties between zebrafish and humans. Furthermore, even with the evolutionary distance, > 80% of human disease-associated genes have orthologs in the zebrafish genome (3). The embryonic and larval characteristics of zebrafish include small size, ex utero development, optical transparency, and rapidity of organogenesis. In combination with the high fecundity of adult zebrafish and a relatively simple and cost-effective animal husbandry, this enables large-scale in vivo investigations. The zebrafish genome has been completely sequenced, thus facilitating genetic and genomic analysis and manipulation (3, 4). For instance, reverse genetics allow for precise investigation of associated phenotypes, by e.g., transient gene knockdown using antisense morpholino oligos or by genome-editing technologies like the CRISPR/Cas9 system (5, 6).

Due to the simplicity of the pronephros that can be readily studied in embryonic and larval stages, the zebrafish is an applicable experimental model system for the analysis of renal development and disease (7). The pronephros, as the earliest nephric stage, contains two nephrons sharing numerous genetic, structural, and functional aspects with the mammalian nephron (8). Phenotypic changes upon genetic alterations can be easily analyzed within intact live animals (9). Large-scale mutagenesis screens have identified various mutants affecting kidney development allowing the exploration of genetic and molecular mechanisms underlying pronephros development and function (10, 11). Moreover, reverse genetics approaches enable researchers to specifically alter orthologous genetic elements potentially associated with human disease. To date, major fields of research where such zebrafish models are being employed include glomerular (i.e., podocytopathies) and cystic renal disorders (i.e., ciliopathies).

Phenotype-based screening for drug discovery applications is increasingly employed in biomedical and pharmaceutical research. In contrast to target-based screening, phenotype-based approaches do not require exact knowledge of the therapeutic target (7). In addition, whole-organism in vivo approaches have the advantage that they can unravel toxic and other side-effects of drugs at a very early stage of the study. Over the last years, due to the versatility and power of the model, the zebrafish has emerged as the main vertebrate model system for high-throughput and high-content chemical screening experiments and large-scale phenotypic scoring (12, 13). Clear and scalable readouts for in vivo large-scale experiments can be readily established and a plethora of mutant and transgenic models expressing fluorescent proteins driven by tissue-specific promoters is available. In combination with automation technologies and dedicated sample handling workflows, this has led to various biomedical screening assays in fields such as genetics (14, 15), toxicology (16, 17), immunology and infection biology (18, 19), cardiovascular research (20, 21), drug discovery and safety (12, 13, 22, 23), personalized medicine (24), non-coding-genome analysis (25) as well as behavioral analysis (26, 27). Notably, several compounds that were identified in the zebrafish model have made it to preclinical and clinical trials, including new substance classes and repurposed drugs (12). For instance, in a chemical genetic screen testing 2.480 compounds, prostaglandin E2 (PGE2) was identified as an evolutionarily conserved regulator of hematopoietic stem cell (HSC) number in zebrafish embryos (28). Based on these results, a chemical derivative of PGE2 (Prohema), has been developed with the aim of improving the efficiency of HSC transplants using umbilical cord blood. Prohema has meanwhile advanced to Phase II clinical trials.

Genetic kidney diseases in zebrafish

Genetic kidney diseases can affect all parts of the kidney and its functions. To date, mutations in more than 150 genes have been identified that cause genetic human kidney diseases such as alterations of kidney development or specific glomerular and tubular diseases (29). Nephrogenesis in vertebrates is an intricate process that includes the successive formation of up to three kidneys depending on the species position in the phylogenetic tree: the pronephros, mesonephros, and metanephros (30). There is an increasing complexity with each successive kidney developing, however the structure and function of the basic renal units, the nephrons, remains largely unvaried across vertebrates (31).

In zebrafish, major vertebrate organ systems form within a few days after fertilization (32). The zebrafish pronephros is functional by 48 hpf (hours post fertilization) accomplishing the functions of blood filtration and osmoregulation (11, 33). It consists of two nephrons with a fused glomerulus at the midline (Figure 1A) (24). The tubular system consisting of the pronephric proximal and distal tubule and pronephric duct contains segment-specific conserved structural and physiological properties and spatio-temporal gene-expression patterns that are homologous to the human kidney (Figure 1B) (8, 35, 36). The zebrafish glomerulus is endowed with podocytes with extended and interdigitating foot processes and fenestrated endothelial cells forming a functional glomerular filtration barrier, analogous to the metanephric glomerulus of higher vertebrates (11).

Figure 1

The zebrafish pronephros: anatomical position and segmental organization. (A) Brightfield dorsal view of a 2 day post fertilization (dpf) zebrafish embryo (upper panel). The rectangle in the anterior trunk indicates the location of the proximal pronephric structures with a fused glomerulus at the midline that connects to the segmented pronephric tubules as labeled in the Tg(wt1b:egfp) zebrafish line by GFP expression (lower panel). (B) Schematic illustration of a zebrafish pronephros showing segmental organization of each nephron into glomerulus (G), neck (N), proximal convoluted tubule (PCT), proximal straight tubule (PST), distal early (DE), corpuscle of Stannius (CS), distal late (DL), and pronephric duct (PD) that fuse to the cloaca (C). Adapted from Wingert and Davidson (34).

At present, zebrafish is predominantly used as a genetic model for normal and abnormal kidney development, for hereditary glomerulopathies (i.e., podocytopathies) and for the study of ciliopathy-associated human cystic kidney diseases. These encompass polycystic kidney diseases and diseases of the nephronophthisis/medullary cystic kidney disease complex including more complex ciliopathies such as Joubert Syndrome, Meckel-Gruber Syndrome, and Bardet-Biedl-Syndrome (37). In this review, we focus on glomerulopathies and cystic kidney diseases.

Zebrafish as a model for human hereditary glomerulopathies

Zebrafish models can recapitulate human genetic glomerulopathies, i.e., a variety of podocytopathies that clinically often manifest by steroid-resistant nephrotic syndrome (SRNS) due to podocyte foot process effacement. SRNS is mostly therapy-resistant and leads to end-stage renal disease (ESRD) within a few years of onset. A growing number of SRNS-causing mutations have been identified. For example, mutations in NPHS1, encoding Nephrin, cause congenital nephrotic syndrome of the Finnish type. Morpholino knockdown of nphs1 in zebrafish results in edema and loss of slit-diaphragms with abnormal podocyte foot processes (38, 39). Mutations in WT1 (Wilms' tumor gene 1) have been associated with syndromic disorders such as Denys-Drash syndrome and Frasier syndrome, but also with diffuse mesangial sclerosis and early-onset isolated nephrotic syndrome (40). In zebrafish, knockdown of wt1a results in defects in podocyte development leading to glomerular injury and nephrosis (41). Mutations in NPHS2, Podocin, are the most relevant cause of autosomal-recessive SRNS of childhood. Zebrafish nphs2 morphants display pronephric glomerular hypoplasia with pericardial edema and ultrastructural glomerular damage of the filtration barrier (38, 39). Mutations in PLCE1 (Nephrocystin-3, NPHS3) have been identified in patients with SRNS and disease onset in the first year of life with a rapid progression to ESRD (42, 43). Zebrafish plce1 morphants display an impairment of the kidney filtration barrier as measured by tubular uptake of filtered 500 kDa fluorescent dextran, accompanied by edema, and severe disorganization of slit diaphragms (43). Other rare human mutations that were mimicked in zebrafish include ADCK4 (AarF domain containing kinase 4 gene) (44), KANK1 (kidney ankyrin repeat-containing protein 1), KANK2, KANK4 (45), CRB2 (Crumbs homolog 2) (46), NUP107 (Nuclear Pore Complex Subunit 107) (47), and ARHGDIA (48).

Whereas in many studies of genetic glomerulopathies the zebrafish has been used to model human disease-associated phenotypes, disruption of the glomerular filtration barrier can only be visualized by ultrastructural techniques like super-resolution or electron microscopy that are not compatible with large-scale chemical screens. Edema formation in zebrafish embryos can indirectly report glomerular barrier impairments; however, despite being easily observed in brightfield microscopy it is not exclusively linked to renal impairment (49). This restricts its value as a phenotypic readout parameter in chemical kidney screens. Functional assessment of glomerular filtration and barrier integrity can be achieved by monitoring the temporal reporter activity after microinjection of fluorescently labeled inulin (50) or dextrans of different molecular weight (51–53) into the vascular system; however, this method is laborious and incompatible with extensive screening experiments. Additionally, a transgenic zebrafish that expresses GFP (green fluorescent protein)-tagged vitamin D-binding protein (VDBP), which acts as a tracer for proteinuria, has been reported (54) and may serve as an attractive alternative for high-content and high-throughput screening (53).

Zebrafish as a model for human cystic kidney diseases

Cystic diseases of the kidney are frequent monogenic disorders in humans (55, 56), with primary cilia dysfunction being the unifying cellular mechanism leading to most if not all cystic kidney diseases (56–58). Mutations in a variety of genes encoding the primary cilia/centrosome complex cause ciliopathies often associated with the development of renal cysts, in both human and zebrafish (59–64). However, it must be noted that cilia in the zebrafish pronephros are motile, whereas human renal cilia are thought to be non-motile (65, 66), suggesting a potential contribution of lack of fluid dynamics to cyst formation in the zebrafish model.

Mutations in polycystin-1 and polycystin-2 are responsible for autosomal dominant polycystic kidney disease (ADPKD), the most common human congenital renal disorder (67). In zebrafish, polycystin-2 morpholino knockdown or mutation of orthologous pkd genes induces kidney cysts, hydrocephalus, left/right asymmetry defects, and strong dorsal axis curvature (63, 68, 69). Autosomal recessive polycystic kidney disease (ARPKD) usually manifests perinatally or in childhood. In addition to PKHD1 (polycystic kidney and hepatic disease 1), mutations in DZIP1L (DAZ interacting zinc finger protein 1 like) have recently been associated with ARPKD (70, 71). Other than dzip1l morphants and dzip1l CRISPR mutants, Lu et al. injected the dzip1l translation-blocking morpholino into Tg(wt1b:egfp) transgenic embryos expressing GFP under Wilms' tumor suppressor (wt1b) promoter in the pronephros (72), allowing for in vivo fluorescence imaging of the cystic pronephros.

Nephronophthisis (NPHP), an autosomal recessive cystic kidney disease, represents the most frequent genetic cause of ESRD in the first three decades of life (55). Nephronophthisis can be accompanied by anomalies in other organs, such as cerebellar vermis hypoplasia, laterality defects, intellectual disability, shortening of bones, retinal degeneration, and hepatobiliary disease (56). These features are represented in a variety of syndromes, including Senior–Løken syndrome, Joubert syndrome, Bardet–Biedl syndrome, and Jeune asphyxiating thoracic dystrophy (73, 74). To date, NPHP-causing mutations have been identified in more than 20 genes (56). Morpholino knockdown in zebrafish has been performed for nphp2 to 6 (75–83) and nphp11 (84). Zebrafish mutants for arl13b/sco (85), ahi1 (86), and cc2d2a (87) develop pronephric cysts to varying degrees and serve as models for Joubert syndrome-related disease.

Intraflagellar transport (IFT) proteins are essential for the development and maintenance of motile and sensory cilia and localize to the cilium, basal body, and/or centrosome (88). Several zebrafish ift mutants demonstrating renal cysts were identified by forward genetic screens (63, 89). IFT80 mutations underlie a subset of Jeune asphyxiating thoracic dystrophy cases, of which 20% are associated with kidney abnormalities including renal cysts (90). Ift80 morphants show a dose-dependent phenotype with strong body curvature, large kidney cyst, and pericardial edema. IFT172 mutations were initially reported in Jeune and Mainzer-Saldino syndromes, but have also been observed in patients with Bardet-Biedl syndrome (91, 92). Zebrafish mutants and morphants for ift172 resemble these phenotypes, with renal cysts readily detectable in brightfield images (91). In our work, we have shown that morpholino knockdown of ift80 and ift172 in Tg(wt1b:egfp) with fluorescently labeled kidney structures allow for visualization of pronephric cysts, providing a model system for large-scale chemical screening studies to identify chemical modifiers of cyst formation (93).

Chemical screening in zebrafish—technical aspects and considerations

The zebrafish can cost-effectively bridge the gap between high-throughput experimentation in cellular models lacking physiological context and low-throughput models such as rodents that are closer to human biology (12, 13, 23, 94). The optically transparent embryos and larvae fit into wells of microtiter plates rendering them amenable for automated microscopy applications using existing screening technologies (12, 13, 17, 23). Due to accessibility to genetic manipulation (37, 95), a plethora of zebrafish transgenic and mutant lines has been generated (www.zfin.org, www.ezrc.kit.edu, www.zebrafish.org) (96), complemented by transient labeling, knockdown and genome-editing techniques (97). This provides an extraordinary rich toolkit to model and visualize the biological processes underlying development and disease. Several transgenic lines labeling pronephric structures such as podocytes (98, 99), tubules, and ducts (98, 100–106) or both (72) have been established. In combination with genetic alteration, either by transient gene knockdown using antisense morpholino oligos or by genome-editing technologies, these lines can mimic e.g., ARPKD (71), ADPKD (107), and other cystic kidney diseases (93, 108) and enable in vivo fluorescence microscopy of the diseased pronephros. Although controversies exist regarding the use of morphants generated by morpholino knockdown (109–113), they remain a valuable tool in altering target gene expression, provided that all mandatory control experiments to validate the observed disease-associated phenotypes have been conducted (110, 111).

To our knowledge, large-scale screening experiments evaluating fluorescent pronephric phenotypes in models of genetic kidney diseases have not been demonstrated. In a chemical modifier screen using a custom library of 115 compounds in pkd2 and ift172 mutants displaying renal cysts (114), morphological parameters such as body axis curvature and/or laterality defects were scored. Histone deacetylase inhibitors Tricostatin A and valproic acid attenuated these phenotypes, and cyst size-reducing effects were confirmed in secondary assays. Additionally, a chemical screen of ~2,000 compounds identified histone deacetylase inhibitors to expand the pool of embryonic renal progenitor cells (115), a mechanism presumably involved in regeneration following acute kidney injury.

In combination with automated microscopy (116), zebrafish disease model systems allow performing large-scale phenotypic whole organism screening assays (117, 118). Phenotypic readouts encompass survival rates, overall morphology, physiological parameters, cell- and tissue-specific phenotypes, reporter gene expression patterns, and behavioral phenotypes (25, 27, 93, 119–123). Phenotypic screens within the context of a live vertebrate provide significant advantage over classical target-based in vitro assays as they do not require a priori knowledge about biochemical pathways affected, thus allowing unbiased identification of drug candidates or toxicological effects of substances. Furthermore, they intrinsically involve potential contributions of biodistribution, metabolism, and pharmacokinetics (1, 12, 13, 17, 95, 124).

Prior to large-scale zebrafish experiments careful planning and preparatory experiments are required to generate robust protocols and tailor conditions toward desired readouts in screening compatible disease models. In brief, handling of thousands of embryos causes logistical challenges, thus protocols for animal husbandry, micromanipulation, embryo culture, and treatment (e.g., anesthesia), and sample handling must be established. When image-based assays are carried out the usage of pigmentation mutant strains (e.g., casper line) or chemical treatment (e.g., 1-phenyl-2-thiourea, PTU) to block pigment formation is often necessary to adequately visualize internal structures. Moreover, to minimize false positive and negative results, non-specific developmental toxicity or off-target effects, control experiments must be carried out to titrate required compound concentration ranges as well as the treatment period during embryonic development. Importantly, controls should also be continuously carried out as a reference readout to benchmark observed phenotypic effects and normalize experimental variation. Finally, image acquisition routines must be balanced with analysis needs to ensure effective and robust scoring of phenotypic alterations.

Despite its amenability to large-scale experimentation, the full exploitation of the zebrafish model in screening assays is often hampered. While sample manipulation can be scaled efficiently, large-scale imaging, and phenotypic scoring remains challenging as available screening methodologies are usually optimized for in vitro assays (Figure 2A) (126). In comparison to cellular models, zebrafish embryos are large three-dimensional objects of complex morphology leading to random orientation of embryos within wells of microtiter plates (126, 127). This can obscure the view on target structures and leads to the generation of non-standardized image data. Therefore, novel sample preparations or automation strategies are needed, as it is unfeasible to upscale classical zebrafish mounting techniques. In our work, we developed orientation tools allowing the generation of agarose cavities within wells of microtiter plates for consistent positioning and orientation of zebrafish embryos (Figure 2B) (25, 93, 125). This enables the automated acquisition of consistent views of 48–96 hpf zebrafish larvae in large-scale screening scenarios (Figure 2C). For instance, we employed that methodology for imaging of embryonic kidneys in automated large-scale microscopy assays to score for morphological alterations of the pronephros upon compound exposure (93, 125, 128) (Westhoff et al. unpublished data) or capture phenotypic changes in cystic kidney disease models (129) (Pandey et al., unpublished data). Other more complex technical solutions employ microfluidic systems that combine automated detection and rotational orientation within glass capillaries followed by microscopic imaging (127, 130).

Figure 2

Overview of screening workflows for organ specific phenotypic screening in zebrafish. Shown are examples from our screening work that illustrate the automatic acquisition of higher resolution datasets of embryonic kidneys in zebrafish embryos. (A) Experimental manipulation of embryos prior to mounting and automated imaging such as microinjection or compound treatment. (B) Mounting of zebrafish embryos in agarose coated microtiter plates generated using 3D printed orientation tools. Agarose layers contain cavities allowing for consistent alignment and orientation of specimen. (C) Automated acquisition of standardized views (e.g., dorsal) of zebrafish embryos arrayed in microtiter plates. (D) Automated acquisition of multidimensional image datasets using smart imaging techniques. Pronephric areas of the Tg(wt1b:egfp) zebrafish transgenic line are detected in low resolution datasets using image processing tools and are subsequently imaged at higher resolution. The hair cross indicates the detected position and the bounding box the field of view in subsequent higher resolution imaging. Scale bars indicate 600 μm (left panel) or 150 μm (right panel). (E) Detailed visualization of kidney regions enabling scoring of kidney phenotypes. Shown are wildtype (first row) or cystic (other rows) kidneys of 72 hpf Tg(wt1b:egfp) embryos. (F) Automated quantitative analysis and phenotypic scoring using image processing techniques. Heatmap shows quantitative measurements of cystic areas as shown in (E). Figure panels are taken or modified from Westhoff et al. (93), Wittbrodt et al. (125), Pandey et al. (unpublished), and www.acquifer.de.

To date, the vast majority of zebrafish screens employ low magnification to capture the entire zebrafish embryo body, followed by subsequent analysis of phenotypic changes (13, 131). However, this significantly attenuates the power of using a fully developed live vertebrate embryo as in vivo visualization of morphological, physiological or genetic events on the cell- or tissue-specific level is hampered. The widespread usage of low resolution assays is largely due to the impracticality of positioning the regions of interest (ROI) within the small field of views of higher magnification objectives in combination with fixed scan-field coordinates of automated microscopes. Additionally, the spatiotemporal location of ROIs within the embryo body might be variable or unpredictable. To overcome that limitation, technologies have been developed that allow the automatic centering of ROIs in front of objective lenses of microscopes. Microfluidic systems (127, 130) can fulfill that requirement but usually require rather complex setups and are potentially challenged when overall morphological changes occur, or developmental and disease-associated processes are observed in time-lapse experiments. Several automated microscopy solutions for microplate-based screening have been reported that allow to teach or detect the position of ROIs followed by automatic centering and multidimensional image acquisition (126, 132–134). This can be an expert operator manually selecting target structures for subsequent automated imaging or more advanced methods employing automatic detection by image processing. These automated smart imaging approaches are based on interfacing the imaging device with external software tools that automatically detect coordinates of features of interest and send back machine commands containing instructions for re-centering, higher resolution imaging, or tracking of target structures. While several solutions have been reported, their application is usually restricted to cellular models and they are often characterized by a high level of complexity requiring expert knowledge in image processing, general programming and hard- and software interfacing. Therefore, to enable a widespread usage of such toolsets novel developments are needed that provide a simplified access to biomedical researchers. To this extent, we have developed a robotic microscopy platform (www.acquifer.de) with a smart imaging interface that allows to manually select ROIs, or to use any image processing software to send back human-readable script commands to the imaging device. We utilize this technique to e.g., acquire high-resolution datasets of cystic kidneys in a genetic zebrafish disease model to screen for modifiers of cystogenesis (Figures 2D–F) (129) (Pandey et al., unpublished data).

Due to the wide variety of potential zebrafish screening assays and thus phenotypic readouts there are myriad of potential quantitative descriptors that can be extracted from image-based datasets. This can include fluorescence intensity, morphological descriptors, or dynamic parameters and ranges from simple whole embryo signal intensity to spatiotemporal activity of fluorescent reporters or tissue dynamics and beyond. For instance, advanced segmentation techniques were used in the Tg(cdh17:egfp) zebrafish line to detect the fluorescently labeled tubular cells of the kidney (105). A full discussion of analysis strategies and development of automated image processing pipelines is beyond the scope of this manuscript (135). However, as post-acquisition analysis strategies are of vital importance for the success of any screening assay, the design of scoring pipelines needs careful consideration at early stages of experimental planning.

Conclusions

Phenotype-based, cost-effective whole-organism chemical screening in zebrafish offers a variety of advantages including the identification of disease-modifying drugs without knowledge of a validated target, the potential to identify compounds with polypharmacological efficacy, and the simultaneous assessment of compound efficacy, toxicity, biodistribution, and pharmacokinetics within a vertebrate model system. While a growing number of genes are being identified to cause human kidney diseases, therapeutic options to combat these diseases are often absent. Ideally, the use of genetically modified zebrafish mimicking human genetic disorders in conjunction with kidney-specific transgenic reporter lines or in conjunction with fluorescently-labeled functional reporter lines (or other secondary readouts), permit the implementation of chemical screening for disease-modifying substances in the field of genetic kidney diseases.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work, and approved it for publication.

Conflict of interest statement

GP and JG are employees of DITABIS AG, Pforzheim, Germany. The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.