Insights into electrosensory organ development, physiology and evolution from a lateral line-enriched transcriptome
Abstract
The anamniote lateral line system, comprising mechanosensory neuromasts and electrosensory ampullary organs, is a useful model for investigating the developmental and evolutionary diversification of different organs and cell types. Zebrafish neuromast development is increasingly well understood, but neither zebrafish nor Xenopus is electroreceptive and our molecular understanding of ampullary organ development is rudimentary. We have used RNA-seq to generate a lateral line-enriched gene-set from late-larval paddlefish (Polyodon spathula). Validation of a subset reveals expression in developing ampullary organs of transcription factor genes critical for hair cell development, and genes essential for glutamate release at hair cell ribbon synapses, suggesting close developmental, physiological and evolutionary links between non-teleost electroreceptors and hair cells. We identify an ampullary organ-specific proneural transcription factor, and candidates for the voltage-sensing L-type Cav channel and rectifying Kv channel predicted from skate (cartilaginous fish) ampullary organ electrophysiology. Overall, our results illuminate ampullary organ development, physiology and evolution.
https://doi.org/10.7554/eLife.24197.001Introduction
The lateral line system of fishes and aquatic amphibians is a good model for studying the diversification of different organs and cell types, both in development and evolution. In jawed vertebrates, this sensory system ancestrally includes mechanosensory neuromasts and electrosensory ampullary organs (‘ampullae of Lorenzini’), both of which develop - together with their afferent neurons - from individual embryonic lateral line placodes (Northcutt et al., 1995; Modrell et al., 2011a; Gillis et al., 2012). (The jawless lampreys have neuromasts and electroreceptors, but the latter are collected in ‘end buds’ at the surface, rather than recessed in ampullary organs, and their embryonic origin is unknown; the lateral line system of hagfishes, which lack electroreceptors altogether, is secondarily reduced; Braun, 1996; Braun and Northcutt, 1997.) The electrosensory division of the lateral line system was lost independently in the lineages leading to extant neopterygian fishes (gars, bowfin and teleosts) and to anuran amphibians (neither of the major anamniote lab models, i.e., the teleost zebrafish and the frog Xenopus, has electroreceptors). However, electrosensory lateral line organs evolved independently at least twice within the teleosts, most likely from neuromast hair cells (Bullock et al., 1983; Northcutt, 1986; Bodznick, 1989; Alves-Gomes, 2001; Bodznick and Montgomery, 2005; Kawasaki, 2009; Baker et al., 2013). The entire lateral line system was lost in amniotes, with the transition to life on land.
The loss of the electrosensory division of the lateral line system in different vertebrate lineages shows that ampullary organ development must be genetically separable from neuromast development. Indeed, even within the same lateral line placode-derived sensory ridge, neuromasts form first, along the center of the ridge, while ampullary organs form later, on the flanks (Schlosser, 2002; Northcutt, 2005a; Piotrowski and Baker, 2014). Neuromasts and ampullary organs are morphologically distinct: neuromasts contain mechanosensory hair cells, plus supporting cells that secrete a gelatinous cupula; ampullary organs comprise a sensory epithelium of electroreceptor and supporting cells located at the base of a duct filled with conductive jelly, leading to a surface pore (Northcutt, 1986; Jørgensen, 2005; Baker et al., 2013). Ampullary electroreceptor cells generally have an apical primary cilium but either no or few apical microvilli, which are not organized into the ‘hair bundle’ (stair-case array) that characterizes hair cells (Northcutt, 1986; Jørgensen, 2005; Baker et al., 2013). The molecular mechanisms underlying neuromast formation from the migrating posterior lateral line primordium in zebrafish have been intensively studied (Chitnis et al., 2012; Piotrowski and Baker, 2014; Thomas et al., 2015), but our molecular understanding of ampullary organ development is very limited.
Like hair cells (and also retinal and pineal photoreceptors, and retinal bipolar neurons), all vertebrate electroreceptors have electron-dense pre-synaptic bodies surrounded by a large pool of synaptic vesicles (Northcutt, 1986; Jørgensen, 2005). Such ‘ribbon synapses’ respond to graded signals and are capable of sustained neurotransmitter release (Matthews and Fuchs, 2010; Pangršič et al., 2012; Safieddine et al., 2012; Nicolson, 2015; Wichmann and Moser, 2015; Moser and Starr, 2016). While the specific proteins involved in ribbon synapse function are increasingly understood in retinal photoreceptors and hair cells (Matthews and Fuchs, 2010; Pangršič et al., 2012; Safieddine et al., 2012; Nicolson, 2015; Wichmann and Moser, 2015; Moser and Starr, 2016), all that is known about neurotransmission at non-teleost electroreceptor ribbon synapses is from work in dissected skate (cartilaginous fish) ampullary organs, which showed that activation of L-type voltage-gated calcium channels results in the release of a ‘glutamate-like’ neurotransmitter (Bennett and Obara, 1986).
Similarly, our only detailed understanding of non-teleost ampullary organ physiology until very recently had come from current- and voltage-clamp approaches to study epithelial currents in dissected single ampullary organ preparations from skates (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005). These revealed that L-type voltage-gated calcium channels are required both for voltage-sensing in the apical (lumenal, i.e., exterior-facing) electroreceptor membrane, and for neurotransmitter release basally. The basal membrane is repolarized by voltage-gated potassium channels and calcium-dependent chloride channels, while the apical membrane is repolarized by the calcium-gated potassium channel BK (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005; King et al., 2016). BK was recently cloned directly from skate ampullary organs (King et al., 2016), a little over 40 years after its properties were initially discovered using the same preparation (Clusin et al., 1975; Clusin and Bennett, 1977a; Clusin and Bennett, 1977b). While the current manuscript was under review, whole-cell patch-clamp experiments on dissociated electroreceptors from adult skates, together with transcriptome profiling, revealed that the L-type voltage-gated calcium channel Cav1.3 mediates the low-threshold voltage-activated inward current, and works together with BK to mediate electroreceptor membrane oscillations (Bellono et al., 2017). Other than this recent exciting advance (Bellono et al., 2017), the specific ion channels and subunits involved in electroreceptor function have not been identified in any vertebrate.
In short, our understanding of the specific molecular basis of both vertebrate electroreceptor development and physiology is still rudimentary. The few candidate gene approaches reported thus far have identified some transcription factors, and a few other genes, expressed in both ampullary organs and neuromasts in larval axolotl (Metscher et al., 1997; Modrell and Baker, 2012), paddlefish (Modrell et al., 2011a, 2011b; Butts et al., 2014), shark and skate (Freitas et al., 2006; Gillis et al., 2012). However, a candidate gene approach will not identify genes important specifically for ampullary organ development or function. Here, we report an unbiased transcriptomic approach in paddlefish (a non-teleost chondrostean fish) to identify such genes, which has yielded novel and wide-ranging insights into the development, physiology and evolution of non-teleost ampullary organs.
Results
We generated transcriptomes at stage 46 (the onset of independent feeding; Bemis and Grande, 1992) from pooled paddlefish opercula (gill-flaps), which are covered in ampullary organs plus some neuromasts, versus fins, which have a generally similar tissue composition but no lateral line organs at all (Modrell et al., 2011a, 2011b). Differential expression analysis yielded 490 genes, excluding duplicates, enriched at least 1.85-fold (log2fold 0.89) in operculum versus fin tissue, hereafter designated ‘lateral line-enriched’ (Supplementary file 1). Of these genes, 112 were uncharacterized loci or could only be assigned to a protein family, while a further 44 were described as being ‘like’ a specific gene, leaving 334 assigned genes (Supplementary file 1). Figure 1 shows a molecular function analysis using Gene Ontology (GO) terms, where available.
Conservation in ampullary organs of transcription factors critical for mechanosensory hair cell development
The lateral line-enriched dataset (Supplementary file 1) includes several transcription factor genes whose expression we had previously reported in both ampullary organs and neuromasts in larval paddlefish, suggesting that the differential expression analysis had been successful. These were: the homeodomain transcription factor genes Six1 and Six2 (Modrell et al., 2011a); the HMG-domain SoxB1 class transcription factor gene Sox3 (Modrell et al., 2011b); and the basic helix-loop-helix (bHLH) transcription factor gene Atoh1, whose lateral line organ expression at stages 40–45 was noted, in passing, in a study on cerebellum development (Butts et al., 2014).
Atoh1 is essential for hair cell formation (Bermingham et al., 1999). Paddlefish Atoh1 is 18.6-fold lateral line-enriched (Supplementary file 1), and expressed by stage 33 within the otic and preopercular neuromast canal lines (Figure 2A). It is expressed within most canal lines and the migrating posterior lateral line primordium by stage 37 (Figure 2B), and in the developing ampullary organ fields by stage 39 (Figure 2C), where it persists at stage 46 (Figure 2D).
In mouse cochlear explants, Six1 and its co-factor Eya1, acting cooperatively with SoxB1 subfamily member Sox2, are sufficient to induce Atoh1 (Ahmed et al., 2012). Like Six1, Eya1 is also expressed in paddlefish ampullary organs (Modrell et al., 2011a), while Sox2 is 3.3-fold lateral line-enriched (Supplementary file 1). Immunostaining with a cross-reactive anti-Sox2 antibody revealed Sox2 expression in developing ampullary organs as well as neuromasts (Figure 2E–H), consistent with Sox2 mRNA expression by in situ hybridization (not shown). Sox2 represses Atoh1 in supporting cells in the mouse cochlea while, conversely, Atoh1 represses Sox2 in developing hair cells (Dabdoub et al., 2008). By stage 46, paddlefish Atoh1 expression is hard to detect in canal neuromasts - perhaps partly due to strong expression in scattered overlying epidermal cells, presumably Merkel cells (Maricich et al., 2009; Whitear, 1989) - but is still strong in ampullary organs (Figure 2I). Sox2 immunoreactivity remains strong in both neuromasts and ampullary organs at stage 46 (Figure 2J), where it is most likely restricted to supporting cells, given its punctate staining pattern.
Of the 66 genes in the lateral line-enriched dataset that encode known transcription factors (three other transcription factor genes in the dataset were assigned only to families) (Supplementary file 1), our top candidate was Pou4f3 (Brn3c, DFNA15), which is 27.5-fold lateral line-enriched (Supplementary file 1). Pou4f3 is a confirmed Atoh1 target in hair cells (Masuda et al., 2011; Ikeda et al., 2015) and required for hearing (Costaridis et al., 1996). In mouse cochlear explants, Six1 and Eya1 - acting independently of Atoh1 - are sufficient to induce Pou4f3, which promotes hair cell differentiation (Ahmed et al., 2012). Like Atoh1, paddlefish Pou4f3 is expressed in ampullary organs, as well as neuromasts (Figure 2K,L).
All three SoxB1 subfamily genes are putative Atoh1 targets in the postnatal mouse cerebellum (Klisch et al., 2011) and present in the lateral line-enriched dataset. Sox3 is 5.2-fold lateral line-enriched (Supplementary file 1); as noted above, we previously reported its expression in paddlefish ampullary organs and neuromasts (Modrell et al., 2011b). Sox1 is 8.6-fold lateral line-enriched (Supplementary file 1) and expressed in ampullary organs from their eruption (Figure 2M–P). It is transiently expressed in the migrating posterior lateral line primordium (Figure 2M,N) and, at later stages, expression is also seen in neuromast canal lines (Figure 2O,P).
Taken together, these data support a high degree of conservation in the transcriptional network regulating the development of both lateral line organ types, including the likely requirement of Atoh1 for electroreceptor as well as hair cell formation.
The proneural transcription factor gene Neurod4 is expressed in ampullary organs but not neuromasts
Another highly lateral line-enriched transcription factor gene is Neurod4 (Ath3, NeuroM) (18.9-fold enriched; Supplementary file 1), a putative cerebellar Atoh1 target (Klisch et al., 2011). In mice, this atonal-related proneural bHLH transcription factor gene is expressed in the brain, spinal cord, retina, trigeminal ganglia and dorsal root ganglia (Takebayashi et al., 1997). It is required for the survival of cerebellar granule cell precursors (Tomita et al., 2000), and cooperates with other bHLH and/or homeodomain transcription factors to regulate the development of retinal bipolar and amacrine cells (see Hatakeyama and Kageyama, 2004), and trigeminal and facial branchiomotor neurons (Ohsawa et al., 2005). In zebrafish, Neurod4 is important for olfactory neuron development (Madelaine et al., 2011); it is transiently expressed during chicken otic neurogenesis, but not hair cell formation (Bell et al., 2008). Although the related gene Neurod1 seems to be important for zebrafish neuromast hair cell differentiation (Sarrazin et al., 2006), we previously showed that paddlefish Neurod1 is not expressed in lateral line organs, only in cranial sensory ganglia (Modrell et al., 2011b). Neurod4 is the only Neurod family member in the lateral line-enriched dataset. At early stages, it is expressed in the brain, olfactory epithelium, eyes and trigeminal ganglion (Figure 3A). The earliest lateral line expression of Neurod4 is in developing ampullary organ fields around stage 39 (Figure 3B), where it persists (Figure 3C–F). Expression is never observed in developing neuromast lines. No other gene has previously been reported to show differential expression in ampullary organs and neuromasts. Hence, Neurod4 may be important for specifying ampullary organ and/or electroreceptor fate.
We went on to examine the expression of transcription factor genes with known links to Neurod4 in other cell types. In trigeminal neurons, Pou4f1 (Brn3a) directly represses Neurod4 (Lanier et al., 2007). Pou4f1 is a putative Atoh1 target in the postnatal mouse cerebellum (Klisch et al., 2011), and is 20.6-fold lateral line-enriched (Supplementary file 1). Expression was seen in both ampullary organs and neuromasts (Figure 3G), suggesting that additional factors besides Pou4f1 are likely involved in controlling the differential expression of Neurod4 in ampullary organs versus neuromasts.
In the spinal cord, Neurod4 cooperatively interacts with the LIM homeodomain transcription factor Lhx3 to specify motor neurons (Lee and Pfaff, 2003). Lhx3 is expressed in all inner ear hair cells, although its expression is differentially regulated by Pou4f3 in cochlear versus vestibular hair cells (Hertzano et al., 2007); it is also a putative cerebellar Atoh1 target gene (Klisch et al., 2011). Lhx3 is 16.3-fold lateral line-enriched (Supplementary file 1), and proved to be expressed in both ampullary organs and neuromasts (Figure 3H). Hence, it is possible that Lhx3 could interact with Neurod4 in paddlefish ampullary organs to specify electroreceptors, and with another partner in neuromasts to specify hair cells.
Finally, we examined the zinc finger transcription factor gene Myt1, another putative cerebellar Atoh1 target gene (Klisch et al., 2011) that is 6.1-fold lateral line-enriched (Supplementary file 1). We selected Myt1 because its expression is upregulated in Xenopus embryos by Neurod4 (Perron et al., 1999; Hardwick and Philpott, 2015), and it synergizes with Neurod4 to induce neuronal differentiation in this system (Perron et al., 1999). Furthermore, Myt1 transcripts are reportedly enriched in both cochlear and vestibular hair cells in postnatal mice (Elkon et al., 2015). We found expression of paddlefish Myt1 in both ampullary organs and neuromasts (Figure 3I).
Overall, our results suggest that the transcription factor networks underlying hair cell development, which center around Atoh1, are likely to be active in paddlefish ampullary organs as well as neuromasts. This suggests significant conservation between the molecular mechanisms underlying hair cell and electroreceptor development. Furthermore, our unbiased RNA-seq approach has also identified the first transcription factor gene expressed in ampullary organs but not neuromasts, Neurod4. Given the importance of members of the proneural bHLH Neurod family for specifying cell fate, we suggest that Neurod4 may be involved in specifying ampullary organ and/or electroreceptor fate in paddlefish.
Hair cell ribbon synapse genes are also expressed in ampullary organs
Some of the most highly lateral line-enriched genes in our paddlefish dataset are required for synaptic transmission in hair cells, which occurs at specialized ‘ribbon synapses’ characterized by electron-dense, pre-synaptic structures called ‘synaptic ribbons’ (Matthews and Fuchs, 2010; Pangršič et al., 2012; Safieddine et al., 2012; Nicolson, 2015; Wichmann and Moser, 2015; Moser and Starr, 2016). These tether glutamate-filled synaptic vesicles and stabilize L-type voltage-gated calcium channels at the plasma membrane (Cav1.3 in hair cells; Cav1.4, in retinal photoreceptors; Joiner and Lee, 2015), enabling rapid and sustained glutamate release in response to activation of these calcium channels by membrane depolarization (Matthews and Fuchs, 2010; Pangršič et al., 2012; Safieddine et al., 2012; Nicolson, 2015; Wichmann and Moser, 2015; Moser and Starr, 2016). Electroreceptors also have synaptic ribbons of varying morphology (Northcutt, 1986; Bodznick and Montgomery, 2005): in paddlefish, they were described as synaptic ‘sheets’ (Jorgensen et al., 1972). In dissected skate ampullary organ preparations, activation of L-type voltage-gated calcium channels in the basal electroreceptor membrane results in release of a ‘glutamate-like’ neurotransmitter (Bennett and Obara, 1986).
In hair cells, glutamate is loaded into synaptic vesicles by the vesicular glutamate transporter Vglut3, which is encoded by Slc17a8 (DFNA25) and essential for hair cell synaptic transmission in mouse (Ruel et al., 2008; Seal et al., 2008) and zebrafish (Obholzer et al., 2008). This is unusual: Vglut1 or Vglut2 are used at glutamatergic synapses in the central nervous system, and at photoreceptor and bipolar cell ribbon synapses (see e.g. Zanazzi and Matthews, 2009; Pangršič et al., 2012). Slc17a8 is one of the most highly enriched genes in our paddlefish dataset (30.9-fold enriched; Supplementary file 1), and proved to be expressed in both ampullary organs and neuromasts (Figure 4A). Hence, synaptic vesicles in electroreceptors are likely to be loaded by the same vesicular glutamate transporter as hair cells. Furthermore, this also provides independent evidence that the afferent neurotransmitter released by non-teleost electroreceptors is indeed glutamate, as in hair cells (and photoreceptors), as suggested by electrophysiology experiments on dissected skate ampullary organs (Bennett and Obara, 1986).
Hair cells are thought to be unique in depending on the multi-C2 domain transmembrane protein otoferlin for synaptic vesicle exocytosis (Yasunaga et al., 1999; Roux et al., 2006; Pangršič et al., 2010; Chatterjee et al., 2015; Strenzke et al., 2016; Vogl et al., 2016), rather than neuronal SNAREs (Nouvian et al., 2011). This contrasts not only with conventional synapses but also with all other ribbon synapses (see e.g. Zanazzi and Matthews, 2009; Mercer and Thoreson, 2011). Otoferlin is a type II ferlin (Lek et al., 2012) encoded by Otof (DFNAB6, DFNAB9), which is 20.9-fold lateral line-enriched (Supplementary file 1). Otof is expressed in both ampullary organs and neuromasts in paddlefish (Figure 4B). This suggests that synaptic vesicle exocytosis at the electroreceptor ribbon synapse, just as at the hair cell ribbon synapse, is otoferlin-dependent.
The L-type voltage-gated calcium channel whose opening triggers synaptic vesicle exocytosis in hair cells is Cav1.3 (Kollmar et al., 1997; Platzer et al., 2000; Brandt et al., 2003; Michna et al., 2003; Dou et al., 2004; Brandt et al., 2005; Baig et al., 2011). This contrasts with retinal photoreceptors, which do express Cav1.3, but rely on Cav1.4 for calcium influx (Matthews and Fuchs, 2010; Joiner and Lee, 2015). The pore-forming (alpha) subunit of Cav1.3 is encoded by Cacna1d, which is required for hearing (Platzer et al., 2000; Dou et al., 2004; Baig et al., 2011), and also for hair cell function in zebrafish (Nicolson et al., 1998; Sidi et al., 2004), where it is expressed in both neuromasts and the inner ear (Sidi et al., 2004). (The two zebrafish cacna1d genes show differential expression: cacna1da is expressed in hair cells plus retinal photoreceptors, while cacna1db is expressed in retinal and pineal photoreceptors, but not hair cells; Sidi et al., 2004.) Paddlefish Cacna1d is 19.3-fold lateral line-enriched (Supplementary file 1) and expressed in both ampullary organs and neuromasts (Figure 4C). Hence, Cav1.3 channels in the basal membrane are likely to mediate glutamate release from both paddlefish electroreceptors and hair cells. Given the homology of non-teleost ampullary organs (Bullock et al., 1983; Northcutt, 1986, Northcutt, 1992; Braun, 1996; New, 1997; Baker et al., 2013), this also suggests that the L-type voltage-gated calcium channels involved in neurotransmitter release from skate ampullary organs are likely to be Cav1.3 channels.
The abundance and function of Cav1.3 channels in inner ear hair cells is regulated by the auxiliary beta subunit Cavβ2, which is required for hearing (Neef et al., 2009). Cavβ2 is encoded by Cacnb2, which is the only other voltage-gated Ca2+ channel subunit gene in the lateral line-enriched dataset (2.7-fold enriched; Supplementary file 1). Although neuromast expression of cacnb2a and cacnb2b was not reported in zebrafish (Zhou et al., 2008), paddlefish Cacnb2 is expressed in both ampullary organs and neuromasts, with seemingly stronger expression in ampullary organs (Figure 4D). This suggests that Cavβ2 may be the auxiliary beta-subunit for Cav1.3 channels in electroreceptors, as well as hair cells.
Furthermore, given that the voltage sensors in skate electroreceptors are L-type voltage-gated calcium channels in the apical membrane (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005), recently demonstrated to be Cav1.3 channels (Bellono et al., 2017), and the fact that an apical calcium conductance is required for the firing of sturgeon ampullary organ afferents (Teeter et al., 1980), it seems likely that apically-located Cav1.3 channels, with Cavβ2 auxiliary beta subunits, act as the voltage-sensing channels, as well as mediating neurotransmitter release basally.
Finally, we report the expression in both ampullary organs and neuromasts of Rims2 (Rim2) and Ctbp2 (Ribeye), which encode proteins associated with synaptic ribbons in photoreceptors as well as hair cells (Matthews and Fuchs, 2010; Pangršič et al., 2012; Safieddine et al., 2012; Nicolson, 2015; Wichmann and Moser, 2015) (Figure 4E,F). Rims2 encodes Rab3-interacting molecules 2α and β, which are required in cochlear inner hair cells for Cav1.3 channel recruitment to the active zone membrane beneath the synaptic ribbon (Jung et al., 2015). Rims2 is 15.7-fold lateral line-enriched (Supplementary file 1) and expressed in all lateral line organs (Figure 4E). Hence, Rims2 may also be involved in Cav1.3 channel recruitment in electroreceptors. Ribeye, the only known synaptic ribbon-specific protein, is the main structural component of synaptic ribbons in both photoreceptors and hair cells, encoded by usage of an alternative start site for the transcription factor gene Ctbp2 that generates an N-terminal A-domain unique to Ribeye (Matthews and Fuchs, 2010; Nicolson, 2015; Wichmann and Moser, 2015). Ribeye is important in zebrafish neuromast hair cells for Cav1.3 channel recruitment to synaptic ribbons and stabilizing synaptic contacts with afferent neurons (Sheets et al., 2011; Lv et al., 2016). Recently, deletion in mice of the exon encoding the A-domain showed that Ribeye is essential in the retina (the ear was not examined) both for ribbon formation per se, and for rapid and sustained neurotransmitter release (Maxeiner et al., 2016). Ctbp2 is not in the lateral line-enriched dataset, but was present in the combined transcriptome, so was easily cloned. As expected, a riboprobe that exclusively recognizes the Ribeye-specific A-domain sequence of paddlefish Ctbp2 reveals expression in both ampullary organs and neuromasts (Figure 4F), suggesting that Ribeye is likely to be a key component of synaptic ribbons in electroreceptors (and hair cells), where it may also be important for Cav1.3 channel recruitment.
Taken together, these data suggest that the mechanisms of neurotransmission at the ribbon synapse in paddlefish electroreceptors are essentially identical to those at the hair cell ribbon synapse, involving otoferlin-dependent exocytosis of synaptic vesicles, loaded with glutamate by Vglut3, in response to the activation of basal Cav1.3 channels. Since non-teleost ampullary organs are homologous (Bullock et al., 1983; Northcutt, 1986, Northcutt, 1992; Braun, 1996; New, 1997; Baker et al., 2013), we predict that these mechanisms will be conserved across all other non-teleost ampullary organs.
Differential expression of beta-parvalbumin genes in ampullary organs and neuromasts
The paddlefish lateral line-enriched dataset contains two parvalbumin (Pvalb) genes: one, annotated as being related to zebrafish ‘pvalb8’, is enriched 22.8-fold; the other, annotated as being related to zebrafish ‘pvalb3’, is enriched 2.1-fold (Supplementary file 1). Parvalbumins are cytosolic EF-hand Ca2+-buffering proteins (Schwaller, 2010). Ca2+ is essential for multiple aspects of hair cell function (Lenzi and Roberts, 1994; Mammano et al., 2007; Ceriani and Mammano, 2012) and different hair cell subtypes are distinguished by different complements of EF-hand Ca2+ buffers, including different parvalbumin family members. Mammals have a single alpha-parvalbumin, encoded by Pvalb, and a single beta-parvalbumin, oncomodulin, encoded by Ocm (Schwaller, 2010). Alpha-parvalbumin is restricted to cochlear inner hair cells, while oncomodulin is restricted to cochlear outer hair cells, and is also expressed in vestibular hair cells (Sakaguchi et al., 1998; Yang et al., 2004; Simmons et al., 2010; Pangršič et al., 2015; Tong et al., 2016). Non-mammalian species have variable numbers of parvalbumin genes (nine in zebrafish; Friedberg, 2005), and the nomenclature for both genes and proteins is inconsistent and confusing. This led us to undertake a phylogenetic analysis of selected vertebrate parvalbumin proteins, including the two predicted paddlefish proteins (Figure 5A). This revealed three distinct clades: one comprising the alpha-parvalbumins (including zebrafish ‘Pvalb6’ and ‘Pvalb7’; Friedberg, 2005), and two containing beta-parvalbumins (Figure 5A). One of the beta-parvalbumin clades includes the oncomodulins (i.e., mammalian beta-parvalbumins), chicken Pvalb3/thymic CPV3 (encoded by Ocm) (Hapak et al., 1994), bullfrog Pvalb3 (Heller et al., 2002) and zebrafish ‘Pvalb8’ and ‘Pvalb9’, which were originally named Pvalb3a and Pvalb3b owing to their similarity to chicken Pvalb3/CPV3 and oncomodulin (Hsiao et al., 2002; Friedberg, 2005) (Figure 5A). The second beta-parvalbumin clade includes chicken thymic Pvalb (avian thymic hormone, encoded by Ocm2) (Brewer et al., 1989), zebrafish ‘Pvalb1-5’ (Friedberg, 2005) and various other beta-parvalbumins (Figure 5A).
Following our phylogenetic analysis, it was clear that both lateral line-enriched Pvalb genes encode beta-parvalbumins. We named the most highly lateral line-enriched gene Pvalbβ1/Ocm (22.8-fold enriched; ‘pvalb8’ in Supplementary file 1), because the predicted protein groups in the beta-parvalbumin clade containing mammalian oncomodulins (Figure 5A). Pvalbβ1/Ocm is expressed in both ampullary organs and neuromasts (Figure 5B). Interestingly, ‘Pvalb8’ was reported to be the most highly expressed transcript in skate ampullary organs: the authors suggest that, following the voltage-dependent Ca2+ influx via Cav1.3 that depolarizes the electroreceptor and activates BK, this parvalbumin could bind Ca2+, thus blocking BK-mediated hyperpolarization and enabling further oscillations (Bellono et al., 2017).
We named the less highly lateral line-enriched gene Pvalbβ2 (2.1-fold enriched; ‘pvalb3’ in Supplementary file 1), as it falls into the second beta-parvalbumin clade (Figure 5A). At all stages examined, Pvalbβ2 proved to be expressed in ampullary organs but not neuromasts (Figure 5C). Therefore, in addition to Neurod4, we have identified another transcript expressed by ampullary organs but not neuromasts. We suggest that this beta-parvalbumin is likely to be involved in ampullary organ-specific aspects of Ca2+ regulation.
Identification of ampullary organ-specific voltage-gated potassium channel subunits
The oscillatory character of the basal membrane voltage of skate electroreceptors depends on voltage-gated potassium channels, which contribute to repolarization of the basal membrane (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005). Noisy voltage oscillations have been recorded from ampullary organ canals in adult paddlefish in vivo (Neiman and Russell, 2004). Hence, we were very interested to find two highly lateral line-enriched shaker-related voltage-gated potassium channel subunit genes in our paddlefish dataset: Kcna5, encoding the pore-forming alpha subunit of Kv1.5 (12.4-fold enriched) and Kcnab3, encoding the beta subunit Kvβ3 (21.4-fold enriched). In contrast to the L-type voltage-gated calcium channel subunit genes Cacna1d and Cacnb2, which are expressed in both ampullary organs and neuromasts (Figure 4C,D), Kcna5 (Figure 6A–D) and Kcnab3 (Figure 6E–H) are only expressed in ampullary organs. Expression of both genes is first seen around stages 38–39, slightly after the eruption of the first ampullary organs at the surface at stage 37 (Modrell et al., 2011a).
Voltage-gated potassium channels are tetramers of four alpha (pore-forming) subunits, each containing six transmembrane segments, of which S4 contains a high density of positively charged residues and is the main transmembrane voltage-sensing component (Barros et al., 2012). Intriguingly, paddlefish Kv1.5 has a number of amino acid substitutions at otherwise highly conserved positions in S4, as well as in the S3-S4 linker, at the top of the S5 helix, and in the channel pore (Figure 6I,J). How the sum of these substitutions affects paddlefish Kv1.5 channel behavior must await studies of channel expression.
Discussion
Here, we took advantage of the abundance of ampullary organs (plus some neuromasts) on the operculum (gill-flap) of late-larval Mississippi paddlefish, and the absence of lateral line organs on the fins (Modrell et al., 2011a, 2011b), to generate a dataset of over 400 identified genes whose transcripts are enriched at least 1.8-fold in operculum versus fin tissue (i.e., lateral line-enriched). The dataset is not exhaustive: it does not include some genes whose expression we have validated in developing ampullary organs and neuromasts in paddlefish, whether via previous candidate gene approaches (e.g. Eya family members and Six family members other than Six1 and Six2; Modrell et al., 2011a), or some genes cloned in this study from transcripts present in the combined operculum plus fin transcriptome (e.g. Ctpb2, encoding the synaptic ribbon protein Ribeye). Nevertheless, this unbiased dataset provides an important foundation for investigating the molecular basis of ampullary organ development. Validation of a selection of genes from the dataset revealed significant molecular conservation between developing paddlefish ampullary organs and neuromasts, both in their expression of transcription factor genes critical for hair cell development, and also of genes essential for transmission specifically at the hair cell ribbon synapse. For the first time in any vertebrate, we also identify genes expressed in ampullary organs but not neuromasts, including a transcription factor, a beta-parvalbumin and voltage-gated potassium channel subunits consistent with predictions from skate ampullary organ electrophysiology.
Atoh1 and Neurod4 are likely to be critical for ampullary organ/electroreceptor development
Key elements of the transcription factor network underlying hair cell development, centering on the bHLH transcription factor Atoh1 (Cai and Groves, 2015; Jahan et al., 2015; Costa et al., 2017), seem to be conserved in developing paddlefish ampullary organs, including expression of Six1, Eya1 (Modrell et al., 2011a), Sox2, Atoh1 itself (this study and Butts et al., 2014) and Pou4f3 (Brn3c). Hence, Atoh1 and Pou4f3 are likely to be critical for the specification and differentiation of electroreceptors, as well as hair cells. This further highlights the importance of developmental context for Atoh1 activity, since it is also required for the specification of mechanosensory Merkel cells and proprioceptive neurons, cerebellar granule cells and intestinal secretory cells (Cai and Groves, 2015; Jahan et al., 2015; Costa et al., 2017). It also raises the question of which transcription factors are involved in the differentiation of ampullary organs/electroreceptors versus neuromasts/hair cells, whether acting downstream of Atoh1, or in parallel with it.
Of the 16 transcription factor genes whose expression we have reported to date in developing paddlefish lateral line organs, all except one are expressed in both ampullary organs and neuromasts: these are Six1, Six2, Six4, Eya1, Eya2, Eya3, Eya4 (Modrell et al., 2011a); the three SoxB1 class genes Sox1, Sox2 (this study) and Sox3 (Modrell et al., 2011b); Atoh1 (this study and Butts et al., 2014), Pou4f3, Pou4f1, Lhx3 and Myt1 (this study). The single exception is the proneural bHLH transcription factor gene Neurod4, a putative cerebellar Atoh1 target (Klisch et al., 2011) present in the lateral line-enriched dataset, which proved to be expressed in developing ampullary organs but not neuromasts. The related family member Neurod1 (which is expressed in paddlefish cranial sensory ganglia, but not lateral line organs; Modrell et al., 2011b) suppresses a hair cell fate in mouse otic neurons, and is important for the specification of outer versus inner hair cells in the cochlea (Jahan et al., 2010), while the even more closely related Neurod6 was identified as enriched in cochlear but not vestibular hair cells (Elkon et al., 2015). Hence, different NeuroD family members may be involved in specifying different hair cell subtypes. Taken together, we propose that Neurod4 is likely to be critical in paddlefish for the specification of ampullary organs/electroreceptors within lateral line placode-derived sensory ridges.
The function of these transcription factors during lateral line development could in principle be tested using the RNA-guided nuclease system CRISPR/Cas9, which has successfully been used to generate biallelic mutations efficiently in, for example, F0-injected axolotl (Fei et al., 2014; Flowers et al., 2014) and lamprey (Square et al., 2015). However, the very restricted annual spawning season for paddlefish, and the availability of 1-cell-stage embryos for injection only in a commercial fishery, rather than a laboratory setting, raise significant technical and logistical obstacles to optimizing CRISPR/Cas9 for this particular species. Another non-teleost chondrostean fish, the sturgeon Acipenser ruthenus (sterlet), has a much longer spawning season, with 1-cell-stage embryos readily available for microinjection in research facilities (e.g. Saito et al., 2014). Therefore, we plan to purse functional experiments on electroreceptor development using this species in the future.
Electroreceptor synaptic transmission mechanisms are conserved with hair cells
Although it was clear from electron microscopy that non-teleost electroreceptors in all species examined have ribbon-type synapses, given the presence of electron-dense presynaptic bodies (‘bars’, ‘sheets’ or ‘spheres’, depending on the species), surrounded by synaptic vesicles (Jørgensen, 2005), nothing was known at the molecular level about transmission mechanisms at electroreceptor synapses, except that activation of L-type voltage-gated calcium channels results in the release of a ‘glutamate-like’ neurotransmitter (Bennett and Obara, 1986). Our data suggest that the electroreceptor ribbon synapse is glutamatergic, and functions in the same way as the hair cell ribbon synapse (Safieddine et al., 2012; Nicolson, 2015; Wichmann and Moser, 2015; Moser and Starr, 2016), with glutamate being loaded into synaptic vesicles by Vglut3, and otoferlin-dependent exocytosis being triggered by the activation of Cav1.3 channels, potentially including the auxiliary beta subunit Cavβ2, as in cochlear inner hair cells (Neef et al., 2009). In contrast, retinal and pineal photoreceptors express Vglut1 and Vglut2 and neuronal SNAREs, and retinal photoreceptors depend on Cav1.4 (see e.g. Zanazzi and Matthews, 2009; Matthews and Fuchs, 2010; Mercer and Thoreson, 2011; Joiner and Lee, 2015). These physiological similarities between electroreceptors and hair cells are consistent with the expression in developing ampullary organs of key hair cell transcription factor genes, some of which presumably control the expression of ribbon synapse-associated genes in both electroreceptors and hair cells.
The expression of Cav and Kv channel subunit genes in paddlefish ampullary organs is consistent with predictions from skate ampullary organ electrophysiology
Our understanding of non-teleost ampullary organ physiology is primarily based on current- and voltage-clamp studies of epithelial currents in dissected single ampullary organ preparations from skates (cartilaginous fishes) (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005). These revealed that L-type voltage-gated calcium channels in the apical (lumenal, i.e., exterior-facing) electroreceptor membrane open in response to low-frequency cathodal electric fields (i.e., lumen-negative at the apical surface of the electroreceptor, relative to the interior of the animal). Ca2+ entry depolarizes the apical membrane and subsequently the basal membrane, which activates L-type voltage-gated calcium channels in the basal membrane. Ca2+ entry through these channels further depolarizes the basal membrane, resulting in synaptic vesicle exocytosis and neurotransmitter release. Basal Ca2+ entry also activates basal voltage-gated potassium channels and calcium-dependent chloride channels, thus repolarizing the basal membrane and de-activating the basal voltage-gated calcium channels (this gives the basal membrane voltage its oscillatory character). The apical membrane is repolarized after sufficient Ca2+ enters through the apical voltage-gated calcium channels to trigger the apical calcium-gated potassium channel BK. This terminates the depolarization and oscillations of the basal membrane. Importantly, the apical L-type voltage-gated calcium channels are partially activated in the absence of stimulus by a bias current across the sensory epithelium (provided by sodium/calcium exchangers and sodium/potassium pumps), so neurotransmitter is released steadily at rest, while the activity of the apical BK channel keeps the electroreceptor near threshold at rest and also during prolonged stimulus. Low-frequency cathodal stimuli accelerate the resting discharge of the afferent nerve, while anodal stimuli decelerate and eventually inhibit discharge (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005). Also in skate, Cav1.3 has just been identified as the L-type voltage-gated calcium channel that mediates the low-threshold voltage-dependent inward current, and works with BK to mediate electroreceptor membrane oscillations (Bellono et al., 2017).
Our data on voltage-gated ion channel expression in larval paddlefish ampullary organs are relevant both for the voltage-sensing calcium channels and for the voltage-gated potassium channel involved in rectifying the basal membrane during oscillations. An apical calcium conductance is required for ampullary organ afferent nerve firing in another non-teleost chondrostean fish, the sturgeon Scaphirhynchus platorynchus (Teeter et al., 1980). The expression in larval paddlefish ampullary organs of Cacna1d, encoding the pore-forming alpha subunit of Cav1.3 channels, and Cacnb2, encoding the auxiliary beta subunit Cavβ2, is consistent with the hypothesis that Cav1.3 channels are candidates for the apical voltage-sensing L-type voltage-gated calcium channels, as well as for the basal L-type voltage-gated calcium channels whose activation triggers otoferlin-dependent synaptic vesicle exocytosis and glutamate release at the ribbon synapse (see previous section). The identification of Cav1.3 channels as the voltage-sensor in skate electroreceptors (Bellono et al., 2017) further supports this hypothesis. Furthermore, given that Cacnb2 seems to be expressed at higher levels in ampullary organs than in neuromasts in paddlefish, it is possible that the voltage-sensing channels include the auxiliary beta subunit Cavβ2. (In contrast, deep-sequencing data from skate ampullary organs suggest that Cavβ1, encoded by Cacnb1, might be the auxiliary beta subunit in this species; Bellono et al., 2017.) Testing these hypotheses will require subcellular localization using paddlefish Cav channel subunit-specific antibodies and/or the optimization of CRISPR/Cas9 in conjunction with larval electroreceptor recordings, in this or related species.
We also identified expression in paddlefish ampullary organs, but not neuromasts, of Kcna5, encoding the pore-forming subunit of Kv1.5 channels, and Kcnab3, encoding an auxiliary beta-subunit. In the cardiac atrium, Kv1.5 channels are thought to conduct the ‘ultra-rapid delayed rectifier’ current (IKur), with relatively slow inactivation properties, which contributes to repolarization (Schmitt et al., 2014; Wettwer and Terlau, 2014). Furthermore, when human KCNA5 and KCNAB3 are heterologously co-expressed in cell culture, they form a novel Kv channel that, upon depolarization, mediates very fast-inactivating (A-type) outward currents (Leicher et al., 1998). Hence, it is possible that Kv1.5 and Kvβ3 in paddlefish ampullary organs together constitute a rapidly-inactivating voltage-gated potassium channel that conducts an ultra-rapid delayed rectifier current. These properties would plausibly be optimal for rapid oscillation of the basal electroreceptor membrane. Noisy voltage oscillations have been recorded from ampullary organ canals in adult paddlefish in vivo (Neiman and Russell, 2004). Importantly, the only evidence until now for the existence of a rectifying Kv channel in non-teleost ampullary organs has been from current- and voltage-clamp approaches in dissected skate (cartilaginous fish) ampullary organs (Bennett and Obara, 1986; Lu and Fishman, 1995; Bodznick and Montgomery, 2005). Our gene expression data suggest that this channel also exists in a chondrostean fish, supporting the hypothesis that non-teleost electroreceptor physiology is conserved (at least across jawed vertebrates), and identify Kv1.5 and Kvβ3 as the candidate channel.
Amino acid substitutions may alter Kv1.5 channel properties in paddlefish
The sequence of paddlefish Kcna5 predicts several amino acid substitutions in highly conserved regions of Kv1.5 that could affect its properties. The S4 segment of voltage-gated ion channels is the voltage sensor, and is highly conserved (Barros et al., 2012). The first positively-charged amino acid in the voltage sensor is normally arginine: since it is the first arginine in S4, it is referred to as the R1 position. In paddlefish Kv1.5, this has been replaced by a threonine (T357). Experimental substitution of methionine at this site in Drosophila Shaker (R362M; Aggarwal and MacKinnon, 1996) decreases the gating charge on the channel and the slope of the conductance voltage curve (Aggarwal and MacKinnon, 1996; Elliott et al., 2012). Aside from its role as part of the voltage sensor, R1 forms a seal with hydrophobic residues in S2 that keeps ions from flowing around the S4 segment: substitution of R1 in Drosophila Shaker with histidine (R362H) allows protons to leak past S4 (Starace and Bezanilla, 2004). Substitution with some other amino acids (alanine, cysteine, serine, valine) allows monovalent cations to leak past S4 when the cell is hyperpolarized (omega current) (Tombola et al., 2005). It is not known if a threonine at this position would encourage an omega current, but an omega current could act like a slow depolarizing leak current, possibly contributing to the spontaneous activity of ampullary electroreceptors.
The other amino acid substitutions in paddlefish Kv1.5 are at conserved positions in the S3-S4 linker, at the top of the S5 helix, and in the channel pore. A variety of studies indicate that the end of the S3/S4 linker and the top of S4 contact the top of S2, S3 and S5/pore in the closed or open state (Kanevsky and Aldrich, 1999; Elliott et al., 2004; Soler-Llavina et al., 2006; Henrion et al., 2009; Lin et al., 2011; Elliott et al., 2012). Generally speaking, amino acid substitutions at some of these amino acid positions cause a rightward shift in the conductance-voltage curve, meaning that the channel would open at more depolarized voltages. Characterization of the effects on Kv1.5 of these various amino acid substitutions will require future paddlefish channel expression in a heterologous system and mutagenesis studies: for example, characterizing the properties of paddlefish Kv1.5, as compared with paddlefish Kv1.5 in which these amino acids (in different combinations) have been substituted with conserved residues, and/or substituting the paddlefish-specific amino acids into, for example, human Kv1.5.
Insights into electroreceptor evolution
As noted above, of the 16 transcription factor genes whose expression we have reported to date in developing paddlefish lateral line organs (this study; Modrell et al., 2011a, 2011b; Butts et al., 2014), all except one (Neurod4) are expressed in ampullary organs as well as neuromasts, including transcription factor genes that are essential for hair cell development, such as Atoh1 and Pou4f3 (Brn3c). Furthermore, genes that are required for synaptic transmission at the hair cell (but not photoreceptor) ribbon synapse are also expressed in larval paddlefish ampullary organs. This suggests very close developmental and, most likely, evolutionary links between hair cells and electroreceptors. This was not necessarily to be expected: electroreceptors could have been more similar to other cell types with ribbon synapses, i.e., retinal or pineal photoreceptors, or retinal bipolar cells (indeed, Neurod4 is also important for bipolar neuron development; Hatakeyama and Kageyama, 2004).
The selective pressure for the evolution of electroreceptors, which enable the detection of living animals in water, is likely related to the transition from filter-feeding to predation in the lineage leading to vertebrates (Gans and Northcutt, 1983; Northcutt and Gans, 1983; Northcutt, 2005b). Various lines of evidence support the homology of all non-teleost electroreceptors, including those of lampreys: their stimulation by weak, low-frequency cathodal fields (and inhibition by strong cathodal fields); their innervation by pre-otic anterior lateral afferents projecting to the dorsal octavolateral nucleus in the hindbrain; and, in jawed vertebrates, their demonstrated embryonic origin, together with neuromasts, from lateral line placodes (Bullock et al., 1983; Northcutt, 1992; Braun, 1996; New, 1997; Baker et al., 2013). Given this, the striking similarities in both development and physiology that we have identified in paddlefish ampullary organs are consistent with electroreceptor evolution in the vertebrate ancestor either via the diversification of hair cells that had already evolved from ancestral mechanoreceptor cells (Duncan and Fritzsch, 2012; Fritzsch and Straka, 2014), or via the diversification of an ancestral mechanoreceptor cell - with an apical primary cilium (lost in lamprey electroreceptors), microvilli, and synaptic ribbons (Northcutt, 1986; Bodznick, 1989; New, 1997), and using otoferlin, Vglut3 and Cav1.3 channels for synaptic transmission - into both hair cells and electroreceptors.
Summary and perspective
Our unbiased RNA-seq approach has shed new light on molecular mechanisms underlying ampullary organ development and physiology in a non-teleost chondrostean fish. The expression in ampullary organs, as well as neuromasts, of key hair cell transcription factor genes, such as Atoh1, Sox2 and Pou4f3, and genes encoding proteins specifically required for glutamate release at the hair cell ribbon synapse, such as otoferlin, Vglut3 and Cav1.3 channels, supports close developmental and physiological, hence evolutionary, relationships between electroreceptors and hair cells. Our identification of the first-reported ampullary organ-specific genes, including the proneural transcription factor gene Neurod4, a beta-parvalbumin gene, and voltage-gated potassium channel subunit genes, provides novel insight into the potential molecular basis of ampullary organ-specific development and physiology. We identify Cav1.3 (with Cavβ2) as a candidate for the apical voltage-sensing channel (indeed, Cav1.3 was recently shown to be the voltage-sensing channel in skate electroreceptors; Bellono et al., 2017), and Kv1.5 (with unusual amino acid substitutions that may affect its properties), together with Kvβ3, as a candidate for the basal membrane rectifying Kv channel predicted from skate ampullary organ electrophysiology. As noted above, the paddlefish has a very limited annual spawning period, so it is technically difficult to optimize e.g. CRISPR/Cas9 for targeted mutation. However, the continuing rapid advances in CRISPR/Cas9 technology, plus the development of other experimentally tractable non-teleost model systems, should eventually make it possible to test specific gene function in this and related species. Overall, our analysis has provided wide-ranging insights into molecular aspects of ampullary organ development and physiology, and an essential framework for future comparative work to determine both the function of these genes, and the extent to which they are conserved across other non-teleost electroreceptive vertebrate groups.
Materials and methods
Tissue and RNA isolation and Illumina sequencing
Request a detailed protocolP. spathula embryos were purchased over multiple spawning seasons from Osage Catfisheries Inc. (Osage Beach, MO, USA) and staged according to Bemis and Grande (1992). All experiments were performed in accordance with the approved institutional guidelines and regulations of the Institutional Animal Care and Use Committee of Kennesaw State University (approved protocol #12–001). Stage 46 yolk-sac larvae were preserved in RNALater (Ambion, Thermo Fisher Scientific Inc., Waltham, MA, USA) overnight at 4°C. Excess solution was removed and samples were stored at −80°C until processed. Opercular (lateral line organ-enriched) and fin (no lateral line organs) tissues were manually dissected and pooled from three different sets of 6–7 specimens each, yielding three biological replicates. RNA was extracted using Trizol reagent (Ambion), according to the manufacturer’s protocol. RNA concentration was assessed using a Nanodrop N1000 spectrophotometer and integrity using an Agilent 2100 Bioanalyzer (Cambridge Genomic Services). Only samples with an RNA integrity number (RIN) greater than nine were used for next-generation sequencing. Illumina RNA-sequencing library preparation and sequencing were performed by The Centre for Applied Genomics, The Hospital for Sick Children, Toronto, Canada. Libraries were prepared following the standard Illumina RNA Library Prep kit and sequenced on an Illumina HiSeq 2500, using Illumina v3 chemistry, following the multiplex paired-end protocol (2 × 100 bases).
Assembly and analysis of transcriptome
Read QC and trimming
Request a detailed protocolReads were subjected to various quality controls, including filtering of high-quality reads based on the score value given in fastq files (FastQC version 0.10.1; http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), removal of reads containing primer/adaptor sequences and trimming of read length using Trimmomatic-0.30 (Bolger et al., 2014).
De novo assembly of the transcriptome
Request a detailed protocolReads were de novo assembled using Velvet version 1.2.10 (Zerbino and Birney, 2008) and Oases version 0.2.08 (Schulz et al., 2012). Velvet was run using different k-mer lengths k27–k77 in k10 increments along with other default parameters. Oases was run using the same k-mer range. Results from these assemblies were merged, again using Velvet and Oases k-mer of k27. All assemblies were performed on a server with 64 cores and 512 Gb of RAM.
Obtaining transcript counts
Request a detailed protocolReads were mapped back to the transcriptome using Bowtie2 version 2–2.1.0 (Langmead and Salzberg, 2012). As reads were uncorrected, some cleaning of the SAM files was needed to remove PCR duplicates and to remove long-stretch (>85%) poly-A sequences. Transcript counts for each sample were obtained with HTseq-count (version 0.5.4p3) (Anders et al., 2015). A locus-to-transcript mapping file was used to collapse related transcripts and obtain locus-level counts. Output was used as input for statistical calculations.
Differential Expression Analysis
Request a detailed protocolThe BioConductor package DESeq (Anders and Huber, 2010) was used for differential expression analysis. A p-value of <0.1 was considered significant after adjustment (multiple testing using Benjamini-Hochberg).
BLAST Annotation
Request a detailed protocolWith a dataset of 189,933 contigs (>200 bp) assembled by Velvet/ Oases, transcripts were searched against chordate and invertebrate protein datasets (obtained from SwissProt and NCBI databases) using NCBI’s Basic Local Alignment Search Tool BLASTX (McGinnis and Madden, 2004), with an expected (E)-value cut-off of ≤1E−05 to reveal sequence conservation. For results above threshold, the UniProt protein record was obtained for the top BLAST hit against each transcript locus.
Functional Annotation and Enrichment Analysis
Request a detailed protocolGene ontology (GO) and protein domain annotations were extracted from relevant UniProt records. Enrichment analysis was performed using hypergeometric mean and Bonferonni multiple testing correction. Enriched genes were classified using GO terms according to molecular function using the web-based resource PANTHER (Mi et al., 2016).
RNA-seq data have been deposited in the NCBI Gene Expression Omnibus (GEO) database under accession code GSE92470.
Phylogenetic analysis of parvalbumins
Request a detailed protocolAmino acid sequences were downloaded from GenBank and aligned using the online version of MUSCLE (Edgar, 2004) from the EMBL-EBI server (Li et al., 2015). The following Parvalbumin (Pvalb) protein sequences were used: for Callorhinchus milii, Pvalb, thymic CPV3-like (AFP11760), Pvalb-alpha-like (AFP11872); for Danio rerio, Pvalb1 (AAH71552), Pvalb2 (AAH93135), Pvalb3 (AAH46001), Pvalb4 (AAH72551), Pvalb5 (AAH92666), Pvalb6 (NP_991136), Pvalb7 (NP_991137), Pvalb8 (NP_891982), Pvalb9 (NP_891983); for Gallus gallus, Pvalb3/thymic CPV3 (AAA17518), Pvalb-alpha (CAX32963), Pvalb thymic/Ocm2 (NP_001007478); for Homo sapiens, Pvalb-alpha (NP_001302461), Ocm (AAH69468); for Mus musculus, Pvalb-alpha (NP_038673), Ocm (NP_149028); for Rana catesbeiana, Pvalb3 (AAL09922), Pvalb-alpha (BAC55948), Pvalb-beta (AC051685); for Rattus norvegicus, Pvalb-alpha (NP_071944), Ocm (P02631); for Salmo salar, Pvalb-alpha (NP_001167235), Pvalb thymic/CPV3 (ACM09534), Pvalb-beta1 (NP_001117190), Pvalb-beta2 (NP_001117189).
Phylogenetic analyses were carried out using a Bayesian framework with the parallel version of MrBayes 3.2.6 (Huelsenbeck and Ronquist, 2001; Ronquist and Huelsenbeck, 2003). One major issue with analyzing highly divergent multi-gene families is finding suitable outgroups to root the tree appropriately. Instead of using an outgroup, a relaxed molecular clock with independent gamma rates was used to infer the position of the root. A ‘mixed’ substitution model prior for the amino acid sequences was used to allow the program to explore and sample across substitution models. To estimate posterior probabilities of all parameters, two Metropolis-coupled Markov chain Monte Carlo (MCMCMC) runs of 100 million generations were performed, sampling every 10,000 generations and discarding the first 25% as burn-in. The resulting annotated consensus tree was subsequently edited in FigTree (http://tree.bio.ed.ac.uk/software/figtree/).
Alignments of Shaker-related voltage-gated potassium channel amino acid sequences
Request a detailed protocolAmino acid sequences for Shaker-related voltage-gated potassium channels were aligned using Seaview (Gouy et al., 2010) with default settings. The following protein sequences were used: Aplysia californica Shaker-like (NP_001191634.1), Callorhinchus milii Kv1.5 (XP_007902374.1); Danio rerio Kv1.5 (XP_005171764.1); Drosophila melanogaster Shaker (CAA29917.1); Gallus gallus Kv1.5 (XP_015147812.1); Homo sapiens Kv1.5 (NP_002225.2), Kv1.1 (NP_000208.2) and Kv1.4 (NP_002224.1); Latimeria chalumnae Kv1.5 (XP_006013184.1); Lepisosteus oculatus Kv1.5 (XP_006633553.2); Nematostella vectensis Shaker4 (AFY09706.1); and Xenopus tropicalis Kv1.5 (XP_004912856.1).
Gene cloning, in situ hybridization and immunohistochemistry
Request a detailed protocolTotal RNA from stage 44–46 embryos or stage 46 opercular tissue was isolated using Trizol (Invitrogen, Carlsbad, CA), as per the manufacturer’s protocol. cDNA was made using the Superscript III First Strand Synthesis kit (Invitrogen, Thermo Fisher Scientific). Gene-specific primers were used under standard PCR conditions to amplify gene fragments, which were cloned into the pDrive vector (Qiagen, Manchester, UK) and individual clones verified by sequencing (Department of Biochemistry Sequencing Facility, University of Cambridge, UK). GenBank accession numbers for cloned paddlefish cDNA fragments are as follows: Cacna1d KY781950, Cacnb2 KY781951, Ctbp2 KY781952, Kcna5 KY781953, Kcnab3 KY781954, Lhx3 KY781955, Myt1 KY781956, Neurod4 KY781957, Otof KY781958, Pou4f1 KY781959, Pou4f3 KY781960, Pvalbβ1/Ocm KY781961, Pvalbβ2 KY781962, Rims2 KY781963, Slc17a8 KY781964, Sox1 KY781965 and Sox2 KY781966. Anti-sense RNA probes were synthesized using T7 or SP6 polymerases (Promega, Southampton, UK) and digoxigenin-labeled dUTPs (Roche, Basel, Switzerland).
Whole-mount in situ hybridization and immunohistochemistry was performed as described (Modrell et al., 2011a). The primary antibody against Sox2 was ab92494 (rabbit, 1:200–1:400; Abcam, Cambridge, UK; reported to work in skate: http://www.abcam.com/sox2-antibody-epr3131-ab92494.html/reviews/37352). A horseradish peroxidase-conjugated goat anti-rabbit secondary (Jackson ImmunoResearch Laboratories, Inc., West Grove, PA, USA) was used at 1:600. Each RNA probe/antibody was tested a minimum of two times, using at least five embryos per experimental trial.
Data availability
-
Data from: Insights into electrosensory organ development, physiology and evolution from a lateral line-enriched transcriptomePublicly available at the NCBI Gene Expression Omnibus (accession no: GSE92470).
References
-
The evolution and development of vertebrate lateral line electroreceptorsJournal of Experimental Biology 216:2515–2522.https://doi.org/10.1242/jeb.082362
-
Cytoplasmic domains and voltage-dependent potassium channel gatingFrontiers in Pharmacology 3:49.https://doi.org/10.3389/fphar.2012.00049
-
Spatial and temporal segregation of auditory and vestibular neurons in the otic placodeDevelopmental Biology 322:109–120.https://doi.org/10.1016/j.ydbio.2008.07.011
-
BookIonic mechanisms and pharmacology of electroreceptorsIn: Bullock TH, Heiligenberg W, editors. Electroreception. New York: Wiley. pp. 157–181.
-
BookComparisons between electrosensory and mechanosensory lateral line systemsIn: Coombs S, Görner P, Münz H, editors. The Mechanosensory Lateral Line. Neurobiology and Evolution. New York: Springer-Verlag. pp. 655–678.
-
BookThe physiology of low-frequency electrosensory systemsIn: Bullock TH, Hopkins CD, Popper AN, Fay RR, editors. Electroreception. New York: Springer. pp. 132–153.
-
Trimmomatic: a flexible trimmer for Illumina sequence dataBioinformatics 30:2114–2120.https://doi.org/10.1093/bioinformatics/btu170
-
Cav1.3 channels are essential for development and presynaptic activity of cochlear inner hair cellsJournal of Neuroscience 26:10832–10840.
-
Few Cav 1.3 channels regulate the exocytosis of a synaptic vesicle at the hair cell ribbon synapseJournal of Neuroscience 25:11577–11585.https://doi.org/10.1523/JNEUROSCI.3411-05.2005
-
The sensory biology of the living jawless fishes: a phylogenetic assessmentBrain, Behavior and Evolution 48:262–276.https://doi.org/10.1159/000113205
-
Avian thymic hormone (ATH) is a parvalbuminBiochemical and Biophysical Research Communications 160:1155–1161.https://doi.org/10.1016/S0006-291X(89)80124-X
-
The role of atonal factors in mechanosensory cell specification and functionMolecular Neurobiology 52:1315–1329.https://doi.org/10.1007/s12035-014-8925-0
-
Calcium signaling in the cochlea - Molecular mechanisms and physiopathological implicationsCell Communication and Signaling 10:20.https://doi.org/10.1186/1478-811X-10-20
-
Building the posterior lateral line system in zebrafishDevelopmental Neurobiology 72:234–255.https://doi.org/10.1002/dneu.20962
-
Calcium-activated conductance in skate electroreceptors: current clamp experimentsThe Journal of General Physiology 69:121–143.https://doi.org/10.1085/jgp.69.2.121
-
Calcium-activated conductance in skate electroreceptors: voltage clamp experimentsThe Journal of General Physiology 69:145–182.https://doi.org/10.1085/jgp.69.2.145
-
Atoh1 in sensory hair cell development: constraints and cofactorsSeminars in Cell and Developmental Biology.https://doi.org/10.1016/j.semcdb.2016.10.003
-
Null mutation of alpha1D Ca2+ channel gene results in deafness but no vestibular defect in miceJournal of the Association for Research in Otolaryngology 5:215–226.https://doi.org/10.1007/s10162-003-4020-3
-
Evolution of sound and balance perception: innovations that aggregate single hair cells into the ear and transform a gravistatic sensor into the organ of CortiThe Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology 295:1760–1774.https://doi.org/10.1002/ar.22573
-
RFX transcription factors are essential for hearing in miceNature Communications 6:8549.https://doi.org/10.1038/ncomms9549
-
Molecular mechanism of voltage sensor movements in a potassium channelThe EMBO Journal 23:4717–4726.https://doi.org/10.1038/sj.emboj.7600484
-
Role of hydrophobic and ionic forces in the movement of S4 of the Shaker potassium channelMolecular Membrane Biology 29:321–332.https://doi.org/10.3109/09687688.2012.710343
-
Developmental origin of shark electrosensory organsEvolution and Development 8:74–80.https://doi.org/10.1111/j.1525-142X.2006.05076.x
-
Parvalbumin isoforms in zebrafishMolecular Biology Reports 32:167–175.https://doi.org/10.1007/s11033-005-2334-4
-
SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree buildingMolecular Biology and Evolution 27:221–224.https://doi.org/10.1093/molbev/msp259
-
Novel avian thymic parvalbumin displays high degree of sequence homology to oncomodulinThe Journal of Biological Chemistry 269:5288–5296.
-
Retinal cell fate determination and bHLH factorsSeminars in Cell and Developmental Biology 15:83–89.https://doi.org/10.1016/j.semcdb.2003.09.005
-
Parvalbumin 3 is an abundant Ca2+ buffer in hair cellsJARO - Journal of the Association for Research in Otolaryngology 3:488–498.https://doi.org/10.1007/s10162-002-2050-x
-
Long QT syndrome-associated mutations in the voltage sensor of IKs channelsCellular Physiology and Biochemistry 24:11–16.https://doi.org/10.1159/000227828
-
Lhx3, a LIM domain transcription factor, is regulated by Pou4f3 in the auditory but not in the vestibular systemEuropean Journal of Neuroscience 25:999–1005.https://doi.org/10.1111/j.1460-9568.2007.05332.x
-
Isolation and expression of two zebrafish homologues of parvalbumin genes related to chicken CPV3 and mammalian oncomodulinMechanisms of Development 119 Suppl 1:S161–S166.https://doi.org/10.1016/S0925-4773(03)00110
-
MRBAYES: Bayesian inference of phylogenetic treesBioinformatics 17:754–755.https://doi.org/10.1093/bioinformatics/17.8.754
-
Voltage-gated Cav1 channels in disorders of vision and hearingCurrent Molecular Pharmacology 8:143–148.https://doi.org/10.2174/1874467208666150507104937
-
BookMorphology of electroreceptive sensory organsIn: Bullock TH, Hopkins CD, Popper AN, Fay RR, editors. Electroreception. New York: Springer. pp. 47–67.
-
The Lorenzinian ampullae of Polyodon spathulaZeitschrift für Zellforschung und Mikroskopische Anatomie 130:362–377.https://doi.org/10.1007/BF00306949
-
Determinants of voltage-dependent gating and open-state stability in the S5 segment of Shaker potassium channelsThe Journal of General Physiology 114:215–242.https://doi.org/10.1085/jgp.114.2.215
-
Evolution of time-coding systems in weakly electric fishesZoological Science 26:587–599.https://doi.org/10.2108/zsj.26.587
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Brn3a target gene recognition in embryonic sensory neuronsDevelopmental Biology 302:703–716.https://doi.org/10.1016/j.ydbio.2006.10.050
-
Coexpression of the KCNA3B gene product with Kv1.5 leads to a novel A-type potassium channelJournal of Biological Chemistry 273:35095–35101.https://doi.org/10.1074/jbc.273.52.35095
-
Calcium signalling in hair cells: multiple roles in a compact cellCurrent Opinion in Neurobiology 4:496–502.https://doi.org/10.1016/0959-4388(94)90049-3
-
The EMBL-EBI bioinformatics web and programmatic tools frameworkNucleic Acids Research 43:W580–W584.https://doi.org/10.1093/nar/gkv279
-
R1 in the Shaker S4 occupies the gating charge transfer center in the resting stateThe Journal of General Physiology 138:155–163.https://doi.org/10.1085/jgp.201110642
-
The diverse roles of ribbon synapses in sensory neurotransmissionNature Reviews Neuroscience 11:812–822.https://doi.org/10.1038/nrn2924
-
BLAST: at the core of a powerful and diverse set of sequence analysis toolsNucleic Acids Research 32:W20–W25.https://doi.org/10.1093/nar/gkh435
-
Homeobox genes in axolotl lateral line placodes and neuromastsDevelopment Genes and Evolution 207:287–295.https://doi.org/10.1007/s004270050116
-
PANTHER version 10: expanded protein families and functions, and analysis toolsNucleic Acids Research 44:D336–D342.https://doi.org/10.1093/nar/gkv1194
-
Cav1.3 (alpha1D) Ca2+ currents in neonatal outer hair cells of miceThe Journal of Physiology 553:747–758.https://doi.org/10.1113/jphysiol.2003.053256
-
Auditory neuropathy--neural and synaptic mechanismsNature Reviews Neurology 12:135–149.https://doi.org/10.1038/nrneurol.2016.10
-
Two distinct types of noisy oscillators in electroreceptors of paddlefishJournal of Neurophysiology 92:492–509.https://doi.org/10.1152/jn.00742.2003
-
The evolution of vertebrate electrosensory systemsBrain, Behavior and Evolution 50:244–252.https://doi.org/10.1159/000113338
-
Ribbon synapses in zebrafish hair cellsHearing Research 330:170–177.https://doi.org/10.1016/j.heares.2015.04.003
-
BookElectroreception in non-teleost bony fishesIn: Bullock TH, Heiligenberg W, editors. Electroreception. New York: Wiley. pp. 257–285.
-
BookThe phylogeny of octavolateralis ontogenies: a reaffirmation of Garstang’s phylogenetic hypothesisIn: Webster DB, Popper AN, Fay RR, editors. The Evolutionary Biology of Hearing. New York: Springer-Verlag. pp. 21–47.
-
Electroreception112–131, Ontogeny of electroreceptors and their neural circuitry, Electroreception, New York, Springer.
-
The New Head Hypothesis revisitedJournal of Experimental Zoology Part B: Molecular and Developmental Evolution 304B:274–297.https://doi.org/10.1002/jez.b.21063
-
The genesis of neural crest and epidermal placodes: a reinterpretation of vertebrate originsThe Quarterly Review of Biology 58:1–28.https://doi.org/10.1086/413055
-
Electroreceptors and mechanosensory lateral line organs arise from single placodes in axolotlsDevelopmental Biology 168:358–373.https://doi.org/10.1006/dbio.1995.1086
-
Vesicular glutamate transporter 3 is required for synaptic transmission in zebrafish hair cellsJournal of Neuroscience 28:2110–2118.https://doi.org/10.1523/JNEUROSCI.5230-07.2008
-
Otoferlin: a multi-C2 domain protein essential for hearingTrends in Neurosciences 35:671–680.https://doi.org/10.1016/j.tins.2012.08.002
-
The development of lateral line placodes: taking a broader viewDevelopmental Biology 389:68–81.https://doi.org/10.1016/j.ydbio.2014.02.016
-
MrBayes 3: Bayesian phylogenetic inference under mixed modelsBioinformatics 19:1572–1574.https://doi.org/10.1093/bioinformatics/btg180
-
The auditory hair cell ribbon synapse: from assembly to functionAnnual Review of Neuroscience 35:509–528.https://doi.org/10.1146/annurev-neuro-061010-113705
-
Oncomodulin is expressed exclusively by outer hair cells in the organ of CortiJournal of Histochemistry and Cytochemistry 46:29–39.https://doi.org/10.1177/002215549804600105
-
Proneural gene requirement for hair cell differentiation in the zebrafish lateral lineDevelopmental Biology 295:534–545.https://doi.org/10.1016/j.ydbio.2006.03.037
-
Cardiac potassium channel subtypes: new roles in repolarization and arrhythmiaPhysiological Reviews 94:609–653.https://doi.org/10.1152/physrev.00022.2013
-
Cytosolic Ca2+ buffersCold Spring Harbor Perspectives in Biology 2:a004051.https://doi.org/10.1101/cshperspect.a004051
-
Oncomodulin identifies different hair cell types in the mammalian inner earThe Journal of Comparative Neurology 518:3785–3802.https://doi.org/10.1002/cne.22424
-
Ampullary electroreceptors in the sturgeon Scaphirhynchus platorynchus (Rafinesque)Journal of Comparative Physiology 138:213–223.https://doi.org/10.1007/BF00657039
-
There and back again: development and regeneration of the zebrafish lateral line systemWiley Interdisciplinary Reviews: Developmental Biology 4:1–16.https://doi.org/10.1002/wdev.160
-
Oncomodulin, an EF-Hand Ca2+ buffer, is critical for maintaining cochlear function in miceJournal of Neuroscience 36:1631–1635.https://doi.org/10.1523/JNEUROSCI.3311-15.2016
-
Pharmacology of voltage-gated potassium channel Kv1.5 - impact on cardiac excitabilityCurrent Opinion in Pharmacology 15:115–121.https://doi.org/10.1016/j.coph.2014.02.001
-
Merkel cells in lower vertebratesArchives of Histology and Cytology 52:415–422.https://doi.org/10.1679/aohc.52.Suppl_415
-
Relating structure and function of inner hair cell ribbon synapsesCell and Tissue Research 361:95–114.https://doi.org/10.1007/s00441-014-2102-7
-
The molecular architecture of ribbon presynaptic terminalsMolecular Neurobiology 39:130–148.https://doi.org/10.1007/s12035-009-8058-z
-
Identification and expression of voltage-gated calcium channel beta subunits in zebrafishDevelopmental Dynamics 237:3842–3852.https://doi.org/10.1002/dvdy.21776
Article and author information
Author details
Funding
Biotechnology and Biological Sciences Research Council (BB/F00818X/1)
- Clare VH Baker
Leverhulme Trust (RPG-383)
- Clare VH Baker
Fisheries Society of the British Isles (Research Grant)
- Melinda S Modrell
National Science Foundation (IOS 1557857)
- Harold H Zakon
National Science Foundation (IOS 1144965)
- Marcus C Davis
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the BBSRC (BB/F00818X/1 to CVHB), the Leverhulme Trust (RPG-383 to CVHB), the Fisheries Society of the British Isles (Research Grant to MSM) and the NSF (IOS 1557857 to HHZ; IOS 1144965 to MCD). Thanks to Rachel Lyne (Cambridge Systems Biology Centre, University of Cambridge) for submitting RNA-seq data to NCBI GEO. Thanks to Peterhouse and the Department of Physiology, Development and Neuroscience at the University of Cambridge for hosting HHZ. We also thank Tatjana Piotrowski and her lab at the Stowers Institute for Medical Research (Kansas City, MO, USA) and Steve and Pete Kahrs and the Kahrs family (Osage Catfisheries, Inc.) for hosting MSM during paddlefish spawning seasons.
Ethics
Animal experimentation: All experiments were performed in accordance with the approved institutional guidelines and regulations of the Institutional Animal Care and Use Committee of Kennesaw State University (approved protocol #12-001).
Copyright
© 2017, Modrell et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,531
- views
-
- 576
- downloads
-
- 40
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Developmental Biology
Eukaryotic cells depend on exocytosis to direct intracellularly synthesized material toward the extracellular space or the plasma membrane, so exocytosis constitutes a basic function for cellular homeostasis and communication between cells. The secretory pathway includes biogenesis of secretory granules (SGs), their maturation and fusion with the plasma membrane (exocytosis), resulting in release of SG content to the extracellular space. The larval salivary gland of Drosophila melanogaster is an excellent model for studying exocytosis. This gland synthesizes mucins that are packaged in SGs that sprout from the trans-Golgi network and then undergo a maturation process that involves homotypic fusion, condensation, and acidification. Finally, mature SGs are directed to the apical domain of the plasma membrane with which they fuse, releasing their content into the gland lumen. The exocyst is a hetero-octameric complex that participates in tethering of vesicles to the plasma membrane during constitutive exocytosis. By precise temperature-dependent gradual activation of the Gal4-UAS expression system, we have induced different levels of silencing of exocyst complex subunits, and identified three temporarily distinctive steps of the regulated exocytic pathway where the exocyst is critically required: SG biogenesis, SG maturation, and SG exocytosis. Our results shed light on previously unidentified functions of the exocyst along the exocytic pathway. We propose that the exocyst acts as a general tethering factor in various steps of this cellular process.
-
- Cell Biology
- Developmental Biology
In most murine species, spermatozoa exhibit a falciform apical hook at the head end. The function of the sperm hook is not yet clearly understood. In this study, we investigate the role of the sperm hook in the migration of spermatozoa through the female reproductive tract in Mus musculus (C57BL/6), using a deep tissue imaging custom-built two-photon microscope. Through live reproductive tract imaging, we found evidence indicating that the sperm hook aids in the attachment of spermatozoa to the epithelium and facilitates interactions between spermatozoa and the epithelium during migration in the uterus and oviduct. We also observed synchronized sperm beating, which resulted from the spontaneous unidirectional rearrangement of spermatozoa in the uterus. Based on live imaging of spermatozoa-epithelium interaction dynamics, we propose that the sperm hook plays a crucial role in successful migration through the female reproductive tract by providing anchor-like mechanical support and facilitating interactions between spermatozoa and the female reproductive tract in the house mouse.