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

    Application of the back-error propagation artificial neural network (BPANN) on genetic variants in the PPAR-γ and RXR-α gene and risk of metabolic syndrome in a Chinese Han population

    2014-04-18 12:37:03XuZhoKngXuHuiShiJinluoChengJinhuYnqinGoQinLiXinhuYeYingLuXiofngYuJunDuWenongDuQingYeLingZhou
    THE JOURNAL OF BIOMEDICAL RESEARCH 2014年2期

    Xu Zho, Kng Xu, Hui Shi, Jinluo Cheng, Jinhu M, Ynqin Go, Qin Li, Xinhu Ye, Ying Lu, Xiofng Yu, Jun Du, Wenong Du, Qing Ye, Ling Zhou,

    aDepartment of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, Jiangsu 211166, China;

    bDepartment of Endocrinology, the Affiliated Changzhou Second Hospital, Nanjing Medical University, Changzhou, Jiangsu 213003, China;

    cDepartment of Endocrinology, the First Affiliated Hospital, Nanjing Medical University, Nanjing, Jiangsu 210029, China;

    dDepartment of Endocrinology, the Third Affiliated Hospital, Nanjing Medical University, Yizheng, Jiangsu 211400, China.

    Application of the back-error propagation artificial neural network (BPANN) on genetic variants in the PPAR-γ and RXR-α gene and risk of metabolic syndrome in a Chinese Han population

    Xu Zhaoa,△, Kang Xua,△, Hui Shia, Jinluo Chengb, Jianhua Mac, Yanqin Gaod, Qian Lic, Xinhua Yeb, Ying Lua, Xiaofang Yua, Juan Dua, Wencong Dua, Qing Yea, Ling Zhoua,

    aDepartment of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, Jiangsu 211166, China;

    bDepartment of Endocrinology, the Affiliated Changzhou Second Hospital, Nanjing Medical University, Changzhou, Jiangsu 213003, China;

    cDepartment of Endocrinology, the First Affiliated Hospital, Nanjing Medical University, Nanjing, Jiangsu 210029, China;

    dDepartment of Endocrinology, the Third Affiliated Hospital, Nanjing Medical University, Yizheng, Jiangsu 211400, China.

    This study was aimed to explore the associations between the combined effects of several polymorphisms in the PPAR-γand RXR-αgene and environmental factors with the risk of metabolic syndrome by back-error propagation artificial neural network (BPANN). We established the model based on data gathered from metabolic syndrome patients (n = 1012) and normal controls (n = 1069) by BPANN. Mean impact value (MIV) for each input variable was calculated and the sequence of factors was sorted according to their absolute MIVs. Generalized multifactor dimensionality reduction (GMDR) confirmed a joint effect of PPAR-γand RXR-αbased on the results from BPANN. By BPANN analysis, the sequences according to the importance of metabolic syndrome risk factors were in the order of body mass index (BMI), serum adiponectin, rs4240711, gender, rs4842194, family history of type 2 diabetes, rs2920502, physical activity, alcohol drinking, rs3856806, family history of hypertension, rs1045570, rs6537944, age, rs17817276, family history of hyperlipidemia, smoking, rs1801282 and rs3132291. However, no polymorphism was statistically significant in multiple logistic regression analysis. After controlling for environmental factors, A1, A2, B1and B2(rs4240711, rs4842194, rs2920502 and rs3856806) models were the best models (cross-validation consistency 10/10, P = 0.0107) with the GMDR method. In conclusion, the interaction of the PPAR-γand RXR-αgene could play a role in susceptibility to metabolic syndrome. A more realistic model is obtained by using BPANN to screen out determinants of diseases of multiple etiologies like metabolic syndrome.

    back-error propagation artificial neural network (BPANN), metabolic syndrome, peroxisome proliferators activated receptor-γ(PPAR) gene, retinoid X receptor-α(RXR-α) gene, adiponectin

    INTRODUCTION

    Metabolic syndrome is characterized by the simultaneous presence of interrelated metabolic risk factors, such as visceral obesity, high blood pressure and carbohydrate and lipid metabolism abnormalities[1], and it leads to an increased risk for type 2 diabetes mellitus and cardiovascular diseases[2,3]. Although the criteria used to diagnose the metabolic syndrome may be different, several studies have demonstrated that the prevalence of metabolic syndrome has increased dramatically in the past two decades[4-6]. Metabolic syndrome is subjected to multifactorial influences including both environmental and genetic factors. Current researches focus more on the cause of one or certain metabolic syndrome components, but rarely pay attention to the etiology of metabolic syndrome.

    Reportedly, adiponectin is an adipose tissue-derived cytokine with anti-inflammatory, anti-atherogenic and cardioprotective properties. Increasing evidence suggested that adiponectin, linked to central obesity and insulin resistance, is a key contributor to the development of metabolic syndrome[7]. The ADIPOQ gene encoding adiponectin is located on chromosome 3q27 and this chromosome region had been mapped as a susceptibility locus for metabolic syndrome and coronary heart disease by genome-wide scans[8]. There are specific functional peroxisome proliferator response elements (PPREs) in the promoter region of adiponectin. PPAR-γbinds to PPREs as a heterodimer with members of the retinoid X receptor (RXR) nuclear receptor subfamily[9], and increases adiponectin promoter activity in cells from humans[10]. The expression of adiponectin in the adipose tissue is maintained and induced by direct binding of endogenous or exogenous PPAR-γ/RXR heterodimer to the PPRE in the adiponectin promoter[10]. Therefore, the PPAR-γand RXR-αgenes in the adiponectin pathway are candidates for obesity and metabolic syndrome.

    The back-error propagation artificial neural network (BPANN) is a computer-based algorithm that is trained to recognize and categorize complex patterns. BPANN as a new approach to etiology research avoids the limitations of case control study and logistic regression without the distribution form and independence of variables. Our group has previously found the association of single nucleotide polymorphisms (SNPs) in thePPAR-γand RXR-αgenes with metabolic syndrome or type 2 diabetes risks in different populations[11,12]. In the present study, we further explored the application characteristics of BPANN in studying the combined effects of genetic variants in the PPAR-γand RXR-αgene and metabolic syndrome risks in a Chinese Han population. We compared the results of BPANN with those of traditional logistic regression analysis to better understand and use artificial neural network in the study of diseases of multiple etiologies.

    SUBJECTS AND METHODS

    Subjects

    A total of 2,081 Chinese patients, including 1,012 metabolic syndrome and 1,069 non-metabolic syndrome controls, were recruited for metabolic studies with their informed consent. The protocol was approved by the Research Ethics Committee of Nanjing Medical University. All participants were genetically unrelated ethic Han Chinese. The cases were consecutively recruited from the inpatient or outpatient departments of three affiliated hospitals of Nanjing Medical University (the Second Affiliated Changzhou Hospital, the Third Affiliated Hospital and the First Affiliated Hospital) between March 2008 and August 2010, without any restriction on age and sex (430 males and 582 females with a mean age of 55.35±10.62 years). Age (±5 years) and sex-matched non-metabolic syndrome controls who underwent routine annual health examinations within the same geographical area and the period (478 males and 591 females with a mean age of 55.78±13.10 years) were also recruited for the study. In order to collect demographic data and information in environmental exposure history, each participant was interviewed face-toface using a standard questionnaire. After the interview, 5 mL venous blood sample was collected from each participant. The level of physical activity was defined as walking or riding≥ 15 minutes/day and/or lifting or carrying heavy objects at work daily and/or doing sports or physical exercise > 2 hours/week. Tobacco smokers were defined as patients who smoked at least one cigarette per day for over 1 year. Alcohol drinkers were defined as those who had the sum of milliliters of alcohol per week from wine, beer, cider or spirits.

    The new (2009) criteria of metabolic syndrome was based on a joint interim statement of the International Diabetes Federation (IDF); National Heart, Lung and Blood Institute (NHLBI); American Heart Association (AHA); World Heart Federation; International Atherosclerosis Society and International Association, which define metabolic syndrome as the presence of three or more of the following features: triglycerides (TG)≥ 150 mg/dL (1.7 mmol/L) (or drug treatment for elevated triglycerides), high-density lipoprotein cholesterol (HDL-C) < 40 mg/dL (1.0 mmol/L) in mates and < 50 mg/dL (1.3 mmol/L) in females (or drug treatment for reduced HDL-C), fasting plasma glucose(FPG)≥ 100 mg/dL (or drug treatment for elevated glucose), systolic blood pressure (SBP)≥130 and/or diastolic blood pressure (DBP)≥85 mmHg (or antihypertensive drug treatment in a patient with a history of hypertension), and waist circumstances (WC)≥85 cm in men and 80 cm in women (current recommended waist circumstance thresholds for abdominal obesity for people in China)[13].

    Measurements

    Weight and height were measured by trained personnel, and body mass index (BMI, in kg/m2) was calculated. Blood pressure was measured on the right arm, with the patient in a sitting position and after a minimum 10-min rest, using a standard mercury sphygmomanometer. After an overnight fast, venous blood samples were drawn and promptly centrifuged and the plasma was stored at -20°C. Serum adiponectin, FPG and TC, HDL-C, LDL-C and TG were analyzed using human adiponectin ELISA kit (RapidBio Co., Calabasas, CA, USA), glucose oxidase method and enzymatic colorimetric method (Au5400; Olympus, Japan), respectively. DNA was extracted from blood samples by using a phenol-chloroform technique. All measurements were conducted by the manufacturers' protocols.

    SNP selection and genotyping

    In the present study, we first used NCBI dbSNPs database (http://www.ncbi.nlm.nih.gov/), the public HapMap SNP database (http://www.hapmap.org/) and previously reported literatures to identify potentially functional SNPs of the PPAR-γand RXR-αgene. Three potentially functional SNPs of the PPAR-γgene (rs2920502, rs3856806 and rs1801282) and four of the RXR-αgene (rs1045570, rs3132291, rs4240711 and rs4842194) with minor allele frequency (MAF)≥0.05 in the Chinese Han population were identified, including four at the coding region and three at the 3' untranslated region. Then, we used public HapMap SNP database to identify PPAR-γand RXR-αgene tagging SNPs by using tagger with greedy algorithm. Finally, two tagging SNPs (rs17817276 and rs38566806) of the PPAR-γ gene and two of the RXR-α gene (rs1045570 and rs6537944) were selected based on pair-wise tagging (r2≥0.80, MAF ≥0.05) by using genotype data from unrelated HapMap CHB individuals. Genomic sequences were obtained from the HapMap database. Primer version 5.0 and polymerase chain reaction (PCR) primer-introduced restriction analysis were used to design the nine primer sets (Table 1).

    Genetic analyses were performed on genomic DNA extracted from leucocytes of venous blood. We used the TaqMan allelic discrimination assay to genotype the polymorphisms on the platform of 7900HT Realtime PCR System (Applied Biosystems, Foster City, CA). The information on assay conditions and the primers and probes are available upon request. Two negative controls were included in each 384-well reaction plate and individual genotype identification were determined by SDS software 2.0 (ABI). Moreo-ver, to confirm genotyping results, 10% of samples were randomly selected to repeat the procedure and the results were 100% concordant.

    Table 1 Primers and annealing temperature used for RET sequencing

    Statistical analysis

    Hardy-Weinberg equilibrium was assessed within controls using the goodness-of-fit χ2test. The distribution of the general characteristics between metabolic syndrome and non-metabolic syndrome controls was compared by using two-sided Chi-square test and/ or Student's t test. The significance level was set at P <0.05. Both univariate and multiple logistic regression analyses were performed to estimate crude and adjusted standardized partial regression coefficient (β), odds ratios (ORs) and 95% confidence intervals (CIs) for the association between genotypes and metabolic syndrome risks. All statistical analyses were performed by SPSS software (SPSS version 13.0, Chicago, IL, USA). Generalized multifactor dimensionality reduction (GMDR) software (version 1.0.1) (www. healthsystem.virginia.edu/internet/Addiction-Genomics/) was applied for detecting gene-gene and geneenvironment interactions[14]. BPANN prediction model was established and each input neurons' mean impact value (MIV) was calculated using MATLAB7.0 software (MathWorks, Natick, MA)[15].

    RESULTS

    Univariate logistic regression analysis of associations between risk factors and metabolic syndrome

    The clinical characteristics of the 2081 subjects are shown in Table 2. The genotype distributions of all the SNPs satisfied Hardy-Weinberg equilibrium (all P > 0.05 in controls). The associations between the 19 variables, including four of the PPAR-γgene, five of the RXR-αgene, serum adiponectin concentration and other environmental factors, was estimated by binary logistic regression analysis. Table 3 shows the sequences of the impact of various factors on metabolic syndrome risk. The top ten were family history of hyperlipidemia or type 2 diabetes, physical activity, family history of hypertension, BMI, alcohol drinking, gender, rs6537944, rs1801282 and rs4842194.

    Multiple logistic regression analysis of associations between risk factors and metabolic syndrome

    Multiple logistic regression analysis was performed to determine the association between the 19 variables and metabolic syndrome with the stepwise regressive method (the removal probability was 0.1). Table 4 illustrates nine factors in the best model and only six factors were statistically significant, including family history of type 2 diabetes or hyperlipidemia, physical activity, gender, alcohol drinking, BMI, rs4240711, rs2920502 and serum adiponectin.

    BPANN multiple analysis of associations between risk factors and metabolic syndrome

    We used 19 factors as input variables and metabolic syndrome diagnosis as output variables to establish model with all available samples. The transfer function was logsig function. Learning rate and training error were 0.1 and 0.01, respectively. Training steps were set to a maximum of 1000 steps. After the completion of training, the MIV was obtained. Table 5summarizes all related factors, consisting of BMI, serum adiponectin, rs4240711, gender, rs4842194, family history of diabetes, rs2920502, physical activity, alcohol drinking, rs3856806, family history hypertension, rs1045570, rs6537944, age, rs17817276, family history of hyperlipidemia, smoking, rs1801282 and rs3132291 in sequence of the absolute value of MIV.

    Table 2 Basic characteristics of the case and control groups

    Table 3 Univariate logistic regression analysis results

    GMDR for the combined effect of the PPAR-γand RXR-αgene

    The SNPs rs4240711, rs4842194, rs2920502 and rs3856806, which were obtained in the top ten in BPANN multiple analysis, were named as A1, A2, B1and B2in turn. GMDR that evaluated the combined effect of the four SNPs detected A1, A2, B1and B2(rs4240711, rs4842194, rs2920502 and rs3856806) model as the best model (Cross-validation consistency 10/10, P = 0.0447) (Table 6). After controlling for age, gender, smoking, alcohol drinking and physical activity, the results showed that A1, A2, B1and B2model still were the best model from GMDR (Crossvalidation consistency 10/10, P = 0.0107) (Table 7). However, we failed to obtain a significant model when all the nineteen SNPs were added into the analysis (Cross-validation consistency 10/10, P = 0.1719) (data not shown).

    DISCUSSION

    Metabolic syndrome, like virtually all human diseases, results from the interactions between geneticand environmental factors. Environmental factors such as obesity, low levels of physical activity and inappropriate dietary habits are strong determinants of metabolic syndrome[16,17]. Traditional statistical methods due to their own limitation do not show the real relationship between gene and environmental elements with the risk of metabolic syndrome. Binary logistic regression analysis, for example, was used to analyze the association between the 19 variables and metabolic syndrome. Table 3 showed that the locant, according to standardized partial regression coefficient (β) of each variable, of top six environmental factors was hyperlipidemia family history, family history of type 2 diabetes, physical activity, smoking, family history of hypertension and BMI. Relative to environmental factors, the contribution of gene polymorphisms to the susceptibility and causation of diseases was hidden and weaker. Binary logistic regression limited to analyzing the combined action of elements did not provide the evidence for gene-environment interaction.

    Table 4 Multiple logistic regression analysis results

    Table 5 Input variables and sorting of mean influence values (MIV)

    By using the maximum likelihood method, only six variables were statistically significant, including family history of type 2 diabetes and hyperlipidemia, physical activity, gender, BMI and serum adiponectin. The result of multiple logistic regression analysis did not detect the role that PPAR-γgene and RXR-αgene polymorphisms played in the development of metabolic syndrome, but BPANN analysis had proved the association[12]. Adiponectin has an obvious protective effect on metabolic syndrome and its cardiovascular complications[18], but the contribution of adiponectin was weakened in the model, because of the complex interaction between variables.

    Artificial neural networks are a non-linear pattern recognition algorithm consisting of a group of processing units that simulate the function of neurons in human brain and reveal relationships among the input data that cannot always be recognized by conventional models[19]. In recent years, artificial neural networks have been increasingly used in complex medical decision-making, such as the diagnosis of various diseases, investigating the predictive values of disease risk factors and analyzing the complex relationships between gene-environment[20,21]. BPANN has not only been the core of feedforward neural networks, but also embodied the essence of artificial neural networks[15]. The units in a BPANN are highly interconnected by weighted links, very similar to neural synapses. In the constant learning process in which the former feedback of error was used to modify corresponding weights and threshold value, BPANN adjusted the weights of links between neurons in order to associate input data with correct output (such as disease diagnosis). BPANN could fully show all the relationship between factors in a simulation. Therefore, we further used BPANN to identify the real association between PPAR-γand RXR-αgene polymorphisms with susceptibility to metabolic syndrome in southernHan Chinese and the key risk factors for metabolic syndrome. Compared to logistic regression analysis, BPANN showed a better understanding of using artificial neural network in epidemiology. Compared with the results of univariate and multiple logistic regression analyses, the seating arrangements of BMI and serum adiponection concentration were obviously raised in BPANN analysis. The sequence of serum adiponection concentration was raised to the second factor only after BMI from the last sequence in multiple logistic regression analysis. Central obesity is considered a pivotal component in metabolic syndrome. Even in subjects without obesity, a higher BMI tends to correlate with a higher number of positive metabolic syndrome components[22]. It has been reported that adiponectin concentration was an important predictor for the risk of metabolic syndrome[23]. Our result further validated that adiponectin had a great influence on the pathogenesis of metabolic syndrome.

    Table 6 The GMDR models for PPAR-γand RXR-αgene interaction on metabolic syndrome

    Table 7 The GMDR models for PPAR-γand RXR-αgene-environment interaction on metabolic syndrome

    PPAR-γrs2920502, rs3856806 and RXR-αrs4240711, rs4842194 had no association with metabolic syndrome in univariate and multiple logistic regression analyses, but the seating arrangements of the four gene polymorphisms were in front in BPANN analysis, which implied an association between the four genetic variants in PPAR-γand RXR-αand metabolic syndrome. Currently, the studies investigating the associations between rs2920502, rs4240711 and rs4842194 polymorphisms and the risk of metabolic syndrome have been reported rarely. In a previous study[12], we also found that those carrying rs2920502CG and CG/GG genotype had a significantly increased risk of metabolic syndrome and rs4240711GG and AG/GG, rs4842194 CC and CT/ CC genotypes were all associated with prominent protective effects for metabolic syndrome. Up to now, numerous studies have focused on the association between C1431T variant (rs3856806) of PPAR-γand the risk of metabolic syndrome, type 2 diabetes and obesity in several populations[24-26], but the conclusions are conflicting. Li et al.[24]reported that polymorphism C1431T of exon 6 of PPAR-γwas associated with metabolic syndrome risks in a Chinese population study of 423 cases with metabolic syndrome and families without metabolic syndrome. However, in another study of 792 Han Chinese in Beijing, Dongxia et al.[27]failed to show the association, but a significant association between rs3856806 and insulin resistance. In a population-based study of 1910 subjects by our group, divided according to the IDF (2005) criteria, we also did not find differences of rs3856806 among metabolic syndrome patients and non-metabolic syndrome controls by logistic regression analysis[12]. In Meirhaeghe's French population study[17], the C1431T variant was only connected with the other three polymorphisms, including P3-681C > G, P2-689C > T, Pro12Ala, and the GTGC haplotype was associated with metabolic syndrome. It may be that the individual C1431T's role in the pathogenesis of metabolic syndrome is not obvious and the combined effect with other polymorphisms exists, but the traditional statistical methods do not show.

    No polymorphisms in the two genes had been found to associate with metabolic syndrome in univariate and multiple logistic regression analyses, but the seating arrangements of the four gene polymorphisms were in front in BPANN analysis. The finding of BPANN analysis and logistic regression analysis suggested that the interactions between four SNPs may be associated with the risk of metabolic syndrome. We used the GMDR method to indentify the combined effect of the four SNPs (rs2920502, rs3856806, rs4240711 and rs4842194). The GMDR is an extended version of the multifactor dimensionality reduction approach, which allows high-dimensional interactions of multiple factors to be simultaneously retrieved. The GMDR method is a non-parametric and genetic model-free approach that efficiently identifies higher-order interactions between genes and/or geneenvironmental factors with both dichotomous and continuous phenotypes in various population-based study designs. The main idea of GMDR is to reduce multi-dimensional genotypes into one-dimensional binary attributes by pooling genotypes of multiple SNPs using a well-defined classifier[14]. Furthermore, GMDR can reduce the complexity substantially and permit adjustment for discrete and quantitative covariates[14]. GMDR has been performed to successfully identifythe combinations of SNPs that significantly influenced complex diseases[28]. The model of A1, A2, B1 and B2 (rs4240711, rs4842194, rs2920502 and rs3856806) were the best model (Cross-validation consistency 10/10, P = 0.0447). It showed that four SNPs loci interact with one another and suggested interactions between the two genes. After controlling age, gender, smoking, alcohol drinking and physical activity, the results showed that A1, A2, B1and B2(rs4240711, rs4842194, rs2920502 and rs3856806) model were still the best model as the highest Test Balanced Accuracy (0.5901) (Cross-validation consistency 10/10, P = 0.0107). After adjusting age, gender, smoking, alcohol drinking and physical activity, the combined effect of PPAR-γgene and RXR-αgene still existed. PPAR-γ acts as a nuclear receptor-transcription factor by forming a heterodimer with RXR. PPAR-γ/RXR heterodimer directly binds to the functional PPAR-responsive element (PPRE) and increases human adiponectin promoter activity in cells. The PPAR-γgene and RXR-αgene are related to the adiponectinsignal transduction pathway. The two genes may play a significant role in insulin resistance, which is a key feature of metabolic syndrome and type 2 diabetes by affecting adiponectin secretion levels.

    BPANN is designed to detect patterns in input data, which may match output data even if the nature of such patterns is not known a priori; thus, all richer relations between factors can be simulated and provided at a time than the ordinary model[29]. Simultaneously, BPANN not requiring the distribution form and independence of variables, also can handle the problem of collinearity better. Compared with the traditional analysis methods, BPANN may be better suited to predict outcomes when the relationships between the variables are complex, multidimensional and nonlinear in complex biological systems[30]. Neural network models have the ability to detect all possible interactions between predictive variables and provide a model consistent with practical situations. Even so, BPANN analysis was a "black box" and had limited ability to explicitly identify possible causal relationships[31]. The model builder of logistic regression was able to select variables which were most strongly predictive of an outcome based on the magnitude of the standardized partial regression coefficients (β) and the associated odds ratios.

    In this article, we put four SNPs, which was selected by BPANN analysis, into model analyses by GMDR software. We obtained an optimal model with statistical significance. The result suggested an interaction between the RXR-αgene and RXR-αgene. However, we failed to obtain an optimal model when all the 19 SNPs were included in the analysis without screening. Therefore, we conclude that BPANN can be used to select influence factors commonly, especially for early screening of genetic factors. We can further analyze gene interactions by using the results of screening, according to the ranking list of MIV as a relative stable reference.

    [1] Sygnowska E, Piwonska A, Waskiewicz A, Broda G. Socioeconomic factors and the risk of metabolic syndrome in the adult Polish population: the WOBASZ study. Kardiol Pol 2012: 70: 718-27.

    [2] Laaksonen DE, Lakka HM, Niskanen LK, Kaplan GA, Salonen JT, Lakka TA. Metabolic syndrome and development of diabetes mellitus: application and validation of recently suggested definitions of the metabolic syndrome in a prospective cohort study. Am J Epidemiol 2002; 156: 1070-7.

    [3] Onat A, Ceyhan K, Basar O, Erer B, Toprak S, Sansoy V. Metabolic syndrome: major impact on coronary risk in a population with low cholesterol levels--a prospective and cross-sectional evaluation. Atherosclerosis 2002: 165: 285-92.

    [4] Ford E, Giles W, Mokdad A. Mokdad, Increasing prevalence of the metabolic syndrome among u.s. Adults Diabetes Care 2004; 27: 2444-9.

    [5] Nestel P, Lyu R, Low LP, Sheu WH, Nitiyanant W, Saito I, et al. Metabolic syndrome: recent prevalence in East and Southeast Asian populations. Asia Pac J Clin Nutr 2007; 16: 362-7.

    [6] Athyros VG, Bouloukos VI, Pehlivanidis AN, Papageorgiou AA, Dionysopoulou SG, Symeonidis AN, et al. The prevalence of the metabolic syndrome in Greece: the Metabolic syndrome-Greece Multicentre Study. Diabetes Obes Metab 2005; 7: 397-405.

    [7] Okamoto Y, Kihara S, Funahashi T, Matsuzawa Y, Libby P. Adiponectin: a key adipocytokine in metabolic syndrome. Clin Sci (Lond) 2006; 110: 267-78.

    [8] Francke S, Manraj M, Lacquemant C, Lecoeur C, Lepretre F, Passa P, et al. A genome-wide scan for coronary heart disease suggests in Indo-Mauritians a susceptibility locus on chromosome 16p13 and replicates linkage with the metabolic syndrome on 3q27. Hum Mol Genet 2001; 10: 2751-65.

    [9] Lehmann JM, Moore LB, Smith-Oliver TA, Wilkison WO, Willson TM, Kliewer SA. An antidiabetic thiazolidinedione is a high affinity ligand for peroxisome proliferator-activated receptor gamma (PPAR gamma) J Biol Chem 1995; 270: 12953-6.

    [10] Iwaki M, Matsuda M, Maeda N, Funahashi T, Matsuzawa Y, Makishima M, et al. Induction of adiponectin, a fat-derived antidiabetic and antiatherogenic factor, by nuclear receptors. Diabetes 2003; 52: 1655-63.

    [11] Lu Y, Ye X, Cao Y, Li Q, Yu X, Cheng J, et al. Genetic variants in peroxisome proliferator-activated receptorgamma and retinoid X receptor-alpha gene and type2 diabetes risk, a case-control study of a Chinese Han population. Diabetes Technol Ther 2011; 13: 157-64.

    [12] Shi H, Yu X, Li Q, Ye X, Gao Y, Ma J, et al. Association between PPAR-γ and RXR-α gene polymorphism and metabolic syndrome risk: a case-control study of a Chinese Han population. Arch Med Res 2012; 43: 233-42.

    [13] Alberti KG, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, et al. Harmonizing the metabolic syndrome: a joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation 2009; 120: 1640-5.

    [14] Lou XY, Chen GB, Yan L, Ma JZ, Zhu J, Elston RC, et al. A generalized combinatorial approach for detecting gene-by-gene and gene-by-environment interactions with application to nicotine dependence. Am J Hum Genet 2007; 80: 1125-37.

    [15] Mukherjee A, Deshpande JM. Application of artificial neural networks in structural design expert s- ystems. Comput Struct 1995; 54: 367-75.

    [16] Bragt MC, Popeijus HE. Peroxisome proliferator-activated receptors and the metabolic syndrome. Physiol Behav 2008; 94: 187-97.

    [17] Meirhaeghe A, Cottel D, Amouyel P, Dallongeville J. Dallongeville, Association between peroxisome proliferator-activated receptor gamma haplotypes and the metabolic syndrome in French men and women. Diabetes 2005; 54: 3043-8.

    [18] Guo W, Che Z, Xu A, Zhou Z. Angiopoietin-like protein 3 and adiponectin levels in patients with metabolic syndrome. Zhong Nan Da Xue Xue Bao Yi Xue Ban (China) 2010; 35: 203-8.

    [19] Bartosch-Harlid A, Andersson R. Diabetes mellitus in pancreatic cancer and the need for diagnosis of asymptomatic disease. Pancreatology 2010; 10: 423-8.

    [20] Coppedè F, Grossi E, Migheli F, Migliore L. Polymorphisms in folate-metabolizing genes, chromosome damage, and risk of Down syndrome in Italian women: identification of key factors using artificial neural networks. BMC Med. Genomics 2010; 3: 42.

    [21] Shi H, Lu Y, Du J, Du W, Ye X, Yu X, et al. Application of Back Propagation Artificial Neural Network on Genetic Variants in Adiponectin ADIPOQ, Peroxisome Proliferator-Activated Receptor-gamma, and Retinoid X Receptor-alpha Genes and Type 2 Diabetes Risk in a Chinese Han Population. Diabetes Technol Ther 2011; 14: 293-300.

    [22] Lee IT, Chiu Y, Hwu C, He C, Chiang F, Lin Y, et al. Central obesity is important but not essential component of the metabolic syndrome for predicting diabetes mellitus in a hypertensive family-based cohort. Results from the Stanford Asia-Pacific Program for Hypertension and Insulin Resistance (SAPPHIRe) Taiwan follow-up study. Cardiovasc Diabetol 2012; 11: 43.

    [23] Zhuo Q, Wang ZQ, Fu P, Piao JH, Tian Y, Xu J, et al. Association between adiponectin and metabolic syndrome in older adults from major cities of China. Biomed. Environ Sci 2010; 23: 53-61.

    [24] Yang LL, Hua Q, Liu RK, Yang Z. Association between two common polymorphisms of PPARgamma gene and metabolic syndrome families in a Chinese population. Arch. Med Res 2009; 40: 89-96.

    [25] Jaziri R, Lobbens S, Aubert R, Pean F, Lahmidi S, Vaxillaire M, et al. The PPARG Pro12Ala polymorphism is associated with a decreased risk of developing hyperglycemia over 6 years and combines with the effect of the APM1 G-11391A single nucleotide polymorphism: the Data From an Epidemiological Study on the Insulin Resistance Syndrome (DESIR) study. Diabetes 2006; 55: 1157-62.

    [26] Prakash J, Srivastava N, Awasthi S, Agarwal C, Natu S, Rajpal N, et al. Association of PPAR-gamma gene polymorphisms with obesity and obesity-associated phenotypes in north indian population. Am J Hum Biol 2012; 24: 454-9.

    [27] Dongxia L, Qi H, Lisong L, Jincheng G. Association of peroxisome proliferator-activated receptorgamma gene Pro12Ala and C161T polymorphisms with metabolic syndrome. Circ J 2008; 72: 551-7.

    [28] Suyama Y, Matsuda C, Isomura M, Hamano T, Karino K, Yamasaki M, et al. Effects of six functional SNPs on the urinary 8-isoprostane level in a general Japanese population; Shimane COHRE Study. Dis Markers 2011; 30: 291-8.

    [29] Curtis D. Comparison of artificial neural network analysis with other multimarker methods for detecting genetic association. BMC Genet 2007; 8: 49.

    [30] Eftekhar B, Mohammad K, Ardebili HE, Ghodsi M, Ketabchi E. Comparison of artificial neural network and logistic regression models for prediction of mortality in head trauma based on initial clinical data. BMC Med Inform Decis Mak 2005; 5: 3.

    [31] Jack V. T. Advantages and disadvantages of using artificial neural networks versus logistic regression for predicting medical outcomes. J Clin Epidemiol 1996; 49: 1225-31.

    Received 03 October 2012, Revised 06 November 2012, Accepted 16 November 2012, Epub 20 March 2013

    The study was supported by the National Natural Science Foundation of China Grant No. 30771858, Jiangsu Provincial Natural Science Foundation Grant No. BK2007229 and Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

    △These authors contributed equally to this study.

    The authors reported no conflict of interests.

    10.7555/JBR.27.20120061

    国产高清激情床上av| 国产亚洲av嫩草精品影院| 亚洲一区高清亚洲精品| 99热这里只有是精品50| 久久国产精品影院| 色吧在线观看| 久久亚洲真实| 色综合欧美亚洲国产小说| 午夜日韩欧美国产| 久久久久久九九精品二区国产| 免费电影在线观看免费观看| 男插女下体视频免费在线播放| 亚洲欧美精品综合久久99| 18+在线观看网站| 12—13女人毛片做爰片一| 国产私拍福利视频在线观看| 欧美激情在线99| 97超级碰碰碰精品色视频在线观看| 国产激情欧美一区二区| 每晚都被弄得嗷嗷叫到高潮| 淫秽高清视频在线观看| 操出白浆在线播放| 精品一区二区三区av网在线观看| 成人特级黄色片久久久久久久| 男人舔女人下体高潮全视频| 嫩草影院入口| 亚洲不卡免费看| 欧美日韩瑟瑟在线播放| 久久精品夜夜夜夜夜久久蜜豆| 日本免费a在线| 久久久久性生活片| 色综合婷婷激情| 老熟妇仑乱视频hdxx| 五月玫瑰六月丁香| 精品熟女少妇八av免费久了| 色综合婷婷激情| 大型黄色视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 女人被狂操c到高潮| 午夜日韩欧美国产| 人人妻人人澡欧美一区二区| 亚洲欧美一区二区三区黑人| 欧美中文综合在线视频| 久久国产精品影院| 欧美黑人欧美精品刺激| 麻豆成人午夜福利视频| 午夜免费男女啪啪视频观看 | a级毛片a级免费在线| av女优亚洲男人天堂| 午夜精品一区二区三区免费看| 变态另类成人亚洲欧美熟女| 国产精品久久久人人做人人爽| 免费搜索国产男女视频| 欧美另类亚洲清纯唯美| 久久国产精品影院| 久久精品影院6| 午夜老司机福利剧场| 亚洲国产精品久久男人天堂| 免费无遮挡裸体视频| 91av网一区二区| 听说在线观看完整版免费高清| 三级毛片av免费| 国产av麻豆久久久久久久| 国产极品精品免费视频能看的| 丰满的人妻完整版| 好男人在线观看高清免费视频| 久久久久久久亚洲中文字幕 | 国产欧美日韩一区二区三| 一边摸一边抽搐一进一小说| 淫妇啪啪啪对白视频| 国产视频一区二区在线看| www.色视频.com| 一区二区三区免费毛片| 欧美最黄视频在线播放免费| 久久精品人妻少妇| 中文字幕熟女人妻在线| 亚洲第一欧美日韩一区二区三区| 亚洲色图av天堂| 精品国产三级普通话版| 国产精华一区二区三区| 99热这里只有是精品50| 久久久久亚洲av毛片大全| 90打野战视频偷拍视频| 国产视频内射| 成人国产综合亚洲| 亚洲精品色激情综合| 男人和女人高潮做爰伦理| 国内精品久久久久精免费| 老汉色∧v一级毛片| 天堂√8在线中文| 别揉我奶头~嗯~啊~动态视频| 啦啦啦韩国在线观看视频| 国产高潮美女av| 国产三级黄色录像| 欧美成人免费av一区二区三区| 国内精品一区二区在线观看| 国产91精品成人一区二区三区| 十八禁人妻一区二区| 九色国产91popny在线| 日韩中文字幕欧美一区二区| 色综合站精品国产| 国产v大片淫在线免费观看| 国产精品乱码一区二三区的特点| 国产蜜桃级精品一区二区三区| 一边摸一边抽搐一进一小说| 一级毛片高清免费大全| 一级a爱片免费观看的视频| 女人十人毛片免费观看3o分钟| 亚洲男人的天堂狠狠| 欧美最新免费一区二区三区 | 色综合站精品国产| 中出人妻视频一区二区| 小说图片视频综合网站| 亚洲精品粉嫩美女一区| 成年女人毛片免费观看观看9| 欧美日韩瑟瑟在线播放| 久久中文看片网| 久久九九热精品免费| 免费观看精品视频网站| 日韩精品青青久久久久久| 黄色视频,在线免费观看| 亚洲精品456在线播放app | 亚洲一区二区三区不卡视频| 一a级毛片在线观看| 中亚洲国语对白在线视频| 一本精品99久久精品77| 国产精品久久久久久亚洲av鲁大| 亚洲国产精品成人综合色| 在线播放国产精品三级| 91字幕亚洲| av在线天堂中文字幕| 特大巨黑吊av在线直播| 精品一区二区三区视频在线观看免费| 亚洲精品粉嫩美女一区| 成人亚洲精品av一区二区| 国产亚洲精品久久久久久毛片| 十八禁网站免费在线| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 国产精品,欧美在线| 最近在线观看免费完整版| 波多野结衣巨乳人妻| 欧美日本亚洲视频在线播放| 嫁个100分男人电影在线观看| 看免费av毛片| 亚洲精品成人久久久久久| 中亚洲国语对白在线视频| 午夜精品一区二区三区免费看| 神马国产精品三级电影在线观看| 天天躁日日操中文字幕| 天堂av国产一区二区熟女人妻| 麻豆成人午夜福利视频| 久久久久国产精品人妻aⅴ院| 国产精品,欧美在线| 亚洲人与动物交配视频| 免费无遮挡裸体视频| 波多野结衣高清作品| 少妇人妻精品综合一区二区 | 在线a可以看的网站| 长腿黑丝高跟| 亚洲av不卡在线观看| 国产老妇女一区| 少妇的逼好多水| 夜夜夜夜夜久久久久| 色老头精品视频在线观看| 国产欧美日韩一区二区三| 欧美日韩一级在线毛片| 国产精品一区二区三区四区免费观看 | 免费大片18禁| 午夜精品一区二区三区免费看| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 夜夜躁狠狠躁天天躁| 日本五十路高清| 欧美在线一区亚洲| xxx96com| 亚洲精品成人久久久久久| 禁无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 国产免费av片在线观看野外av| 99在线视频只有这里精品首页| 又黄又粗又硬又大视频| 国产精品一区二区三区四区久久| 99久国产av精品| 九九在线视频观看精品| 超碰av人人做人人爽久久 | 国产成人系列免费观看| 熟女人妻精品中文字幕| 国产精品1区2区在线观看.| 天堂√8在线中文| 久久草成人影院| 国产视频一区二区在线看| 中文字幕熟女人妻在线| 国产色婷婷99| 亚洲片人在线观看| 久久久久性生活片| 人妻夜夜爽99麻豆av| 欧美一区二区亚洲| 免费观看精品视频网站| 搡女人真爽免费视频火全软件 | 免费av不卡在线播放| 99在线人妻在线中文字幕| 欧美黄色淫秽网站| 日本免费一区二区三区高清不卡| 欧美xxxx黑人xx丫x性爽| 日日摸夜夜添夜夜添小说| 亚洲精品成人久久久久久| 1024手机看黄色片| 欧美乱色亚洲激情| 丰满乱子伦码专区| 国产色爽女视频免费观看| 精品人妻1区二区| www.www免费av| 麻豆一二三区av精品| 男人舔女人下体高潮全视频| 18禁黄网站禁片午夜丰满| 国产精品98久久久久久宅男小说| 久久久久久大精品| 国产精品 欧美亚洲| ponron亚洲| 女人十人毛片免费观看3o分钟| 成人欧美大片| 美女高潮喷水抽搐中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 人妻夜夜爽99麻豆av| 欧美激情在线99| 亚洲av中文字字幕乱码综合| 色吧在线观看| 久久草成人影院| 日韩欧美 国产精品| 两个人视频免费观看高清| 亚洲欧美日韩高清在线视频| 丁香欧美五月| 精品无人区乱码1区二区| 国产成+人综合+亚洲专区| 国产精品一区二区三区四区久久| 欧美成人免费av一区二区三区| 91久久精品电影网| 久久久久九九精品影院| 麻豆一二三区av精品| 国产欧美日韩一区二区三| 成年免费大片在线观看| 伊人久久大香线蕉亚洲五| 女生性感内裤真人,穿戴方法视频| 夜夜夜夜夜久久久久| 国产99白浆流出| 日本免费a在线| 免费观看精品视频网站| 日韩欧美三级三区| 日韩欧美在线二视频| 国产精品久久久久久亚洲av鲁大| 在线观看日韩欧美| 亚洲中文字幕日韩| 欧美xxxx黑人xx丫x性爽| 香蕉久久夜色| 久久精品综合一区二区三区| 一本精品99久久精品77| 在线观看一区二区三区| 国产69精品久久久久777片| 小说图片视频综合网站| 欧美黑人巨大hd| 欧美日韩国产亚洲二区| 国产精品三级大全| 精品久久久久久久久久免费视频| 美女免费视频网站| 无人区码免费观看不卡| 欧美日韩黄片免| 两个人视频免费观看高清| 午夜免费男女啪啪视频观看 | av专区在线播放| av国产免费在线观看| 热99re8久久精品国产| 中文亚洲av片在线观看爽| 午夜亚洲福利在线播放| 亚洲午夜理论影院| 亚洲精品色激情综合| 熟女少妇亚洲综合色aaa.| 亚洲真实伦在线观看| 亚洲最大成人手机在线| 成年版毛片免费区| 嫩草影院精品99| 亚洲在线自拍视频| 欧美午夜高清在线| 亚洲久久久久久中文字幕| 亚洲乱码一区二区免费版| 好男人在线观看高清免费视频| 婷婷精品国产亚洲av| 亚洲性夜色夜夜综合| 夜夜躁狠狠躁天天躁| 国产蜜桃级精品一区二区三区| 国产单亲对白刺激| 亚洲真实伦在线观看| 别揉我奶头~嗯~啊~动态视频| 女同久久另类99精品国产91| 一个人观看的视频www高清免费观看| 搡老熟女国产l中国老女人| 变态另类丝袜制服| 亚洲人成网站高清观看| 久久久久久久久大av| 成人亚洲精品av一区二区| 国产男靠女视频免费网站| 精品熟女少妇八av免费久了| 国产精品综合久久久久久久免费| 男人和女人高潮做爰伦理| 在线a可以看的网站| 亚洲精品一区av在线观看| 男人和女人高潮做爰伦理| 久久国产精品影院| 中国美女看黄片| 日本熟妇午夜| 日本免费a在线| 国产亚洲精品综合一区在线观看| 亚洲国产欧美人成| 欧美最黄视频在线播放免费| 九色成人免费人妻av| 国产高清三级在线| 宅男免费午夜| 日韩欧美国产在线观看| 亚洲精品成人久久久久久| 日本撒尿小便嘘嘘汇集6| 国产乱人视频| 久99久视频精品免费| 国产精品久久久久久久电影 | 91字幕亚洲| 最近在线观看免费完整版| 全区人妻精品视频| 久久性视频一级片| 夜夜看夜夜爽夜夜摸| 国产精品 国内视频| 国产免费男女视频| 狠狠狠狠99中文字幕| 91九色精品人成在线观看| 亚洲国产欧美人成| 日韩大尺度精品在线看网址| 丰满人妻一区二区三区视频av | 无遮挡黄片免费观看| 热99re8久久精品国产| 日本黄大片高清| 一区二区三区免费毛片| 亚洲精品一卡2卡三卡4卡5卡| 俄罗斯特黄特色一大片| 国产亚洲精品av在线| 美女高潮喷水抽搐中文字幕| 男女那种视频在线观看| 国产精品久久久人人做人人爽| 色av中文字幕| 丁香欧美五月| 久久久国产成人精品二区| 波多野结衣高清无吗| 久久午夜亚洲精品久久| 亚洲人成网站在线播放欧美日韩| 国产精品 国内视频| 久久精品国产亚洲av涩爱 | www.999成人在线观看| svipshipincom国产片| 国产午夜精品久久久久久一区二区三区 | 禁无遮挡网站| 亚洲国产高清在线一区二区三| 免费人成视频x8x8入口观看| 免费看光身美女| 精品久久久久久久毛片微露脸| 国产精品野战在线观看| 久久精品综合一区二区三区| 丁香六月欧美| 九色成人免费人妻av| 欧美一级a爱片免费观看看| 色吧在线观看| 男女那种视频在线观看| 亚洲18禁久久av| 在线播放国产精品三级| 久久久久久久久大av| 欧美国产日韩亚洲一区| 丰满的人妻完整版| xxx96com| 亚洲精品国产精品久久久不卡| 日日摸夜夜添夜夜添小说| 成年女人永久免费观看视频| 午夜免费男女啪啪视频观看 | 亚洲专区国产一区二区| 中国美女看黄片| 国产精品亚洲一级av第二区| 久久久色成人| 91久久精品电影网| 老司机午夜十八禁免费视频| 亚洲国产精品sss在线观看| 首页视频小说图片口味搜索| 国产精品香港三级国产av潘金莲| 亚洲精品乱码久久久v下载方式 | 真人做人爱边吃奶动态| 国产一区二区三区在线臀色熟女| 欧美绝顶高潮抽搐喷水| 99久久无色码亚洲精品果冻| 在线天堂最新版资源| a在线观看视频网站| 欧美区成人在线视频| 一边摸一边抽搐一进一小说| 身体一侧抽搐| 18禁黄网站禁片午夜丰满| 日韩欧美国产一区二区入口| 免费看日本二区| 亚洲精品日韩av片在线观看 | 最好的美女福利视频网| 久久精品国产99精品国产亚洲性色| 久久人人精品亚洲av| 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看| 亚洲在线观看片| 十八禁网站免费在线| 美女被艹到高潮喷水动态| 九九热线精品视视频播放| 久久久久久久久中文| 国内毛片毛片毛片毛片毛片| 国内精品一区二区在线观看| 日韩高清综合在线| 日本黄大片高清| 操出白浆在线播放| 亚洲美女黄片视频| 亚洲不卡免费看| 黄片小视频在线播放| 日韩大尺度精品在线看网址| 噜噜噜噜噜久久久久久91| 91在线观看av| 日日干狠狠操夜夜爽| 国产精品99久久99久久久不卡| 久久天躁狠狠躁夜夜2o2o| 精品不卡国产一区二区三区| 欧美成狂野欧美在线观看| 国产精品综合久久久久久久免费| 成年免费大片在线观看| 久久亚洲精品不卡| 淫妇啪啪啪对白视频| 国产亚洲欧美在线一区二区| 久久精品国产亚洲av涩爱 | 久久久久精品国产欧美久久久| 99热只有精品国产| av专区在线播放| 乱人视频在线观看| 久久久久国产精品人妻aⅴ院| 国产精品久久视频播放| 1000部很黄的大片| 色视频www国产| 国产精品香港三级国产av潘金莲| 黑人欧美特级aaaaaa片| 成年女人永久免费观看视频| 搞女人的毛片| 最近最新免费中文字幕在线| 波多野结衣巨乳人妻| 99精品欧美一区二区三区四区| 脱女人内裤的视频| 最近视频中文字幕2019在线8| 亚洲最大成人中文| 中国美女看黄片| 国产免费一级a男人的天堂| 精品无人区乱码1区二区| 美女免费视频网站| 两个人看的免费小视频| 久久久久久国产a免费观看| 99久久99久久久精品蜜桃| 国产高清有码在线观看视频| 一区二区三区激情视频| bbb黄色大片| 久久人人精品亚洲av| 老司机福利观看| 一进一出好大好爽视频| 麻豆一二三区av精品| 一夜夜www| 在线观看免费视频日本深夜| 在线观看日韩欧美| 美女免费视频网站| 两个人看的免费小视频| 淫秽高清视频在线观看| 亚洲av五月六月丁香网| 国产欧美日韩一区二区精品| 亚洲欧美日韩无卡精品| 黄色丝袜av网址大全| 亚洲精品粉嫩美女一区| 成人高潮视频无遮挡免费网站| 性色avwww在线观看| 每晚都被弄得嗷嗷叫到高潮| 美女高潮的动态| 99久久精品热视频| 精品国内亚洲2022精品成人| 精品免费久久久久久久清纯| 99视频精品全部免费 在线| 国产精华一区二区三区| 国产一级毛片七仙女欲春2| 在线观看日韩欧美| 亚洲精品成人久久久久久| 最近最新中文字幕大全电影3| 丰满乱子伦码专区| 老司机在亚洲福利影院| 真人一进一出gif抽搐免费| 波多野结衣高清无吗| 精品久久久久久成人av| 九色成人免费人妻av| 亚洲一区二区三区不卡视频| 俺也久久电影网| 此物有八面人人有两片| 国产精品 欧美亚洲| 19禁男女啪啪无遮挡网站| 天天一区二区日本电影三级| 俺也久久电影网| 波多野结衣高清无吗| 内地一区二区视频在线| 国产一区二区激情短视频| www日本在线高清视频| 亚洲黑人精品在线| 色哟哟哟哟哟哟| 一本精品99久久精品77| 久久久久久国产a免费观看| 国内精品一区二区在线观看| 国产欧美日韩一区二区精品| 五月伊人婷婷丁香| 免费在线观看日本一区| 天美传媒精品一区二区| 夜夜看夜夜爽夜夜摸| 国产日本99.免费观看| 99久久精品热视频| 亚洲av成人精品一区久久| 偷拍熟女少妇极品色| 久久香蕉精品热| 在线看三级毛片| 色综合亚洲欧美另类图片| 看免费av毛片| 欧美在线一区亚洲| 日韩免费av在线播放| 精品99又大又爽又粗少妇毛片 | 久久久成人免费电影| 小蜜桃在线观看免费完整版高清| 最近最新免费中文字幕在线| 成人国产一区最新在线观看| 熟女少妇亚洲综合色aaa.| 亚洲熟妇熟女久久| 亚洲最大成人中文| 国产真实乱freesex| 91av网一区二区| 国产成人av激情在线播放| 性色av乱码一区二区三区2| 欧美日韩黄片免| 91麻豆av在线| 麻豆国产97在线/欧美| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩高清专用| 丁香欧美五月| 免费人成在线观看视频色| 99国产精品一区二区蜜桃av| 日韩人妻高清精品专区| 亚洲人成网站在线播| 变态另类丝袜制服| 欧美日韩精品网址| 亚洲成人久久爱视频| 婷婷六月久久综合丁香| 亚洲成av人片在线播放无| 久久久久九九精品影院| 亚洲国产精品sss在线观看| 亚洲欧美日韩高清在线视频| 99国产精品一区二区蜜桃av| 亚洲 国产 在线| 99久久成人亚洲精品观看| 90打野战视频偷拍视频| 国产亚洲精品久久久久久毛片| 最新在线观看一区二区三区| 精品乱码久久久久久99久播| 国产成+人综合+亚洲专区| 欧美乱妇无乱码| eeuss影院久久| 天堂av国产一区二区熟女人妻| 国产精品久久视频播放| 男女午夜视频在线观看| 真实男女啪啪啪动态图| 制服丝袜大香蕉在线| 两个人的视频大全免费| 国产久久久一区二区三区| 色综合站精品国产| 亚洲人成网站在线播放欧美日韩| 少妇的丰满在线观看| 亚洲内射少妇av| 免费av不卡在线播放| e午夜精品久久久久久久| 免费电影在线观看免费观看| 一个人看视频在线观看www免费 | www日本黄色视频网| 高潮久久久久久久久久久不卡| 免费人成视频x8x8入口观看| 久久久久亚洲av毛片大全| 日本一本二区三区精品| 欧美激情在线99| 国产精品久久久久久亚洲av鲁大| 男女那种视频在线观看| 亚洲av熟女| 波多野结衣高清无吗| 日本一本二区三区精品| 亚洲aⅴ乱码一区二区在线播放| 国产爱豆传媒在线观看| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 国产乱人视频| xxx96com| 免费搜索国产男女视频| 女人高潮潮喷娇喘18禁视频| 怎么达到女性高潮| 最近最新免费中文字幕在线| 无遮挡黄片免费观看| 99久久精品一区二区三区| a在线观看视频网站| 3wmmmm亚洲av在线观看| 18+在线观看网站| 亚洲av免费高清在线观看| 久久亚洲精品不卡| 国产伦在线观看视频一区| 99精品久久久久人妻精品| 男女做爰动态图高潮gif福利片| 国产亚洲精品av在线| 99在线视频只有这里精品首页| 日韩欧美一区二区三区在线观看|