- Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, MI, United States
Efficient excitation-contraction coupling in ventricular myocytes depends critically on the presence of the t-tubular network. It has been recently demonstrated that cholesterol, a major component of the lipid bilayer, plays an important role in long-term maintenance of the integrity of t-tubular system although mechanistic understanding of underlying processes is essentially lacking. Accordingly, in this study we investigated the contribution of membrane cholesterol to t-tubule remodeling in response to acute hyposmotic stress. Experiments were performed using isolated left ventricular cardiomyocytes from adult mice. Depletion and restoration of membrane cholesterol was achieved by applying methyl-β-cyclodextrin (MβCD) and water soluble cholesterol (WSC), respectively, and t-tubule remodeling in response to acute hyposmotic stress was assessed using fluorescent dextran trapping assay and by measuring t-tubule dependent IK1 tail current (IK1,tail). The amount of dextran trapped in t-tubules sealed in response to stress was significantly increased when compared to control cells, and reintroduction of cholesterol to cells treated with MβCD restored the amount of trapped dextran to control values. Alternatively, application of WSC to normal cells significantly reduced the amount of trapped dextran further suggesting the protective effect of cholesterol. Importantly, modulation of membrane cholesterol (without osmotic stress) led to significant changes in various parameters of IK1,tail strongly suggesting significant but essentially hidden remodeling of t-tubules prior to osmotic stress. Results of this study demonstrate that modulation of the level of membrane cholesterol has significant effects on the susceptibility of cardiac t-tubules to acute hyposmotic stress.
Introduction
Transverse-axial t-tubular system (TATS) is a complex network of invaginations of the surface membrane necessary for efficient excitation-contraction coupling in cardiac and skeletal muscle cells. In particular, a dense TATS is present in adult ventricular myocytes of likely all mammalian hearts. The TATS becomes significantly remodeled in response to a variety of experimental challenges [e.g., osmotic stress (Kawai et al., 1999; Moench et al., 2013)] and in clinically relevant stress conditions (e.g., heart failure; Guo et al., 2013, for review). It has been well established that a number of t-tubular structural proteins are involved in the maintenance and remodeling of TATS (Kline and Mohler, 2013 for review). However, it is less clear how other membrane components affect TATS structure and remodeling.
T-tubular system membranes are highly enriched in cholesterol (Sumnicht and Sabbadini, 1982) and are stiffer than other biological membranes (Hidalgo, 1985). Cholesterol is an important player in the formation of membrane rafts including caveolae (Simionescu et al., 1983; Rothberg et al., 1990), which can also be found in cardiac t-tubules (Page, 1978). Changes in cellular cholesterol has also been suggested to contribute to electrical remodeling in cardiomyocytes and may have pro- or anti-arrhythmic effects (Coronel, 2017). However, it remains unknown whether in cardiac t-tubules cholesterol exerts its effects through its numerous molecular targets or it can affect their structural stability. Interestingly, patients undergoing statin therapy sometimes experience myalgia associated with skeletal muscle t-tubular dilation and vacuolation (Mohaupt et al., 2009). These skeletal muscle t-tubular defects could be recapitulated in vitro in isolated human and mouse muscle fibers by extracting cholesterol using methyl-β-cyclodextrin (MβCD) (Draeger et al., 2006). Recent work has shown that modification of cholesterol content in cardiomyocytes affects cardiac TATS structure (Zhu et al., 2016) although another report has contradicted this initial finding (Gadeberg et al., 2017). One potential underlying reason contributing to this controversy is that these studies, whether using whole animal or isolated cell models, are focused on the analysis of the final result of long term remodeling of TATS and thus may miss the immediate and largely hidden changes in the TATS properties.
Given the important structural role that cholesterol plays in TATS, we hypothesize that one of the likely early consequences of membrane cholesterol modification could be a change in the susceptibility of TATS to various stresses. In this study we show that manipulation of membrane cholesterol leads to significant changes in the susceptibility of TATS to hyposmotic challenge while having no significant changes in the overall appearance of TATS. Furthermore, cholesterol modulation causes significant changes in the electrophysiological properties of cardiomyocytes reflecting diffusional accessibility of TATS that may indicate changes in the underlying t-tubular structure.
Materials and Methods
Animals
This study was carried out in accordance with the recommendations provided in the Guide for the Care and Use of Laboratory Animals (8th edition; The National Academic Press, Washington, DC, United States). The protocol was approved by the veterinary staff of the University Committee on Use and Care of Animals at the University of Michigan.
Two- to six-months old male and female C57BL/6 mice were included in this study.
Solutions (mM)
All solutions were filtered using a 0.22 μm filter and pH adjusted to 7.35 with NaOH. Osmolarity was measured in previous study (Uchida et al., 2016) using a Vapro 5520 osmometer (Wescor, France; mean ± standard deviation; sample size = 3).
Modified Tyrode’s solution (Tyr; 281 ± 4 mOsm/l): 137 NaCl, 5.4 KCl, 0.5 MgCl2, 0.3 CaCl2, 0.16 NaH2PO4, 3 NaHCO3, 5 HEPES, 10 glucose.
Myocyte storage solution (C solution; 290 ± 3 mOsm/l): 122 NaCl, 5.4 KCl, 4 MgCl2, 0.16 NaH2PO4, 3 NaHCO3, 15 HEPES, 10 glucose, 5 mg/mL of bovine serum albumin, 1.38 mg/mL taurine.
Hyposmotic Tyrode’s solution (0.6 Na; 186 ± 3 mOsm/l): prepared as Tyr but with 60% of NaCl.
Hyposmotic (0.7 Na; 211 ± 2 mOsm/l) solution: prepared by mixing 0.6 Na and Tyr solution in a 3:1 ratio.
Chemicals
HEPES (Calbiochem, United States); KCl, NaHCO3, NaH2PO4 (Mallinckrodt Chemicals, United States); 3 kDa tetramethylrhodamine dextrans in anionic, lysine fixable form (Thermo Fisher Scientific Inc., Waltham, MA, United States). Collagenase (Type 2) (Worthington Biochemical Corp., Lakewood, NJ, United States). Methyl-β-cyclodextrin (MβCD; C4555), Water Soluble Cholesterol (WSC; C4951), Filipin (F9765) and all other chemicals and reagents were purchased from Sigma, St. Louis, MO, United States.
Isolation of Ventricular Myocytes
Myocytes were isolated from the hearts essentially as described in the study by Moench and Lopatin (2014) and used for experiments within 1–8 h post-isolation.
Cholesterol Modulation
Methyl-β-cyclodextrin was dissolved in C solution at concentrations of 4, 8, and 12 mg/mL, corresponding to approximate concentrations of 3, 6, and 9 mM. WSC was dissolved in C solution at 5 mg/mL, corresponding to approximately 50 μM cholesterol. All experiments, unless specifically noted, were carried out at room temperature (RT, ≈19–22°C). One of the main reasons for using RT vs. 37°C was to minimize damaging effects of stronger and faster cholesterol depletion at higher temperatures, as well as to reduce general stress on cardiomyocytes.
Filipin Staining
Quantification of membrane cholesterol was performed using staining with filipin with all procedures carried out at RT. Cardiomyocytes were fixed for 20 min using 2% paraformaldehyde, washed with paraformaldehyde-free solution and incubated with 65 μg/mL filipin for 20 min. After washing out filipin cardiomyocytes were imaged on Nikon TE 20000 microscope using 60× oil immersion objective (NA = 1.4) and CoolSnapEZ camera (Photometrix, Tucson, AZ, United States). Filipin was excited at 340 nm and the emission was collected using 400 nm dichroic mirror and no emission filter in the light path. The data from individual images were corrected for system background, average fluorescence of filipin-free cardiomyocytes (<1–3% of the useful signal) and time-dependent decline in filipin fluorescence due to its washout from the cells (∼30% per hour).
Electrophysiological Measurements
Ionic currents were recorded in the whole-cell configuration essentially as described in a previous study (Uchida et al., 2016). In brief (Figures 2A,B), accumulation of K+ in TATS was induced by applying 400 ms depolarizing step to +50 mV in order to activate voltage-dependent K+ currents. Repolarization back to the holding potential of -75 mV leads to appearance of inward-going IK1 current, known as IK1,tail, originating due to transiently increased concentration of K+ in TATS. The density of IK1 itself was quantified by measuring the peak of IK1 in response to 400 ms voltage step to -120 mV prior to the following voltage ramp (Figure 2A). The ramp data were not used for any analysis and are presented here just for illustration purpose. Initial fast decline of IK1 from its peak value (black dot in Figure 2A) is due to depletion of t-tubular K+ but not due to change in membrane potential which follows later as a ramp. In order to minimize the variation of various parameters due to the highly variable size of cardiomyocytes the data were normalized to the cell size (cross-sectional, XY, area).
Patched cardiomyocytes were imaged with a MD500 microscope eyepiece camera and the AmScope 3.7 software (AmScope, Irvine, CA, United States). The cross-sectional area was calculated by manually outlining the cell border using ImageJ1.
Time Lapse Cell Width Measurements
Control and MβCD treated (4 mg/mL in C solution for 1 h) cardiomyocytes were plated onto a RC-20 perfusion chamber (Warner Instruments, Hamden, CT, United States). Only cells that settled in the center of the perfusion lane were imaged. Images were obtained at 10 s intervals. After a short (∼5 min) period of washout of C solution with Tyr, the perfusion was switched to 0.6 Na solution for 7 min.
The kinetics of solution exchange in the center of the perfusion chamber was measured by recording the change in fluorescence when switching from distilled water to distilled water containing 1:200 diluted 3 kDa dextran. After switching solutions, there is a delay of 3–5 s depending on the location in the bath chamber. The time course of fluorescence change was characterized by a time constant of <500 ms when fit with a single exponential function.
In response to hyposmotic stress, cardiomyocytes behave as nearly perfect osmometers and increase intracellular volume by expanding in width rather than in length (Drewnowska and Baumgarten, 1991). Furthermore, the cell depth to cell width ratio remains constant during hyposmotic swelling (Ogura et al., 2002), indicating that changes in the cell width serve as a good measure of cell volume changes. Experimentally, cardiomyocyte cell width was calculated from time-lapse image stacks using a combination of custom ImageJ and MATLAB scripts.
Di-8-ANEPPS Labeling of Cardiomyocytes
Stock solution of di-8-ANEPPS was prepared in DMSO at a concentration of 8.4 mM. The stock di-8-ANEPPS was mixed with 20% pluronic acid in 1:1 ratio and the resulting mixture was diluted (1:600) in C solution to prepare the working solution as previously described (Moench et al., 2013). Right before imaging, cardiomyocytes were incubated in the working solution for 15 min followed by washout with C solution without di-8-ANEPPS.
Dextran Trapping Assay
Dextran trapping assay was performed essentially as described in earlier publications (Moench et al., 2013; Uchida et al., 2016). In brief, 3 kDa dextran was added to a suspension of isolated ventricular myocytes during the swelling phase in hyposmotic 0.6 Na solution in order to fill t-tubules with this fluorescent marker. Cells were then returned to Tyr solution, still containing dextran, and finally extracellular dextran was washed out using normal Tyr. The cells were further washed and stored in C solution on ice prior to confocal imaging. Control myocytes were treated identically except that they were exposed to Tyr solution instead of 0.6 Na solution.
Confocal Imaging
Confocal imaging was performed in Microscopy and Image Analysis Laboratory (University of Michigan, Ann Arbor, MI, United States) on an Olympus FV-500 microscope using 60× 1.4 NA oil objective. Images of myocytes were manually outlined and mean intracellular fluorescence of trapped dextran per unit area calculated using ImageJ. Further data analysis, e.g., correction for background fluorescence, was performed in Microsoft Excel.
Skeletonization of TATS
Skeletonization of TATS was performed using approach similar to that described in the paper by Guo and Song (2014). Cardiomyocytes were labeled with di-8-ANEPPS and imaged on confocal microscope using 60× 1.4 NA oil objective and 68 nm pixel size. Image analysis was performed using ImageJ and custom macro to automate various steps in the procedure. Images were rotated to bring cardiomyocytes to the same (horizontal) orientation, smoothed using 1 pixel Gaussian Blur. Application of “Auto Threshold” using Huang method followed by “Analyze Particles” to exclude all objects smaller than the size of cardiomyocyte produces the outline of the cell. The outlined cell is then filled with black color and the resulting object uniformly reduced in size using “Erode” function thus creating a mask of cell interior (i.e., cell border excluded). Application of cell interior mask to the original smoothed image followed by “Auto Local Threshold” using Otsu method with radius = 20 provides a binary image of TATS. The images were then skeletonized using Skeletonize (2D/3D) plug-in. Separation of axial ant transverse t-tubules was carried out by application of appropriate morphological filters using MorphoLibJ plug-ins. Specifically, radial and axial t-tubule segments were isolated using 3 pixels long lines (structural element) at 90 and 45 degrees (for radial) and 0 and 135 degrees (for axial), respectively. The use of 45 and 135 degrees lines allows for assigning slanted t-tubular segments to either radial or axial group, and because of left/right symmetry of the cell the effects of application of asymmetrical lines (45 vs. 135 degrees) cancel out. It should be noted that segments smaller than 3 pixels will be lost during the above procedure but this does not affect the ratio of radial vs. axial t-tubules. The density of TATS was calculated as the % of t-tubule pixels relative to the total number of pixels in the cell mask.
Statistics
The data (mean ± standard error) in each experimental series are from at least two heart preparations. Statistical significance was determined using a one-way ANOVA with Bonferroni correction or two-sample t-test assuming equal or unequal variances (whichever is appropriate) and considered significant if p < 0.05. With some data one-way ANOVA could not be applied (e.g., due to large disparity in variance) and the data were then analyzed using two-sample t-tests. In figures, ∗, ∗∗, ∗∗∗ and #, ##, ### correspond to ANOVA or t-test with p-values of 0.05, 0.01, and 0.001, respectively.
Results
Effects of Membrane Cholesterol Modulation on the Overall Organization of TATS
We first tested the tolerance of cardiomyocytes to various concentrations of MβCD at room temperature (Supplementary Figure 1). There were no easily observable changes in TATS appearance at 3 or 6 mM MβCD but at 9 mM MβCD cells displayed clear deterioration of TATS and increased death rate.
In particular, we also found that the treatment of cardiomyocytes with 1 mM MβCD at 37°C (one of the commonly used conditions) is more damaging than that performed with 3 mM MβCD at room temperature (same exposure time; 1 h). Importantly, with 1 mM MβCD at 37°C only few cardiomyocytes survived the following standard hyposmotic detubulation while cells treated with 1 or 3 mM MβCD at room temperature did not show any overt changes in mortality following detubulation.
Accordingly, the following experiments were performed using 3 mM MβCD to minimize the detrimental consequences of cholesterol depletion in order to help unmask the underlying reasons behind its action.
As expected, the data in Figures 1A–C show that application of either 3 mM MβCD or WSC (5 mg/mL) has little effect on the overall appearance of TATS. In particular, the regularity of TATS estimated as the amplitude of the first harmonic of fluorescence spectra of membrane bound di-8-ANEPPS dye was not affected by MβCD and WSC (Figures 1B,C). However, quite small (∼3.4%) but statistically significant reduction in the wavelength of the first harmonic corresponding to the sarcomeric length of cardiomyocytes was observed in MβCD-treated cells (Figure 1D). TATS density was estimated in two different ways. First, Figure 1E shows that the intensity of intracellular (t-tubular) fluorescence of di-8-ANEPPS was not affected by MβCD but was increased by ∼16% in cardiomyocytes treated with WSC. Second, skeletonization of the TATS shows no significant effects of the drugs relative to control, although small, but statistically significant difference, can be observed when one would compare the effects of MβCD and WSC (Figure 1F). Additional analysis also shows that the balance between radial and axial t-tubules is not affected as well (Figure 1G). Consistent with relatively low effective concentration of MβCD the fluorescence of filipin, a cholesterol specific agent, was reduced only by ∼12% compared to that in control cardiomyocytes (see section “Discussion” and Supplementary Figure 2 on limitations of filipin-based assay). Importantly, the magnitude of reduction in filipin fluorescence observed with treatment using 3 mM MβCD for 1 h at RT (∼12% reduction; Figure 1H) is significantly smaller than the reduction observed with 1 mM MβCD for 1 h at 37°C (∼50% reduction; Supplementary Figure 2), suggesting that cholesterol extraction treatment in our study is milder than that used in similar studies (e.g., Gadeberg et al., 2017).
FIGURE 1. Effects of MβCD and WSC on TATS appearance (A; left column). Representative confocal images of control (Ctrl), 3 mM MβCD- and 5 mg/mL WSC-treated ventricular myocytes stained with di-8-ANEPPS. Scale bars: 20 μm. (B; middle column). Corresponding amplitude spectra of fluorescence (normalized to the amplitude of zero order component (A0). (C) Normalized amplitude (A1) and (D) wavelength (λ1) of the first harmonic of fluorescence. (E) Density of t-tubular fluorescence (TATSfluo) normalized to that in Ctrl myocytes. (F,G) Densities of various components of TATS derived from skeletonization of confocal images of cardiomyocytes labeled with di-8-ANEPPS. The density is calculated as percentage of pixels belonging to TATS skeleton. n = 20, 20, and 20 for Ctrl, MβCD- and WSC-treated cells, respectively. (H) Effects of MβCD and WSC on the intensity of filipin (cholesterol specific agent) staining. The data are normalized to that obtained in control cardiomyocytes. n = 67, 44 and 59 for Ctrl, MβCD- and WSC-treated cells, respectively.)
Overall, the general appearance of the TATS network remains largely unchanged following cholesterol depletion or supplementation at the indicated concentrations of cholesterol modifying agents.
Electrophysiological Characterization of the Effects of Membrane Cholesterol Modulation
Many ionic currents originate from ion channels concentrated in TATS. However, the IK1,tail current (IK1,tail) is a unique current that is in large degree dependent on the lumenal K+ within TATS (Clark et al., 2001; Cheng et al., 2011), which in turn are dependent on their fine geometrical structure (Uchida and Lopatin, 2018). Figures 2A–C explains the origin of IK1,tail and highlights some important features of this current related to TATS structure. In brief, during a prolonged depolarization, outward potassium currents cause potassium ions to flow into the t-tubular lumen and accumulate due to restricted diffusion to the extracellular space (Figure 2C, top). Upon repolarization, the accumulated potassium depletes by (1) flowing into the cell through IK1 channels manifesting as the IK1,tail and by (2) diffusing out of the t-tubule lumen to the extracellular space (Figure 2C, bottom). It has been shown that under conditions similar to that used in our experiments (e.g., holding membrane potential -75 mV) depletion of accumulated potassium by diffusion constitutes a large portion of total efflux and thus IK1,tail can serve as a quantitative measure of diffusional component of TATS (Clark et al., 2001).
FIGURE 2. Electrophysiological properties of cholesterol depleted and supplemented cardiomyocytes. (A) Cartoon trace of IK1 current in response to a slow ramp protocol. (B) Cartoon trace of outward K+ currents in response to a 400 ms step depolarization to +50 mV and inward IK1,tail in response to a step repolarization back to –75 mV. The amplitude of IK1,tail and its kinetics (characterized by time constants τ1 and τ2) are measures of diffusional properties of TATS. (C) Cartoon diagram depicting the origin of the IK1,tail. Top: outward potassium currents cause the accumulation of potassium within the t-tubular lumen during depolarization. Bottom: upon repolarization, the accumulated potassium can either flow into the cell through IK1 channels producing the IK1,tail or diffuse out of the t-tubular lumen to the extracellular space. (D–G) Quantification of the membrane capacitance normalized to the cell cross-sectional area (D), peak IK1 current normalized to the cell cross-sectional area (E), steady-state IK,end (F), and IK1,tail amplitude normalized to the steady-state IK,end current. n = 35 control, 40 MβCD, and 16 WSC-treated cells. (H,I) Quantification of the slow (H) and fast (I) time constants of IK1,tail decline. n = 37 control, 39 MβCD-, and 16 WSC-treated cells.
We first assessed the effects of cholesterol modulation on membrane capacitance. The averaged values of membrane capacitance (177.1 ± 9.3 pF, 168.2 ± 7.7 pF, and 173.5 ± 9.9 pF) and cell size (cross-sectional area; 3821.4 ± 153.6 μm2, 4005.1 ± 164.0 μm2, and 3388.5 ± 159.1 μm2) were not different (p = NS by ANOVA) between Ctrl, MβCD, and WSC groups, respectively. In order to better estimate potential changes in specific membrane capacitance, cell size should be taken into account. As shown in Figure 2D, treatment with MβCD significantly decreases normalized (to cell cross-sectional area) membrane capacitance. Conversely, treatment with WSC leads to an increase in normalized membrane capacitance although with no statistical significance (p = 0.06). Since cholesterol modulation significantly affects specific membrane capacitance, membrane currents were normalized to cell cross sectional area. Importantly, we found that IK1 is unaffected by cholesterol modulation (Figure 2E). In contrast, the outward potassium current at the end of the depolarizing pulse (IK,end) was significantly increased by MβCD treatment compared to that in control cells. Somewhat surprisingly, IK,end in WSC treated cells was essentially unchanged (Figure 2F).
Since the magnitude of potassium accumulation in TATS is dependent on the magnitude of outward potassium current, IK1,tail were normalized to IK,end (Cheng et al., 2011). Consistent with the data in Figure 1, which shows that the overall appearance of TATS is essentially unaffected by cholesterol modulation, no significant differences in the amplitude of the normalized IK1,tail were observed as well (Figure 2G). However, the kinetics of IK1,tail decline, which primarily reflects the rate of potassium diffusion out of the t-tubule lumen, was significantly affected (Figures 2H,I). Specifically, cells treated with MβCD display IK1,tail that decline faster than in control cells as both τ1 and τ2 of the two-exponential fit were decreased. Conversely, WSC treated cells display a slower decline in IK1,tail with both τ1 and τ2 significantly increased. This data suggests significant changes in fine (sub-microscopic) geometrical structure of TATS which may occur due to modulation of membrane cholesterol.
Effects of Modulation of Membrane Cholesterol on the Magnitude of Detubulation
Despite there being no or minor observable effects of cholesterol modifying agents on the overall architecture of TATS (Figure 1) the electrophysiological data above show significant, largely hidden but likely important changes in TATS structure strongly suggesting that there might be other consequences of cholesterol modulation. Specifically, we tested whether susceptibility of TATS to osmotic challenge, known to strongly affect the integrity of TATS (Moench et al., 2013), is affected by MβCD or WSC. In this regard, the magnitude of stress-induced sealing of TATS measured by the amount of trapped extracellular dextran serves as a useful and quantitative measure of TATS stability (Moench et al., 2013). Hyposmotic detubulation with 0.6 Na stress, however, leads to nearly complete detubulation in control cardiomyocytes. In order to maximize the putative observable effect, a milder osmotic stress (0.7 Na) was used to expand the dynamic range of dextran trapping. The data in Figures 3A,B, along with that in Supplementary Figures 3, 4, highlight one of the central findings of this work: depletion of cholesterol by MβCD significantly compromises the resistance of TATS to sealing in response to hyposmotic shock. Importantly, the effect is fully reversible by re-introducing cholesterol with WSC. Significant (>twofold) protective effect of membrane cholesterol is also confirmed in experiments using “full strength” hyposmotic shock with 0.6 Na solution where WSC is applied to normal cardiomyocytes (not previously treated with MβCD; Figures 3C,D).
FIGURE 3. Effects of MβCD and WSC on the amount of dextran trapped in sealed t-tubules. (A) Representative confocal images of ventricular myocytes (cell border zone indicated by white outline) highlighting the magnitude of t-tubular sealing in response to hyposmotic shock under various conditions. Detubulation with 0.6 Na solution (1) served as a reference, (2) no treatment, (3) MβCD for 1 h at RT, and (4) 3 mM MβCD for 1 h at RT followed by washout and application of 5 mg/mL WSC for 45 min at RT, all in response to shock using 0.7 Na solution. Scale bar: 20 μm. (B) Quantification of the data in panel (A). n = 68, 58, 55, and 40 cells, respectively. (C) Representative confocal images of control (top) and WSC-treated (5 mg/mL for 45 min) ventricular myocytes detubulated with standard 0.6 Na. (D) Quantification of dextran trapping of the data in pane (C). n = 35 cells each.
The data in Supplementary Figure 3 show that cardiomyocytes treated with 1 mM MβCD at 37°C display significantly greater dextran trapping than cardiomyocytes treated with 1 or 3 mM MβCD at room temperature for the same duration. It also follows from the data in Supplementary Figure 3 that treatment with 3 mM MβCD followed by detubulation with 0.7 Na solution places the response in a “dynamic” range allowing for observation of both increases and decreases in dextran trapping.
Notably, the magnitude of dextran trapping in MβCD treated cells inversely correlates with the remaining filipin fluorescence (Supplementary Figure 4), consistent with the notion that membrane cholesterol plays a critical role in TATS susceptibility to osmotic stress.
Consistent with the above findings, the magnitude of reduction in both normalized membrane capacitance and normalized IK1,tail amplitude (due to sealing of TATS) in response to 0.7 Na hyposmotic stress is more pronounced in MβCD treated cells compared to that in control cells (Figure 4). Overall, the above results suggest that MβCD treated cells are more susceptible to t-tubular remodeling in response to hyposmotic stress.
FIGURE 4. Electrophysiological effects of MβCD after 0.7 Na detubulation. (A) Quantification of cell area normalized membrane capacitance of control and MβCD treated cells after detubulation with 0.7 Na. n = 13 cells each. (B) Quantification of normalized IK1,tail in control and MβCD treated cells after detubulation with 0.7 Na. n = 12 and 11 cells, respectively.
Cell Swelling Is Unaffected by Cholesterol Depletion
Previous work suggested that cholesterol depletion may accelerate cardiomyocyte swelling in response to osmotic stress (Kozera et al., 2009) and thus may affect the magnitude of changes in TATS following osmotic detubulation. Accordingly, the kinetics of cell swelling in 0.6 Na solution was measured in control and MβCD treated cardiomyocytes. As shown in Figure 5, in both control and MβCD treated cardiomyocytes cell width increases in a nearly identical manner upon perfusion with 0.6 Na solution. The change in cell width can be fitted with a single exponential function and the time constant characterizing the kinetics of cell swelling is not significantly different between those measured in control and MβCD treated cells. Furthermore, there is no significant difference in the magnitude of cell swelling in 0.6 Na solution. These data suggest that cholesterol depletion has no effect on hyposmotic swelling and argues against an osmotic effect underlying the differences in dextran trapping between control and MβCD treated cells.
FIGURE 5. Effect of MβCD treatment on hyposmotic cell swelling. Traces of average changes in relative cell width of control (black) and MβCD treated (gray) cardiomyocytes as a function of time. Vertical line denotes the start of perfusion with hyposmotic 0.6 Na solution. Inset: (A) Quantification of the kinetics of cell swelling. A single exponential function was fitted to the time course of cell width during 0.6 Na perfusion and the time constant was averaged for each group (n = 7 and 9 for control and MβCD treated cells, respectively). (B) Quantification of the maximum change in cell width.
Discussion
Role of Cholesterol in T-Tubular Membranes
The importance of cholesterol in t-tubular membranes is highlighted in recent papers. In particular, cholesterol depletion using MβCD was shown to affect TATS integrity in a concentration and time-dependent manner (Zhu et al., 2016). This effect of MβCD seems controversial as more recent work demonstrated that cholesterol depletion does not appear to disrupt TATS organization (Gadeberg et al., 2017). However, a closer look at the specific experimental conditions in the mentioned studies can reconcile the findings. Specifically, in the study by Gadeberg et al. (2017) threefold lower concentration of MβCD was used compared to the lowest concentration employed in the study by Zhu et al (1 mM vs. 3 mM) (Zhu et al., 2016). The data in our study confirm this view and essentially eliminate that formal discrepancy. Although we performed experiments at RT (vs. 37°C in mentioned papers) a different range of MβCD concentrations was tested.
The data in Figure 1 and Supplementary Figure 1 show that, indeed, at RT and 3 mM MβCD concentration essentially no or only minor changes in the overall structure of TATS can be observed while at 6 mM and higher concentrations the disruptive effect of cholesterol depletion becomes apparent. It should be noted here, that the amount of membrane cholesterol measured using filipin approach is likely overestimated, and thus the effect of MβCD on membrane cholesterol is underestimated (Figure 1H), likely due to membrane permeability of filipin leading to staining intracellular pools of cholesterol (see Supplementary Figure 2 for further detail). Alternatively, the amount of membrane cholesterol may be overestimated in WSC-treated cardiomyocytes likely due to membrane binding of WSC aggregates or micelles.
The two parameters which we found to be affected by modulation of membrane cholesterol were sarcomeric length (Figure 1D; ∼3.4%; MβCD vs. control) and the density of TATS estimated using either di-8-ANEPPS labeling (Figure 1E; ∼19%; control vs. WSC) or skeletonizing algorithm (Figure 1F; ∼10%; MβCD vs. WSC). The change in the sarcomeric length is very small, and the observation became possible in part due to relatively high precision of the measurements employing Fourier transformation. The effect is likely linked to changes in the concentration of resting intracellular Ca2+, however, it was impractical to proceed with further investigation of this phenomenon keeping in mind the magnitude of the changes. Measurements of TATS density using di-8-ANEPPS can be significantly confounded by unpredictable effects of cholesterol modulation on the binding of the dye (likely more than due to minor changes in the resting membrane potential of cardiomyocytes), and therefore the data should be interpreted with caution. In contrast, skeletonization approach using local thresholding algorithms is not sensitive to the variation in the intensity of di-8-ANEPPS. Changes in TATS density measured with this approach are very small and could only be reliably detected by comparing two extreme treatments (MβCD vs. WSC). We speculate that the observed effects are not due to true changes in TATS density but are rather apparent (i.e., algorithm dependent) changes due to submicroscopic transformations in the morphology of individual t-tubules. Overall, quantitative data show that the appearance of TATS is essentially unaffected by treatment with MβCD and WSC.
Importantly, however, one has to make a clear distinction between the overall TATS integrity and its stability. In this regard, one of the central findings of our study is that although the TATS integrity appears essentially unaffected following cholesterol modulation with 3 mM MβCD, the stability of the TATS in response to stress is greatly compromised. In particular, the data in Figure 3 demonstrate that cholesterol depletion promotes dextran trapping following hyposmotic detubulation, while cholesterol supplementation protects against osmotic detubulation. These data suggest that cholesterol plays a significant role in stabilizing t-tubules against osmotic stress. Importantly, as shown in Figure 5, modulation of membrane cholesterol did not affect the magnitude and kinetics of osmotic swelling suggesting that the observed results cannot be simply explained by changes in mechanical forces caused by osmotic stress. Furthermore, this data may suggest that Zhu et al. observed the combined effects of MβCD treatment and long-term culture, a condition known to cause t-tubule loss (Louch et al., 2004; Zhu et al., 2016; Guo et al., 2018).
Effects of Cholesterol Modulation on Electrophysiological Properties of Cardiomyocytes
It is commonly expected that cholesterol modulation would have significant effects on membrane capacitance. However, the available data are quite controversial, even in studies using artificial membranes. For example, Ohki (1969) found that capacitance of phosphatidylcholine bilayers increases with addition of cholesterol while in other studies an opposite effect was observed (Naumowicz et al., 2005; Budvytyte et al., 2013), and the effect of cholesterol on membrane capacitance may even show biphasic relationship in lipid monolayers (Pasek et al., 2008). Assuming that increasing membrane cholesterol monotonically decreases membrane capacitance and that cholesterol modulation does not alter the total membrane surface area, one would expect that treatment with MβCD would result in increased membrane capacitance. In contrast, Gadeberg et al. (2017) have recently found that in mouse cardiomyocytes depletion of membrane cholesterol does not alter the membrane capacitance. In even greater contrast, we observe the opposite effect of cholesterol on membrane capacitance (Figure 2). As mentioned earlier, differences in experimental conditions may underlie quantitative differences between the studies. In particular, we used higher concentration of MβCD but lower (RT) temperature which altogether likely led to a stronger modulation of membrane cholesterol. Overall, the data highlight the complexity of cholesterol action on membrane capacitance, which likely involves cholesterol-dependent caveolae that contribute significantly to membrane area (Levin and Page, 1980). In this regard, EM images of cardiomyocytes treated with MβCD display significantly fewer caveolae than control cells (Kozera et al., 2009). However, it remains unknown whether these caveolae are internalized thus decreasing total membrane area or whether the proteins supporting the caveolae structure merely disassemble thus retaining the lipids in the surface membrane. If these caveolae are removed from the total surface membrane, then the membrane capacitance would be expected to decrease, consistent with our findings. Unfortunately, the effect of cholesterol supplementation on the number of caveolae remains unknown and, therefore, it is unclear whether additional cholesterol increases the total surface membrane area as suggested by the membrane capacitance measurements. It should be noted that direct correlation between the intensity of di-8-ANEPPS fluorescence and membrane area cannot be made since di-8-ANEPPS fluorescence is strongly dependent on the presence of membrane cholesterol (Gross et al., 1994).
In contrast to none or relatively small effects of cholesterol modulation on membrane capacitance, there are numerous reports on significant roles of this lipid in the activity of various ion channels (Dart, 2010). In this study we focused on two types of potassium currents highly useful in determining the diffusional properties of TATS: cardiac outward rectifier and inward rectifier potassium currents (Figure 2). We have previously shown that the magnitude of K+ accumulation in TATS maximizes at about 400 ms after membrane depolarization and is primarily determined by IK,end which is, in turn, is carried primarily by several members of Kv1 and Kv2 subfamilies of voltage-gated K+ channels underlying IK,slow (Nerbonne, 2000). Consistent with previous reports we find that depletion of membrane cholesterol leads to significant increase in late IK,end (Balse et al., 2009), although cholesterol enrichment with WSC did not have any significant effect (Figure 2). Accordingly, because of sensitivity of IK,end to changes in membrane cholesterol, for a meaningful interpretation of the data the magnitude of IK1,tail should be normalized to the magnitude of IK,end (leading to INK1,tail; Figure 2).
As mentioned in Figure 2, IK1,tail dissipates due to two fluxes: movement of potassium into the cell through IK1 channels carried by Kir2 subfamily (Anumonwo and Lopatin, 2010) and diffusion of potassium out of t-tubular lumen. Surprisingly, we found that in mouse cardiomyocytes IK1 is essentially insensitive to cholesterol modulation (Figure 2E). This is different from the effects observed in endothelial cells endogenously expressing inward-rectifier K+ channels (Romanenko et al., 2002) or in cells exogenously expressing Kir2.1 channels (Romanenko et al., 2004) where treatment with MβCD or cholesterol loaded MβCD causes robust changes in the current density. Due to insensitivity of cardiac IK1 to membrane cholesterol there was no need to further normalize the INK1,tail to IK1 amplitude.
Consistent with no effects of cholesterol modulation on the overall structure of TATS we found no significant changes in the amplitude of INK1,tail upon treatment with MβCD or WSC (Figure 2). It should be noted, however, that the interpretation of the INK1,tail amplitude with regard to TATS structure is not that straightforward, in particular, because it depends on both the density of outward K+ current and t-tubular diffusion. In contrast, kinetics of IK1,tail (Figure 2C) depends in a large degree on the diffusional properties of TATS (Clark et al., 2001). In this regard, we recently demonstrated that t-tubular dilations and constrictions, which generally can be viewed as submicroscopic irregularities of t-tubular lumen shape, have a significant impact on the diffusion of molecules within TATS (Uchida and Lopatin, 2018). T-tubule dilations are common features of normal cardiac TATS (Savio-Galimberti et al., 2008; Pinali et al., 2013) and a recent study has reported that t-tubular dilations are associated with regions where Cav3 and RyR2 colocalize (Wong et al., 2013). The loss of these sub-microscopic structures resulting in more uniformly shaped t-tubules would be expected to accelerate the diffusion rate while introduction of more constrictions/dilations will restrict diffusion and slow diffusion rate. The data in Figures 2H,I show that modulation of membrane cholesterol leads to significant changes in diffusional properties of TATS, consistent with the notion that disruption of cholesterol domains [e.g., t-tubular caveolae (Levin and Page, 1980; Burton et al., 2017)] by MβCD results in reduced irregularity of t-tubule diameters and thus faster diffusion of potassium while supplementation with additional cholesterol by WSC does the opposite. An alternative mechanism may involve the effect of cholesterol on t-tubular cBIN1-microfolds which were suggested to contribute significantly to diffusional properties of TATS (Hong et al., 2014). Quantitative analysis of the changes in cBIN1-microfolds, however, would require optical super-resolution or electron microscopy imaging in future projects.
Overall, the results of this study show that in mouse ventricular myocytes modulation of membrane cholesterol leads to significant changes in susceptibility of TATS to acute hyposmotic stress. These findings suggest that long term effects of membrane cholesterol on the integrity of TATS may be explained, at least in part, by its effects of the stability of TATS.
Author Contributions
AN, KU, and AL designed the research, performed the experiments, analyzed the data, and wrote the manuscript. IM performed the experiments and analyzed the data. All the authors contributed to the final version of the manuscript.
Funding
This work was supported by the National Institutes of Health National Heart, Lung, and Blood Institute (Grant HL127023; AL) and AHA Predoctoral Fellowship (17PRE33350049; KU).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank Yasmine Elghoul for assistance with isolation of cardiomyocytes.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.01516/full#supplementary-material
Footnotes
References
Anumonwo, J. M., and Lopatin, A. N. (2010). Cardiac strong inward rectifier potassium channels. J. Mol. Cell. Cardiol. 48, 45–54. doi: 10.1016/j.yjmcc.2009.08.013
Balse, E., El-Haou, S., Dillanian, G., Dauphin, A., Eldstrom, J., Fedida, D., et al. (2009). Cholesterol modulates the recruitment of Kv1.5 channels from Rab11-associated recycling endosome in native atrial myocytes. Proc. Natl. Acad. Sci. U.S.A. 106, 14681–14686. doi: 10.1073/pnas.0902809106
Budvytyte, R., Mickevicius, M., Vanderah, D. J., Heinrich, F., and Valincius, G. (2013). Modification of tethered bilayers by phospholipid exchange with vesicles. Langmuir 29, 4320–4327. doi: 10.1021/la304613a
Burton, R. A. B., Rog-Zielinska, E. A., Corbett, A. D., Peyronnet, R., Bodi, I., Fink, M., et al. (2017). Caveolae in rabbit ventricular myocytes: distribution and dynamic diminution after cell isolation. Biophys. J. 113, 1047–1059. doi: 10.1016/j.bpj.2017.07.026
Cheng, L., Wang, F., and Lopatin, A. N. (2011). Metabolic stress in isolated mouse ventricular myocytes leads to remodeling of t-tubules. Am. J. Physiol. Heart Circ. Physiol. 301, H1984–H1995. doi: 10.1152/ajpheart.00304.2011
Clark, R. B., Tremblay, A., Melnyk, P., Allen, B. G., Giles, W. R., and Fiset, C. (2001). T-tubule localization of the inward-rectifier K(+) channel in mouse ventricular myocytes: a role in K(+) accumulation. J. Physiol. 537, 979–992. doi: 10.1111/j.1469-7793.2001.00979.x
Coronel, R. (2017). The pro- or antiarrhythmic actions of polyunsaturated fatty acids and of cholesterol. Pharmacol. Ther. 176, 40–47. doi: 10.1016/j.pharmthera.2017.02.004
Dart, C. (2010). Lipid microdomains and the regulation of ion channel function. J. Physiol. 588, 3169–3178. doi: 10.1113/jphysiol.2010.191585
Draeger, A., Monastyrskaya, K., Mohaupt, M., Hoppeler, H., Savolainen, H., Allemann, C., et al. (2006). Statin therapy induces ultrastructural damage in skeletal muscle in patients without myalgia. J. Pathol. 210, 94–102. doi: 10.1002/path.2018
Drewnowska, K., and Baumgarten, C. M. (1991). Regulation of cellular volume in rabbit ventricular myocytes: bumetanide, chlorothiazide, and ouabain. Am. J. Physiol. 260, C122–C131. doi: 10.1152/ajpcell.1991.260.1.C122
Gadeberg, H. C., Kong, C. H. T., Bryant, S. M., James, A. F., and Orchard, C. H. (2017). Cholesterol depletion does not alter the capacitance or Ca handling of the surface or t-tubule membranes in mouse ventricular myocytes. Physiol. Rep. 5:e13500. doi: 10.14814/phy2.13500
Gross, E., Bedlack, R. S. Jr., and Loew, L. M. (1994). Dual-wavelength ratiometric fluorescence measurement of the membrane dipole potential. Biophys. J. 67, 208–216. doi: 10.1016/S0006-3495(94)80471-0
Guo, A., Chen, R., Wang, Y., Huang, C. K., Chen, B., Kutschke, W., et al. (2018). Transient activation of PKC results in long-lasting detrimental effects on systolic [Ca(2+)]i in cardiomyocytes by altering actin cytoskeletal dynamics and T-tubule integrity. J. Mol. Cell. Cardiol. 115, 104–114. doi: 10.1016/j.yjmcc.2018.01.003
Guo, A., and Song, L. S. (2014). AutoTT: automated detection and analysis of t-tubule architecture in cardiomyocytes. Biophys. J. 106, 2729–2736. doi: 10.1016/j.bpj.2014.05.013
Guo, A., Zhang, C., Wei, S., Chen, B., and Song, L. S. (2013). Emerging mechanisms of T-tubule remodelling in heart failure. Cardiovasc. Res. 98, 204–215. doi: 10.1093/cvr/cvt020
Hidalgo, C. (1985). Lipid phase of transverse tubule membranes from skeletal muscle. an electron paramagnetic resonance study. Biophys. J. 47, 757–764. doi: 10.1016/S0006-3495(85)83978-3
Hong, T., Yang, H., Zhang, S. S., Cho, H. C., Kalashnikova, M., Sun, B., et al. (2014). Cardiac BIN1 folds T-tubule membrane, controlling ion flux and limiting arrhythmia. Nat. Med. 20, 624–632. doi: 10.1038/nm.3543
Kawai, M., Hussain, M., and Orchard, C. H. (1999). Excitation-contraction coupling in rat ventricular myocytes after formamide-induced detubulation. Am. J. Physiol. 277, H603–H609. doi: 10.1152/ajpheart.1999.277.2.H603
Kline, C. F., and Mohler, P. J. (2013). Evolving form to fit function: cardiomyocyte intercalated disc and transverse-tubule membranes. Curr. Top. Membr. 72, 121–158. doi: 10.1016/B978-0-12-417027-8.00004-0
Kozera, L., White, E., and Calaghan, S. (2009). Caveolae act as membrane reserves which limit mechanosensitive I(Cl,swell) channel activation during swelling in the rat ventricular myocyte. PLoS One 4:e8312. doi: 10.1371/journal.pone.0008312
Levin, K. R., and Page, E. (1980). Quantitative studies on plasmalemmal folds and caveolae of rabbit ventricular myocardial cells. Circ. Res. 46, 244–255. doi: 10.1161/01.RES.46.2.244
Louch, W. E., Bito, V., Heinzel, F. R., Macianskiene, R., Vanhaecke, J., Flameng, W., et al. (2004). Reduced synchrony of Ca2 + release with loss of T-tubules-a comparison to Ca2 + release in human failing cardiomyocytes. Cardiovasc. Res. 62, 63–73. doi: 10.1016/j.cardiores.2003.12.031
Moench, I., and Lopatin, A. N. (2014). Ca homeostasis in sealed t-tubules of mouse ventricular myocytes. J. Mol. Cell. Cardiol. 72, 374–383. doi: 10.1016/j.yjmcc.2014.04.011
Moench, I., Meekhof, K. E., Cheng, L. F., and Lopatin, A. N. (2013). Resolution of hypo-osmotic stress in isolated mouse ventricular myocytes causes sealing of t-tubules. Exp. Physiol. 98, 1164–1177. doi: 10.1113/expphysiol.2013.072470
Mohaupt, M. G., Karas, R. H., Babiychuk, E. B., Sanchez-Freire, V., Monastyrskaya, K., Iyer, L., et al. (2009). Association between statin-associated myopathy and skeletal muscle damage. CMAJ 181, E11–E18. doi: 10.1503/cmaj.081785
Naumowicz, M., Petelska, A. D., and Figaszewski, Z. A. (2005). Impedance analysis of phosphatidylcholine-cholesterol system in bilayer lipid membranes. Electrochim. Acta 50, 2155–2161. doi: 10.1016/j.electacta.2004.09.023
Nerbonne, J. M. (2000). Molecular basis of functional voltage-gated K + channel diversity in the mammalian myocardium. J. Physiol. 525, 285–298. doi: 10.1111/j.1469-7793.2000.t01-1-00285.x
Ogura, T., Imanishi, S., and Shibamoto, T. (2002). Osmometric and water-transporting properties of guinea pig cardiac myocytes. Jpn. J. Physiol. 52, 333–342. doi: 10.2170/jjphysiol.52.333
Ohki, S. (1969). The electrical capacitance of phospholipid membranes. Biophys. J. 9, 1195–1205. doi: 10.1016/S0006-3495(69)86445-3
Page, E. (1978). Quantitative ultrastructural analysis in cardiac membrane physiology. Am. J. Physiol. 235, C147–C158. doi: 10.1152/ajpcell.1978.235.5.C147
Pasek, M., Brette, F., Nelson, A., Pearce, C., Qaiser, A., Christe, G., et al. (2008). Quantification of t-tubule area and protein distribution in rat cardiac ventricular myocytes. Prog. Biophys. Mol. Biol. 96, 244–257. doi: 10.1016/j.pbiomolbio.2007.07.016
Pinali, C., Bennett, H., Davenport, J. B., Trafford, A. W., and Kitmitto, A. (2013). 3-D reconstruction of the cardiac sarcoplasmic reticulum reveals a continuous network linking t-tubules: this organization is perturbed in heart failure. Circ. Res. 113, 1219–1230. doi: 10.1161/CIRCRESAHA.113.301348
Romanenko, V., Fang, Y., Byfield, F., Travis, A. J., Vandenberg, C., Rothblat, G. H., et al. (2004). Cholesterol sensitivity and lipid raft targeting of kir 2.1 channels. Biophys. J. 87, 3850–3861. doi: 10.1529/biophysj.104.043273
Romanenko, V. G., Rothblat, G. H., and Levitan, I. (2002). Modulation of endothelial inward-rectifier K + current by optical isomers of cholesterol. Biophys. J. 83, 3211–3222. doi: 10.1016/S0006-3495(02)75323-X
Rothberg, K. G., Ying, Y. S., Kamen, B. A., and Anderson, R. G. (1990). Cholesterol controls the clustering of the glycophospholipid-anchored membrane receptor for 5-methyltetrahydrofolate. J. Cell Biol. 111, 2931–2938. doi: 10.1083/jcb.111.6.2931
Savio-Galimberti, E., Frank, J., Inoue, M., Goldhaber, J. I., Cannell, M. B., Bridge, J. H., et al. (2008). Novel features of the rabbit transverse tubular system revealed by quantitative analysis of three-dimensional reconstructions from confocal images. Biophys. J. 95, 2053–2062. doi: 10.1529/biophysj.108.130617
Simionescu, N., Lupu, F., and Simionescu, M. (1983). Rings of membrane sterols surround the openings of vesicles and fenestrae, in capillary endothelium. J. Cell Biol. 97, 1592–1600. doi: 10.1083/jcb.97.5.1592
Sumnicht, G. E., and Sabbadini, R. A. (1982). Lipid composition of transverse tubular membranes from normal and dystrophic skeletal muscle. Arch. Biochem. Biophys. 215, 628–637. doi: 10.1016/0003-9861(82)90124-2
Uchida, K., and Lopatin, A. N. (2018). Diffusional and electrical properties of t-tubules are governed by their constrictions and dilations. Biophys. J. 114, 437–449. doi: 10.1016/j.bpj.2017.11.3742
Uchida, K., Moench, I., Tamkus, G., and Lopatin, A. N. (2016). Small membrane permeable molecules protect against osmotically induced sealing of t-tubules in mouse ventricular myocytes. Am. J. Physiol. Heart Circ. Physiol. 311, H229–H238. doi: 10.1152/ajpheart.00836.2015
Wong, J., Baddeley, D., Bushong, E. A., Yu, Z., Ellisman, M. H., Hoshijima, M., et al. (2013). Nanoscale distribution of ryanodine receptors and caveolin-3 in mouse ventricular myocytes: dilation of T-tubules near junctions. Biophys. J. 104, L22–L24. doi: 10.1016/j.bpj.2013.02.059
Keywords: mouse ventricular myocytes, t-tubule, cholesterol, osmotic stress, potassium currents
Citation: Nikouee A, Uchida K, Moench I and Lopatin AN (2018) Cholesterol Protects Against Acute Stress-Induced T-Tubule Remodeling in Mouse Ventricular Myocytes. Front. Physiol. 9:1516. doi: 10.3389/fphys.2018.01516
Received: 30 May 2018; Accepted: 09 October 2018;
Published: 12 November 2018.
Edited by:
Alexey V. Glukhov, University of Wisconsin System, United StatesReviewed by:
Andrew F. James, University of Bristol, United KingdomTingTing Hong, Cedars-Sinai Medical Center, United States
Copyright © 2018 Nikouee, Uchida, Moench and Lopatin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Anatoli N. Lopatin, alopatin@umich.edu
†These authors have contributed equally to this work