Abstract
Pathogenic variants in HCN1 are associated with a range of epilepsy syndromes including a developmental and epileptic encephalopathy. The recurrent de novo HCN1 pathogenic variant (M305L) results in a cation leak, allowing the flux of excitatory ions at potentials where the wild-type channels are closed. The Hcn1M294L mouse recapitulates patient seizure and behavioral phenotypes. As HCN1 channels are highly expressed in rod and cone photoreceptor inner segments, where they shape the light response, mutated channels are likely to impact visual function. Electroretinogram (ERG) recordings from male and female mice Hcn1M294L mice revealed a significant decrease in the photoreceptor sensitivity to light, as well as attenuated bipolar cell (P2) and retinal ganglion cell responses. Hcn1M294L mice also showed attenuated ERG responses to flickering lights. ERG abnormalities are consistent with the response recorded from a single female human subject. There was no impact of the variant on the structure or expression of the Hcn1 protein in the retina. In silico modeling of photoreceptors revealed that the mutated HCN1 channel dramatically reduced light-induced hyperpolarization, resulting in more Ca2+ flux during the response when compared with the wild-type situation. We propose that the light-induced change in glutamate release from photoreceptors during a stimulus will be diminished, significantly blunting the dynamic range of this response. Our data highlight the importance of HCN1 channels to retinal function and suggest that patients with HCN1 pathogenic variants are likely to have a dramatically reduced sensitivity to light and a limited ability to process temporal information.
SIGNIFICANCE STATEMENT Pathogenic variants in HCN1 are emerging as an important cause of catastrophic epilepsy. HCN1 channels are ubiquitously expressed throughout the body, including the retina. Electroretinogram recordings from a mouse model of HCN1 genetic epilepsy showed a marked decrease in the photoreceptor sensitivity to light and a reduced ability to respond to high rates of light flicker. No morphologic deficits were noted. Simulation data suggest that the mutated HCN1 channel blunts light-induced hyperpolarization and consequently limits the dynamic range of this response. Our results provide insights into the role HCN1 channels play in retinal function as well as highlighting the need to consider retinal dysfunction in disease caused by HCN1 variants. The characteristic changes in the electroretinogram open the possibility of using this tool as a biomarker for this HCN1 epilepsy variant and to facilitate development of treatments.
Introduction
Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are found in the heart, and in the CNS and peripheral nervous system. HCN1-4 encode four HCN channel subtypes (Biel et al., 2009), which control neuronal activity through a variety of mechanisms that include the modulation of rhythmic activity, dendritic integration, transmission of synaptic responses, and synaptic plasticity (Biel et al., 2009). Pathogenic variation in HCN1 has emerged as an important cause of epilepsy (Nava et al., 2014; Bonzanni et al., 2018; Marini et al., 2018; Porro et al., 2021). De novo HCN1 pathogenic variants typically cause developmental and epileptic encephalopathies (DEEs), the most severe group of epilepsies characterized by drug-resistant seizures and developmental impairment (Nava et al., 2014; Marini et al., 2018). Despite their localization throughout the CNS and peripheral nervous system (Moosmang et al., 2001), little is known about sensory deficits associated with HCN1 pathogenic variants.
Recently, we reported a child with HCN1 DEE who had difficulty discriminating colors (McKenzie et al., 2022). Although all four HCN channel isoforms/subtypes are expressed in the retina, HCN1 is the most widely expressed isoform being localized to rod and cone inner segments (Moosmang et al., 2001; Demontis et al., 2002). Rod photoreceptors (Phs) express only HCN1 (Moosmang et al., 2001; Müller et al., 2003), whereas cone photoreceptors also express HCN2 and HCN3, although at lower levels (Müller et al., 2003). Additionally, HCN1 channels are found on cone bipolar cells (Müller et al., 2003) as well as amacrine and retinal ganglion cells (Müller et al., 2003).
HCN1 channels have an established role in photoreceptors, generating the hyperpolarization-activated cation current Ih (Moosmang et al., 2001; Demontis et al., 2002; Knop et al., 2008; Tanimoto et al., 2012). HCN1 channels are activated during hyperpolarization in bright light and depolarize the cell toward the dark membrane potential (Vm), making the light response more transient (Baylor et al., 1984; Hestrin, 1987; Barrow and Wu, 2009). Their slow activation plays a role in defining the photoreceptor membrane potential and therefore the activity of other voltage sensitive channels (Della Santina et al., 2012). HCN1 channels may also modulate bipolar cell synaptic transmission (Chevaleyre and Castillo, 2002; Müller et al., 2003; Della Santina et al., 2010), particularly those bipolar cells that show high temporal tuning (Euler et al., 2014) and contribute to motion detection (Yonehara et al., 2013). Thus, response shaping throughout the retina by HCN1 channels has a role in facilitating temporal resolution of visual stimuli. In line with a role for HCN1 in retinal function, genetic deletion of Hcn1 in mice results in prolongation of the scotopic and photopic electroretinogram (ERG) responses (Seeliger et al., 2011; Sothilingam et al., 2016). At mesopic light levels, rod HCN1 channels shape the temporal characteristics of cone responses (Lankford et al., 2022). Moreover, deletion of Hcn1 results in a reduction of the dynamic range, with saturation of the downstream retinal network during bright light stimulation (Knop et al., 2008; Seeliger et al., 2011; Sothilingam et al., 2016). Whether retinal dysfunction also occurs in patients carrying HCN1 pathogenic variants remains unclear.
Here we investigated the impact of an HCN1 pathogenic variant on phototransduction and downstream retinal processing, including temporal resolving capacity using ERG on the Hcn1M294L mouse model of DEE. The Hcn1M294L mouse carries the homolog of the de novo HCN1 M305L variant identified in two unrelated patients with DEEs (Marini et al., 2018; Bleakley et al., 2021). The Hcn1M294L mouse recapitulates many of the features of the human disease, including spontaneous seizures, hyperactivity, and cognitive impairment (Bleakley et al., 2021). This mutation, located in the S5 transmembrane domain, disrupts the voltage dependence, rendering the HCN1 channel constitutively open and resulting in aberrant cation “leak” (Marini et al., 2018; Bleakley et al., 2021; Hung et al., 2021). Neuronal depolarization caused by the cation leak is proposed to underlie excitability leading to seizures (Bleakley et al., 2021).
Materials and Methods
Ethics and animal welfare
Experiments were performed in accordance with the Prevention of Cruelty to Animals Act, 1986, under the guidelines of the National Health and Medical Research Council (NHMRC) of Australia Code of Practice for the Care and Use of Animals for Experimental Purposes. All experiments were approved by the Animal Ethics Committee at the Florey Institute of Neuroscience and Mental Health before commencement. Animals were monitored in line with protocols approved by this committee. Anesthesia and analgesia were used where appropriate. Mice were acclimatized to experimental rooms for at least 30 min before experimentation. At the conclusion of experimentation, mice were killed by cervical dislocation or decapitation following deep anesthesia.
Human procedures and ethics
The parents of the patient provided written informed consent for research participation. The research was approved by the Human Research Ethics Committee of Austin Health (H2007/02,961). Electroretinogram recordings were made from a patient carrying the HCN1M305L variant, with pupil dilation measured using lower eye lid electrodes with skin electrodes (Kriss and Russell-Eggitt, 1992). With the room lights on, a Ganzfeld stimulator (ColorDome, Diagnosys) was used to generate a single moderate-intensity flash (∼5 cd.s/m2) to produce a photopic cone ERG followed by a 30 Hz flicker stimulus of the same intensity.
Generation of the Hcn1M294L mouse model
A guide RNA target site close to amino acid 294 of the mouse Hcn1 gene was identified using the UCSC Genome Browser [https://genome.ucsc.edu/ (accessed June 22, 2021)]. The guide RNA sequence used was 5′-GAAGCAGCATCATGCCAATG-3′. CRISPR RNA (crRNA; Integrated DNA Technologies) was annealed with transactivating crRNA (tracrRNA) to form a functional crRNA:tracrRNA guide RNA duplex. Cas9 nuclease (Integrated DNA Technologies) was incubated with the guide RNAs to form a ribonucleoprotein (RNP) complex. A single-stranded repair oligo containing the desired amino acid change (Integrated DNA Technologies) was used as a repair template. The C147-HCN1-M294L repair template sequence was 5′-TGCTGACAGATATTCCACATGACCTATGACCTCGCCAGTGCTGTGGTGAGGATCTTCAACCTCATCGGAATGTTGCTGCTTCTGTGCCACTGGGATGGCTGTCTTCAGTTCCTGGTTCCCCTGCTGCAGGACTTCCCA-3′. Cas9 nuclease (30 ng/µl), the RNPs (30 ng/µl,) and the single-stranded DNA homology-directed repair template (30 ng/µl) were microinjected into the pronucleus of C57BL/6J zygotes at the pronuclei stage. Injected zygotes were transferred into the uterus of pseudopregnant F1 females.
Animal housing and husbandry
Mice were bred on site at the Florey Institute of Neuroscience and Mental Health. Male heterozygous Hcn1M294L mice were mated with wild-type C57BL/6J female mice, yielding either Hcn1+/+ [wild-type (WT)] or Hcn1M294L/+ offspring. Mice were housed in standard rodent cages with access to dry pellet food and water ad libitum. Room lighting was maintained at <50 lux with a 12 h dark/light cycle (lights on at 8:00 A.M.) to minimize the potential for light-induced retinal injury. For all experiments, both male and female mice were used, with littermates used wherever possible.
General procedures
Ten Hcn1M294L mice and 10 Hcn1WT mice were used. Animals underwent noninvasive assessment of retinal structure and function. Before the procedure, animals were dark adapted overnight, then weighed and anesthetized with an intraperitoneal injection of ketamine (80 mg/kg) and xylazine (10 mg/kg; Troy Laboratory). The mixture was diluted in sterile injectable saline (1:10) to aid with hydration and ease of administration (10 μl/g).
Topical anesthesia and pupil mydriasis were achieved with single drops of proxymetacaine 0.5% (Alcaine, Alcon) and tropicamide 0.5% (Mydriacyl, Alcon). Corneal hydration was maintained with lubricating eye drops during ERG (Celluvisc, Allergan) and imaging (Systane, Alcon). During both ERG and optical coherence tomograph (OCT) assessment, body temperature was maintained at 37°C to prevent hypothermia. At the end of in vivo assessment, anesthetized animals were killed by cervical dislocation and eyes were collected for immunohistochemistry.
Electroretinography
Mice were dark adapted overnight before ERG. Experiments were conducted in a lightproof room, and all preparation was undertaken with only the aid of a dim red headlight to minimize light exposure. Following the induction of anesthesia and mydriasis, animals were lightly secured to a water-heated platform to minimize movement and breathing artifacts. Custom-made silver chloride active electrodes (A&E Metal Merchants) were placed on the apex of the corneas, with the reference ring-shaped electrode placed around the equator of the eye. A grounding needle electrode (Grass) was then inserted subcutaneously into the tail, before positioning the Ganzfeld bowl at eye level. Animals were set up in a Faraday cage to block out electromagnetic noise. Raw signals were processed via a preamplifier (model P511 AC Amplifier, Grass) and a main amplifier (model ML785 Powerlab 8SP, ADInstruments) with a bandpass filter setting of 0.1–1000 Hz. Signals were recorded using Scope software (ADInstruments) and exported for offline analysis. Analysis of line noise characteristics showed that there was no difference between the two strains (Fig. 3D inset).
Scotopic responses.
Following electrode setting a further 10 min was allowed to ensure that the eyes had completely recovered from the red light used during setup. ERGs were then elicited using a range of luminous energies (from −5.53 to 2.07 log cd.s/m2, 0.1–4 ms) starting from the dimmest, where 20 responses were recorded and averaged. As the light intensity increased, the number of repeated flashes was reduced and the time between subsequent flashes was progressively lengthened to avoid adaption. This wide range of stimuli allowed for the isolation of retinal ganglion cell-mediated responses [scotopic threshold response (STR)] at the dimmest light levels, as well as inner retinal inhibitory circuits [oscillatory potentials (OPs)], P2, and photoreceptor (P3) responses at moderate and high scotopic light levels.
Photopic responses.
Following the collection of scotopic responses, animals were light adapted for 10 min to a white background of 140 cd/m2. Once signals had stabilized, a sequence of luminous energies ranging from 0.76 to 2.67 log cd.s/m2 elicited the cone pathway-driven a-wave, b-wave, and OPs. At each intensity, 10 signals were recorded with an interstimulus interval of 2 s. To assess both transient and sustained responses, a longer flash (2.67 log cd.s/m2, 256 ms) was delivered.
Temporal resolution.
To quantify the temporal processing capacity of the eye, we elicited responses using progressively higher temporal frequencies with a fixed light intensity (3.99 log cd.s/m2). Flickering stimuli ranged from 5 to 45 Hz (in 5 Hz steps) with a 50% duty cycle without any background light.
Spectral response.
Spectral response characteristics of the eye were assessed using responses to flicker lights (10 Hz, 50% duty cycle) of varying peak wavelength, generated by a xenon light source with a monochromator (150 W Polychrome V, Till Photonics). Wavelengths ranged from 420 to 680 nm in 20 nm steps (bandwidth, 10 nm).
Photoreceptor P3 analysis.
As previously described (Bui et al., 2005), to analyze the scotopic photoreceptoral response (P3), the first corneal electronegative component of the ERG waveform (the a-wave) was modeled as a function of time (in milliseconds) and intensity (in log cd.s/m2) using a delayed Gaussian function (Eq. 1; Lamb and Pugh, 1992) to quantify the amplitude (i.e., number of receptors and their outer segment length (RmP3; in microvolts) and the sensitivity [amplification of the phototransduction cascade (S); in cd−1.m2/s3]. The delay term (td; in seconds) describes both biochemical and other recording latencies. The model was fit to the two brightest intensities (1.55 and 2.07 log cd.s/m2), which elicit photoreceptor (a-wave) responses of saturated amplitude and free of b-wave intrusion. The model was optimized by changing all parameters to minimize the root mean square (rms) error term (Solver module in Excel, Microsoft). The photopic a-wave was modeled using the same equation fit as an ensemble to the brightest three light levels (2.07, 2.37, and 2.67 log cd.s/m2), as follows:
Bipolar cell P2 and amacrine cell oscillatory potential analysis.
For both scotopic and photopic waveforms, the bipolar cell-dominated P2 component of the ERG was quantified by first subtracting the family of P3 responses from their respective intensities of the raw ERG waveform. The derived waveform contains the underlying P2 and the OPs. The OPs provide information about inner retinal inhibitory pathways involving amacrine cells (Wachtmeister, 1998). The OPs and P2 were isolated by bandpass (−3 dB at 50 and 280 Hz) and low-pass filtering (−3 dB at 50 Hz), respectively. OPs were then quantified by the amplitude from baseline of their peak response at the brightest light level (scotopic, 2.07 log cd.s/m2; photopic, 2.67 log cd.s/m2).
The response amplitude (P2max; in microvolts) and sensitivity (1/K, log cd.s/m2) of bipolar cells was then quantified by fitting the P2 amplitude as a function of all stimulus intensities (i; log cd.s/m2) using a saturating hyperbolic function (Eq. 2). The function was fit to the range of either scotopic or photopic responses by minimizing the rms error term (Excel), as follows:
Ganglion cell scotopic threshold response.
Ganglion cell responses are known to dominate the rodent ERG waveform at light levels near absolute threshold (Saszik et al., 2002; Bui and Fortune, 2004). The amplitude of the positive lobe of the STR (pSTR) was taken at −5.01 and −4.90 log cd.s/m2.
Optical coherence tomography
Following ERG recordings, animals underwent in vivo retinal structure imaging using spectral-domain OCT (Spectralis SD-OCT2, Heidelberg Engineering). A drop of ocular lubricant gel (Genteal Tears, Alcon) and a coverslip were used to rehydrate and remove any opacities that developed from general anesthesia used during ERG recordings. Once clear optics were achieved (10–15 min), excess eye gel was removed, and a single drop of eye lubricant (Systane, Alcon) was applied to improve tear film optics for OCT imaging. The center of cornea was aligned to the lens of OCT via a three-dimensional adjustable platform. Images were acquired with a volumetric scan pattern (7.6 × 6.3 × 1.9 mm) centered over the optic nerve head. Each volume scan consisted of 121 vertical B-scans (five repeats), and each B-scan was made up of 768 A-scans (3.87 µm axial, 9.8 µm lateral resolution).
Analysis used the automated segmentation algorithm of the manufacturer software (HEYEX, Heidelberg Engineering), which was checked for errors. Thicknesses were averaged across an annulus from 3 to 6 mm (i.e., Early Treatment Diabetic Retinopathy Study outer ring) for the retinal nerve fiber layer (RNFL), ganglion cell inner plexiform layer (IPL), inner nuclear layer (INL), outer plexiform layer (OPL), outer nuclear layer (ONL), photoreceptor (Ph), as well as total retinal thickness (TRT).
Immunohistochemistry
Following in vivo assessment, animals were killed, and eyes were enucleated and immersion fixed in 4% paraformaldehyde prepared in 0.1 m PBS, pH 7.4, for 15 min at room temperature, and cryoprotected in 30% sucrose at 4°C overnight. Eyes were then embedded in Tissue-Tek O.C.T. Compound (Miles), and frozen at −20°C. Serial sections (thickness, 14 µm) were collected on Superfrost Plus slides. For immunostaining, sections were permeabilized in PBS containing 0.1% Triton X-100 for 10 min, followed by 1 h incubation in CAS-Block (catalog #008120, Thermo Fisher Scientific) to block nonspecific labeling. Sections were incubated O/N at 4°C with polyclonal rabbit anti-HCN1 (1 : 100; catalog #APC-056 Alomone Labs), washed (0.1 m PBS, 3 × 5 min) and incubated with secondary antibody conjugated to Alexa Fluor-488 for 2 h at room temperature (1 : 500; Thermo Fisher Scientific). Sections were visualized with a confocal microscope (model LSM780, Zeiss), and image files were processed with Zen Blue software (Carl Zeiss).
Photoreceptor modeling and simulations
A passive electrical model (Table 1) for a single photoreceptor cell was based on previously published models comprising parallel conductances, membrane capacitance, and photocurrent source (Kamiyama et al., 1996; Liu and Kourennyi, 2004; Publio et al., 2006). Current through each conductance was defined by Equation 3, as follows:
Modeling parameters
In Equation 3,
Statistical analysis
All statistical analyses were performed using Prism software (version 9; GraphPad Software). Data are reported and plotted as the mean ± SEM. Other data were first analyzed using a Shapiro–Wilk test for normality. Statistical significance was then determined using a two-tailed Student's t test (paired or unpaired, as appropriate) for normally distributed data, or a Mann–Whitney U test for non-normally distributed data. For parameters collected across a range of stimulus conditions a two-way repeated-measures (RM) ANOVA was used to compare the two mouse strains, with Sidak post hoc test used at specific stimulus conditions. Statistical significance was set at p < 0.05.
Results
Retinal function and structure in Hcn1M294L mice
Figure 1A shows the average scotopic ERG waveform at selected intensities in both Hcn1WT and Hcn1M294L mice. At the brightest intensity (2.07 log cd.s/m2), there was a small deficit in the photoreceptoral a-wave, with moderate intensities highlighting a bipolar cell b-wave attenuation (−2.75 log cd.s/m2). At the dimmest light levels, the ganglion cell pSTR was also reduced in size (−5.01 and −4.90 log cd.s/m2). Modeling of the a-wave leading edge revealed no difference in phototransduction amplitude (Fig. 1B; t(14) = 1.39, p = 0.19, t test) but significantly reduced sensitivity to light (Fig. 1C; t(14) = 2.67, p = 0.018, t test) in Hcn1M294L mice compared with their WT counterparts. Downstream of photoreceptors, rod bipolar cell responses were both significantly reduced in amplitude (Fig. 1D; t(14) = 3.78, p = 0.002, t test) and less sensitive to light (Fig. 1E; t(14) = 4.86, p = 0.0003, t test). Additionally, scotopic ganglion cell pSTRs (Fig. 1F; t(14) = 4.80, p = 0.0003, t test) were significantly reduced, whereas amacrine cell OPs were not affected (Fig. 1G; t(14) = 0.65, p = 0.53, t test).
Deficits in scotopic and photopic electroretinogram responses in Hcn1M294L mice. A, Group averaged Hcn1WT mice (WT, n = 8, blue) and Hcn1M294L mice (M294L, n = 8, orange) scotopic waveforms at selected light levels. Photoreceptor (a-wave), bipolar cell-mediated (b-wave) and amacrine cell (Oscillatory potential, OPs) components are analyzed at brighter light levels, whereas the ganglion cell pSTR was isolated at very dim light levels (−4.9 and −5.01 log cd.s/m2). Rod photoreceptor response were analyzed by modeling the phototransduction cascade. B, Group average (±SEM) scotopic phototransduction amplitude (RmP3) for Hcn1WT (n = 8, blue) and Hcn1M294L (n = 8, orange) mice. Individual eyes are shown. C, Scotopic phototransduction amplification (log S). D, Scotopic bipolar cell amplitude (Vmax). E, Scotopic bipolar cell sensitivity (K). F, pSTR. G, Scotopic amacrine cell amplitude (OPs). H, Group-averaged photopic waveforms at selected light levels. I, Photopic phototransduction amplitude. J, Photopic phototransduction amplification. K, Photopic bipolar cell amplitude. L, Photopic bipolar cell sensitivity. M, Amplitude of transient b-wave measured after light onset. N, Amplitude of sustained b-wave measured just before stimulus offset (256 ms). Additional analysis of ON and OFF components can be found in Figure 2. For all panels: ns, not significant; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.
Hcn1M294L mice that showed photopic ERG responses elicited with short flashes on a bright-light background, have a smaller b-wave that was less sustained than Hcn1WT mice (Fig. 1H). This effect is readily appreciated with the longer flash, whereby Hcn1M294L mice also showed a transient component but only a very small, sustained b-wave tail. Modeling the photopic a-wave (Fig. 1H, inset) showed that cone photoreceptors in Hcn1M294L mice produced response amplitudes similar to those in Hcn1WT controls (Fig. 1I; t(14) = 0.69, p = 0.49, t test), but were significantly less sensitive to light (Fig. 1J; t(14) = 3.29, p = 0.005, t test). These photoreceptoral deficits resulted in similar cone bipolar cell response amplitudes (Fig. 1K; t(14) = 1.95, p = 0.07, t test), which were also less sensitive to light (Fig. 1L; t(14) = 4.20, p = 0.0009, t test). Analysis of the longer flash showed that both transient (Fig. 1M; t(14) = 2.45, p = 0.028, t test) and sustained (Fig. 1N; t(14) = 6.97, p < 0.0001, t test) components were attenuated in Hcn1M294L mice compared with WT controls. Figure 2 shows that in Hcn1M294L mice the OFF component was more attenuated in amplitude (Fig. 2B–D) and was faster (Fig. 2E–G) in comparison with the ON component.
Comparison of ON and OFF electroretinogram responses in Hcn1M294L mice. A, A longer flash (256 ms) was used to elicit the murine ERG response from ON and OFF components. Amplitudes and implicit times were taken after light onset, as well as following light offset. OFF amplitudes were taken as the difference between the amplitudes at light offset (256 ms) and the OFF trough; also, the implicit time of the OFF trough. B, Group-averaged Hcn1WT mice (WT, n = 8, blue) and Hcn1M294L mice (M294L, n = 8, orange) ON amplitude (t(14) = 2.5, p = 0.03). C, OFF amplitude (t(14) = 6.9, p < 0.0001). D, Amplitudes of ON and OFF components are compared by expressing the amplitudes in Hcn1M294L mice relative to Hcn1WT mice (t(14) = 4.8, p = 0.0003). E, ON-peak time (t(14) = 7.1, p < 0.0001). F, OFF-trough time (t(14) = 4.1, p = 0.0012). G, Implicit times of ON and OFF components are compared by expressing implicit time for Hcn1M294L mice relative to Hcn1WT mice (t(14) = 4.2, p = 0.0009). *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.
Given previous reports of color vision difficulties in human HCN1 DEE (McKenzie et al., 2022), we sought to examine responses elicited by colored stimuli to identify the potential functional absence of photoreceptors. Additionally, we assessed the temporal characteristics of Hcn1M294L retinal responses using full-field flickering stimuli. Figure 3A shows that waveforms elicited using 10 Hz flickering stimuli of various peak colors were smaller in Hcn1M294L mice compared with Hcn1WT mice. Responses to various colored stimuli were smaller in Hcn1M294L mice (Fig. 3B; interaction F(1,112) = 490.0, p < 0.0001, RM-ANOVA), particularly for 420–600 nm stimuli (Sidak multiple-comparisons test, p < 0.01). However, when normalized to the maximal amplitude, the spectral response function was similar between the two groups (Fig. 3C). Response amplitudes for flicker rates from 5 to 45 Hz differed between Hcn1M294L and Hcn1WT mice (Fig. 3D; interaction: F(1,126) = 107.0, p < 0.001, RM-ANOVA), and post hoc comparison between groups identified significant attenuation in Hcn1M294L mice for frequencies from 5 to 20 Hz (Sidak multiple-comparisons test, p < 0.05).
Deficits in temporal responses in Hcn1M294L mice. A, Representative Hcn1WT mice (WT, blue) and Hcn1M294L mice (M294L, orange) waveforms at all wavelengths tested. B, Group average (±SEM) response to stimuli of different colors for Hcn1WT mice (n = 5, blue) and Hcn1M294L mice (n = 5, orange). C, Normalized spectral response functions for Hcn1WT mice (n = 5, blue) and Hcn1M294L mice (n = 5, orange). D, Group average (±SEM) responses for Hcn1WT mice (n = 8, blue) and Hcn1M294L mice (n = 8, orange) as a function of flicker frequency. *p < 0.05 post hoc test between groups. Quantification of noise amplitude at 50 Hz for all flicker frequencies can be found in the inset. Group average (±SD) noise amplitude obtained from spectral analysis of the ERG signal at 50 Hz in Hcn1WT (WT, n = 8, blue) and Hcn1M294L mice (M294L, n = 8, orange) collected across the range of temporal frequencies used. Two-way repeated-measures ANOVA showed that there was no statistically significant interaction between genotype and frequency (F(1,8) = 0.09, p = 0.99) or a difference in genotype (F(1,8) = 1.65, p = 0.19) or frequency (F(1,8) = 0.54, p = 0.83).
Retinal structure was assessed using in vivo optical coherence tomography imaging (Fig. 4A). There was no significant difference in layer thickness between Hcn1M294L and Hcn1WT mice, as follows: RNFL (Fig. 4B: t(18) = 0.067, p = 0.95, t test); GCL (Fig. 4C: t(18) = 1.79, p = 0.09, t test); IPL (Fig. 4D: t(18) = 1.08, p = 0.29, t test); INL (Fig. 4E: t(18) = 0.95, p = 036, t test); OPL (Fig. 4F: t(18) = 1.60, p = 0.13, t test); ONL (Fig. 4G: t(18) = 1.36, p = 0.19, t test); PR (Fig. 4H: t(18) = 0.43, p = 0.67, t test); and TRT (Fig. 4I: t(18) = 1.46, p = 0.16, t test).
Normal gross retinal structure in Hcn1M294L mice. A, Example of optical coherence tomography image of the mouse retina centered at the optic nerve, with the segmentation location of the retinal layers shown. B, Group averaged Hcn1WT mice (WT, n = 10, blue) and Hcn1M294L mice (M294L, n = 10, orange) RNFL thickness. C, Ganglion cell layer (GCL). D, inner plexiform layer (IPL). E, inner nuclear layer (INL). F, outer plexiform layer (OPL). G, outer nuclear layer (ONL). H, Photoreceptor (PR) thickness. I, total retinal thickness (TRT). J–L, Retinal cross section from Hcn1WT mice stained for HCN1 channels (green) and counterstained for cell nuclei with DAPI (blue). M–O, Retinal cross sections from Hcn1M294L mice. ns, Not significant.
Immunohistochemistry showed that there were no gross retinal abnormalities in the Hcn1M294L mouse retina (Fig. 4M–O) compared with Hcn1WT mice (Fig. 4J–L) and no difference in the HCN1 channel location, which is highly expressed in the photoreceptors as well as in the outer and inner plexiform layers.
ERG responses in a female patient with HCN1 DEE
Table 2 compares the available parameters for the right and left eyes with normal data (Bradshaw et al., 2004; Meredith et al., 2004) collected using similar stimulus settings and electrodes (Kriss et al., 1992; Kriss and Russell-Eggitt, 1992; Kriss, 1994; Bradshaw et al., 2004; Marmoy et al., 2022). While the available normative ranges are for slightly older children (4–14 years), previous studies have suggested that the major ERG components are similar in amplitude and implicit time to adult ERGs by 6–10 months of age (Fulton et al., 1989). This comparison suggests that ERG responses for the patient fall below the most conservative 95% confidence limit for b-wave amplitude and timing, as well as for 30 Hz flicker amplitude. While the implicit time of the a-wave in the left eye was delayed, this was not the case for the right eye. The mouse ERG data reported above broadly recapitulate these clinical observations.
Clinical data from a patient with HCN DEE compared with normal data collected using similar stimulus settings and skin electrodes
Computer simulation modeling of the photoreceptor responses
To gain further insight into the cellular mechanism underlying the altered ERG response observed for Hcn1M294L mice, we performed simulations using a model for a single photoreceptor. The model was based on previous studies (see Materials and Methods) and incorporated a leak, the nonspecific potassium conductance, and an L-type calcium conductance as well as a hyperpolarization-activated (Ih) component. For the Ih-related conductance, we used a 10-state allosteric model (Altomare et al., 2001; Hung et al., 2021) to predict the channel open probability. This was reduced to a four-state model to account for the altered kinetics of Hcn1M294L channels. We adjusted the original parameter set in Hung et al. (2021) to account for the faster mouse kinetics using brain slice voltage-clamp data as a template (Bleakley et al., 2021). As shown in Figure 5 (Table 1) under simulated voltage-clamp conditions, the hyperpolarization-activated current for Hcn1M294L mice shows a time-dependent activation component superimposed on an instantaneous current step, as we have previously described (Bleakley et al., 2021; Hung et al., 2021). The steady-state current–voltage relationships (Fig. 5C) from these data underscored the leaky behavior of Hcn1M294L channels and lacked the inward rectification characteristic of Hcn1WT channels over the physiologically relevant voltage range. The effect of HCN1 channels on simulated photoreceptor membrane potential is illustrated in Figure 5D. We assumed equal HCN1 maximum conductance for the WT and M294L cases. This was based on similar immunohistochemical expression patterns in the retina (Fig. 4) and similar maximal currents measured in layer 5 pyramidal neurons (Bleakley et al., 2021).
A, B, Simulated voltage-clamp records for Hcn1WT mice (A) and Hcn1M294L mice (B) showing HCN1-related current in response to a series of voltage steps from −110 to −45 mV in 5 mV intervals from a −50 mV holding potential and returning to −100 mV. The pulse to −110 mV has been highlighted for clarity. This pulse protocol matches that used in Bleakley et al. (2021, their Fig. 5A) from mouse brain slice whole-cell patch recordings. C, The corresponding steady-state current–voltage relationships normalized to −100 mV show a clear absence of inward rectification for M294L (orange) compared with the WT (blue). D–F, Simulated response of the generic photoreceptor to a bright flash (see text and inset) for membrane potential (Vm) for Hcn1WT mice (blue) and Hcn1M294L mice (orange) and in the absence of HCN1 (black; D), HCN-related current (Ih; E), and voltage-dependent ICa (F) for WT (blue traces) and M294L (orange traces), and in the absence of HCN1 (black traces). Note the timescale offset to allow steady state to be reached before the flash onset. Details on modeling parameters can be found in Table 1.
Under dark-adapted conditions preceding the flash, Hcn1WT and Hcn1M294L channels would not contribute significantly to the net membrane conductance as the predicted membrane potentials were −28 mV (Hcn1WT) and −32 mV (Hcn1M294L), respectively: Hcn1WT channels would be closed, whereas for Hcn1M294L channels, the predicted Vm corresponds to their reversal potential (Fig. 5B,D). With Hcn1WT in the model, the flash induced an instantaneous hyperpolarization, and as Hcn1WT channel-activated hyperpolarization was suppressed throughout the time course of the simulated photocurrent, giving rise to a pronounced “nose” in the Vm response (Barrow and Wu, 2009; Della Santina et al., 2010). In contrast, for Hcn1M294L, hyperpolarization was significantly weaker during the flash and flatter in the subsequent plateau phase as the mutant channels conduct a larger inward current (Fig. 5D). As a consequence of the effect of HCN1 channel activity in shaping the time course of Vm, the corresponding time course of the L-type Ca current (ICa) was also affected (Fig. 5F). At the flash onset, our model predicted that ICa was fully suppressed for Hcn1WT as Vm would be sufficiently hyperpolarized to close these channels and subsequently they only begin to open toward the end of the simulated photoreceptor current time course as Vm recovers. In contrast, for the Hcn1M294L, our simulations predicted only partial suppression of ICa at flash onset with an inward ICa of ∼0.7 pA during the same period. We expect that the presence of a small, standing ICa during and following a bright-flash ICa would have consequences for the suppression of glutamate release at the synapse between photoreceptor and bipolar cells.
Discussion
In this study, we show that an HCN1 pathogenic variant has a significant impact on retinal function. Most significantly, photoreceptoral responses under both scotopic and photopic conditions were significantly less sensitive to light. Furthermore, postphotoreceptor responses were smaller, less sensitive to light and showed an impaired capacity to resolve stimuli of high temporal frequency. These mouse results are in line with ERG recordings from a patient with HCN DEE because of the HCN1 M305L pathogenic variant. Visual deficits driven by retinal dysfunction are therefore highly likely to be a feature of HCN1 DEE caused by pathogenic variants, resulting in cation leak.
HCN1 channels are abundantly expressed in rod and cone photoreceptor inner segments (Müller et al., 2003; Ivanova and Müller, 2006; Della Santina et al., 2010) and are activated after a short delay to drive the membrane potential back to the resting state (Knop et al., 2008; Della Santina et al., 2010; Seeliger et al., 2011). This makes the photoreceptoral response faster, creating the characteristic nose profile (Barrow and Wu, 2009; Della Santina et al., 2012). Previous studies show that HCN1 channel knockout (KO) does not affect the ERG a-wave and only manifests as a prolongation of the b-wave in mice (Seeliger et al., 2011; Sothilingam et al., 2016). This is consistent with the fact that HCN1 channels are only open after light has closed cyclic nucleotide-gated channels (Moosmang et al., 2001; Demontis et al., 2002; Knop et al., 2008; Tanimoto et al., 2012). As such, HCN1 channels are not expected to affect the a-wave or light sensitivity (Della Santina et al., 2012). However, we show here that the Hcn1M294L variant causes a reduced a-wave sensitivity to light as well as a downstream reduction in b-wave sensitivity to light. A key difference between Hcn1 knock-out and Hcn1M294L knock-in mice is that the knock-in missense variant would impact cellular excitability. By being constitutively open, abnormal inward cation flux through the mutant channels is likely to impact passive cell properties. This idea is supported by previous reports that constitutively open HCN1 channels cause a significant change in resting membrane potential in cortical and hippocampal neurons (Bleakley et al., 2021). Furthermore, a standing “cation leak” will shape light-induced changes in membrane potential. Consistent with this, our computer simulations (Fig. 5D) suggest that light-induced hyperpolarization was weaker at flash onset because of the presence of a standing leak conductance. Our simulations also indicate that following the initial hyperpolarization, Vm shows a flatter plateau phase (less prominent nose) as the mutant channels conduct a larger inward current at these voltages when compared with WT channels, and the time-dependent activation component constitutes a smaller proportion of total Ih (Fig. 5D,E). Reduced light-induced hyperpolarization along with impaired shaping of the photoreceptoral response might account for slower scotopic (Fig. 1C) and photopic (Fig. 1J) a-wave responses (i.e., less sensitivity to light).
Our simulations show that in the complete absence of HCN1, light-induced hyperpolarization is larger, but has a flatter plateau phase. This is consistent with patch–clamp photoreceptor recordings from Hcn1–/– as well as Hcn1WT mice in the presence of the channel blocker ivabradine (Della Santina et al., 2012). A similar flatter plateau phase is also described for the HCN1 pathogenic variant, a consequence of the significantly attenuated voltage dependence of the channels. The lack of this dynamic change in Vm during a light stimulus may account for the fact that deficits in temporal processing are seen in both the mutant and knock-out cases (Seeliger et al., 2011; Sothilingam et al., 2016).
Seeliger et al. (2011) show in Hcn1−/− mice that if the rods are silenced, then cone responses improve. Thus, at mesopic light levels HCN1 channels on rods influence the cone response (Lankford et al., 2022). Here we show, using photopic conditions, which saturate rod photoreceptors, that there is a robust deficit in cone-mediated responses both in response to short flashes as well as with high temporal frequency flickering stimuli. This shows that in the Hcn1M294L mouse HCN1 channels on both rod and cone receptors are impacted. Indeed, by fitting the a-wave leading edge, we found that the reduction in light sensitivity was greater in the photopic ERG compared with the scotopic mixed response (−0.23 ± 0.05 vs −0.08 ± 0.03 log units, t(14) = 2.39, p = 0.03, t test).
Our simulation modeling suggests that light-induced changes in Ca2+ concentration in the simulated rod response is reduced, which will result in the glutamate release during a light flash being attenuated, potentially explaining the downstream reduction of the bipolar b-wave amplitude and sensitivity under scotopic conditions and the reduced sensitivity under photopic conditions (Fig. 1D,E,L). It also is likely that the combination of slower and more prolonged photoreceptoral responses (Fig. 4E) along with reduced glutamatergic neurotransmission between rod photoreceptors and bipolar cells accounts for the impairment of flicker fusion frequency in Hcn1M294L mutant mice (Fig. 2). However, a differential effect of the mutation on ON and OFF responses, as observed in Figure 2, and as has been previously reported following ivabradine blockade of HCN1 channels (Bemme et al., 2017), may also contribute to the impairment of temporal responses in Hcn1M294L mutant mice.
It is noteworthy that in Hcn1−/− mice the rod b-wave is prolonged (Seeliger et al., 2011; Sothilingam et al., 2016), whereas in Hcn1M294L mutant mice the ERG b-wave is more truncated, which is particularly apparent in the photopic ERG (Fig. 1H). Despite these difference in the single flash ERG, we observed reductions in the flicker fusion frequency (Fig. 3B) that are consistent with the deficit seen in Hcn1−/− mice (Seeliger et al., 2011; Sothilingam et al., 2016). Based on our modeling, we speculate that in Hcn1−/− mice a larger and more prolonged light-induced photoreceptoral hyperpolarization produces a larger and more sustained b-wave. Furthermore, the absence of temporal shaping of photoreceptoral response results in impaired temporal responses (reduced flicker fusion frequency). In the case of the Hcn1M294L mutant mouse, a smaller light-induced hyperpolarization results in a smaller b-wave, and at the same time the loss of temporal shaping leads to the same flicker fusion deficit as in Hcn1−/− mice.
Within the scotopic pathway, we also noted that the reduction in bipolar cell b-wave light sensitivity (−0.42 ± 0.25 log units) was much greater than would be predicted based on the attenuation of the a-wave sensitivity (−0.08 ± 0.09 log units) to light. Furthermore, the reduction in the ganglion cell-mediated pSTR amplitude is more than would be expected based on attenuation of bipolar cell b-wave amplitude. This raises the possibility that within the rod pathway abnormal HCN1 channel activity on inner retinal neurons may contribute to the phenotype. Specifically, HCN1 channels are also found on bipolar cells and are concentrated at synaptic terminals and play a role in modulating bipolar cell synaptic transmission (Chevaleyre and Castillo, 2002; Müller et al., 2003; Della Santina et al., 2010). This is true particularly for bipolar cells that show high temporal tuning (Euler et al., 2014) and contribute to motion detection pathways (Yonehara et al., 2013). Recently, Hellmer et al. (2020) show that HCN currents are important for shaping of subgroups of bipolar cells and that such current can be modulated by a range of receptors. Whether the deficit in temporal processing reflects change to photoreceptoral as well as inner retinal HCN channels requires further investigation. Moreover, studies quantifying synaptic terminals mediating communication between photoreceptor and bipolar cells would be important to fully understand the nature of the functional deficits seen here. Furthermore, it may be likely that development and refinement of retinal circuitry is also impacted. The potential for developmental delays has been addressed through comprehensive characterization demonstrating that Hcn1M294L show clear behavioral deficits, some of which are indicative of developmental delays (Bleakley et al., 2021).
One may have expected a proportional decrease in the OPs, given the attenuation of photoreceptoral and bipolar cell-mediated responses seen in the Hcn1M294L mutant mouse, as is the case in Hcn-/- mice (Knop et al., 2008). However, the impact of HCN1 cation leak may vary depending on retinal location. For instance, HCN1 channels may modulate OPs via their expression in the IPL on ON–bipolar cell axon terminals, as well in the OFF layer of the IPL on amacrine cells (Müller et al., 2003; Knop et al., 2008). Furthermore, HCN channels can modulate GABAergic neurotransmission (Cai et al., 2022), which can impact bipolar cell response shaping. Should aberrant HCN channels modulate the balance of GABAergic signaling from amacrine to bipolar cells, this can lead to changes in OPs, and in some cases can lead to paradoxically larger OPs (Kapousta-Bruneau, 2000; McCall et al., 2002), despite upstream a-wave and b-wave attenuation.
The clinical ERG assessment of light-adapted response in a patient with the same HCN1 pathogenic variant, which leads to constitutively open HCN1 channels, shows a pattern of smaller b-wave and 30 Hz flicker responses, which is recapitulated in the Hcn1M294L mouse. While there has been a previous report of color vision deficits (McKenzie et al., 2022), this is the first report of reduced temporal processing in HCN1 DEE. Critically, these data indicate that individuals with a cation-leak HCN1 pathogenic variant are likely to have impaired vision, including reduced sensitivity to light and significant problems with processing visual stimuli of moderate to high temporal frequencies. Furthermore, altered dynamics of the Ca2+ concentration in photoreceptors is also likely to impact light adaptation (Chen et al., 2010). Such findings may have clinical implications for patients, given that our eyes are constantly sampling information of images projected onto the retina. Information gathered periodically is integrated so that objects around us appear to be stable or move smoothly. Thus, good temporal processing is important to all visual activities. Further understanding of the precise visual deficit will help clinicians to design visual environments that optimize visual function.
In summary, this is the first study to examine the visual impact of a pathogenic variant in HCN1 that leads to a cation leak. Together, our data identify the impact of HCN1 pathogenic variants on vision and provide further insights into the important role that HCN1 channels play. Finally, our study highlights the potential usefulness of clinical measures of comorbid states, such as visual dysfunction, which may provide useful biomarkers for monitoring and assessing treatment efficacy in relevant forms of epilepsy.
Footnotes
↵*D.Z., P.P.-G., and L.E.B. are joint first authors.
↵†C.A.R. and B.V.B. are joint senior authors.
This research was supported by National Health and Medical Research Council (NHMRC) Program Grant 10915693 to C.A.R. and I.E.S. I.E.S. is also the recipient of NHMRC Practitioner Fellowship 1104831 and Investigator Grant 1172897. This work was made possible through the Victorian State Government Operational Infrastructure Support and Australian Government NHMRC Independent Research Institute Infrastructure Support Scheme. We thank the patient and her family for participating in our research.
The authors declare no competing financial interests.
- Correspondence should be addressed to Bang V. Bui at bvb{at}unimelb.edu.au