Abstract
Vertebrate hearing organs manifest cellular asymmetries across the radial axis that underlie afferent versus efferent circuits between the inner ear and the brain. Therefore, understanding the molecular control of patterning across this axis has important functional implications. Radial axis patterning begins before the cells become postmitotic and is likely linked to the onset of asymmetric expression of secreted factors adjacent to the sensory primordium. This study explores one such asymmetrically expressed gene, Wnt9a, which becomes restricted to the neural edge of the avian auditory organ, the basilar papilla, by embryonic day 5 (E5). Radial patterning is disrupted when Wnt9a is overexpressed throughout the prosensory domain beginning on E3. Sexes were pooled for analysis and sex differences were not studied. Analysis of gene expression and afferent innervation on E6 suggests that ectopic Wnt9a expands the neural-side fate, possibly by re-specifying the abneural fate. RNA sequencing reveals quantitative changes, not only in Wnt-pathway genes, but also in genes involved in axon guidance and cytoskeletal remodeling. By E18, these early patterning effects are manifest as profound changes in cell fates [short hair cells (HCs) are missing], ribbon synapse numbers, outward ionic currents, and efferent innervation. These observations suggest that Wnt9a may be one of the molecules responsible for breaking symmetry across the radial axis of the avian auditory organ. Indirectly, Wnt9a can regulate the mature phenotype whereby afferent axons predominantly innervate neural-side tall HCs, resulting in more ribbon synapses per HC compared with abneural-side short HCs with few ribbons and large efferent synapses.
SIGNIFICANCE STATEMENT Wnts are a class of secreted factors that are best known for stimulating cell division in development and cancer. However, in certain contexts during development, Wnt-expressing cells can direct neighboring cells to take on specific fates. This study suggests that the Wnt9a ligand may play such a role in the developing hearing organ of the bird cochlea. This was shown through patterning defects that occur in response to the overexpression of Wnt9a. This manipulation increased one type of sensory hair cell (tall HCs) at the expense of another (short HCs) that is usually located furthest from the Wnt9a source. The extraneous tall HCs that replaced short HCs showed some physiological properties and neuronal connections consistent with a fate switch.
Introduction
Vertebrate hearing organs show tremendous variability over the course of evolution (Manley, 2017) while also sharing mechanosensory hair cells (HCs) as their primary receptors. Despite their variability, significant organizational homologies remain. The auditory organs of birds and mammals show similar specializations across the radial axis (Fuchs, 1992). The bird basilar papilla (BP) is populated by tall HCs (THCs) and short HCs (SHCs), although some gradation in size is apparent. THCs and the inner HCs (IHCs) of mammals are each positioned on the “neural” side of the organ, so called because of its proximity to the auditory ganglion (AG). These HCs make many afferent synapses with the peripheral processes of the auditory/spiral ganglion neurons, thereby alerting the brain to sound stimuli. SHCs and the outer HCs (OHCs) of mammals each reside on the “abneural” side and primarily receive efferent synaptic input from the brain (Tanaka and Smith, 1978; Fischer, 1992). In addition, small efferent terminals interact with THCs and some afferents contact SHCs.
The acquisition of different cell fates and connectivity across the organ likely requires the action of secreted ligands during embryogenesis. Asymmetric expression of Wnts and other signaling molecules (Fgfs, Bmps) influences radial patterning of the mammalian cochlea (Groves and Fekete, 2012). Likewise, in the BP on embryonic days 5 and 6 (E5–E6), before the first prosensory cells cease dividing, Wnt7b is expressed by abneural sensory cells, whereas Wnt7a and Wnt9a are present in a neural-side nonsensory domain (Sienknecht and Fekete, 2008). Wnt responsiveness and Wnt effects on cochlear proliferation have been explored in cultured chicken cochleas (Jacques et al., 2014). In addition to boosting proliferation, the regulation of cell patterning by Wnts (and Bmp) signaling was demonstrated using mouse cochlea organ cultures: the neural (IHC) compartment was enlarged and specification of the abneural (OHC) compartment was inhibited by Wnt activators (Munnamalai and Fekete, 2016). Both approaches involved pharmacological manipulations of the canonical Wnt-β-catenin signaling pathway. Potential problems with this approach are teasing apart off-target effects, activation of genes downstream of several potential Wnt ligands, bypassing the roles of Wnt receptors and secreted Wnt inhibitors, and a recognition that the involvement of β-catenin does not necessarily equate to a requirement for a Wnt ligand because β-catenin can interplay with other pathways (Itasaki and Hoppler, 2010; Collu et al., 2014). Morever, pharmacological approaches preclude being able to pinpoint the effects of an individual Wnt ligand. Finally, the organ cultures did not include the ganglion or brainstem and thus could not address changes in afferent and efferent innervation when confronted with alterations in cell types within the sensory epithelium. The current study was performed in vivo to address some of these questions.
One way to assess the bioactivity of a specific ligand is to misexpress a candidate molecule before specific cellular phenotypes manifest and look for phenotypic transformations. This approach was used to reveal a role for Bmp7 and retinoic acid in patterning HC size and sterociliarly morphology along the longitudinal (tonotopic) axis of the BP (Mann et al., 2014; Thiede et al., 2014). Here, we explore the molecular basis of patterning the organ's orthogonal axis, focusing on Wnt9a. We hypothesized that an asymmetric source of Wnt9a might endow neighboring sensory cells with neural-side identities. To test this idea, we overexpressed Wnt9a in the otocyst using retrovirus-mediated gene transfer. This created a wider BP with greatly expanded neural-side identity based on gene expression, neurite projections, and HC anatomy and physiology. THCs differentiated across the entire organ, whereas SHCs were missing. In summary, Wnt9a has the appropriate spatiotemporal expression in the normal avian cochlea to influence its radial patterning and, indeed, when overexpressed, it can change the local identity, cell fates, and innervation of the auditory sensory organ.
Materials and Methods
Virus injection.
Avian retrovirus encoding full-length chicken Wnt9a (formerly Wnt14) was generated by transfecting UMNSAH/DF-1 chicken embryo fibroblasts (ATCC #CRL-12203) with pRCAS(A)/Wnt9a plasmid DNA (Hartmann and Tabin, 2001). Growth and titration of virus on UMNSAH/DF-1 cells were performed as described previously (Morgan and Fekete, 1996). Specific pathogen-free G. gallus eggs (Charles River Laboratories) were windowed on E2, incubated at 38°C, and injected with concentrated RCAS(A)-derived viruses (1–2 × 109 infectious units/ml) into the E2 otic mesenchyme or E3 otic vesicle, as described previously (Stevens et al., 2003). Embryos were incubated in ovo until harvest. No sex determination was performed and all male and female embryos were pooled together. A total of 268 Wnt9a-overexpression (Wnt9a-OE) embryos and 18 RCAS(A) embryos from 28 batches survived the injection procedure to be analyzed histologically between E6 and E19. In all experiments, the left uninjected control BP was compared with the right Wnt9a-injected BP unless specified otherwise. An additional nine Wnt9a-OE embryos and five RCAS(A) embryos from six injection batches were subjected to electrophysiological analysis on E18–E20. For RNA sequencing, 6 injection batches were used to obtain 144 right BPs 3 d after injection of virus into the right E3 otocyst, with half Wnt9a-OE injected and half RCAS(A) injected.
Histology and transmission electron microscopy.
Embryos were fixed in 4% paraformaldehyde in PBS either by immersion or by transcardiac perfusion. For cryosectioning, BPs were processed through graded sucrose solutions and frozen in either 7.5% sucrose/15% gelatin/PBS or TFM tissue freezing medium (TBS). Other BPs were processed as whole mounts after fixation. For immunolabeling studies, primary antibodies from the Developmental Studies Hybridoma Bank (DSHB) were as follows: otoferlin/HCS-1 (antibody registry-AB_10804296; DSHB #HCS-1; 1:50; Goodyear et al., 2010), 3A10 (AB_531874; #3A10; 1:50; Sienknecht and Fekete, 2008), and AMV-3C2 (AB_528098; #AMV-3C2; 1:40; Potts et al., 1987). Ribbon synapses were immunolabeled with anti-CtBP2 (BD Transduction Laboratories #612044; 1:300; Khimich et al., 2005; Wahlin et al., 2008). Samples were processed for diaminobenzidine (DAB) colorimetric staining according to manufacturer's protocol (Vector Laboratories), or with fluorescent secondary antibodies (1:500; Invitrogen). For TEM studies, an E18 embryo was harvested and the inner ear epithelium exposed, after which the head was immersed in 2% paraformaldehyde/2% glutaraldehyde in 0.1 m cacodylate buffer (CB), pH 7.4, overnight at 4°C. Isolated BPs were washed 3 times in 0.1 m CB, postfixed with 1% OsO4/0.1 m CB for 1 h before washing in CB and water. Samples were dehydrated in an ascending ethanol series and once in propylene oxide (PO). Samples were infiltrated with sequential PO:Epon solutions (3:1, 1:1, 1:3 overnight) and then 100% Epon for 6 h. Semithin sections were stained with toluidine blue. Thin sections were collected and imaged on an FEI/Philips CM-10 transmission electron microscope using an accelerating voltage of 80 kV.
HC morphometric analysis.
Dissected BPs of a series of embryos ranging from E10 to posthatch day 35 were fixed in 4% paraformaldehyde/1% glutaraldehyde in PBS and postfixed in 2% OsO4 in PBS before Durcupan (EMS) embedding. A series of BP cross-sections (5 μm) were stained with toluidine blue (Sigma-Aldrich). HC height and widths were measured in radial cross-sections at 30%, 50%, 70%, and 90% locations from the proximal end using the Soft Imaging System (SIS; Olympus). A third-order polynomial line was fitted to visualize the trend of HC measurement distributions (see Fig. 1).
HC counts were performed on single 15 μm cryostat sections that were immunostained for HCS-1 and axons (3A10); HCS-1-positive cells from all focal planes were included in the counts when a clearly identified nucleus surrounded by HCS-1 immunoreactivity was observed in the HC layer.
Efferent labeling.
NeuroVue red-plus nylon filters coated with lipophilic dye (Molecular Targeting Technologies) were used to trace the BP efferents. After perfusion–fixation on E17–E19, the cerebellum was removed to reveal the brainstem. Filter strips of ∼0.5 × 4 mm size were inserted bilaterally below the vestibular nuclei with their long axis aligned with the anterior posterior dimension and their short axis angled at 35–45° from the midline (dorsal nearer to the midline than ventral). This targets the cell bodies of the octavolateral efferents (Rasmussen, 1946; Whitehead and Morest, 1981; Zidanic and Fuchs, 1995). Because contralaterally projecting efferent axons would be cut near the dorsal midline when the filters are inserted and pushed ventrally, this method is expected to label only ipsilateral efferents. The heads were incubated in 4% paraformaldehyde at 37°C for 3–4 weeks, after which the AGs and BPs were isolated and imaged.
Statistical analysis.
Tests were selected after confirming normal distributions of the data. Statistical significance was determined based on the test used and p-values are listed in the figure legends. A blinding protocol was not used in these analyses. Statistical tests and estimates of variance were performed using the following statistical software: Prism (GraphPad), JMP (SAS Institute), and SAS (SAS Institute). To quantify HCs and CtBP2 foci, high resolution z-stack images were acquired with the 60× objective on the neural and abneural edges at the 33–40% and 60–66% positions along the longitudinal axis of the BP with 0% at the proximal end. Each rectangular sample area was parallel to the long edge of the BP offset from the edge by the width of one sample area (∼60 μm) and encompassed 10,350 μm2, which is more than twice the length of the areas shown in Figure 3B. CtBP2 foci were counted manually through z-stacks visualized with the Nikon NIS elements image processing software. The number of ribbon foci per HC was quantified based on the superimposition of CtBP2-positive puncta with HCS-1 cytoplasmic labeling of individual HCs. Statistical significance was determined with a two-tailed t test with Bonferroni adjustment (α = 0.05/12 = 0.0042).
BP width measurements were quantified on ImageJ at each percentage position along the longitudinal axis and significance was determined with a one-tailed paired t test with Bonferroni adjustment (α = 0.05/8 = 0.006). Widths were measured by outlining the natural curvature of apical edge of the BPs (i.e., the sensory region only) from the neural to abneural side.
PH3+ cells were counted (within the sensory region only) from every third section through the BP and analyzed based on stage of development. A generalized linear mixed model with a least-squares method to fit the linear regression and a random embryo effect was used to evaluate side (left/right; represents treatment condition), stage of development, and their interactions. The analysis takes into account repeated measures on embryos (independence between embryos but dependence within an embryo). Post hoc Tukey's multiple-comparisons analysis was used to determine significance.
Whole-cell voltage-clamp electrophysiology.
HC physiology was assessed by recordings from semi-intact BPs from embryos between E18 and E20, as described previously (Li et al., 2009). Temporal bones were extracted into extracellular fluid (ECF) consisting of the following (in mm): 154 NaCl, 6 KCl, 5 CaCl2, 2 MgCl2, and 5 HEPES, pH 7.4, and the BPs were placed into ECF with 0.01% protease (type XXIV; Sigma-Aldrich) for 1–5 min. Segments of the sensory epithelium, ∼1 mm in length, were excised from a region between 50% and 75% distance from the basal end of the BP and mounted in a recording chamber containing ECF such that recordings could be made from HCs lining the neural or abneural edge. Current–voltage relationships were generated for four experimental groups: neural and abneural HCs from Wnt9a-OE and left BPs. Recordings from Wnt9a-OE BPs were limited to only those exhibiting morphological evidence of Wnt9a activity (i.e., increased width or presence of ectopic otoconia), although sampling regions were selected from longitudinal positions devoid of otoconia. Data were acquired with a Multiclamp 700B and Digidata 1440A under pCLAMP software control (Molecular Devices). Glass electrodes were filled with an internal solution containing the following (in mm): 112 KCl, 0.1 CaCl2, 2 MgCl2, 11 EGTA, 10 HEPES, 5 Na2ATP, 15 4-aminopyridine, and 8 glucose, buffered to pH 7.2 with KOH. To isolate calcium-sensitive currents, preparations were locally perfused with ECF constructed without CaCl2. All recordings were made at 21–24°C. Leak currents were subtracted offline. Voltages were not corrected for junction potentials.
In situ hybridization.
For whole mounts, embryos were harvested on E7 (stage 30–31) and BPs were dissected in sterile HBSS (Sigma-Aldrich). Both the left and the right BPs were dissected and processed together for in situ hybridization. They were fixed with 4% paraformaldehyde in PBS and subjected to an ethanol dehydration series of 50%, 70%, 90%, and 100% ethanol, followed by brief washes with xylene and then 100% ethanol washes. Digoxigenin (DIG)-labeled antisense probes were synthesized with DIG-RNA labeling mix (Roche). Whole mounts were hybridized with 1 μg antisense probes to Wnt9a or Ser1, transcribed from plasmid templates, for 2 d at 70°C. Samples were incubated overnight at 4°C with alkaline phosphatase (AP)-conjugated anti-DIG antibody (Roche) and washed with Tris buffer, pH 9.5. The AP reaction was developed with nitroblue tetrazolium/5-bromo-4chloro-3-indolyl-phosphate or BM purple (Roche) and postfixed with 4% paraformaldehyde. In situ hybridization protocols and probes for Wnt9a, Ser1, Fzd10, Sfrp2, Wnt7b, and Wnt5b were described previously (Sienknecht and Fekete, 2008, 2009). The Sema3D probe was made using published primers (Jin et al., 2006).
RNA sequencing.
The right otocyst was injected with RCAS(A) parent virus for controls and with RCAS(A)-Wnt9a for experimental subjects. Both groups were killed on E6 (stage 28–30) and the infected right BP was separated from the lagena macula, saccular macula, and the bulk of the AG. Freshly isolated BPs were pooled and homogenized in QIAzol. Three biological samples consisting of a pool of 24 BPs per condition were prepared for sequencing. Total RNA was extracted using the miRNeasy kit (Qiagen). Quality control for total RNA samples with an RNA Integrity Number (RIN) ≥ 8 were validated on an Agilent Technologies Bioanalyzer. cDNA libraries were constructed using the TruSeq Stranded kit (Illumina). The transcriptome was sequenced on an Illumina HiSeq2500 platform with >60 million paired reads per sample. Fragments were aligned to the Galgal4 Ensembl G.gallus genome database version 83.4 using Tophat2 (Kim et al., 2013). DESeq2 (Love et al., 2014) and EdgeR were used to identify differentially expressed genes with a false discovery rate ≤ 0.05 (Robinson et al., 2010). Ingenuity pathway analysis was performed on the Edge R analysis to identify canonical pathways that were overrepresented among the differentially expressed genes. Pathways were considered to be significantly overrepresented if p < 0.05, calculated using Fisher's exact test. The full dataset of differentially expressed genes can be downloaded from the Gene Expression Omnibus database (GEO accession #GSE95295).
Confocal imaging and image processing.
Images of in situ hybridization, dye-labeled AGs, and DAB-stained samples were acquired on a Nikon Eclipse E800 epifluorescence microscope equipped with differential interference contrast optics and a Spot digital camera. Immunolabeled or NeuroVue-labeled whole-mount BPs were imaged on a Nikon Eclipse 90i confocal microscope. Z-stacks of low-magnification images were acquired with a 20×, 0.75 numerical aperture (NA) objective. Z-stacks of high-magnification images were acquired with a 60×, 1.4 NA oil objective. Figures were assembled with Adobe Illustrator and/or Photoshop CS6.
Results
Forced expression of Wnt9a alters radial patterning of the BP
Young HCs are uniformly tall and narrow before diverging into THCs and SHCs by the third week in ovo. Between E12 and E18, abneural HCs become progressively shorter and wider, whereas the neural HCs maintain their shape (Fig. 1) (Fischer, 1992). Long before HC morphologies diverge, radial asymmetries within the BP are present by E6 in cell birthdates (Katayama and Corwin, 1989) and gene expression patterns (Sienknecht and Fekete, 2008). This suggests that secreted factors may orchestrate radial asymmetries on or before E6. We first considered Wnt9a as a candidate ligand for this function based on the timing and distribution of its transcripts and the simultaneous presence of six different Wnt receptors. Wnt9a is localized to the neural half of the BP on E4.5 and, on E5, it withdraws from the sensory organ and is confined to a neural-side nonsensory domain that will become the homogene cells (Sienknecht and Fekete, 2008). It extends along the entire longitudinal axis on E6 (Fig. 2A, arrowhead) and remains there through E10 (data not shown). Therefore, the prosensory domain is poised for Wnt responsiveness before HCs and their associated supporting cells are born.
To test whether exposure to Wnt9a could influence the acquisition of distinct radial cell fates, we used the RCAS(A) retroviral vector to overexpress Wnt9a in the developing otocyst by viral injections on E3. Three days after injection, ectopic Wnt9a transcripts were widely distributed across the BP with a labeling pattern resembling the RCAS viral gag epitope recognized by the 3C2 antibody compared with left controls (Fig. 2B–D). An initial batch of embryos received RCAS(A) parent virus inoculum into the right periotic mesenchyme on E2 (stage 11.5–14) (Hamburger and Hamilton, 1951), followed by injection of RCAS(A)-Wnt9a virus into the lumen of the right otocyst the next day (stage 16–18). This approach should confine Wnt9a-bioactive virus primarily to the otic epithelium (Stevens et al., 2003). Embryos surviving to E18 (7/23) varied between stages 39.5 and 45, indicating developmental delays in some specimens. BPs with Wnt9a overexpression (Wnt9a-OE) were more heavily innervated compared with the left sides (Fig. 2E,F). Strikingly, differentiated HCs appeared to be tall and narrow across the width of the BP after Wnt9a-OE and this was consistent throughout the entire longitudinal axis. This manipulation also produced another qualitatively distinct phenotype in a subset of samples: ectopic vestibular sensory patches. When present, ectopic otoconia were located along the extreme abneural edge of the BP. In one sample, the vestibular identity of a patch underlying these ectopic otoconia was confirmed with an anti-gm2 antibody that labels vestibular but not auditory supporting cells (Fig. 2G,H′; Goodyear et al., 1994). Pre-infection of the mesenchyme with parent virus reduced the incidence of ectopic otoconia from 18% to 7% in Wnt9a-OE BPs. Because vestibular conversion was a relatively modest effect and was reported previously in response to overexpression of activated β-catenin (Stevens et al., 2003), this report is focused instead on the cell fate changes arising within the bulk of the BP that retains an auditory identity.
By E18, Wnt9a-OE BPs were significantly wider throughout all longitudinal positions, with an average increase in width of 36% (Fig. 2I). These wider BPs accommodate a larger number of cells. HC counts were taken from single sections at the 50% longitudinal position through the same BPs quantified in Figure 2I. More HCs were present in the right Wnt9a-OE BPs compared with the uninjected left BPs with a right/left ratio of 1.43 ± 0.25 (mean±SD; n = 7). This difference was statistically significant (one-tailed paired t test, p < 0.001).
HCs differentiate with an increased number of presynaptic ribbons in response to Wnt9a-OE
Because all HCs start out tall and narrow, the apparent absence of the SHC phenotype might simply reflect a developmental delay in response to Wnt9a-OE. To further assess HC differentiation by transmission electron microscopy, we selected a specimen pair in which the right (Wnt9a-OE) BP was obviously wider than the left. Semi-thin cross-sections through the BP at approximately the 50% longitudinal position confirmed the presence of SHCs only on the left side. Nearby thin sections taken from the same blocks showed well developed HCs across the width of both BPs (Fig. 3A). HCs contained electron-dense cytoplasm, a centrally located nucleus, enrichment of mitochondria above the nucleus, and an elongated stereociliary bundle that terminated in a defined cuticular plate. Control SHCs received an efferent terminal. Ribbon foci were identified near the plasma membrane of both neural and abneural HCs. However, a ribbon found in an abneural Wnt9a-OE HC faced a supporting cell rather than an afferent terminal, suggesting the presence of a hemisynapse (Fig. 3A, abneural, right Wnt9a-OE). On the abneural half of infected organs, HCs were indistinguishable from the THCs on the neural side of either the same (right) BP or the opposite (left) BP. Because differentiated THC-like cells were present abneurally and the SHC phenotype was missing, it is possible that high Wnt9a exposure imposed a neural-side fate throughout the organ, with the THC morphology being one obvious readout of this identity shift.
THCs are more densely packed and normally form significantly more presynaptic ribbons than SHCs. These features extended to the abneurally located THC-like cells in Wnt9a-OE BPs. Both neural THCs and abneural THC-like cells contained excess ribbons, identified by immunostaining for CtBP2 (Fig. 3B). Figure 3B′ indicates schematically the sampling regions used for quantitative analysis of ribbon numbers from the BPs of three embryos. Two of three embryos sampled at the 66% position from the base (distal) showed statistically significant increases in the number of ribbons per HC on both the neural and abneural sides (Fig. 3C); this includes Embryo 1 with a portion of the distal sample regions shown in Figure 3B. Likewise, at a more proximal (33–40%) position, all three embryos showed increases in the number of ribbons per HC on the neural side, and two of three embryos showed increases on the abneural side (Fig. 3C). In total, only 3 of the 12 left–right comparisons did not show changes in ribbon numbers in reponse to Wnt9a-OE. Because there was no systematic localization of these unresponsive regions, they may reflect insufficient viral infection; this could not be verified due to CtBP2 antibody incompatibility with the 3C2 antibody. The observation that neural-side HCs had more ribbons than normal is notable. This suggests either that Wnt9a-OE created a “super” THC phenotype on the neural side that was perhaps linked to the increased afferents or that synapse maturation was altered. Qualitatively, the numerous ribbon foci that emerged in the Wnt9a-OE condition appeared smaller than in the uninjected controls.
Wnt9a-OE alters outward currents in abneural HCs
Whole-cell patch-clamp recordings of cells from E18 BPs showed that THCs (neural) and SHCs (abneural) exhibited characteristic electrical features in controls. THCs expressed two major outward potassium currents, including a slow delayed rectifier current and a fast, calcium-sensitive potassium current that appears after E18 (Fuchs and Sokolowski, 1990). In contrast, SHCs exhibit an inactivating A-type potassium current. We observed that outward potassium currents from HCs along the neural edge were rapidly activating regardless of treatment condition. Transient A-type currents were found in control abneural HCs. Interestingly, both of these electrical signatures were found among abneural HCs upon Wnt9a-OE (Fig. 4A). Overall, the magnitude and voltage sensitivity of steady-state currents of abneural HCs from Wnt9a-OE BPs were more similar to THCs of either treatment condition (Fig. 4B, left). Steady-state current amplitudes elicited by voltage steps to +20 mV were significantly different across treatment groups (Table 1). Specifically, the average current amplitude at this voltage step was 3.8-fold greater in Wnt9a-OE abneural HCs than in control abneural SHCs, but not significantly different from the control or Wnt9a-OE neural HCs. For the fast-activating potassium current (IKf), the amplitudes in Wnt9a-OE abneural HCs increased by 1.7-fold compared with control abneural SHCs (difference not significant). Interestingly, the fast current amplitudes in Wnt9a-OE abneural THC-like cells were also not significantly different from control neural THCs (Fig. 4B, right; Table 2). Cells were exposed to extracellular solution with and without CaCl2 to further examine acquisition of IKf because this current arises from BK-type potassium channels that are activated by the influx of Ca2+ from voltage-dependent calcium channels (Jiang et al., 1997; Li et al., 2009). Removal of extracellular calcium reduced the amplitude of the outward current and slowed the activation time constant for HCs in Wnt9a-OE BPs located along both the neural and abneural aspects of the BP (Fig. 4C). Therefore, abneural HCs appear to demonstrate intermediate characteristics between normal THCs and SHCs.
Wnt9a-OE alters efferent neurite patterning across the radial axis of the BP
To observe efferent innervation, axons were anterogradely labeled by the insertion of lipophilic dyes near their cell bodies (Whitehead and Morest, 1981) on E18 (stages 39.5–45). AGs and BPs were separated for imaging. In controls, efferents traveled longitudinally as a coherent bundle along the abneural edge of the ganglion, split into several major branches, and then turned to penetrate the BP (Fig. 5A; Zidanic, 2002). Efferents were mostly confined to the proximal half of the BP as expected (Fig. 5B). Small-caliber fibers with numerous fine boutons projected into the neural side. Thicker and straighter radial fibers and bundles advanced abneurally (Fig. 5D, arrow). Because our analysis was confined to embryonic stages, we did not observe the large, cup-shaped terminals onto SHCs described for 3-week-old chicks (Zidanic, 2002). To correlate efferent trajectories with HC types, we postprocessed for CtBP2. There was an abrupt change in ribbon density that correlated with the pentration of axons into the BP at the so-called habenula perforata; this obviously segregated the normal BP into two distinct domains: neural and abneural (Fig. 5C,D).
Wnt9a-OE often reduced the efferent bundle within the SAG and yielded fewer secondary branches that did not extend as far in the neural direction (Fig. 5E). Right BPs showed 3 categories of phenotypic alteration: no change (n = 6), moderate (n = 8), and extreme (n = 7). The habenula perforata in the moderate phenotype was disorganized, broader, and shifted abneurally (Fig. 5F). Higher magnification showed that the BP had fewer radial bundles that were thinner and some of these were misdirected orthogonally (Fig. 5H, arrowhead). Ribbon foci were increased across the radial axis with one marked difference on the neural side: the puncta were smaller (Fig. 5G,H). Compared with the left control sides, BPs with the extreme phenotype had few to no efferents within a widened BP (Fig. 5I–M).
Wnt9a-OE induces changes in gene expression by E6
To begin to understand the mechanisms underlying the changes in cell fates seen near the end of embryogenesis resulting from ectopic Wnt9a delivery, we looked at earlier stages of development. The earliest detected phenotypic change was broader innervation on E6 within 3 d after viral injections (Fig. 6A). This suggests that, by E6, the afferents are already responding to some change in the environment of the Wnt9a-OE BP that now renders the abneural side an attractive target. To reveal changes in gene expression at this time point, RNA sequencing was performed on right BPs isolated from E6 Wnt9a-OE and RCAS(A)-injected animals. A total of 500 differentially regulated genes was found (296 upregulated and 204 downregulated; p < 0.05). Wnt9a transcripts increased significantly, validating the efficacy of gene delivery.
RNA sequencing showed that several axon guidance molecules were among the genes with significant changes in expression on E6. For example, the expression of Sema3d, a repulsive guidance cue, was decreased (Table 3). In adjacent sections, we examined its spatial localization by in situ hybridization; Sema3d was expressed across the prosensory epithelium in control BPs, whereas Sema3d was abolished in Wnt9a-OE BPs (Fig. 6A). We previously found no evidence for a direct outgrowth-promoting effect of Wnt9a on E4 otic neurons using an established in vitro assay (Fantetti, 2011). Nevertheless, several genes linked to cell motility showed level changes, such as RhoA, calcium signaling molecules (Wipf1, Rhpn1, Arpc4, Actc1, Rac3, Pak3, Casq2, Iptr3), and axon guidance molecules (Bdnf, Nrtk2, Sema3d; Table 3). Fgf19 and Bmp7, which promote statoacoustic ganglion outgrowth (Fantetti and Fekete, 2012), were upregulated (Table 3).
RNA sequencing also revealed several Wnt-related genes that were expressed differentially, including Fzd10, Sfrp2, and Wnt7b (Table 3). In situ hybridization performed on E6 showed that expression of Fzd10, an abneural transcript, was abolished. Sfrp2 is most strongly expressed in the center of the BP and in the future homogene cell region (where Wnt9a is expressed), but was reduced only in the sensory domain. Ser1/Jag1 is a known Wnt target gene, so it serves as a positive control. Ser1 is normally expressed in a gradient across the BP, being stronger on the neural side. Unexpectedly, it was unchanged by Wnt9a-OE in the E6 RNA-sequencing dataset and was only modestly increased by in situ hybridization on E6 (Fig. 6A). An overall increase in Ser1 transcript levels was more obvious on E7 (Fig. 6B). In whole mounts, the signal strength for Ser1 on the Wnt9a-OE abneural side appeared comparable to the neural side of the left BP with isolated patches even stronger (arrow). Nonetheless, a neural-high to abneural-low gradient persisted in the context of strong Wnt9a-OE, suggesting that at least some neural–abneural patterning information was still present in the treated samples on E7.
In situ hybridization was also performed on E10 to determine whether the gene expression changes were maintained at a later stage, but before HC types were morphologically distinguishable (Fig. 6C). The abneural/central genes, Wnt7b, Fzd10, and Sfrp2 were all abolished in the sensory domain in Wnt9a-OE samples, whereas nonsensory expression of Sfrp2 and Wnt5b was unchanged. Ser1 expression was increased, but still displayed a gradient that was highest neurally. Although RNA sequencing showed that Wnt7b transcripts were present on E6 and were downregulated by Wnt9a-OE, they could not be detected using in situ hybridization at earlier stages; however, they were clearly reduced on the treated sides on E10.
Wnt9a-OE increases proliferation during the same time window of differential gene expression
We found that the superfluous cells populating Wnt9a-OE BPs resulted from both enhanced proliferation and a prolonged period of cell division (Fig. 7A). When RCAS(A)-Wnt9a was injected into E3 otocysts (without an E2 preinjection of parent virus) and analyzed on E6–E10, there was a delay of ∼2–3 d in HC differentiation (Fig. 7B). This may be a reflection of the general delay in cell cycle withdrawal, allowing additional time for cell fate reprogramming by Wnt9a.
Discussion
The THC-only phenotype resulting from Wnt9a-OE can be interpreted in light of the asymmetric source of Wnt9a mRNA along the neural edge of normal BPs. Via either a short-range signal or dilution of Wnt signaling propogated directionally through dividing cells (Farin et al., 2016), Wnt9a may function to subdivide the BP into neural versus abneural fates. This signaling would directly upregulate Wnt target genes (e.g., Ser1), as was observed by E6 in the BP based on neural–abneural asymmetries in gene expression. The reduction of both central (Sfrp2) and abneural genes (Fzd10 and Wnt7b) on E6, in response to the misexpression of Wnt9a, may reflect the loss of signaling by an abneural inducer or it may simply result from the absence of abneural fates. By E18, the entire organ was populated by THCs and THC-like cells. There is prior evidence for related morphological and transcriptional changes to Wnt-activated embryonic mammalian cochleas: pharmacological activation of the Wnt pathway promotes the specification of the medial domain containing IHCs while suppressing formation of the lateral-domain-containing OHCs (Munnamalai and Fekete, 2016).
Together, results from mouse and chicken suggest some evolutionary conservation for Wnts to control the size and identity of the neural/medial compartment of the cochlea while suppressing the abneural/lateral compartment. This raises the question of how the organs respond to the absence or reduction of a Wnt signal during the critical period when radial asymmetry is being established. Embryonic mouse cochleas cultured in the presence of a Wnt inhibitor (IWR-1) are smaller due to a 40% reduction in the medial compartment (Munnamalai and Fekete, 2016). In chicken, E7.5 BPs cultured with IWR-1 are also smaller (Jacques et al., 2014), although radial identities (THC/SHC fates) were not evaluated and would have required longer culture times to differentiate HC types. Because a pharmacological approach to blocking Wnt signaling is insufficient to reveal the role of a particular Wnt ligand, we attempted to specifically knock down Wnt9a using short-hairpin RNAs encoded in the introns of viral or plasmid vectors. However, we were unable to detect a reduction in Wnt9a. We also examined E18.5 Wnt9a-KO mouse cochleas (gift from Christine Hartmann), but saw no disruptions in cochlear patterning. Low levels of expression of Wnt9a in the embryonic mouse cochlea (Shen et al., 2015) could explain why it is apparently not required for radial patterning in this species. Therefore, birds and mammals have likely diverged as to which Wnt ligands influence radial axis patterning.
The challenge in deciphering whether enhanced Wnt9a signaling was directly responsible for an abneural-to-neural conversion is that Wnt signaling can also stimulate proliferation in developing auditory organs (Jacques et al., 2014; Munnamalai and Fekete, 2016), thus potentially expanding the neural-side population. Indeed, we detected increased proliferation with Wnt9a-OE. However, increasing the size of the organ fails to explain the absence of SHCs. During normal development of the BP, cell cycle exit progresses from neural to abneural and from proximal to distal (Katayama and Corwin, 1989). These directional gradients were unaffected by Wnt9a-OE: this is shown for the radial axis in Figure 7A (dividing cells are mostly confined to the abneural half of the BP, similar to untreated BPs) and was also observed for the longitudinal axis because residual mitotic cells in the BP on E10 were confined to the apex (data not shown, n = 3 samples). It is not known what enforces cell cycle arrest, but it occurs even in the face of enhanced Wnt activity. Several abneural marker genes that are missing in the context of Wnt9a-OE are normally restricted to the abneural side even before those cells pull out of division. This indicates that the neural–abneural fate decision is being made by cycling progenitors rather than postmitotic cells, thus mitigating the relevance of Wnt's proliferative effect. This leads us to posit that Wnt9a promotes neural-side identity directly and/or blocks the acquisition of abneural fates.
We showed previously that overexpression of constitutively active β-catenin in the BP generated ectopic vestibular patches at various locations across its radial axis (Stevens et al., 2003). Wnt9a-OE also stimulated vestibular fates, but only on the extreme abneural side of the BP. The spatial constraints on vestibular conversion in the context of ligand overexpression, but not when intracellular Wnt signaling is enhanced, may be related to the combined spatial expression of Fzd receptors and Wnt inhibitors (Sienknecht and Fekete, 2008, 2009) that will affect the bioavailability of Wnt ligands. Moreover, whereas robust Wnt signaling may first specify vestibular fates, it could be repurposed for radial patterning after commitment of the BP to an auditory fate.
IKf (BK) is a distinguishing feature of THCs and appears after E18, reaching mature levels ∼E19, only days before hatching. Similarly, IA, a transient outward potassium current, is characteristic of SHCs (Fuchs and Sokolowski, 1990; Tan et al., 2013). Interestingly, Wnt9a-OE resulted in an intermediate phenotype for some abneural HCs that exhibited features of both THCs and SHCs. Perhaps this reflects a delay in BK acquisition because we recorded at a stage when BK channels begin to be expressed (E18–E20). Alternatively, Wnt9a-OE may alter the complement of channels appropriate to the tonotopic position. BK macroscopic conduction and activation/deactivation kinetics vary systematically along the longitudinal axis (Pantelias et al., 2001; Duncan and Fuchs, 2003), so differences in calcium sensitivity and amplitude could mask underlying radial fate changes. Conversely, ∼25% of the Wnt9a-infected abneural HCs exhibited A-type currents, clearly demonstrating a SHC phenotype and suggesting that there is a mixed response among the Wnt9a-OE BPs. Indeed, labeling of ribbons in HCs also suggests that not all abneural domains were fully converted in all animals.
Among the changes found when Wnt9a was overexpressed was an increase in THC numbers and, based on ultrastructural morphology, a complete absence of SHCs by E17. Both sides of the BP showed an increased number of ribbons per HC. This implies that excess Wnt9a may convert neural-side THCs into “super” THCs, although effects on ribbon synapse maturation could also underlie these changes. Manipulation of intracellular calcium has been shown to affect the size of ribbon foci in HCs (Sheets et al., 2012) and two differentially expressed genes (Casq2 and Itpr3) may be responsible for increasing intracellular calcium, resulting in smaller ribbons in Wnt9a-OE BPs.
In the normal BP, spatial targeting of afferent axons predominantly to the THC side (Fischer, 1992) is already evident on E6. Upon Wnt9a-OE, the axons spread abneurally by E6. The radial spread of axons is sustained on E18 and mirrored by presynaptic changes (ribbon counts), although the latter may be an independent phenomenon. Increased abneural projections may reflect Wnt9a-induced changes in molecules that stimulate auditory axon outgrowth: FGFs, BMP7, or survival neurotrophins (Garrido et al., 1998; Fantetti and Fekete, 2012). Such genes that are elevated upon Wnt9a-OE include Fgf3, Fgf19, and Bmp7. However, RNA-sequencing data showed a downregulation of both Bdnf and its receptor, Nrtk2, upon Wnt9a-OE. Therefore, this outgrowth of afferent axons is independent of BDNF but perhaps dependent on FGFs (Appler and Goodrich, 2011), BMPs, or other growth factors.
Several of the genes differentially regulated by Wnt9a-OE are involved in axon guidance, calcium signaling, F-actin cytoskeletal remodeling, or other cellular processes that correlate with axon outgrowth. Sema3d downregulation was notable because, in the mammalian cochlea, type I afferents that explore the OHC region undergo retraction and become restricted to the neural side due in part to Sema3f-mediated repulsion (Coate et al., 2015; Druckenbrod and Goodrich, 2015). Consistent with this, afferents spread into abneural territory when Sema3d expression was abolished by Wnt9a-OE. However, we found that Sema3d transcripts are normally present across the width of the BP on E6 even at locations where afferents are present; this is inconsistent with a repulsive effect on afferents at this age. Several genes that are components of RhoA-mediated axonal steering were also downregulated. These could be downstream of a loss in Semaphorin-mediated growth cone repulsion. Alternatively, these genes may be regulated directly by the Wnt pathway in afferent neurons, leading to select neurites losing the ability to steer away from the abneural half. Another possibility is that some of these F-actin-remodeling proteins are required for afferent axon pruning. In the chicken, excess branches are pruned as synapses mature (Whitehead and Morest, 1985a). In mammals, the retraction of afferents from immature OHCs coincides with the arrival of efferent neurites (Pujol, 1986).
At present, we are lacking a mechanistic understanding of why efferent trajectories are abnormal after Wnt9a-OE. The proposed abneural-to-neural fate switch may alter an efferent guidance cue or the efferents may arrive on schedule but subsequently retract in the absence of their normal SHC target cells. During normal development, processes that innervate the abneural edge traverse enlarged extracellular spaces, which provide a pathway for axons to reach SHCs (Whitehead and Morest, 1985a, 1985b). Further investigation will be necessary to determine whether overcrowding the target tissue with cells and afferents resists the penetration of efferents into the epithelium.
In summary, this study shows that ectopic Wnt9a-mediated signaling not only increases cell proliferation, but also influences the specification of radial identities in the BP at a time when many HCs and supporting cells are still proliferating. This imposition of regional identity is associated with changes in spatial gene expression and can affect axonal projections and later differentiation of HCs into two distinct types. This is the first study to identify a specific Wnt ligand that can have such a profound effect on radial identity in the developing cochlea.
Footnotes
This work was supported by the Collaborative Core for Cancer Bioinformatics, the Purdue Statistical Consultant Services (Department of Statistics), the Purdue University Center for Cancer Research [National Institutes of Health (NIH) Grant P30CA023168], the Indiana University Simon Cancer Center (NIH Grant P30CA082709), the Walther Foundation, the Purdue Bioinformatics Core, the National Institute on Deafness and Other Communication Disorders (Grant R01DC002756 to D.M.F. and Grant P30DC005188 to R.K.D.), and the Deutsche Forschungsgemeinschaft (Grant SPP1608 SI823/3-1 to U.J.S.). We thank Christine Hartmann for providing mouse Wnt9a-KO cochleas and Lisa Goodrich for critiquing an earlier version of this manuscript and Qing Zhou, Pin-Chao Liao, and Kyle Patton for technical assistance.
The authors declare no competing financial interests.
- Correspondence should be addressed to Donna M. Fekete, Department of Biological Sciences, Purdue University, 915 W. State St., West Lafayette, IN 47907-2054. dfekete{at}purdue.edu