Involvement of the zebrafish trrap gene in craniofacial development

Trrap (transformation/transcription domain-associated protein) is a component shared by several histone acetyltransferase (HAT) complexes and participates in transcriptional regulation and DNA repair; however, the developmental functions of Trrap in vertebrates are not fully understood. Recently, it has been reported that human patients with genetic mutations in the TRRAP gene show various symptoms, including facial dysmorphisms, microcephaly and global developmental delay. To investigate the physiological functions of Trrap, we established trrap gene-knockout zebrafish and examined loss-of-function phenotypes in the mutants. The trrap zebrafish mutants exhibited smaller eyes and heads than the wild-type zebrafish. The size of the ventral pharyngeal arches was reduced and the mineralization of teeth was impaired in the trrap mutants. Whole-mount in situ hybridization analysis revealed that dlx3 expression was narrowly restricted in the developing ventral pharyngeal arches, while dlx2b expression was diminished in the trrap mutants. These results suggest that trrap zebrafish mutants are useful model organisms for a human disorder associated with genetic mutations in the human TRRAP gene.

Trrap is a common component of various histone acetyltransferase (HAT) complexes and participates in transcription and DNA repair by recruiting HAT complexes to chromatin 1 ; however, the physiological functions of TRRAP are not fully understood. Trrap possesses FAT (FRAP, ATM, TRRAP), PIKK-TRRAP (pseudokinase domain of TRRAP) and FATC (FRAP, ATM, TRRAP C-terminal) domains. The kinase domain of PIKK-TRRAP lacks catalytic activity and the functions of the FAT and FATC domains are not fully understood. Trrap-knockout mice exhibit peri-implantation lethality due to blocked proliferation of blastocytes 2 , and conditional disruption of the Trrap gene achieved by crossing Trrap-floxed mice and Nestin-Cre mice causes premature differentiation of neural progenitors in the mouse brain 3 . Disruption of the mouse Trrap gene in embryonic stem cells (ESCs) causes unscheduled differentiation 4 , suggesting contributions of Trrap to self-renewal and appropriate differentiation. These results suggest that Trrap plays important roles in the generation of various organs during early vertebrate embryogenesis.
In a recent study, patients with missense variants in the human TRRAP gene predominantly presented with facial dysmorphisms (19/24 individuals), global developmental delay (24/24 individuals) and intellectual disability (17/20 individuals) 5 . Other phenotypes, such as microcephaly (7/24 individuals), hearing impairment (3/24 individuals) and visual impairment (4/24 individuals), were also observed in the patients, who presented various clinical spectra associated with TRRAP pathogenic missense variants. RNA sequencing analysis revealed that skin fibroblasts from patients exhibit significant expression differences in several genes, suggesting the involvement of TRRAP in the transcriptional regulation of various genes. Another group independently showed that the human TRRAP gene is responsible for autosomal dominant nonsyndromic hearing loss (ADNSHL) 6 . The p.Arg171Cys variant of the human TRRAP gene cosegregated with hearing loss in ADNSHL. The authors demonstrated that knockdown using an antisense morpholino targeting the trrap start codon or knockout via a 7-bp deletion near the trrap start codon (frameshift mutation) caused a decreased number of posterior lateral line (PLL) neuromasts and an impaired acoustic startle response (ASR) in zebrafish, but they did not describe other morphological abnormalities in their trrap mutants. Because patients with human TRRAP missense mutations exhibit multiple malformations, the physiological and developmental functions of Trrap remain unclear.
Zebrafish are useful animal models with which to investigate the physiological and pathological functions of genes responsible for human disorders 7 . Zebrafish mutants for a human disease can be established to analyze the processes of morphological abnormalities during early embryogenesis. Importantly, organogenesis, including craniofacial morphogenesis, is well conserved between zebrafish and humans 8 . In this study, we established a zebrafish trrap mutant and investigated loss-of-function phenotypes in the heads of the mutant fish. Our findings   5 . Therefore, we examined the morphological phenotypes in the heads of the trrap-mutant zebrafish. We observed that the diameter of the eye in the trrap mutants (trrap −/− : n = 8) at 3 days post-fertilization (dpf) was slightly smaller than that in the wild-type larvae containing the wild-type allele (trrap +/+ : n = 2; trrap +/− : n = 6) ( Fig. 2). Histological analysis using toluidine blue staining revealed that lamination of the retina to produce three nuclear layers (the retinal ganglion cell layer and the inner and outer nuclear layers) and two plexiform layers (the inner and outer plexiform layer) progressed normally in both the wild-type (trrap +/+ : n = 2; trrap +/− : n = 6) and the mutant zebrafish (trrap −/− : n = 8) (Fig. 2). The area of the head was measured excluding the eyes. The size of the head was smaller in the trrap mutants (trrap −/− : n = 8) than in the wild-type zebrafish (trrap +/+ : n = 1; trrap +/− : n = 7), while the distance between the eyes was comparable (Fig. 3).
Next, we examined the morphologies of the pharyngeal arches using Alcian blue to stain sulfated and carboxylated acid mucopolysaccharides 9 . We found that the length of ceratohyal cartilage was reduced in the mutants at 5 dpf ( Fig. 4 and Table 1). The angle of the paired ceratohyals in the mutants was larger than that in the wild-type fish, whereas the morphologies of the ethmoid plate and ceratobranchials appeared to be normal in the mutants ( Fig. 4 and Table 1).
Furthermore, we examined bone formation in the head using Alizarin red to stain calcium deposits in tissues 9 . We found that ossification hypoplasia was present in the head but not in cleithrum in the trrap mutants (trrap −/− : n = 14) compared to the wild-type fish (trrap +/+ : n = 5; trrap +/− : n = 11) (Fig. 4). Notably, wild-type fish possessed mineralized teeth on ceratobranchial 5, whereas the trrap mutants exhibited impairment of tooth mineralization. We observed similar defects in eyes, head, ventral pharyngeal arches and tooth development (Supplemental Fig. S3), when other trrap crRNAs (trrap-crRNA3 and trrap-crRNA4) were injected with tracrRNA and Cas9 into zebrafish embryos in the one-cell stage. These results suggest that the trrap gene is required for appropriate craniofacial development, including the development of the eyes, head, ventral pharyngeal arches and teeth.
Expression of pharyngeal arch and tooth marker genes in the trrap mutant. Because the trrap mutants exhibited craniofacial abnormalities, including pharyngeal arch and tooth defects, we investigated the expression of pharyngeal arch (dlx2a, dlx3 and nkx2.3) and tooth (dlx2b and pitx2) marker genes by WISH. The dlx2a is expressed in the pharyngeal arch ectomesenchyme 10 and the dlx3 is expressed in the ventral mesenchymal cells 11 , while the nkx2.3 is expressed in the lateral pharyngeal endoderm 12 . The expression patterns of pharyngeal arch genes in the trrap mutants were similar to those of wild-type zebrafish at 24 hpf (Supplemental Fig. S4). The expression of dlx2a, dlx3 and nkx2.3 at 72 hpf was narrowly restricted in the pharyngeal arches   www.nature.com/scientificreports/ of the mutants compared to those of the wild-type zebrafish (Fig. 4). The dlx2b and pitx2 genes at 72 hpf were weakly expressed in the developing teeth of the mutants (Fig. 5). The expression levels of neural genes, elavl3/ huC (neuron), gfap (radial glial and neural progenitors), olig2 (primary motor neuron and oligodendrocyte) and gli2a (central nervous system and pharyngeal arch), in the central nervous system (CNS) at 54 hpf were comparable between the wild-type and the trrap-mutant zebrafish (Supplemental Fig. S5). These results suggest that the expression of pharyngeal arch and tooth marker genes was restricted and diminished in the trrap mutant, respectively.

Discussion
Recent accumulating evidence demonstrates that human patients with genetic mutations in the TRRAP gene have various symptoms, including facial dysmorphisms, microcephaly, global developmental delay and intellectual disability 5 . Because the trrap-mutant zebrafish exhibited tooth hypoplasia, small eyes and head, and short ventral pharyngeal arches (Figs. 2, 3, 4), we proposed that the functional impairment of the human TRRAP and zebrafish trrap genes causes craniofacial abnormalities, including microcephaly. Recently, Xia W. et al. reported that the p.Arg171Cys variant near the T-terminus of the human TRRAP cosegregated with hearing loss 6 . They established trrap-knockout zebrafish that contained 7-bp deletion near the start codon, leading to a frameshift at amino acid 7 and a premature stop codon after eight amino acids. Their trrap-mutant zebrafish exhibited decreased numbers of posterior lateral line (PLL) neuromasts and an impaired acoustic startle response (ASR), however, they did not mention other morphological abnormalities in the mutants. Cogné et al. independently reported that most patients with missense TRRAP variants exhibited primarily facial dysmorphisms, global developmental delay and intellectual disability 5 . Hearing impairment (3/24 individuals) and visual impairment (4/24 individuals) were observed in a small proportion of the patients. www.nature.com/scientificreports/ They also found that the cluster of human TRRAP mutations containing 13 variants carried by 24 patients was located between codons 1031 and 1159. In this study, we established a zebrafish trrap mutant that presumably generated the 1002 N-terminal amino acids. We found that our trrap-mutant zebrafish exhibited several defects in craniofacial development, but the mutants had normal numbers of PLL neuromasts. It is not clear why the two different trrap mutant alleles show different phenotypes in zebrafish. Because their trrap mutant zebrafish has a premature stop codon near the start codon 6 , there is a possibility of the presence of trrap transcripts from the second methionine codon (18-3841 amino acids). In the case of tif1γ/moonshine (mon) mutant zebrafish, the hypomorphic phenotype in the mon m262 allele is interpreted by the presence of N-terminal truncated Tif1γ protein from another methionine downstream of a premature stop codon 13 . The patients with different missense mutations in human TRRAP gene exhibit various symptoms, including facial dysmorphisms and microcephaly, therefore, further analysis is required to clarify what kind molecules associate with the uncharacterized domain of TRRAP (1031-1159 amino acids). The expression of the zebrafish trrap gene was detected in the developing head, including in the eyes and pharyngeal arches, during early development (Fig. 1). We found that the eye diameter was slightly reduced in the trrap mutants at 3 dpf, whereas lamination of the retina to produce three nuclear layers and two plexiform layers occurred in the mutants (Fig. 2). The trrap mutants exhibited small heads compared to the wild-type fish, while the distance between the eyes was comparable (Fig. 3). The angle of the paired ceratohyal cartilage was larger than that of the wild-type ceratohyal cartilage (Fig. 4). We observed that the mineralization of teeth on ceratobranchial 5 was inhibited in the trrap mutants, while the wild-type fish had normally mineralized teeth. We confirmed that trrap-crispants injected with trrap-crRNA3, trrap-crRNA4, tracrRNA and Cas9 developed tooth hypoplasia, small eyes and head, and short ventral pharyngeal arches. Thus, these morphological abnormalities in the head may be associated with the craniofacial malformations of patients with TRRAP mutations 5 .
The contribution of Trrap to craniofacial development may be dependent on the activity of the HAT complex. Because the inhibition of histone deacetylation is involved in neural differentiation 2 , the ability to make the chromatin structure acceptable for transcriptional activation is an important aspect of targeted gene activation. One possible explanation is that the craniofacial defects in the trrap mutants were due to insufficient maintenance of pharyngeal arch gene expression and insufficient activation of tooth genes. We found that the expression patterns of pharyngeal arch markers (dlx2a, dlx3 and nkx2.3) at 24 hpf were similar to those in the wild-type, but these expression patterns were narrowly restricted at 72 hpf (Fig. 5, Supplemental Fig. S4). The expression of tooth markers (dlx2b and pitx2) was diminished in the trrap mutants (Fig. 6). On the other hand, the expression patterns of neural genes, such as elavl3/huC, gfap, olig2 and gli2a, were comparable in the CNS in the wild-type and mutant zebrafish. Possible targeted genes regulated by Trrap may be required for appropriate craniofacial development of structures including the head, eye, pharyngeal arches and teeth. We cannot exclude the possibility that the transcriptional difference might have been caused by dysfunction other than impaired HAT activity. Further analysis will be required to clarify how the Trrap protein contributes to craniofacial development in vertebrates. Because trrap-mutant zebrafish partly recapitulate craniofacial malformations, including

Zebrafish maintenance and ethics statement.
One-year-old adult heterozygous trrap-mutant zebrafish were maintained in a controlled aquatic facility with purified water by a reverse osmosis system with the following conditions: 14/10 h light/dark photoperiod, 28.5 °C (± 1 °C), pH 7.0 (± 1) and conductivity 450 mS/ cm. Zebrafish embryos were obtained from mating of adult heterozygous trrap-mutant zebrafish. The collected embryos were washed with E3 water, and fertilized embryos were selected with an optical microscope (Olympus SZ61). All animal experiments were performed in accordance with institutional and national guidelines and regulations. The study was carried out in compliance with the ARRIVE guidelines 14  Alcian blue staining and Alizarin red staining. Fixed larvae at 5 dpf were incubated overnight with 4% paraformaldehyde. After three washes with phosphate-buffered saline (PBS) containing 0.1% Tween 20 (PBS-T), the larvae were dehydrated with ethanol and incubated overnight with Alcian blue (0.02%) in 30% acetic acid and 70% ethanol 17 . The larvae were rehydrated with ethanol and incubated with 2% KOH and subsequently incubated with bleaching buffer (1% H 2 O 2 and 1% KOH) for 1 h. The larvae were washed three times with 2% KOH. Alizarin red staining was performed as described previously 17 . Larvae at 10 dpf were incubated with 4% paraformaldehyde in PBS and washed twice with PBS-T. The larvae were incubated with bleaching buffer for 30 min and washed with PBS-T twice. After 1 h of incubation with 1 mg/ml Alizarin red in 0.5% KOH at room temperature, the larvae were washed with 0.5% KOH.
Histological analysis. Fixed larvae were incubated in ethanol and embedded using a Technovit 8100 kit (Kulzer). The embedded larvae were sectioned at 5 μm on a Leica RM2125 microtome and mounted on slides. The larvae were stained with toluidine blue (0.1%) after sectioning 18 . Whole-mount in situ hybridization (WISH). The expression of trrap, dlx2a 19 , dlx2b, dlx3 11 , nkx2.3 12 , pitx2, elavl3 20 , gfap, olig2 21 and gli2a was examined by WISH as previously described 22 . Zebrafish embryos and larvae hybridized with the digoxygenin (DIG)-labeled antisense RNA probe were incubated with an alkaline phosphatase-conjugated anti-DIG antibody. The embryos and larvae were incubated with BM Purple (Roche) as the substrate to visualize the RNA probe recognized by the anti-DIG antibody. After three washes with PBST, the embryos and larvae were incubated with 4% paraformaldehyde.
Lateral line neuromast labeling. Larvae at 54 hpf were fixed in 4% paraformaldehyde for 3 h at room temperature and washed with PBS-T three times. The larvae were incubated in alkaline phosphatase buffer containing NBT and BCIP (Nacalai Tesque) for 30 min.