• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Neurodynamics analysis of cochlear hair cell activity

    2020-03-27 03:43:14WeifengRongRubinWngJinhiZhngWnzengKong

    Weifeng Rong, Rubin Wng,b,*, Jinhi Zhng, Wnzeng Kong

    a School of Computer Science, Hangzhou Dianzi University, Hangzhou 310018, China

    b Institute for Cognitive Neurodynamics, East China University of Science and Technology, Shanghai 200237, China

    Keywords:Hair cells Sound frequencies Membrane potential Power Outer hair cell amplification

    ABSTRACT There have been many studies on the effect of cochlea basal membrane movement on the resolution of different frequencies and intensities. However, these studies did not take into account the influence of power and energy consumption of the hair cells in the process of the electromotility movement, as well as the neurodynamic mechanism that produced this effect. This makes previous studies unable to fully clarify the function of outer hair cells (OHCs) and the mechanism of sound amplification. To this end, we introduce the gate conductance characteristics of the hair cells in the mechanical process of increasing frequency selectivity. The research finds that the low attenuation of OHCs membrane potential and the high gain in OHC power and energy consumption caused that OHC amplification is driven by electromotility. The research results show that the amplification of the OHCs is driven by low attenuation of membrane potential and high gain of power and energy consumption. This conclusion profoundly reveals the physiological mechanism of the electromotility movement.

    1 Introduction

    The auditory nervous system converts external mechanical signals into neural signals, exhibiting ingenious signal processing mechanism. The frequency selection in the basilar membrane of the inner ears reflects its unique spatial encoding,wherein the sound frequencies are transformed into resonance at different positions in the basilar membrane [1-4]. In addition,hair cells, which are the receptors in the inner ear, transform the mechanical signals of basilar membrane oscillations into neural signals, which are critical for the transduction of auditory information [1, 2, 5].

    The mechanical vibration of sound waves is then transferred into the inner ear after collection in the outer ear and amplification in the middle ear, leading to perilymph vibration in the inner ear [2, 6]. This vibration in turn causes basilar membrane oscillation, resulting in the movement of the organ of Corti as well as tectorial membrane. As the basilar and tectorial membranes move up and down around different axes, it causes shearing motion of the reticular lamina on the organ of Corti and the underlying surface of the tectorial membrane [1, 7, 8]. This shearing motion leads to deflection of stereocilia of the outer hair cells(OHCs) that connect with tectorial membrane in a specific direction. On the other hand, the stereocilia present in the inner hair cells do not connect with the tectorial membrane, and its deflection is caused by the flow of endolymph between the reticular lamina and the tectorial membrane. The deflections of the stereocilia modify the opening of mechanoelectrical transduction(MET) channels located on the stereocilia.

    The stereocilia that exist on the top of the hair cells are infiltrated in the endolymph. However, the body of hair cells is infiltrated into the perilymph [9]. The endolymph in the scala media contained high concentration of K+and low concentration of Na+, which is similar to that in the intracellular fluid. In contrast,the perilymph in the scala vestibule and scala tympani comprises of low concentration of K+and high concentration of Na+[1, 3, 10]. If the potential of perilymph is 0 mV, then the potential of the endolymph in the scala media is as high as 80 mV, which is known as the endocochlear potential (EP). After the MET channels on stereocilia of hair cells are opened, the potential difference formed by the EP and negative potentials in the hair cells led to the flow of K+in the endolymph into hair cells to form receptor potential in the lateral membrane of hair cells [2, 11, 12].

    The MET of hair cells occurs on the stereocilia, and the MET channels are present on the tip of the stereocilia. During the basilar membrane oscillations, the basilar membrane moves towards the scala vestibule, leading to shearing motion between the hair cells and the tectorial membrane. This subsequently causes deflection of stereocilia towards the exciting direction,and turns near the tall edge. This process increases the opening of MET channels, where more ions (mainly K+) enter the cell and depolarize the hair cells. If the basilar membrane moves towards the scala tympani, then the stereocilia deflects towards the inhibiting directions, deviating from the taller edge. This subsequently decreases the opening of MET channels, and repolarize the hair cells [9, 13, 14].

    Therefore, basilar membrane oscillation causes deflection of stereocilia towards different directions, modulating the current flowing into the hair cells. The current flows into the hair cells through the apical membrane of the hair cells from the endolymph and flows out to the perilymph from the lateral membrane, and this is known as the transducer current. The voltage produced by the current on the lateral membrane is termed as the receptor potential [3, 15].

    Through MET of hair cells, a mechanical signal of basilar membrane oscillation is then converted into a neural signal in hair cells. Hence, the auditory information can be transmitted through the auditory nerve to the auditory center of the brain,producing auditory perception [1, 16, 17].

    The electromotility of the OHCs was found by Brownell et al.[18]. The electromotility of the OHCs converts the membrane potential to change the length of the OHCs, and the membrane potential depolarization causes the OHCs to shorten, while the membrane potential hyperpolarization leads to the elongation of OHCs [19]. The amplification of sound is driven by the electromotility of OHCs in the cochlea [18-20], which then depends on the lateral membrane potential of OHCs [20]. This amplification effect can be explored by studying the lateral membrane potential in OHCs.

    In order to deeply understand the neurodynamic mechanism of hair cell activity during the cochlear basilar membrane oscillation, a hair cell model based on Hodgkin-Huxley equation was established in this study. Instead of MET channel conductance proposed by Howard and Hudspeth [21], MET channel conductance proposed by Mountain and Cody [22] was introduced . Numerical simulations were performed to analyze the neurodynamics of hair cell membrane potential, power and energy consumption by stimulating them at different frequencies.

    The novelty of this paper is to calculate the power and energy consumption of hair cells during the process of movement,and to reveal the physiological mechanism of the electromotility,which has not been reported so far.

    2 Mathematical model

    2.1 Basilar membrane model

    A model that maps sound frequencies to the best frequency positions of the basilar membrane was required.

    In the present study, a basilar membrane model was utilized[23], and its simplified form was as follows

    wherein, ρ represents the density of lymph fluid, h indicates the height of the cochlear duct, x is the distance from the stapes,S(x)is the stiffness coefficient of the basilar membrane per unit,M(x)represents the mass of the basilar membrane per unit,β(x)=is the damping of basilar membrane, ζ is the damping coefficient per unit, α is the motility factor of the OHCs,representing the outer hair cell force of basilar membrane displacement per unit, γ is the feed-forward/feed-backward ratio, representing the relative effects of the OHC force, d indicates the tilt distance, δ is the basilar membrane displacement, t1and t2represent the delay in action of OHC, and B represents a constant coefficient.

    2.2 Hair cell model

    The Hodgkin-Huxley model [24] is considered important in many neuron models, and is regarding as one of the best models to match the experiment results. The model was described as follows

    However, the Hodgkin-Huxley model considers only the potassium current (IK), sodium current (INa), leakage current (IL)and external stimulus current (Iext) that do not correspond to the actual current in hair cells.

    In order to make the Hodgkin-Huxley model suitable for hair cells, Catacuzzeno et al. [25, 26] analyzed the ion current during the release process of hair cells and established a hair cell model based on the Hodgkin-Huxley model. After that, Neiman et al.[27, 28] improved the model by adding mechanically-gated ion current, which coincides with the physiological mechanism of the hair cells. However, they did not consider the effects of different sound frequencies on hair cells and did not calculate the power and energy consumption of hair cells, making it impossible to understand the neurodynamic mechanism of hair cell activity during the basilar membrane oscillation.

    This study established a hair cell model at the best frequency position in the basilar membrane based on the Hodgkin-Huxley model. The outward current included inward-rectifier potassium current (IK), delayed rectifier potassium current (IDRK), and calcium-activated potassium current (IKCa). The inward current included cation h-type current (Ih), voltage-gated calcium current (ICa), leakage current (IL), and mechanically-gated ion current (IMET). The specific model was described as follows

    wherein, C was the membrane capacitor of hair cell and V was the lateral membrane potential.

    Inward-rectifier potassium current (IK) included one rapid and one slow activation channel [26]. In the following equations,EK=-95 mV and gK=29 nS.

    Cation h-type cationic current (Ih) had three independent activation channels [25]. In the following equations, Eh=-45 mV,and gh=2.2 nS.

    Delayed rectifier potassium current (IDRK) had two independent activation channels [26]. In the following equations,PDRK=2.4×10-14L/s, [K]in=112 mM, [K]out=2 mM, F was Faraday constant, R was a universal gas constant and T=295.15 K.

    Voltage-gated calcium current (ICa) had three independent activation channels [25]. In the following equations, ECa=42.5 mV, gCa=1.2 nS.

    Calcium-activated potassium current (IKCa) had one activation channel. In the following equations, EKCa=-70 mV, and F was Faraday constant. For conductance, gKCa_IHC=50 nS (IHC means inner hair cell), while gKCa_OHCshowed variation [29, 30],with a logarithmic linear increase from the cochlear apex to the cochlear base. In this paper, gKCa_OHC=5×101.68bmx+1nS, wherein,“bmx” was the basilar membrane proportional distance from the cochlear apex.

    In the leakage current (IL), EL=0 mV, gL=0.3 nS

    IMETwas different from the other currents present on the lateral membrane. It flew into the hair cell from the top of the hair cells via the MET channels on the stereocilia (the main ion was potassium ion). The apical membrane potential of hair cells was subtracted from the endocochlear potential (VEP=80 mV). The K+concentration of the endolymph was similar to that of the intracellular fluid, and therefore EMET=0 mV [1]. Thus, the MET current (IMET) was described as follows

    Herein, the method proposed by Mountain and Cody [22]was used for calculating the MET conductance (gMET), given by

    wherein, gmaxwas the maximum of MET conductance, and x was the deflecting distance of the stereocilia. For gmax, gmax_IHC=11.2 nS, while gmax_OHCwas a variable, with a logarithmic linear increase from the cochlear apex to the cochlear base [29]. In this paper, gmax_OHC=1.12×100.7806bmx+1nS.

    The movement of the stereocilia was sinusoidal [10], and Eq.(13) was thus set up as below

    wherein A was a constant coefficient, and f was sound frequency.

    Since 10% channels on the stereocilia were open when it was at a resting position without deflection (x = 0) [14, 31], the value of A could be calculated.

    Through the above hair cell model, the hair cell membrane potential at different sound frequencies can be obtained.

    2.3 Hair cell power model

    Hair cell power could be calculated according to the Hodgkin-Huxley model as the hair cell model was based on the Hodgkin-Huxley model [32, 33]. The specific method was described as follows.

    The energy contained in a hair cell was described as the sum of the potential energy stored in the membrane capacitor and the stored energy of the equivalent battery generated by Nernst potential of ions, which was presented by the equation below

    The latter seven terms on the right side of Eq. (14) were difficult to calculate directly. But the rate of change for each term was equal to the current through the equivalent battery multiplied by its electromotive force (Nernst potential) [31]. So, the total power P(t) was the derivative of H(t) over time

    Equation (16) was obtained by putting Eq. (4) into Eq. (15)

    wherein, the left side was the total power, while various items on the right side indicated the power of various ion channels,respectively.

    By the above method, the hair cell power at different sound frequencies can be obtained. Then, by integrating Eq. (16) during the stimulation time, the hair cell energy consumption at different sound frequencies during the period can be obtained.

    3 Results and discussion

    3.1 Position response of basilar membrane at different sound frequencies

    In Eq. (1), k was a function of x. Therefore, k could be solved by iteration, which was later placed into Eq. (2) to solve the basilar membrane displacement δ. Thus, the plot of the position response of basilar membrane at any sound frequency could be obtained (Fig. 1).

    The peak in Fig. 1 indicated the best frequency positions of the basilar membrane at different frequencies. This suggested that the resonance position of the characteristic frequency was consistent with that of the physiological characteristics, so that the best frequency position of the basilar membrane stimulated by high-frequency sound was close to the base of the cochlea,while that stimulated by the low-frequency sound was distal from the base of the cochlea.

    Based on Fig. 1, the characteristic frequency and best frequency position could be fitted as shown in Fig. 2.

    Fig. 1. Position response of basilar membrane at different sound frequencies (SPL means sound pressure level)

    Fig. 2. Best frequency positions of characteristic frequencies

    In Fig. 2, the horizontal axis represented the equal distance of the resonance position of characteristic frequency to the top of the cochlea, while the vertical axis was the voice frequency,and the curve expression was as follows

    Figure 2 was highly consistent with the physiological data [34].

    3.2 Hair cell membrane potential at different sound frequencies

    To study the effect of sound frequencies on hair cell membrane potential during the basilar membrane oscillation, the simulation experiments in Figs. 3 and 4 were conducted in this section. Among them, the sound stimuli ranged from 10 ms to 50 ms, and the sound frequencies were divided into three levels f =500 Hz, 1000 Hz, and 2000 Hz.

    Fig. 3. Membrane potentials of IHC (from top to bottom) at 500 Hz,1000 Hz, and 2000 Hz

    Fig. 4. Membrane potentials of OHC (from top to bottom) at 500 Hz, 1000 Hz, and 2000 Hz

    In Figs. 3 and 4, as the sound frequencies increased, the peak values of the hair cell membrane potential were gradually decreased, while the valley values were gradually increased. The overall result suggested that the range of the hair cell membrane potential was gradually lowered. This was consistent with the physiological experimental data [35], and reproduced the detailed characteristics of the hair cell membrane potential during the basilar membrane oscillation, confirming the neurodynamic properties of hair cells. This indicated that the hair cell model established at different sound frequencies during the basilar membrane oscillation in this paper was reasonable.

    Additionally, with the increase of sound frequencies and the change of hair cell membrane potential, the relationship between sound frequency and hair cell membrane potential was extensively studied. The sound frequencies were in the range of 0.1-20 kHz, and the range of hair cell membrane potential was ΔV. The simulation results were shown in Fig. 5.

    In Fig. 5, the solid blue line and the green dotted line represents the range of membrane potential of IHC and OHC at the best frequency position of the basilar membrane at different sound frequencies, respectively. When the sound frequency was increased, both IHC and OHC membrane potentials were attenuated. Before 1000 Hz, OHC was decayed faster than IHC, but both showed less attenuation, and after 1000 Hz, IHC was attenuated faster than OHC.

    According to Fig. 5, the attenuation of IHC membrane potential was calculated to be -16.3516 dB/decade, and the attenuation of OHC membrane potential was lower, which was-8.7562 dB/decade. The low attenuation of OHC membrane potential suggested that the cochlear amplification was caused due to the electromotility of OHCs [29]. This calculation result was consistent with the experimental data in Ref. [29].

    3.3 Hair cell membrane potential at different basilar membrane regions

    Fig. 5. Range of hair cell membrane potential at different frequencies

    The amplitude of basilar membrane was increased as the basilar membrane regions approached the best frequency position [1, 3], which was also confirmed in Fig. 1. At the position that is not the best frequency position, a decrease in the amplitude of basilar membrane led to a decrease in the deflection of stereocilia, thereby indicating a decrease in the A value.

    If suppose, the best frequency position at a specific frequency was x (x was the distance to the stapes), the amplitude of basilar membrane could be calculated as y by the basilar membrane model. If a position that is not the best frequency position x1, the amplitude of the basilar membrane at this position was y1 that could be calculated by the same method. Thus, the deflection of stereocilia at this position wasthereby providing the hair cell membrane potential.

    Since the traveling wave of basilar membrane attenuated rapidly after the peak, and stopped vibration at a position slightly away from that [1, 2], the hair cell membrane potential of the basilar membrane before the peak was only analyzed.

    To study the effect of basilar membrane regions on hair cell membrane potential during the basilar membrane oscillation,the hair cell membrane potentials at 5 mm and 15 mm to the staples, as well as the best frequency position were plotted.Among them, the sound stimuli ranged from 10 ms to 50 ms, and the sound frequencies were 500 Hz.

    In Fig. 6, the peak value of IHC membrane potential was gradually increased as the basilar membrane region approached the best frequency position. However, the valley value showed a very slight change. The changes in the peak value and the valley value led to a gradual increase in the range of hair cell membrane potential. And from Fig. 7, it has been deduced that OHC have the same pattern.

    To analyze the relationship between basilar membrane regions and hair cell membrane potential on it more intuitively,the peak value, valley value, and range at different basilar membrane regions at 500 Hz was plotted.

    The horizontal axis of Fig. 8 indicated that the ratio of basilar membrane amplitude (y1) to the basilar membrane amplitude at the best frequency position (y). y 1/y represented represents the distance of basilar membrane region to the best frequency position (0 meant the farthest and 1 meant the nearest).

    Fig. 6. Membrane potentials of IHC at a 5 mm, b 15 mm, and c best frequency positon at 500 Hz

    Figure 8 showed that the peak value of IHC membrane potential was gradually increased as the basilar membrane region approached the best frequency position; however, the valley value remained almost the same. The overall result suggested that the range of IHC membrane potential was gradually increased. This indicated that the effect of the basilar membrane region on hair cell membrane potential was mainly at the peak value. And in case of OHC, it was similar.

    3.4 Hair cell power and energy consumption at different sound frequencies

    To study the effect of sound frequencies on hair cell power and energy consumption during the basilar membrane oscillation, simulation experiments (Figs. 9 and 10) were conducted in this section. Among them, the sound stimuli ranged from 10 ms to 50 ms, and the sound stimulation frequencies were divided into three levels: f = 500 Hz, 1000 Hz, and 2000 Hz.

    Fig. 7. Membrane potentials of OHC at a 5 mm, b 15 mm, and c best frequency positon at 500 Hz

    Fig. 8. Peak value, valley value, and range (from top to bottom) of IHC membrane potential at different basilar membrane regions

    For IHC (Fig. 9), the peak value, valley value, and the range of the power were slightly changed with increasing sound frequencies. For OHC (Fig. 10), the peak value was gradually increased,while the valley value was gradually decreased with the increasing sound frequencies, leading to a gradual increase in the range. These results showed that the sound frequencies had significant effect on the power of OHC rather than IHC. The numerical simulation results of hair cells in Figs. 10 and 11 are consistent with the experimental data of hair cell physiology under different sound frequencies in Ref. [34].

    Fig. 9. Power of IHC (from top to bottom) at 500 Hz, 1000 Hz, and 2000 Hz.

    Fig. 10. Power of OHC (from top to bottom) at 500 Hz, 1000 Hz, and 2000 Hz (1 pW = 1×10-12 W)

    In addition, the increasing sound frequencies and changing hair cell power were studied to analyze the relationship between sound frequency and hair cell power extensively. The sound frequencies were in the range of 0.1-20 kHz, and the range of hair cell power was Δ P. The simulation results were shown in Fig. 11.

    In Fig. 11, the solid blue line and the green dotted line represents the range of power of IHC and OHC at the best frequency position of the basilar membrane at different sound frequencies, respectively. In Fig. 11, if the sound frequency was increased, the IHC power was attenuated, but the amplitude remained small. OHC power demonstrated a large gain. According to Fig. 11, the attenuation of IHC power was calculated to be-0.0084 dB/decade, and the gain of OHC power was 6.9592 dB/decade.

    Then, the energy consumption of different frequencies during 40 ms stimulation time was calculated by integrating the hair cell power.

    In Fig. 12, the solid blue line and the green dotted line represents the energy consumption of IHC and OHC at the best frequency position of the basilar membrane at different sound frequencies, respectively. According to Fig. 12, when the sound frequency was increased, the energy consumption of IHC was increased, but the amplitude remained small. The energy consumption of OHC was significantly increased. According to Fig.12, the gain of IHC energy consumption was calculated to be 0.1216 dB/decade, and the gain of OHC energy consumption was 6.5259 dB/decade.

    Fig. 11. Range of hair cell power at different frequencies

    Fig. 12. Hair cell energy consumption at different frequencies (1 pJ= 1×10-12 J)

    Compared to IHC, the power and energy consumption of OHC showed a significantly high gain with the increase of sound frequencies, reflecting the electromotility of OHCs. When the sound frequencies are increased, then the oscillation frequency of the basilar membrane was also increased, resulting in increased number of lengthening and shortening OHCs, and the power and energy consumed are naturally increased. The high gain of OHC power and energy consumption also supported, in a certain sense, that the OHC amplification was driven by electromotility.

    4 Conclusion

    To profoundly understand the neurodynamic mechanism of hair cell activity and function during the cochlear basilar membrane oscillation, a hair cell model based on Hodgkin-Huxley equation was established. Numerical simulations were performed to analyze the neurodynamics of hair cell membrane potential, power and energy consumption by stimulating at different frequencies. The results showed that in the range of 0.1-20 kHz, the attenuation of the membrane potential of OHCs was lower than that of IHCs, which was only -8.7562 dB/decade. The effect of the basilar membrane region on the hair cell membrane potential is mainly present at the peak value. The gain in the power and energy consumption of OHCs was much higher than those of IHCs, which were 6.9592 dB/decade and 6.5259 dB/decade, respectively.

    The results showed that the IHCs and OHCs have different functions under different frequencies. The functional differences between them indicated that the IHCs exist as receptive cells, while the OHCs play an important role in sound amplification.

    The low attenuation of OHC membrane potential and the high gain in OHC power and energy consumption supported the view that OHC amplification is driven by electromotility. The differences in the membrane potentials of IHC and OHC reflected different functions, where IHCs work as sensory cells, and OHCs play a role in sound amplification [9, 29]. This novel conclusion profoundly illustrated the physiological mechanism of electrokinesis, which is of great significance for improving the operating principle of the artificial cochlea in future, and provides an important guiding significance for improving the working quality of the artificial cochlea in the field of manufacturing.

    Acknowledgments

    This study was funded by the National Natural Science Foundation of China (Grants 11232005, 11472104, 11872180,61633010, and 61473110).

    九九热线精品视视频播放| 高清在线视频一区二区三区 | 欧美一区二区精品小视频在线| 亚洲成人中文字幕在线播放| 最近视频中文字幕2019在线8| 亚洲av免费高清在线观看| 在线播放国产精品三级| av福利片在线观看| 草草在线视频免费看| 国产黄色小视频在线观看| 国产片特级美女逼逼视频| 尤物成人国产欧美一区二区三区| 精品不卡国产一区二区三区| 久久亚洲精品不卡| 六月丁香七月| 色尼玛亚洲综合影院| 亚洲在久久综合| 亚洲精品久久久久久婷婷小说 | av黄色大香蕉| 在线观看66精品国产| 色噜噜av男人的天堂激情| 国产黄色小视频在线观看| 成年免费大片在线观看| 日韩国内少妇激情av| 国产精品国产三级专区第一集| 亚洲av成人av| 熟女人妻精品中文字幕| 99久久精品一区二区三区| 国内精品美女久久久久久| 我的女老师完整版在线观看| 日本一二三区视频观看| 啦啦啦韩国在线观看视频| 精品熟女少妇av免费看| 色综合亚洲欧美另类图片| 国内揄拍国产精品人妻在线| 亚洲婷婷狠狠爱综合网| 精品人妻熟女av久视频| 午夜久久久久精精品| 国产精品1区2区在线观看.| 在线观看一区二区三区| av国产久精品久网站免费入址| 亚洲人成网站在线播| 亚洲最大成人中文| 热99在线观看视频| 午夜福利视频1000在线观看| 黄色欧美视频在线观看| 日韩av不卡免费在线播放| 丰满乱子伦码专区| 亚洲欧美一区二区三区国产| 看片在线看免费视频| 一区二区三区高清视频在线| 亚洲最大成人中文| 日本与韩国留学比较| 亚洲内射少妇av| 简卡轻食公司| 国产探花在线观看一区二区| 国产极品天堂在线| a级毛片免费高清观看在线播放| 乱系列少妇在线播放| av专区在线播放| 久久久久久九九精品二区国产| 亚洲内射少妇av| 简卡轻食公司| 国产免费男女视频| 日本wwww免费看| 国产成人a∨麻豆精品| 高清毛片免费看| h日本视频在线播放| 精品久久久久久电影网 | 一区二区三区高清视频在线| 欧美性感艳星| 高清毛片免费看| av福利片在线观看| 少妇的逼好多水| 国产精品一区www在线观看| av专区在线播放| 又粗又爽又猛毛片免费看| 欧美性猛交黑人性爽| 少妇熟女欧美另类| 色5月婷婷丁香| 黄片无遮挡物在线观看| 99久久无色码亚洲精品果冻| 午夜精品一区二区三区免费看| 大话2 男鬼变身卡| 激情 狠狠 欧美| 欧美丝袜亚洲另类| 少妇高潮的动态图| 国产精品人妻久久久影院| 好男人在线观看高清免费视频| 精品久久久噜噜| 一区二区三区高清视频在线| 午夜亚洲福利在线播放| 亚洲国产成人一精品久久久| 夜夜看夜夜爽夜夜摸| 国产欧美另类精品又又久久亚洲欧美| 最近的中文字幕免费完整| 国产伦精品一区二区三区四那| 男女啪啪激烈高潮av片| 一级毛片久久久久久久久女| 如何舔出高潮| av.在线天堂| 久久人人爽人人爽人人片va| 久久99热6这里只有精品| 亚洲在久久综合| 男女下面进入的视频免费午夜| 一个人观看的视频www高清免费观看| 亚洲图色成人| 中文亚洲av片在线观看爽| 亚洲精品乱久久久久久| 在线免费观看不下载黄p国产| 日韩在线高清观看一区二区三区| 男女下面进入的视频免费午夜| 亚洲国产最新在线播放| 国产免费又黄又爽又色| 在线观看av片永久免费下载| 校园人妻丝袜中文字幕| 久久久久久久久久久丰满| 在线天堂最新版资源| 亚洲国产欧美在线一区| 国产亚洲一区二区精品| 看十八女毛片水多多多| 97超视频在线观看视频| 精品人妻熟女av久视频| 岛国毛片在线播放| 亚洲av福利一区| 国产麻豆成人av免费视频| 在线播放无遮挡| 日韩中字成人| videossex国产| 欧美一级a爱片免费观看看| 91午夜精品亚洲一区二区三区| 亚洲av电影在线观看一区二区三区 | 久久欧美精品欧美久久欧美| av视频在线观看入口| 国产精品久久电影中文字幕| 国产极品精品免费视频能看的| 国产黄a三级三级三级人| 精品无人区乱码1区二区| 欧美激情国产日韩精品一区| 我的女老师完整版在线观看| 久久99热6这里只有精品| 99久久九九国产精品国产免费| 七月丁香在线播放| 99热这里只有是精品在线观看| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 中文字幕制服av| 在线天堂最新版资源| 麻豆国产97在线/欧美| 色视频www国产| 久久久久久久久久久免费av| 国产一区二区亚洲精品在线观看| 免费黄网站久久成人精品| 亚洲国产高清在线一区二区三| 亚洲内射少妇av| 婷婷色av中文字幕| 人体艺术视频欧美日本| 亚洲av电影在线观看一区二区三区 | 99国产精品一区二区蜜桃av| 免费观看精品视频网站| 欧美日韩精品成人综合77777| 最后的刺客免费高清国语| 成人欧美大片| 久久亚洲精品不卡| 免费在线观看成人毛片| 人体艺术视频欧美日本| 欧美性感艳星| 夜夜爽夜夜爽视频| 久久精品久久久久久噜噜老黄 | 蜜臀久久99精品久久宅男| 欧美zozozo另类| 久久久久久久亚洲中文字幕| 国产真实伦视频高清在线观看| 最后的刺客免费高清国语| 成人午夜精彩视频在线观看| 久久99热这里只有精品18| 久99久视频精品免费| 国产探花极品一区二区| 长腿黑丝高跟| 亚洲精品乱码久久久v下载方式| 1000部很黄的大片| 少妇裸体淫交视频免费看高清| 色尼玛亚洲综合影院| ponron亚洲| 日日干狠狠操夜夜爽| 亚洲在线自拍视频| 免费看a级黄色片| av国产久精品久网站免费入址| 国产伦精品一区二区三区四那| 欧美日韩精品成人综合77777| 午夜福利视频1000在线观看| 久久久久免费精品人妻一区二区| 亚洲国产欧洲综合997久久,| 又爽又黄无遮挡网站| 欧美一区二区精品小视频在线| 免费观看人在逋| 国产又黄又爽又无遮挡在线| 成人三级黄色视频| 亚洲国产色片| 在线a可以看的网站| 国国产精品蜜臀av免费| 级片在线观看| 久久久午夜欧美精品| 六月丁香七月| 国产探花极品一区二区| 久久久久久伊人网av| 热99re8久久精品国产| 看十八女毛片水多多多| 少妇裸体淫交视频免费看高清| 丝袜喷水一区| 亚洲欧美日韩无卡精品| 成人av在线播放网站| 日本熟妇午夜| 色哟哟·www| 成年av动漫网址| 欧美精品一区二区大全| 69人妻影院| 中文字幕久久专区| 久久亚洲精品不卡| 免费看av在线观看网站| 欧美精品一区二区大全| 国模一区二区三区四区视频| 亚洲国产日韩欧美精品在线观看| 中文精品一卡2卡3卡4更新| 成人二区视频| 男人舔奶头视频| av国产久精品久网站免费入址| 在线播放无遮挡| 国产精品国产三级国产av玫瑰| 美女高潮的动态| 天堂中文最新版在线下载 | 久久精品熟女亚洲av麻豆精品 | 超碰av人人做人人爽久久| 极品教师在线视频| av天堂中文字幕网| 国产伦精品一区二区三区视频9| 一区二区三区乱码不卡18| 联通29元200g的流量卡| 国内精品美女久久久久久| 久久久久久久午夜电影| 国产精品国产三级国产专区5o | 一级毛片电影观看 | 大话2 男鬼变身卡| 舔av片在线| 成人亚洲精品av一区二区| 成人午夜精彩视频在线观看| 亚洲精品乱码久久久久久按摩| 男女国产视频网站| av线在线观看网站| 日本午夜av视频| 国产精品美女特级片免费视频播放器| 你懂的网址亚洲精品在线观看 | 最近中文字幕2019免费版| 欧美性感艳星| kizo精华| 成年女人永久免费观看视频| 午夜精品一区二区三区免费看| 日韩精品有码人妻一区| 麻豆精品久久久久久蜜桃| 久久精品91蜜桃| 成人午夜高清在线视频| 久久久久免费精品人妻一区二区| 91久久精品国产一区二区三区| 又爽又黄a免费视频| 深夜a级毛片| 成人性生交大片免费视频hd| 亚洲第一区二区三区不卡| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费| 床上黄色一级片| 久久久久久久久久久丰满| 精品人妻视频免费看| 99在线视频只有这里精品首页| 三级经典国产精品| 日韩av在线大香蕉| 搞女人的毛片| 欧美一区二区国产精品久久精品| 成人综合一区亚洲| 日韩大片免费观看网站 | 蜜臀久久99精品久久宅男| 国产精品嫩草影院av在线观看| 国产精品国产高清国产av| 中文字幕熟女人妻在线| 小蜜桃在线观看免费完整版高清| 高清av免费在线| 国产大屁股一区二区在线视频| 精华霜和精华液先用哪个| 在线免费观看不下载黄p国产| 成年女人看的毛片在线观看| 美女cb高潮喷水在线观看| 国产精品一区二区三区四区免费观看| 日日干狠狠操夜夜爽| 在线播放无遮挡| 有码 亚洲区| 国产一区有黄有色的免费视频 | 床上黄色一级片| 最近最新中文字幕大全电影3| 免费播放大片免费观看视频在线观看 | 国产精品电影一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产又黄又爽又无遮挡在线| 99久久中文字幕三级久久日本| 欧美另类亚洲清纯唯美| av国产免费在线观看| 国产精品国产三级国产专区5o | 欧美成人免费av一区二区三区| 免费看a级黄色片| 日韩 亚洲 欧美在线| 日本黄色视频三级网站网址| 男人的好看免费观看在线视频| 网址你懂的国产日韩在线| 久久99热这里只频精品6学生 | 亚洲内射少妇av| 男女那种视频在线观看| 精品久久久久久久久亚洲| 亚洲四区av| 美女高潮的动态| 三级国产精品片| 欧美+日韩+精品| 亚洲综合色惰| 亚洲成人久久爱视频| 中国美白少妇内射xxxbb| 日韩制服骚丝袜av| 国产精品美女特级片免费视频播放器| 免费av观看视频| 亚洲精品日韩在线中文字幕| 美女xxoo啪啪120秒动态图| 久久久久久久国产电影| 欧美三级亚洲精品| or卡值多少钱| 日本三级黄在线观看| 男的添女的下面高潮视频| 欧美一区二区亚洲| 又爽又黄a免费视频| 国产三级中文精品| 毛片一级片免费看久久久久| 天堂av国产一区二区熟女人妻| 美女大奶头视频| 最近手机中文字幕大全| 日韩在线高清观看一区二区三区| 亚洲伊人久久精品综合 | videos熟女内射| 免费看a级黄色片| 亚洲综合精品二区| 欧美日韩综合久久久久久| 日韩欧美三级三区| 人妻系列 视频| 成年版毛片免费区| 少妇的逼好多水| 中文字幕亚洲精品专区| av在线蜜桃| 欧美zozozo另类| 天堂av国产一区二区熟女人妻| 国内揄拍国产精品人妻在线| 国产人妻一区二区三区在| 看免费成人av毛片| АⅤ资源中文在线天堂| 91精品国产九色| 麻豆久久精品国产亚洲av| 日韩成人av中文字幕在线观看| 国产高清三级在线| 深夜a级毛片| 简卡轻食公司| 人体艺术视频欧美日本| 日韩成人av中文字幕在线观看| 啦啦啦观看免费观看视频高清| 亚洲一级一片aⅴ在线观看| 国产av一区在线观看免费| 精品久久久久久电影网 | 国产亚洲午夜精品一区二区久久 | 国产一区有黄有色的免费视频 | 青青草视频在线视频观看| kizo精华| 久久99热这里只频精品6学生 | 舔av片在线| 国产伦理片在线播放av一区| 日本av手机在线免费观看| 欧美+日韩+精品| 日韩成人av中文字幕在线观看| 搞女人的毛片| 国产精品久久久久久av不卡| 小说图片视频综合网站| 国产精华一区二区三区| 纵有疾风起免费观看全集完整版 | 国产v大片淫在线免费观看| 哪个播放器可以免费观看大片| 亚洲精品456在线播放app| 深夜a级毛片| 晚上一个人看的免费电影| 看非洲黑人一级黄片| 亚洲人成网站在线播| 99久久精品热视频| 男女那种视频在线观看| 亚洲欧洲日产国产| 日日摸夜夜添夜夜添av毛片| 一夜夜www| 日韩,欧美,国产一区二区三区 | 身体一侧抽搐| 日韩 亚洲 欧美在线| 婷婷色麻豆天堂久久 | 国产精品一及| 成人一区二区视频在线观看| 综合色av麻豆| 免费看光身美女| 少妇猛男粗大的猛烈进出视频 | 欧美激情久久久久久爽电影| 丰满少妇做爰视频| 国内精品美女久久久久久| 国产精品一区www在线观看| 日韩av在线免费看完整版不卡| 久久久久性生活片| 精华霜和精华液先用哪个| 九九爱精品视频在线观看| 亚洲av男天堂| 久久精品熟女亚洲av麻豆精品 | 亚洲高清免费不卡视频| 欧美日韩一区二区视频在线观看视频在线 | 边亲边吃奶的免费视频| 精品无人区乱码1区二区| 国产午夜精品久久久久久一区二区三区| 女人被狂操c到高潮| 日韩中字成人| 丝袜喷水一区| 久久久久久大精品| 内射极品少妇av片p| 99国产精品一区二区蜜桃av| 中文亚洲av片在线观看爽| 国产片特级美女逼逼视频| 男人舔奶头视频| 国产单亲对白刺激| 久久这里只有精品中国| 国产高清三级在线| 国产亚洲最大av| 色5月婷婷丁香| 欧美激情在线99| 免费av不卡在线播放| 春色校园在线视频观看| 日韩一本色道免费dvd| 麻豆精品久久久久久蜜桃| 日本午夜av视频| 国产av一区在线观看免费| 亚洲18禁久久av| 精品久久久久久久人妻蜜臀av| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| 国产精品久久视频播放| 边亲边吃奶的免费视频| 免费看av在线观看网站| 久久久久久久久大av| 国产探花在线观看一区二区| 国产精品女同一区二区软件| 高清毛片免费看| 美女国产视频在线观看| 夜夜爽夜夜爽视频| 18禁在线播放成人免费| 老司机影院毛片| 久久久国产成人免费| 久久99精品国语久久久| 亚洲不卡免费看| 国产成年人精品一区二区| av专区在线播放| 综合色丁香网| 欧美最新免费一区二区三区| 亚洲av免费高清在线观看| 欧美又色又爽又黄视频| 免费黄网站久久成人精品| 爱豆传媒免费全集在线观看| 国产精品人妻久久久久久| 久久6这里有精品| 色网站视频免费| 好男人在线观看高清免费视频| 汤姆久久久久久久影院中文字幕 | 国产成人精品一,二区| 国产亚洲午夜精品一区二区久久 | 真实男女啪啪啪动态图| 国产精品爽爽va在线观看网站| 久久久久久久久久成人| 欧美成人一区二区免费高清观看| 国产爱豆传媒在线观看| 91精品国产九色| 国产精品,欧美在线| 一级毛片aaaaaa免费看小| 少妇裸体淫交视频免费看高清| 在线a可以看的网站| 能在线免费观看的黄片| eeuss影院久久| 丝袜喷水一区| 18禁在线播放成人免费| 国产老妇女一区| 熟女人妻精品中文字幕| 欧美bdsm另类| 免费观看精品视频网站| 国产黄片视频在线免费观看| 亚洲欧美精品自产自拍| 成人二区视频| 国产精品人妻久久久影院| 久久亚洲国产成人精品v| 亚洲婷婷狠狠爱综合网| 男女啪啪激烈高潮av片| 欧美成人午夜免费资源| 成年女人永久免费观看视频| 少妇高潮的动态图| 能在线免费观看的黄片| 国产真实伦视频高清在线观看| 国产久久久一区二区三区| 精品久久久久久久久av| 亚洲最大成人中文| 大香蕉97超碰在线| 久久99蜜桃精品久久| 直男gayav资源| 久久久精品94久久精品| 网址你懂的国产日韩在线| 国产精品一区www在线观看| 国产真实伦视频高清在线观看| 在线播放国产精品三级| 噜噜噜噜噜久久久久久91| 三级经典国产精品| 欧美xxxx性猛交bbbb| 久久久久精品久久久久真实原创| 一级毛片aaaaaa免费看小| 日日干狠狠操夜夜爽| 三级国产精品欧美在线观看| 亚洲一区高清亚洲精品| 美女黄网站色视频| 中文资源天堂在线| 国产精华一区二区三区| 少妇人妻一区二区三区视频| 午夜精品一区二区三区免费看| 午夜亚洲福利在线播放| 三级国产精品欧美在线观看| 老司机福利观看| 99在线人妻在线中文字幕| 欧美激情久久久久久爽电影| 国产精品一区www在线观看| 久久精品久久久久久噜噜老黄 | 黄色配什么色好看| 久久久久网色| 少妇被粗大猛烈的视频| 国产激情偷乱视频一区二区| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠久久av| 中文字幕免费在线视频6| 99在线人妻在线中文字幕| 欧美激情久久久久久爽电影| 国产精品美女特级片免费视频播放器| 国产精品av视频在线免费观看| 小蜜桃在线观看免费完整版高清| 一区二区三区免费毛片| 久久99精品国语久久久| 国产精品麻豆人妻色哟哟久久 | 女的被弄到高潮叫床怎么办| 中文精品一卡2卡3卡4更新| 久久久久国产网址| 亚洲国产精品成人久久小说| 精品国产三级普通话版| 国产午夜精品一二区理论片| 一区二区三区乱码不卡18| 国产午夜精品久久久久久一区二区三区| 国产一区有黄有色的免费视频 | 国产精品嫩草影院av在线观看| 精华霜和精华液先用哪个| 欧美潮喷喷水| 好男人视频免费观看在线| 久久99精品国语久久久| 国产日韩欧美在线精品| 99国产精品一区二区蜜桃av| 亚洲精品亚洲一区二区| 乱码一卡2卡4卡精品| 午夜亚洲福利在线播放| 欧美精品一区二区大全| 色尼玛亚洲综合影院| 一级二级三级毛片免费看| 亚洲第一区二区三区不卡| 欧美成人精品欧美一级黄| 亚洲国产欧美人成| 99视频精品全部免费 在线| 又爽又黄无遮挡网站| 变态另类丝袜制服| 三级经典国产精品| 久久99热6这里只有精品| 国产精品一区二区三区四区久久| 精品酒店卫生间| 韩国av在线不卡| 亚洲自偷自拍三级| 欧美不卡视频在线免费观看| 久久精品夜色国产| 97超视频在线观看视频| 欧美激情久久久久久爽电影| 99热6这里只有精品| 日韩三级伦理在线观看| 成人高潮视频无遮挡免费网站| 国产精品一区www在线观看| 99热网站在线观看| 欧美人与善性xxx| 在线免费十八禁| 美女国产视频在线观看| 不卡视频在线观看欧美| 国产麻豆成人av免费视频| 国产精品国产三级国产av玫瑰| 黄色配什么色好看| 日韩一区二区三区影片| 欧美精品国产亚洲| .国产精品久久| 婷婷色av中文字幕| 久久久久久国产a免费观看| 九九在线视频观看精品| 女的被弄到高潮叫床怎么办| 狂野欧美激情性xxxx在线观看| 国产一级毛片七仙女欲春2| 高清毛片免费看| 青青草视频在线视频观看| 国内少妇人妻偷人精品xxx网站| 中文在线观看免费www的网站| 一区二区三区四区激情视频| 亚洲精品久久久久久婷婷小说 |