首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 45 毫秒
1.
A recent genome-wide association study (GWAS) of subjects from Japan and South Korea reported a novel association between the TP63 locus on chromosome 3q28 and risk of lung adenocarcinoma (p = 7.3 × 10(-12)); however, this association did not achieve genome-wide significance (p ≤ 10(-7)) among never-smoking males or females. To determine if this association with lung cancer risk is independent of tobacco use, we genotyped the TP63 SNPs reported by the previous GWAS (rs10937405 and rs4488809) in 3,467 never-smoking female lung cancer cases and 3,787 never-smoking female controls from 10 studies conducted in Taiwan, Mainland China, South Korea, and Singapore. Genetic variation in rs10937405 was associated with risk of lung adenocarcinoma [n = 2,529 cases; p = 7.1 × 10(-8); allelic risk = 0.80, 95% confidence interval (CI) = 0.74-0.87]. There was also evidence of association with squamous cell carcinoma of the lung (n = 302 cases; p = 0.037; allelic risk = 0.82, 95% CI = 0.67-0.99). Our findings provide strong evidence that genetic variation in TP63 is associated with the risk of lung adenocarcinoma among Asian females in the absence of tobacco smoking.  相似文献   

2.
BackgroundObservational studies examining associations between adult height and risk of colorectal, prostate, and lung cancers have generated mixed results. We conducted meta-analyses using data from prospective cohort studies and further carried out Mendelian randomization analyses, using height-associated genetic variants identified in a genome-wide association study (GWAS), to evaluate the association of adult height with these cancers.ConclusionsOur study provides evidence for a potential causal association of adult height with the risk of colorectal and lung cancers and suggests that certain genetic factors and biological pathways affecting adult height may also affect the risk of these cancers.  相似文献   

3.
Recent evidence suggests that inflammation plays a pivotal role in the development of lung cancer. In this study, we used a two-stage approach to investigate associations between genetic variants in inflammation pathways and lung cancer risk based on genome-wide association study (GWAS) data. A total of 7,650 sequence variants from 720 genes relevant to inflammation pathways were identified using keyword and pathway searches from Gene Cards and Gene Ontology databases. In Stage 1, six GWAS datasets from the International Lung Cancer Consortium were pooled (4,441 cases and 5,094 controls of European ancestry), and a hierarchical modeling (HM) approach was used to incorporate prior information for each of the variants into the analysis. The prior matrix was constructed using (1) role of genes in the inflammation and immune pathways; (2) physical properties of the variants including the location of the variants, their conservation scores and amino acid coding; (3) LD with other functional variants and (4) measures of heterogeneity across the studies. HM affected the priority ranking of variants particularly among those having low prior weights, imprecise estimates and/or heterogeneity across studies. In Stage 2, we used an independent NCI lung cancer GWAS study (5,699 cases and 5,818 controls) for in silico replication. We identified one novel variant at the level corrected for multiple comparisons (rs2741354 in EPHX2 at 8q21.1 with p value = 7.4 × 10?6), and confirmed the associations between TERT (rs2736100) and the HLA region and lung cancer risk. HM allows for prior knowledge such as from bioinformatic sources to be incorporated into the analysis systematically, and it represents a complementary analytical approach to the conventional GWAS analysis.  相似文献   

4.
To date, genome-wide association studies have identified thousands of statistically-significant associations between genetic variants, and phenotypes related to a myriad of traits and diseases. A key goal for human-genetics research is to translate these associations into functional mechanisms. Popular gene-set analysis tools, like MAGMA, map variants to genes they might affect, and then integrate genome-wide association study data (that is, variant-level associations for a phenotype) to score genes for association with a phenotype. Gene scores are subsequently used in competitive gene-set analyses to identify biological processes that are enriched for phenotype association. By default, variants are mapped to genes in their proximity. However, many variants that affect phenotypes are thought to act at regulatory elements, which can be hundreds of kilobases away from their target genes. Thus, we explored the idea of augmenting a proximity-based mapping scheme with publicly-available datasets of regulatory interactions. We used MAGMA to analyze genome-wide association study data for ten different phenotypes, and evaluated the effects of augmentation by comparing numbers, and identities, of genes and gene sets detected as statistically significant between mappings. We detected several pitfalls and confounders of such “augmented analyses”, and introduced ways to control for them. Using these controls, we demonstrated that augmentation with datasets of regulatory interactions only occasionally strengthened the enrichment for phenotype association amongst (biologically-relevant) gene sets for different phenotypes. Still, in such cases, genes and regulatory elements responsible for the improvement could be pinpointed. For instance, using brain regulatory-interactions for augmentation, we were able to implicate two acetylcholine receptor subunits involved in post-synaptic chemical transmission, namely CHRNB2 and CHRNE, in schizophrenia. Collectively, our study presents a critical approach for integrating regulatory interactions into gene-set analyses for genome-wide association study data, by introducing various controls to distinguish genuine results from spurious discoveries.  相似文献   

5.
Frozen shoulder is a painful condition that often requires surgery and affects up to 5% of individuals aged 40–60 years. Little is known about the causes of the condition, but diabetes is a strong risk factor. To begin to understand the biological mechanisms involved, we aimed to identify genetic variants associated with frozen shoulder and to use Mendelian randomization to test the causal role of diabetes. We performed a genome-wide association study (GWAS) of frozen shoulder in the UK Biobank using data from 10,104 cases identified from inpatient, surgical and primary care codes. We used data from FinnGen for replication and meta-analysis. We used one-sample and two-sample Mendelian randomization approaches to test for a causal association of diabetes with frozen shoulder. We identified five genome-wide significant loci. The most significant locus (lead SNP rs28971325; OR = 1.20, [95% CI: 1.16–1.24], p = 5x10-29) contained WNT7B. This variant was also associated with Dupuytren’s disease (OR = 2.31 [2.24, 2.39], p<1x10-300) as were a further two of the frozen shoulder associated variants. The Mendelian randomization results provided evidence that type 1 diabetes is a causal risk factor for frozen shoulder (OR = 1.03 [1.02–1.05], p = 3x10-6). There was no evidence that obesity was causally associated with frozen shoulder, suggesting that diabetes influences risk of the condition through glycemic rather than mechanical effects. We have identified genetic loci associated with frozen shoulder. There is a large overlap with Dupuytren’s disease associated loci. Diabetes is a likely causal risk factor. Our results provide evidence of biological mechanisms involved in this common painful condition.  相似文献   

6.
It is a fundamental challenge to discover the association of genetic traits with phenotypic traits. In this study, we aimed to identify possible genetic traits related to horse temperament. Based on previous findings, we selected 71 candidate genes related to temperamental trait and examined them in the human and horse reference genomes (hg38 and equCab2, respectively). We found 16 orthologous genes and, by comparing with the human reference genome, 17 homologous genes in the horse reference genome. We designed probes specific for the 33 horse genes. Using the probes, we built sequencing libraries of the genomic DNA samples from eight aggressive and eight docile horses, and sequenced the constructed libraries using the Illumina Hiseq2500 platform. Through the analysis of the targeted exome sequences, we identified single nucleotide polymorphisms (SNPs) in the genes. SNPs could be served as genetic markers to evaluate aggressive or docile levels of horses. To examine whether any genetic variants are associated with horse temperament, we performed genome-wide association study (GWAS) using the SNP data. GWAS analysis identified ten variants (p-value?<0.05) which could be related to horse temperament. We validated the variants using Sanger sequencing. The most significant variants were found in MAOA (c.1164+41T>C) and AR (c.1047+27G>T) genes with 8.09?×?10?4 p-value. We suggest that the variants might be used to assess horse temperament and to determine superior horses for riding or racing.  相似文献   

7.
Over a decade of genome-wide association studies (GWAS) have led to the finding of extreme polygenicity of complex traits. The phenomenon that “all genes affect every complex trait” complicates Mendelian Randomization (MR) studies, where natural genetic variations are used as instruments to infer the causal effect of heritable risk factors. We reexamine the assumptions of existing MR methods and show how they need to be clarified to allow for pervasive horizontal pleiotropy and heterogeneous effect sizes. We propose a comprehensive framework GRAPPLE to analyze the causal effect of target risk factors with heterogeneous genetic instruments and identify possible pleiotropic patterns from data. By using GWAS summary statistics, GRAPPLE can efficiently use both strong and weak genetic instruments, detect the existence of multiple pleiotropic pathways, determine the causal direction and perform multivariable MR to adjust for confounding risk factors. With GRAPPLE, we analyze the effect of blood lipids, body mass index, and systolic blood pressure on 25 disease outcomes, gaining new information on their causal relationships and potential pleiotropic pathways involved.  相似文献   

8.
Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes.  相似文献   

9.
Aspirin-exacerbated respiratory disease (AERD) is a nonallergic clinical syndrome characterized by a severe decline in forced expiratory volume in one second (FEV1) following the ingestion of non-steroidal anti-inflammatory drugs (NSAIDs) such as aspirin. The effects of genetic variants have not fully explained all of the observed individual differences to an aspirin challenge despite previous attempts to identify AERD-related genes. In the present study, we performed genome-wide association study (GWAS) and targeted association study in Korean asthmatics to identify new genetic factors associated with AERD. A total of 685 asthmatic patients without AERD and 117 subjects with AERD were used for the GWAS of the first stage, and 996 asthmatics without AERD and 142 subjects with AERD were used for a follow-up study. A total of 702 SNPs were genotyped using the GoldenGate assay with the VeraCode microbead. GWAS revealed the top-ranked variants in 3′ regions of the HLA-DPB1 gene. To investigate the detailed genetic effects of an associated region with the risk of AERD, a follow-up targeted association study with the 702 single nucleotide polymorphisms (SNPs) of 14 genes was performed on 802 Korean subjects. In a case–control analysis, HLA-DPB1 rs1042151 (Met105Val) shows the most significant association with the susceptibility of AERD (p = 5.11 × 10?7; OR = 2.40). Moreover, rs1042151 also shows a gene dose for the percent decline of FEV1 after an aspirin challenge (p = 2.82 × 10?7). Our findings show that the HLA-DPB1 gene polymorphism may be the most susceptible genetic factor for the risk of AERD in Korean asthmatics and confirm the importance of HLA-DPB1 in the genetic etiology of AERD.  相似文献   

10.
Genome-wide association studies (GWAS) conducted using commercial single nucleotide polymorphisms (SNP) arrays have proven to be a powerful tool for the detection of common disease susceptibility variants. However, their utility for the detection of lower frequency variants is yet to be practically investigated. Here we describe the application of a rare variant collapsing method to a large genome-wide SNP dataset, the Wellcome Trust Case Control Consortium rheumatoid arthritis (RA) GWAS. We partitioned the data into gene-centric bins and collapsed genotypes of low frequency variants (defined here as MAF ≤0.05) into a single count coupled with univariate analysis. We then prioritised gene regions for further investigation in an independent cohort of 3,355 cases and 2,427 controls based on rare variant signal p value and prior evidence to support involvement in RA. A total of 14,536 gene bins were investigated in the primary analysis and signals mapping to the TNFAIP3 and chr17q24 loci were selected for further investigation. We detected replicating association to low frequency variants in the TNFAIP3 gene (combined p = 6.6 × 10?6). Even though rare variants are not well-represented and can be difficult to genotype in GWAS, our study supports the application of low frequency variant collapsing methods to genome-wide SNP datasets as a means of exploiting data that are routinely ignored.  相似文献   

11.
Following the widespread use of genome-wide association studies (GWAS), focus is turning towards identification of causal variants rather than simply genetic markers of diseases and traits. As a step towards a high-throughput method to identify genome-wide, non-coding, functional regulatory variants, we describe the technique of allele-specific FAIRE, utilising large-scale genotyping technology (FAIRE-gen) to determine allelic effects on chromatin accessibility and regulatory potential. FAIRE-gen was explored using lymphoblastoid cells and the 50,000 SNP Illumina CVD BeadChip. The technique identified an allele-specific regulatory polymorphism within NR1H3 (coding for LXR-α), rs7120118, coinciding with a previously GWAS-identified SNP for HDL-C levels. This finding was confirmed using FAIRE-gen with the 200,000 SNP Illumina Metabochip and verified with the established method of TaqMan allelic discrimination. Examination of this SNP in two prospective Caucasian cohorts comprising 15,000 individuals confirmed the association with HDL-C levels (combined beta = 0.016; p = 0.0006), and analysis of gene expression identified an allelic association with LXR-α expression in heart tissue. Using increasingly comprehensive genotyping chips and distinct tissues for examination, FAIRE-gen has the potential to aid the identification of many causal SNPs associated with disease from GWAS.  相似文献   

12.
Several variations in the nicotinic receptor genes have been identified to be associated with both lung cancer risk and smoking in the genome-wide association (GWA) studies. However, the relationships among these three factors (genetic variants, nicotine dependence, and lung cancer) remain unclear. In an attempt to elucidate these relationships, we applied mediation analysis to quantify the impact of nicotine dependence on the association between the nicotinic receptor genetic variants and lung adenocarcinoma risk. We evaluated 23 single nucleotide polymorphisms (SNPs) in the five nicotinic receptor related genes (CHRNB3, CHRNA6, and CHRNA5/A3/B4) previously reported to be associated with lung cancer risk and smoking behavior and 14 SNPs in the four ‘control’ genes (TERT, CLPTM1L, CYP1A1, and TP53), which were not reported in the smoking GWA studies. A total of 661 lung adenocarcinoma cases and 1,347 controls with a smoking history, obtained from the Environment and Genetics in Lung Cancer Etiology case-control study, were included in the study. Results show that nicotine dependence is a mediator of the association between lung adenocarcinoma and gene variations in the regions of CHRNA5/A3/B4 and accounts for approximately 15% of this relationship. The top two CHRNA3 SNPs associated with the risk for lung adenocarcinoma were rs1051730 and rs12914385 (p-value = 1.9×10−10 and 1.1×10−10, respectively). Also, these two SNPs had significant indirect effects on lung adenocarcinoma risk through nicotine dependence (p = 0.003 and 0.007). Gene variations rs2736100 and rs2853676 in TERT and rs401681 and rs31489 in CLPTM1L had significant direct associations on lung adenocarcinoma without indirect effects through nicotine dependence. Our findings suggest that nicotine dependence plays an important role between genetic variants in the CHRNA5/A3/B4 region, especially CHRNA3, and lung adenocarcinoma. This may provide valuable information for understanding the pathogenesis of lung adenocarcinoma and for conducting personalized smoking cessation interventions.  相似文献   

13.
Large-scale genome-wide association studies (GWAS) have established chromosome 5q31.1 as a susceptibility locus for colorectal cancer (CRC), which was still lack of causal genetic variants. We searched potentially regulatory single nucleotide polymorphisms (SNPs) in the overlap region between linkage disequilibrium (LD) block of 5q31.1 and regulatory elements predicted by histone modifications, then tested their association with CRC via a case-control study. Among three candidate common variants, we found rs17716310 conferred significantly (heterozygous model: OR = 1.273, 95% confidence interval (95%CI) = 1.016–1.595, P = 0.036) and marginally (dominant model: OR = 1.238, 95%CI = 1.000–1.532, P = 0.050) increase risk for CRC in a Chinese population including 695 cases and 709 controls. This variation was suggested to be regulatory altering the activity of enhancer that control PITX1 expression. Using epigenetic information such as chromatin immunoprecipitation-sequencing (ChIP-seq) data might help researchers to interpret the results of GWAS and locate causal variants for diseases in post-GWAS era.  相似文献   

14.
Accurate measurement of the effects of disease status on healthcare costs is important in the pragmatic evaluation of interventions but is complicated by endogeneity bias. Mendelian Randomization, the use of random perturbations in germline genetic variation as instrumental variables, can avoid these limitations. We used a novel Mendelian Randomization analysis to model the causal impact on inpatient hospital costs of liability to six prevalent diseases and health conditions: asthma, eczema, migraine, coronary heart disease, Type 2 diabetes, and depression. We identified genetic variants from replicated genome-wide associations studies and estimated their association with inpatient hospital costs on over 300,000 individuals. There was concordance of findings across varieties of sensitivity analyses, including stratification by sex and methods robust to violations of the exclusion restriction. Results overall were imprecise and we could not rule out large effects of liability to disease on healthcare costs. In particular, genetic liability to coronary heart disease had substantial impacts on costs.  相似文献   

15.
Individuals with fast nicotine metabolism typically smoke more and thus have a greater risk for smoking-induced diseases. Further, the efficacy of smoking cessation pharmacotherapy is dependent on the rate of nicotine metabolism. Our objective was to use nicotine metabolite ratio (NMR), an established biomarker of nicotine metabolism rate, in a genome-wide association study (GWAS) to identify novel genetic variants influencing nicotine metabolism. A heritability estimate of 0.81 (95% CI 0.70–0.88) was obtained for NMR using monozygotic and dizygotic twins of the FinnTwin cohort. We performed a GWAS in cotinine-verified current smokers of three Finnish cohorts (FinnTwin, Young Finns Study, FINRISK2007), followed by a meta-analysis of 1518 subjects, and annotated the genome-wide significant SNPs with methylation quantitative loci (meQTL) analyses. We detected association on 19q13 with 719 SNPs exceeding genome-wide significance within a 4.2 Mb region. The strongest evidence for association emerged for CYP2A6 (min p = 5.77E-86, in intron 4), the main metabolic enzyme for nicotine. Other interesting genes with genome-wide significant signals included CYP2B6, CYP2A7, EGLN2, and NUMBL. Conditional analyses revealed three independent signals on 19q13, all located within or in the immediate vicinity of CYP2A6. A genetic risk score constructed using the independent signals showed association with smoking quantity (p = 0.0019) in two independent Finnish samples. Our meQTL results showed that methylation values of 16 CpG sites within the region are affected by genotypes of the genome-wide significant SNPs, and according to causal inference test, for some of the SNPs the effect on NMR is mediated through methylation. To our knowledge, this is the first GWAS on NMR. Our results enclose three independent novel signals on 19q13.2. The detected CYP2A6 variants explain a strikingly large fraction of variance (up to 31%) in NMR in these study samples. Further, we provide evidence for plausible epigenetic mechanisms influencing NMR.  相似文献   

16.
Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer’s disease (AD). These low estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer’s Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7–11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.  相似文献   

17.
ObjectivesAlthough type 2 diabetes mellitus is a known risk factor for pancreatic cancer, the existence of shared genetic susceptibility is largely unknown. We evaluated whether any reported genetic risk variants of either disease found by genome-wide association studies reciprocally confer susceptibility.MethodsData that were generated in previous genome-wide association studies (GENEVA Type 2 Diabetes; PanScan) were obtained through the National Institutes of Health database of Genotypes and Phenotypes (dbGaP). Using the PanScan datasets, we tested for association of 38 variants within 37 genomic regions known to be susceptibility factors for type 2 diabetes. We further examined whether type 2 diabetes variants predispose to pancreatic cancer risk stratified by diabetes status. Correspondingly, we examined the association of fourteen pancreatic cancer susceptibility variants within eight genomic regions in the GENEVA Type 2 Diabetes dataset.ResultsFour plausible associations of diabetes variants and pancreatic cancer risk were detected at a significance threshold of p = 0.05, and one pancreatic cancer susceptibility variant was associated with diabetes risk at threshold of p = 0.05, but none remained significant after correction for multiple comparisons.ConclusionCurrently identified GWAS susceptibility variants are unlikely to explain the potential shared genetic etiology between Type 2 diabetes and pancreatic cancer.  相似文献   

18.
Genome-wide association studies (GWAS) have successfully identified common genetic variants that contribute to breast cancer risk. Discovering additional variants has become difficult, as power to detect variants of weaker effect with present sample sizes is limited. An alternative approach is to look for variants associated with quantitative traits that in turn affect disease risk. As exposure to high circulating estradiol and testosterone, and low sex hormone-binding globulin (SHBG) levels is implicated in breast cancer etiology, we conducted GWAS analyses of plasma estradiol, testosterone, and SHBG to identify new susceptibility alleles. Cancer Genetic Markers of Susceptibility (CGEMS) data from the Nurses' Health Study (NHS), and Sisters in Breast Cancer Screening data were used to carry out primary meta-analyses among ~1600 postmenopausal women who were not taking postmenopausal hormones at blood draw. We observed a genome-wide significant association between SHBG levels and rs727428 (joint β = -0.126; joint P = 2.09 × 10(-16)), downstream of the SHBG gene. No genome-wide significant associations were observed with estradiol or testosterone levels. Among variants that were suggestively associated with estradiol (P<10(-5)), several were located at the CYP19A1 gene locus. Overall results were similar in secondary meta-analyses that included ~900 NHS current postmenopausal hormone users. No variant associated with estradiol, testosterone, or SHBG at P<10(-5) was associated with postmenopausal breast cancer risk among CGEMS participants. Our results suggest that the small magnitude of difference in hormone levels associated with common genetic variants is likely insufficient to detectably contribute to breast cancer risk.  相似文献   

19.
20.
Colorectal cancer (CRC) is a complex disease, and therefore its development is determined by the combination of both environmental factors and genetic variants. Although genome-wide association studies (GWAS) of SNP variation have conveniently identified 20 genetic variants so far, a significant proportion of the observed heritability is yet to be explained. Common copy-number variants (CNVs) are one of the most important genomic sources of variability, and hence a potential source to explain part of this missing genetic fraction. Therefore, we have performed a GWAS on CNVs to explore the relationship between common structural variation and CRC development. Phase 1 of the GWAS consisted of 881 cases and 667 controls from a Spanish cohort. Copy-number status was validated by quantitative PCR for each of those common CNVs potentially associated with CRC in phase I. Subsequently, SNPs were chosen as proxies for the validated CNVs for phase II replication (1,342 Spanish cases and 1,874 Spanish controls). Four common CNVs were found to be associated with CRC and were further replicated in Phase II. Finally, we found that SNP rs1944682, tagging a 11q11 CNV, was nominally associated with CRC susceptibility (p value = 0.039; OR = 1.122). This locus has been previously related to extreme obesity phenotypes, which could suggest a relationship between body weight and CRC susceptibility.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号