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

    Effect of cognitive training on brain dynamics

    2024-02-29 09:17:50GuiyangLv呂貴陽TianyongXu徐天勇FeiyanChen陳飛燕PingZhu朱萍MiaoWang王淼andGuoguangHe何國光
    Chinese Physics B 2024年2期
    關(guān)鍵詞:王淼國光貴陽

    Guiyang Lv(呂貴陽), Tianyong Xu(徐天勇), Feiyan Chen(陳飛燕),Ping Zhu(朱萍), Miao Wang(王淼), and Guoguang He(何國光)

    School of Physics,Zhejiang University,Hangzhou 310027,China

    Keywords: brian dynamics,functional brain networks,cognitive training,abacus-based mental calculation

    1.Introduction

    The human brain is composed of billions of neurons in various levels of complex structures.Revealing the relationship between complex structures and functions of the brain is one of the current focuses of neuroscience.Su′arezet al.built artificial neural networks with the network structures of brain regions.[1]The networks have shown stronger performance in the memory coding task than any other networks, which indicates that the network structures are essential for the brain to perform tasks efficiently.The work by Kanakaet al.revealed that a randomly connected recurrent neural network(RNN)can carry out memory work by modifying certain connections of the network and a signal is transmitted through the interaction between synapse and external input.It means that the neural signal sequence may emerge through learning from the unstructured network architectures.[2]Therefore,flexible and variable connections of brain regions are necessary for the brain to realize various cognitive tasks.At present,most of the whole-brain dynamics researches have focused on the synchronization or chimera state among brain regions by the interaction between the coupled nodes with the structural brain network.[3–5]It is lack to understand how the structures in brain networks affect the overall performance of the brain,and how the information contained in the network structures is expressed through dynamic activities.

    In terms of cognitive experiments, it is known that the human brain is highly plastic.[6]Cognitive training can lead to the modification of the brain’s functional excitation mode and achieve better task performance.Abacus-based mental calculation(AMC)is one of typical cognitive training method.[7–9]Previous studies have shown that the resting-state functional connections between the brain regions of the subjects would be changed by long-term abacus mental arithmetic training.[10,11]The experiments demonstrated that AMC training can enhance the subjects’ numerical processing efficiency[12]and improve their working memory.[13]Based on task-state fMRI,it was found that the brain excitation patterns and functional structures between brain regions for the experimental group with AMC training are different from those for the control group.[14]The overall performance of the brain is determined by the network structures and node dynamics.[15,16]Up to date,most AMC researches are focused on revealing the correlation between cognitive abilities and functional structures.The influence of functional connections on brain dynamics is rarely reported.How functional connectivity affects the overall dynamics of the brain remains unclear.

    Brain activity not only is based on structural networks,but also depends on factors such as external stimuli and neurotransmitter levels.[5,17]It is found that some brain regions are not directly connected in structure, but they show functionally relations,[18]which means that the functional networks are the embodiment of brain interaction after various factors are coupled.Therefore, we utilize the functional brain networks as the connection structures of brain regions in our work.The resting-state fMRI data of the AMC group and the control group are collected to establish the functional networks and dynamic behaviors of functional networks from the AMC group and the control group are investigated in this work.Indeed, small networks based on simple oscillators are also capable of reproducing the collective dynamics of the brain and are expected to be used for disease-related simulations.[19]For example,whole brain dynamics simulations based on the FHN model are able to demonstrate spontaneous brain activity.[20]Like the FHN model,the Hindmarsh–Rose(HR)model has a wealth of dynamic characteristics from single neuron to neuron clusters and is often used in brain dynamics modeling.[21]The neural network model based on neuron models can reproduce EEG-like oscillatory signals.[22]Mitchellet al.used the same method to reproduce the process of seizures in mice.[23]Chouzouriset al.and Kanget al.used the HR model and FHN model to replace a brain region to study the synchronization and chimera state of brain.[24,25]Ramlowet al.used the FHN model to study unihemispheric sleep in the human brain.[26]A similar approach was used in the early diagnosis of brain tumors.[27]Antonopouloset al.verified that the brain evolves by maximizing its information flow capacity through HR neural networks.[28]Focusing on dynamic modeling based on brain networks,Santoset al.discovered chimera-like states in cat brain networks with the HR model.Our previous study has demonstrated that the dynamics of the resting-state brain can be simulated using the HR model and functional networks.[29]Ansarinasabet al.used the HR model and functional networks based on EEG signals to identify chimera states in both attention deficit hyperactivity disorder (ADHD) patients and control groups.In the presence of emotional stimuli, the ADHD group exhibited premature synchronization,leading to abnormal brain emotional processing function.[30]

    In view of wide applications of the HR model, it is feasible to replace one brain region with the HR model in the brain networks.A complete brain dynamics model can be constructed on the basis of functional connections of the brain and the HR model.Our previous work proved that this model can well simulate the dynamic behavior of the brain at resting state.[29]Using this model, we try to study the dynamic differences between the AMC group and the control group, and reveal the functions of AMC training on brain dynamics.

    The rest of this article is structured as follows.The resting-state functional brain networks are introduced in Second 2.In Section 3, the resting-state brain network dynamics model is described and simulations are conducted.The task-state brain dynamics model is given in Section 4,and the dynamic differences between the AMC group and the control group are pointed out too.The differences in the topological structure of the resting-state functional networks between the AMC group and the control group are analyzed in Section 5.The conclusion is given in Section 6.

    2.Functional brain networks of the resting-state

    The structures of functional networks in the restingstate are similar to those in the task-induced states to some extent,[31,32]and the task-related activation patterns can be explained partly by the functional connections between brain regions in the resting state.[33]In order to study the effect of functional connectivity on brain dynamics, we collect fMRI data of 25 subjects trained by AMC and other 24 subjects as a control group.The data are measured in the resting state.The corresponding functional networks in the resting-state are constructed.All subjects are from urban families,have normal hearing,normal or corrected vision,no neurological or mental disorders,and no special educational assistance requirements.Participants in the experimental group and the control group are randomly selected.The details about the experiment can be found in Ref.[34].

    Table 1.Regions of interest included in the AAL-atlas(45 in each cerebral hemisphere,90 in total,odd/even numerical order for left/right hemisphere).

    Different templates for brain regions would lead to different number of nodes and topology structures in brain networks.[35]We believe that the main dynamic properties would be similar in different templates.In order to make comparisons with previous studies, we used the widely adopted AAL template,[36]which divides the brain into 90 regions as shown in Table 1.The correlations between brain regions are expressed by Pearson correlation.The correlation coefficient matrix of the functional brain network for a subject is denoted byR.Rijis the correlation coefficient between thei-th andj-th brain regions.

    Usually, there are individual differences in functional connections among subjects.[37]To highlight the changes in brain functional connectivity caused by AMC training and to avoid the influence of individual differences,we take an average on functional connections for different groups,that is,two averaged correlation coefficient matrices are obtained for the AMC group and the control group.The correlation coefficient matrices are shown in Fig.1.

    Fig.1.Correlation coefficient matrices of the control group and AMC group.Left: the control group.Right: the AMC group.

    In addition, in constructing a functional brain network,experimental noise should be removed while the connection density between brain regions is in a reasonable range to ensure that there are no isolated nodes.It is therefore necessary to select a thresholdTfor the correlation coefficient matrix.When the absolute value of the correlation coefficientRijis less thanT, it is considered that there is no functional connection between the two brain regions, and theRi jis turned to be 0.For averaged correlation coefficient matrices of the AMC group and the control group,we choose thresholdTas 0.22 based on the following criteria: no isolated nodes and a connection density of approximately 20%in networks.The matrices after thresholding are shown in Fig.2.It can be found that the connection structures of brain networks for the AMC group and the control group are similar,but not identical.

    Fig.2.The correlation coefficient matrices of the control group and the experimental group after threshold processing.Left: the control group.Right: the AMC group.

    3.Resting-state brain dynamics

    3.1.Model of resting-state brain dynamics

    The HR model can exhibit most dynamic properties observed at the level of real neurons,such as regular and irregular spiking and bursting,and at the level of the collective brain behaviors.In order to explore the influence of functional connections between brain regions on the dynamic activities of the brain,the HR model is introduced to functional brain networks of the AMC group and the control group.A dynamics model for a resting-state functional brain network was established by using the HR model instead of brain regions and the correlation coefficients between brain regions as the connection weights.[29]The dynamics model can be expressed as

    wherexirepresents the cortical electrical signal of thei-th brain region, anda,b,c,d,r,s, andxrare model parameters;Ibiasis the bias current, which regulates the intensity of dynamic activity of brain regions;σis the gain coefficient, which regulates the connection strength between brain regions.The parameters are set as follows,a=1.0,b=3.0,c=1.0,d=5.0,s=4.0,r=0.006,xr=-1.56,Ibias=1.6,andσ=1.2.We find that all brain regions in the functional brain networks of the AMC group and the control group exhibit bursts of spikes, as shown in Fig.3(a).Different from a single HR model, the nodes in the networks show obvious entrainment behavior, that is, the fast oscillations of the encoded signals are mixed with lower frequency oscillations,which is an effective means for information transmission between neurons.[38,39]The dynamic behavior shown in Fig.3(a)is consistent with the remote communication mode in the brain,in which low-frequency resonances determine the longdistance information exchange while high-frequency oscillations characterize the local neural activities in the brain.[40]

    We perform Fourier transform on the cortical electrical signal of brain regions.Fourier transform of Fig.3(a)is shown in Fig.3(b).Compared with the resting state EEG signal shown in Fig.3(c),it can be found that the spectrum characteristics of the cortical electrical signal generated by our model is consistent with that of the resting state EEG signals.

    In order to analyze the dynamic characteristics of brain regions,we count the spikes in bursts.In accordance with the patterns of node bursts, the time width of the burst is defined as the interval from the time when the potential,denoted asx,exceeds the threshold of-1.3 to the time when it falls below the threshold of-1.3 (ending the bursting discharge).Additionally, within this bursting discharge interval, a spike was confirmed when the potential change goes from negative to positive.

    The spikes in bursts are almost stable or random in some brain regions,while in other brain regions spikes in bursts look similar periodic, as shown in Figs.4(a)–4(c).In order to distinguish the dynamic behavior of the brain regions, we make spatial transformation on the spikes according to the following equation:

    whereMnrepresents the number of spikes in then-th burst,andNrepresents the total number of bursts.The results of spatial transformation on the spikes are shown in Figs.4(d)–4(f).The dynamic behavior of brain regions can be divided into the stable mode,the chaotic mode,and the periodic mode.

    Fig.3.(a) The signal of PUT(R) of control group.(b) Spectrum distribution of PUT(R)’s signal.(c) Typical resting-state EEG spectrum distribution.

    Fig.4.The numbers of spikes and their spatial transformations of different modes (control group).[(a), (d)] IFGoperc(L), stable mode.[(b),(e)]TPOsup(R),chaotic mode.[(c),(f)]PUT(R),periodic mode.

    3.2.Results of resting-state dynamics

    Our previous work revealed that the brain regions in periodic mode are active ones.[29]Under the same parameter settings,we obtain all activated brain regions of the AMC group and the control group.As shown in Fig.5, the AMC group and the control group have four common activated brain regions under the same parameter settings, namely CAU(L,R)and PUT(L,R)brain areas.CAU and PUT form the basal ganglia,which play an important role in advanced cognitive activities such as learning, memory, reward, motivation, emotion,and romantic interaction.[41,42]Meanwhile, the AMC group has other three activated brain regions which are AMYG(L),and THA(L,R).Further analysis revealed that the activated brain regions in both the AMC group and the control group are located in the basal ganglia.Previous research indicates that the basal ganglia network is similar to the default mode network and is one of the main activation areas in the restingstate brain.[43]More brain regions are inspired in the AMC group.It is confirmed that AMC training induces functional changes in brain activation and such plasticity may be transferable beyond the AMC.[44]Wanget al.,based on task-based fMRI,confirmed similar group difference of the activation regions between the AMC group and the control group,and the difference was closely related to the subjects’performance in then-back task.[34]

    Our dynamical simulations based on the resting-state functional brain network prove that the differences of brain activation between the AMC group and the control group exist.Since the parameter settings of our dynamic model are the same for both the AMC group and the control group,we guess that the group differences on the activation of brain regions are caused by the different connectivity of their functional brain networks.

    Fig.5.(a)Activated brain area of control group,(b)activated brain area of AMC group,plotted with BrainNet Viewer.[45]

    4.Task state brain dynamics

    4.1.Dynamic characteristics of task states

    Transcranial alternating current stimulation (tACS) can induce brain activation or inhibition and has been widely cited in studies related to brain dynamics.[46]For example,tACS at the appropriate frequency induces phosphenes, activating the visual cortex.[47]We therefore add sinusoidal signals to brain networks to simulate the brain dynamics in the task states.The dynamic model of the brain is transferred as follows:

    whereIextis an external sine stimulus, and the other parameters are the same as those used in Section 3.After adding the sinusoidal signals,the activities of most brain regions become periodic mode.The result further demonstrates that the periodic mode is a kind of active performance.We count the excitation rate(the proportion of brain regions in the periodic mode)and the response intensity(the average dynamic intensity of brain regions obtained under Fourier transform at the same cycle as the external sinusoidal signals).The response intensity represents the synchronization of networks.

    As shown in Fig.6(a), when the cycle of the sinusoidal signals is near 150(the time length in the HR model),the dynamics of the brain networks varies significantly.At small cycle, the excitation rates of the AMC group and the control group are similar,and both show high excitation rates.When the cycle is about 150 or above 150,the periodic mode is suppressed, and the AMC group exhibits a higher excitation rate than the control group.The response intensities at different cycles are plotted in Fig.6(b).The response intensity is dependent on the cycle.When the cycle is less than 150, the response intensity displays positive correlation with the cycle,and the response intensity of the AMC group is lower than that of the control group.When the cycle is greater than 150,the response intensities of the two groups are relatively stable,while the response intensity of the AMC is higher than that of the control group.

    Fig.6.(a)The excitation rates.(b)The average response intensities.

    4.2.State space analysis

    The state space is a two-dimensional space shaped by the average time width of the bursts and the average number of spikes in the bursts, and is used to describe the communication behavior between the neuron clusters.[38]We introduce the state space to find out the origin of the dynamic differences between the AMC group and the control group.We mark the dynamic behaviors of all brain regions for the AMC group and the control group with the stimulus of sinusoidal signals at different cycles in the state space,as shown in Fig.7.

    Analyzing the distributions of brain region activities in the state space,one can observe that obvious mode switching exists for both groups as the cycle of sinusoidal signals increases, as shown in Figs.7(a)–7(c).In the state space, the power exponents of the state distributions of the brain activities evolve from positive to negative as the cycle of sinusoidal signals increases.The power law is strong evidence that the brain works near a critical state.[48]Based on multiple datasets,Fonteneleet al.found that the critical exponents point to a universality class different from mean-field directed percolation.[49]Further research by Fosqueet al.revealed that as external stimuli drive the cortex, quasi criticality predicts a departure from criticality along a Widom line with exponents that decrease in absolute value,while still approximately maintaining a dynamical scaling relation.[50]

    Fig.7.The distributions of brain activities in the state space with the stimulus of sinusoidal signals at different cycles:(a)cycle=8,(b)cycle=150,(c)cycle=256.(d)Exponent distribution over cycles.

    Near the transition points between positive and negative exponents,the average excitation rates of both groups decrease rapidly as shown in Fig.6(a).The response intensity represents the synchronization of networks.From Fig.6(b), one can find that the overall synchronization of the brain networks for both groups has positive correlation with the cycle below the transition points and reach the maximum near the transition points.We define the states of brain networks with a high excitation rate and positive power exponent as the overall task states, and others are the critical state with high synchronous oscillation.

    As shown in Fig.6, the AMC group has a higher excitation rate, but a lower response intensity during the overall task states than the control group,which means that the AMC group can do tasks with less concussion and higher efficiency by integrating the functional brain networks.The difference of the power exponents between the AMC group and the control group shown in Fig.7(d) indicates that the AMC group exhibits the overall task states in longer length of cycles,that is,the AMC group has stronger ability to perform task.In addition,the power exponents with different cycles also show a power law, which indicates further that the brain works near the critical state.

    Previous research showed that AMC training can lead to a change in the excitation patterns of the subjects’brain,[9,11]strengthen working memory,[13]and improve the efficiency of numerical processing.[12,51]Evidence shows that these changes may be transferable beyond the AMC.[44]Our simulations based on functional brain networks confirm the differences of brain dynamics between the AMC group and control group,and AMC training can enhance the abilities of subjects to perform tasks.Since we use the same models for the AMC group and control group, our results indicate that the differences of brain dynamics come from the change of functional structures in brain networks by AMC training.

    5.Graphy theoretical analysis of resting-state functional networks

    Because the parameter settings for each node are consistent in our dynamic model, the differences in the overall and local dynamics of the network are attributed to differences in network structure only.Here, we use graph theory metrics to further analyze the differences in the topological structure of the resting-state functional networks between the AMC group and the control group.We calculate the connection density,connection strength, clustering coefficient, and characteristic path length for both networks.The specific definitions of these metrics are as follows:

    Connectivity densityThe ratio of connections between brain regions to the possible number of connections, i.e., the proportion of non-zero elements in the correlation coefficient matrix after thresholding.

    Connection strengthThe connection strength of a node represents the average strength of its connections to all other nodes.The network’s connection strength is the average of the connection strengths of all nodes,i.e.,the mean of the correlation coefficient matrix after thresholding.

    Clustering coefficientThe clustering coefficient of a node is the ratio of the number of triangles formed by the node with other nodes to the total possible number of triangles.The network’s clustering coefficient is the mean of the clustering coefficients of all nodes,representing local efficiency.

    The characteristic path lengthThe shortest path length of a node is the average number of minimum connecting edges required for it to be connected to other nodes.The network’s characteristic path length is the mean of the shortest path lengths for all pairs of nodes,representing global efficiency.

    The topological characteristics of the AMC group and the control group networks are shown in Fig.8.Compared to the control group, the AMC group exhibits higher connection strength and connection density.Furthermore, the AMC group has a higher clustering coefficient and a lower shortest path length, indicating higher local and global efficiency.We believe that this optimization in network connectivity is the reason for the improvement in the dynamics of the AMC group.

    Fig.8.Differences in graph theory metrics between the networks of the control group and the AMC group: (a)connection density,(b)connection strength, (c) clustering coefficient, (d) the characteristic path length.

    6.Conclusion

    The dynamic behaviors of the brain at the resting and task states for the AMC group and the control group are simulated with the dynamic model established based on functional brain networks and the HR model.For the resting state, besides the common activated brain regions, more brain regions are inspired in the AMC group.A stimulus with sinusoidal signals to brain networks is used to simulate the brain dynamics in the task states.The dynamic characteristics are extracted by the excitation rates, the response intensities and the state distributions.Observing these quantities at different cycles of sinusoidal signals,obvious mode switching exists near the cycle of 150 for both groups.The state distributions of brain in the state space are expressed as a power law.The power exponents are dependent on the cycle of the external signals and turn from positive power exponent to negative ones as the cycle of the external signals increases.The transition points are near the cycle of 150.It is confirmed that brain works near a critical state.The AMC group displays a more extended active period in the overall task states,suggesting superior capabilities of task performance.

    Since the dynamic parameters of each brain region are consistent during the experimental process, the differences in dynamics between the AMC group and the control group are attributed to differences in network connectivity only.Analysis reveals that compared to the control group,the AMC group has a higher network clustering coefficient and a smaller characteristic path length, indicating better local and global efficiency.We believe that the optimization in network structure is a potential reason for the superior dynamic performance of the AMC group.Our work proves that the changes in the functional structures of the brain with AMC training make the dynamic behaviors of subjects more efficient.The research method in this work can be used to simulate dynamic behavior of the brain and evaluate its dynamic performance.

    Acknowledgement

    Project supported by the National Natural Science Foundation of China(Grant Nos.62276229 and 32071096).

    猜你喜歡
    王淼國光貴陽
    膽小的迪姆賽
    神秘的外衣
    保衛(wèi)紅房子
    Application of the edge of chaos in combinatorial optimization?
    找朋友
    貴陽分店
    藝術(shù)品(2019年4期)2019-05-30 04:45:40
    尋味貴州——貴陽
    四種解法中到底誰對誰錯
    叫桃的女人
    尋找一個叫桃的女人
    短篇小說(2014年7期)2014-09-20 15:08:53
    久久久水蜜桃国产精品网| 精品国产乱码久久久久久男人| 又大又爽又粗| 高清在线国产一区| 老司机深夜福利视频在线观看| 在线亚洲精品国产二区图片欧美| 亚洲男人天堂网一区| 交换朋友夫妻互换小说| 午夜成年电影在线免费观看| 我要看黄色一级片免费的| 亚洲精品国产区一区二| 一二三四社区在线视频社区8| 视频在线观看一区二区三区| 精品久久蜜臀av无| 亚洲成人手机| 国产日韩欧美亚洲二区| 露出奶头的视频| 久久久国产成人免费| 欧美精品av麻豆av| 亚洲欧美激情在线| 亚洲久久久国产精品| 色精品久久人妻99蜜桃| 久久午夜综合久久蜜桃| 动漫黄色视频在线观看| 国产精品二区激情视频| 亚洲精华国产精华精| 亚洲国产欧美在线一区| 午夜91福利影院| 纯流量卡能插随身wifi吗| 午夜福利,免费看| 视频在线观看一区二区三区| 在线播放国产精品三级| 亚洲欧美日韩另类电影网站| 欧美黄色片欧美黄色片| 波多野结衣一区麻豆| 两个人看的免费小视频| 少妇粗大呻吟视频| 国产一区二区激情短视频| 成年人黄色毛片网站| 欧美日韩亚洲国产一区二区在线观看 | 女同久久另类99精品国产91| 欧美变态另类bdsm刘玥| 桃红色精品国产亚洲av| 国产成人av激情在线播放| 国产精品亚洲一级av第二区| av视频免费观看在线观看| 国产成人精品无人区| 可以免费在线观看a视频的电影网站| 精品国产国语对白av| 三上悠亚av全集在线观看| 天堂动漫精品| 国产欧美日韩一区二区精品| 亚洲欧美日韩高清在线视频 | 热99re8久久精品国产| 免费在线观看日本一区| 好男人电影高清在线观看| av有码第一页| 国产国语露脸激情在线看| 国产免费视频播放在线视频| 伊人久久大香线蕉亚洲五| 成年人免费黄色播放视频| 一区在线观看完整版| 亚洲情色 制服丝袜| 精品国内亚洲2022精品成人 | 亚洲国产欧美日韩在线播放| 精品少妇黑人巨大在线播放| 成人亚洲精品一区在线观看| 国产成人免费观看mmmm| xxxhd国产人妻xxx| 下体分泌物呈黄色| 午夜福利一区二区在线看| 人成视频在线观看免费观看| 国产日韩一区二区三区精品不卡| 12—13女人毛片做爰片一| 99精品欧美一区二区三区四区| a在线观看视频网站| 久久午夜综合久久蜜桃| 狠狠精品人妻久久久久久综合| 99在线人妻在线中文字幕 | xxxhd国产人妻xxx| 亚洲精品在线观看二区| 国产福利在线免费观看视频| 搡老乐熟女国产| 亚洲五月色婷婷综合| 一进一出好大好爽视频| 最黄视频免费看| 丁香六月欧美| 韩国精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 成人国语在线视频| 国产精品久久久久成人av| 国产又爽黄色视频| 国产黄色免费在线视频| 亚洲精品av麻豆狂野| 纯流量卡能插随身wifi吗| 亚洲av日韩精品久久久久久密| 国产国语露脸激情在线看| 亚洲av片天天在线观看| 丝袜喷水一区| 色精品久久人妻99蜜桃| 久久青草综合色| 99精品久久久久人妻精品| av国产精品久久久久影院| 久久久久久久国产电影| 精品免费久久久久久久清纯 | 欧美日韩av久久| 精品人妻在线不人妻| 亚洲精品自拍成人| av国产精品久久久久影院| 欧美+亚洲+日韩+国产| 老司机影院毛片| 色婷婷av一区二区三区视频| 亚洲自偷自拍图片 自拍| 亚洲国产毛片av蜜桃av| 国产成人精品无人区| 成人18禁在线播放| 天堂俺去俺来也www色官网| 91精品三级在线观看| 欧美人与性动交α欧美精品济南到| 国产精品电影一区二区三区 | 国产一区二区三区在线臀色熟女 | 99香蕉大伊视频| 窝窝影院91人妻| 老熟女久久久| 亚洲精华国产精华精| 国产aⅴ精品一区二区三区波| 不卡一级毛片| 亚洲av日韩精品久久久久久密| 国产成人免费观看mmmm| 在线观看免费视频日本深夜| 亚洲精品粉嫩美女一区| 久久久久久亚洲精品国产蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕精品免费在线观看视频| 国产一卡二卡三卡精品| 人人妻人人澡人人看| 欧美乱码精品一区二区三区| 他把我摸到了高潮在线观看 | 丁香六月天网| 国产成人影院久久av| 午夜免费鲁丝| 免费看a级黄色片| 51午夜福利影视在线观看| 黑人猛操日本美女一级片| 一区二区三区激情视频| 精品国产一区二区三区久久久樱花| 色综合欧美亚洲国产小说| 少妇 在线观看| 久久av网站| 18禁黄网站禁片午夜丰满| 91九色精品人成在线观看| 99re在线观看精品视频| 久久性视频一级片| 亚洲av成人一区二区三| 久久久欧美国产精品| 久久午夜亚洲精品久久| 欧美人与性动交α欧美软件| 人人妻人人爽人人添夜夜欢视频| 久久精品国产99精品国产亚洲性色 | 午夜福利欧美成人| 一区在线观看完整版| 天天躁日日躁夜夜躁夜夜| 国产xxxxx性猛交| 中亚洲国语对白在线视频| 日本av手机在线免费观看| 久久久久久久国产电影| 久久免费观看电影| 丝袜美足系列| 国产不卡av网站在线观看| 欧美人与性动交α欧美精品济南到| 欧美乱码精品一区二区三区| 一进一出好大好爽视频| 两人在一起打扑克的视频| 曰老女人黄片| 大香蕉久久网| 国产精品欧美亚洲77777| 亚洲精品在线美女| 国产av精品麻豆| 婷婷丁香在线五月| 亚洲自偷自拍图片 自拍| av天堂在线播放| 男女免费视频国产| 少妇猛男粗大的猛烈进出视频| 亚洲成人免费av在线播放| 久久天躁狠狠躁夜夜2o2o| 日日爽夜夜爽网站| 在线永久观看黄色视频| 夜夜夜夜夜久久久久| 国产精品98久久久久久宅男小说| 日韩视频一区二区在线观看| 久久热在线av| 搡老岳熟女国产| 国产精品偷伦视频观看了| 午夜91福利影院| 热99久久久久精品小说推荐| 久久久精品国产亚洲av高清涩受| 丰满饥渴人妻一区二区三| 精品久久久久久久毛片微露脸| 欧美成人免费av一区二区三区 | 午夜免费成人在线视频| 又大又爽又粗| 99九九在线精品视频| 久久精品aⅴ一区二区三区四区| 69精品国产乱码久久久| 欧美 亚洲 国产 日韩一| 日韩欧美三级三区| 人妻久久中文字幕网| 另类精品久久| 亚洲精华国产精华精| 91精品国产国语对白视频| 久久狼人影院| 老司机午夜十八禁免费视频| 不卡一级毛片| 男女免费视频国产| 久久精品国产综合久久久| 蜜桃在线观看..| 两人在一起打扑克的视频| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人| 我的亚洲天堂| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线| 黄色视频在线播放观看不卡| 丝袜美腿诱惑在线| 91精品国产国语对白视频| 人妻 亚洲 视频| 国产国语露脸激情在线看| 黄色成人免费大全| 大香蕉久久网| 青青草视频在线视频观看| 欧美av亚洲av综合av国产av| 精品国产一区二区久久| 激情视频va一区二区三区| 性色av乱码一区二区三区2| 怎么达到女性高潮| 久热这里只有精品99| 亚洲人成伊人成综合网2020| 蜜桃国产av成人99| 男女边摸边吃奶| 亚洲av成人一区二区三| 19禁男女啪啪无遮挡网站| 国产精品免费大片| 人人妻,人人澡人人爽秒播| 91国产中文字幕| 欧美乱码精品一区二区三区| 五月开心婷婷网| 国产精品欧美亚洲77777| 老熟女久久久| 99香蕉大伊视频| 国产精品99久久99久久久不卡| 性高湖久久久久久久久免费观看| 91老司机精品| 国产精品国产av在线观看| av片东京热男人的天堂| 日韩免费av在线播放| 男女边摸边吃奶| 女人爽到高潮嗷嗷叫在线视频| 日本精品一区二区三区蜜桃| 黑人猛操日本美女一级片| 国产男靠女视频免费网站| 如日韩欧美国产精品一区二区三区| 欧美激情久久久久久爽电影 | 亚洲一区中文字幕在线| 丝袜美腿诱惑在线| 精品国内亚洲2022精品成人 | 午夜福利,免费看| 婷婷丁香在线五月| 波多野结衣av一区二区av| 亚洲天堂av无毛| 成年版毛片免费区| 男女免费视频国产| 亚洲自偷自拍图片 自拍| 国产三级黄色录像| 精品国内亚洲2022精品成人 | 女性生殖器流出的白浆| 啦啦啦中文免费视频观看日本| 一区福利在线观看| 最新在线观看一区二区三区| 韩国精品一区二区三区| 操出白浆在线播放| 亚洲五月色婷婷综合| 天堂中文最新版在线下载| 国产1区2区3区精品| 精品乱码久久久久久99久播| 熟女少妇亚洲综合色aaa.| 国产精品影院久久| 国产精品美女特级片免费视频播放器 | 国产伦人伦偷精品视频| 色播在线永久视频| 十八禁网站免费在线| 男女高潮啪啪啪动态图| 国产成人系列免费观看| 十分钟在线观看高清视频www| 999精品在线视频| av超薄肉色丝袜交足视频| 亚洲自偷自拍图片 自拍| 日韩熟女老妇一区二区性免费视频| 久久热在线av| 亚洲第一青青草原| 成人影院久久| 亚洲视频免费观看视频| 国产精品欧美亚洲77777| 天天躁狠狠躁夜夜躁狠狠躁| 女人精品久久久久毛片| 国产日韩欧美亚洲二区| 夜夜夜夜夜久久久久| 日韩免费av在线播放| 中文欧美无线码| 亚洲精品一卡2卡三卡4卡5卡| 精品免费久久久久久久清纯 | 新久久久久国产一级毛片| 成人黄色视频免费在线看| 亚洲专区国产一区二区| 久久精品91无色码中文字幕| 国产精品98久久久久久宅男小说| a级毛片在线看网站| 不卡av一区二区三区| 午夜福利免费观看在线| 大片免费播放器 马上看| 国产欧美亚洲国产| 国产精品电影一区二区三区 | 成人黄色视频免费在线看| 一级黄色大片毛片| 麻豆av在线久日| 欧美激情久久久久久爽电影 | 丰满饥渴人妻一区二区三| 男女边摸边吃奶| 亚洲午夜精品一区,二区,三区| 极品教师在线免费播放| 成人永久免费在线观看视频 | 成年版毛片免费区| 电影成人av| 久久久久网色| 在线看a的网站| av电影中文网址| 亚洲av电影在线进入| 免费一级毛片在线播放高清视频 | 香蕉久久夜色| 中文字幕人妻丝袜制服| 中文字幕精品免费在线观看视频| 叶爱在线成人免费视频播放| 国产在视频线精品| 久久人妻熟女aⅴ| netflix在线观看网站| 成人免费观看视频高清| 丝袜美足系列| 国产精品av久久久久免费| 真人做人爱边吃奶动态| 色94色欧美一区二区| 国产成人免费无遮挡视频| 啦啦啦 在线观看视频| 淫妇啪啪啪对白视频| a级片在线免费高清观看视频| 午夜成年电影在线免费观看| 久久久欧美国产精品| 男女高潮啪啪啪动态图| a级片在线免费高清观看视频| 高清在线国产一区| 国产伦理片在线播放av一区| 9色porny在线观看| 99在线人妻在线中文字幕 | 国产高清视频在线播放一区| 99热国产这里只有精品6| 免费在线观看影片大全网站| 69精品国产乱码久久久| www.自偷自拍.com| 极品人妻少妇av视频| 嫩草影视91久久| 黄色片一级片一级黄色片| 女人被躁到高潮嗷嗷叫费观| 国产成人av教育| 丝瓜视频免费看黄片| 国产精品美女特级片免费视频播放器 | 中国美女看黄片| 在线观看66精品国产| 麻豆av在线久日| 麻豆成人av在线观看| 1024视频免费在线观看| 色老头精品视频在线观看| 黄色视频不卡| 国产精品欧美亚洲77777| 丝瓜视频免费看黄片| 男女边摸边吃奶| 搡老熟女国产l中国老女人| 亚洲精品美女久久av网站| 怎么达到女性高潮| videosex国产| 一本大道久久a久久精品| 性色av乱码一区二区三区2| 久久精品亚洲熟妇少妇任你| 久久久久久久精品吃奶| 黄色成人免费大全| 久久精品亚洲av国产电影网| 午夜视频精品福利| 亚洲伊人久久精品综合| 天天影视国产精品| 久久久精品94久久精品| 中文字幕人妻丝袜一区二区| 日日夜夜操网爽| 成人免费观看视频高清| 国产在视频线精品| 欧美精品人与动牲交sv欧美| 黑人巨大精品欧美一区二区mp4| 嫁个100分男人电影在线观看| 日韩一区二区三区影片| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品在线观看二区| 精品欧美一区二区三区在线| 考比视频在线观看| 色视频在线一区二区三区| 他把我摸到了高潮在线观看 | 国产一区二区三区综合在线观看| 久久 成人 亚洲| 99re6热这里在线精品视频| 亚洲精品国产色婷婷电影| 久久久国产欧美日韩av| 一进一出抽搐动态| 热re99久久精品国产66热6| 欧美日韩成人在线一区二区| 免费在线观看影片大全网站| 国产色视频综合| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆69| 精品亚洲成国产av| 国产成人啪精品午夜网站| 他把我摸到了高潮在线观看 | 丝袜美腿诱惑在线| 成年版毛片免费区| 成人手机av| 飞空精品影院首页| 最新在线观看一区二区三区| 国产主播在线观看一区二区| 黄片播放在线免费| 亚洲av电影在线进入| 免费在线观看视频国产中文字幕亚洲| 男女边摸边吃奶| 不卡一级毛片| 性少妇av在线| 成年动漫av网址| 成人18禁高潮啪啪吃奶动态图| 天堂动漫精品| 丰满少妇做爰视频| 老鸭窝网址在线观看| 亚洲三区欧美一区| 三级毛片av免费| 精品亚洲成国产av| 男女高潮啪啪啪动态图| 国产一区二区三区综合在线观看| 女同久久另类99精品国产91| 亚洲熟女精品中文字幕| 欧美+亚洲+日韩+国产| 老司机午夜十八禁免费视频| av在线播放免费不卡| 熟女少妇亚洲综合色aaa.| 1024香蕉在线观看| 日日摸夜夜添夜夜添小说| 黄色视频,在线免费观看| 在线十欧美十亚洲十日本专区| 亚洲精品自拍成人| 悠悠久久av| 国产麻豆69| 巨乳人妻的诱惑在线观看| 国产精品 欧美亚洲| 91精品三级在线观看| 99久久99久久久精品蜜桃| 国产福利在线免费观看视频| 久久久久久免费高清国产稀缺| 无人区码免费观看不卡 | 制服人妻中文乱码| 亚洲av日韩精品久久久久久密| 免费在线观看视频国产中文字幕亚洲| 久久久久国内视频| 亚洲全国av大片| 亚洲精品乱久久久久久| 亚洲成人国产一区在线观看| 精品高清国产在线一区| 亚洲一码二码三码区别大吗| 国产亚洲一区二区精品| 性少妇av在线| 欧美变态另类bdsm刘玥| 国产免费视频播放在线视频| 国产成人精品无人区| 丰满饥渴人妻一区二区三| 国产精品秋霞免费鲁丝片| 精品国产一区二区三区四区第35| 黄片播放在线免费| 国内毛片毛片毛片毛片毛片| 麻豆成人av在线观看| 免费在线观看黄色视频的| 搡老岳熟女国产| 国产日韩一区二区三区精品不卡| 国产精品一区二区在线观看99| 美女国产高潮福利片在线看| 亚洲天堂av无毛| 午夜福利视频在线观看免费| 伦理电影免费视频| 欧美日韩中文字幕国产精品一区二区三区 | videos熟女内射| 午夜久久久在线观看| 一本一本久久a久久精品综合妖精| 在线观看免费高清a一片| 国产福利在线免费观看视频| 婷婷成人精品国产| 一本色道久久久久久精品综合| 久久天堂一区二区三区四区| 大片免费播放器 马上看| av网站免费在线观看视频| 超碰97精品在线观看| 国产97色在线日韩免费| 国产精品免费一区二区三区在线 | 精品福利观看| 国产熟女午夜一区二区三区| 丁香欧美五月| 精品国产乱子伦一区二区三区| 欧美激情高清一区二区三区| 一边摸一边做爽爽视频免费| 国产一区二区在线观看av| 人妻 亚洲 视频| 国产免费视频播放在线视频| 欧美精品人与动牲交sv欧美| 亚洲专区中文字幕在线| 欧美另类亚洲清纯唯美| 岛国毛片在线播放| 一边摸一边抽搐一进一出视频| 免费观看a级毛片全部| 国产又爽黄色视频| xxxhd国产人妻xxx| 国产成人精品无人区| 免费在线观看影片大全网站| 一进一出抽搐动态| 国产男靠女视频免费网站| 成人国产av品久久久| 男女边摸边吃奶| 亚洲欧美激情在线| 一区在线观看完整版| 狠狠精品人妻久久久久久综合| 人人妻人人澡人人看| 男女午夜视频在线观看| 精品午夜福利视频在线观看一区 | 黄片大片在线免费观看| 丝袜在线中文字幕| 亚洲国产av新网站| 国产精品久久久av美女十八| www.精华液| 久久中文看片网| 人人妻人人添人人爽欧美一区卜| 99在线人妻在线中文字幕 | 国产精品国产高清国产av | 国产欧美日韩综合在线一区二区| av天堂在线播放| 久久久久国产一级毛片高清牌| 精品一品国产午夜福利视频| 久久影院123| 国产在线视频一区二区| 亚洲成人手机| 中文字幕人妻丝袜制服| 欧美日韩av久久| a在线观看视频网站| 欧美亚洲日本最大视频资源| 国产欧美日韩综合在线一区二区| 亚洲少妇的诱惑av| 国产亚洲欧美精品永久| 蜜桃在线观看..| 精品国产乱子伦一区二区三区| 日韩熟女老妇一区二区性免费视频| 在线观看舔阴道视频| 亚洲国产精品一区二区三区在线| 一级黄色大片毛片| 久9热在线精品视频| 日韩欧美国产一区二区入口| 久久精品国产99精品国产亚洲性色 | 欧美精品亚洲一区二区| 国产黄频视频在线观看| 少妇粗大呻吟视频| 女人精品久久久久毛片| 久久av网站| 欧美国产精品va在线观看不卡| 午夜精品久久久久久毛片777| 午夜福利乱码中文字幕| 男男h啪啪无遮挡| 国产成人影院久久av| tube8黄色片| 亚洲 欧美一区二区三区| 国产精品久久久久久人妻精品电影 | 麻豆av在线久日| 国产男靠女视频免费网站| 91老司机精品| 天天操日日干夜夜撸| 国产深夜福利视频在线观看| 成人18禁高潮啪啪吃奶动态图| 黄频高清免费视频| 纯流量卡能插随身wifi吗| 黄色丝袜av网址大全| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美日韩在线播放| 多毛熟女@视频| 精品久久久久久久毛片微露脸| 激情视频va一区二区三区| av网站在线播放免费| 久久精品人人爽人人爽视色| 亚洲中文字幕日韩| 精品欧美一区二区三区在线| 国产精品.久久久| 在线观看免费视频日本深夜| 亚洲第一欧美日韩一区二区三区 | 国产精品一区二区在线不卡| 久久人妻福利社区极品人妻图片| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 久久中文字幕人妻熟女| 国产亚洲一区二区精品| 熟女少妇亚洲综合色aaa.| 日韩视频一区二区在线观看| 欧美日韩黄片免| 欧美国产精品一级二级三级|