首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
We developed a new statistical framework to find genetic variants associated with extreme longevity. The method, informed GWAS (iGWAS), takes advantage of knowledge from large studies of age-related disease in order to narrow the search for SNPs associated with longevity. To gain support for our approach, we first show there is an overlap between loci involved in disease and loci associated with extreme longevity. These results indicate that several disease variants may be depleted in centenarians versus the general population. Next, we used iGWAS to harness information from 14 meta-analyses of disease and trait GWAS to identify longevity loci in two studies of long-lived humans. In a standard GWAS analysis, only one locus in these studies is significant (APOE/TOMM40) when controlling the false discovery rate (FDR) at 10%. With iGWAS, we identify eight genetic loci to associate significantly with exceptional human longevity at FDR < 10%. We followed up the eight lead SNPs in independent cohorts, and found replication evidence of four loci and suggestive evidence for one more with exceptional longevity. The loci that replicated (FDR < 5%) included APOE/TOMM40 (associated with Alzheimer’s disease), CDKN2B/ANRIL (implicated in the regulation of cellular senescence), ABO (tags the O blood group), and SH2B3/ATXN2 (a signaling gene that extends lifespan in Drosophila and a gene involved in neurological disease). Our results implicate new loci in longevity and reveal a genetic overlap between longevity and age-related diseases and traits, including coronary artery disease and Alzheimer’s disease. iGWAS provides a new analytical strategy for uncovering SNPs that influence extreme longevity, and can be applied more broadly to boost power in other studies of complex phenotypes.  相似文献   

2.
Identifying signatures of selection can provide a straightforward insight into the mechanism of artificial selection and further uncover the causal genes related to the phenotypic variation. Based on Illumina Porcine60KSNP chip data, four complementary methods, Long-Range Haplotype (LRH), Tajima’s D, Cross Population Extend Haplotype Homozygosity Test (XPEHH) and FST, were implemented in this study to detect the selection signatures in the whole genome of one typical Chinese indigenous breed, Rongchang, one Chinese cultivated breed, Songliao, and two western breeds, Landrace and Yorkshire. False Discovery Rate (FDR) was implemented to control the false positive rates. In our study, a total of 159, 127, 179 and 159 candidate selection regions with average length of 0.80 Mb, 0.73 Mb, 0.78 Mb and 0.73 Mb were identified in Landrace, Rongchang, Songliao and Yorkshire, respectively, that span approximately 128.00 Mb, 92.38 Mb, 130.30 Mb and 115.40 Mb and account for approximately 3.74–5.33% of genome across all autosomes. The selection regions of 11.52 Mb shared by Landrace and Yorkshire were the longest when chosen pairs from the pool of the four breeds were examined. The overlaps between Yorkshire and Songliao, approximately 9.20 Mb, were greater than those of Yorkshire and Rongchang. Meanwhile, the overlaps between Landrace and Songliao were greater than those of Landrace and Rongchang but less than those of Songliao and Ronchang. Bioinformatics analysis showed that the genes/QTLs relevant to fertility, coat color, and ear morphology were found in candidate selection regions. Some genes, such as LEMD3, MC1R, KIT, TRHR etc. that were reported under selection, were confirmed in our study, and this analysis also demonstrated the diversity of breeds.  相似文献   

3.
4.
It is recognized that genetic factors contribute to human longevity. Besides the hypothesis of existence of longevity genes, another suggests that a lower frequency of risk alleles decreases the incidence of age-related diseases in the long-lived people. However, the latter finds no support from recent genetic studies. Considering the crucial role of epigenetic modification in gene regulation, we then hypothesize that suppressing disease-related genes in longevity individuals is likely achieved by epigenetic modification, e.g. DNA methylation. To test this hypothesis, we investigated the genome-wide methylation profile in 4 Chinese female centenarians and 4 middle-aged controls using methyl-DNA immunoprecipitation sequencing. 626 differentially methylated regions (DMRs) were observed between both groups. Interestingly, genes with these DMRs were enriched in age-related diseases, including type-2 diabetes, cardiovascular disease, stroke and Alzheimer’s disease. This pattern remains rather stable after including methylomes of two white individuals. Further analyses suggest that the observed DMRs likely have functional roles in regulating disease-associated gene expressions, with some genes [e.g. caspase 3 (CASP3)] being down-regulated whereas the others [i.e. interleukin 1 receptor, type 2 (IL1R2)] up-regulated. Therefore, our study suggests that suppressing the disease-related genes via epigenetic modification is an important contributor to human longevity.  相似文献   

5.
IgM provides a first line of defense during microbial infections. Serum IgM levels are detected routinely in clinical practice. And IgM is a genetically complex trait. We conducted a two-stage genome-wide association study (GWAS) to identify genetic variants affecting serum IgM levels in a Chinese population of 3495, including 1999 unrelated subjects in the first stage and 1496 independent individuals in the second stage. Our data show that a common single nucleotide polymorphism (SNP), rs11552708 located in the TNFSF13 gene was significantly associated with IgM levels (p = 5.00×10−7 in first stage, p = 1.34×10−3 in second stage, and p = 4.22×10−9 when combined). Besides, smoking was identified to be associated with IgM levels in both stages (P<0.05), but there was no significant interaction between smoking and the identified SNP (P>0.05). It is suggested that TNFSF13 may be a susceptibility gene affecting serum IgM levels in Chinese male population.  相似文献   

6.
BACKGROUND: Gastric cancer (GC) is the fifth leading cause of cancer-related deaths worldwide. As an effective and easily performed method, microscopy-based Lauren classification has been widely accepted by gastrointestinal surgeons and pathologists for GC subtyping, but molecular characteristics of different Lauren subtypes were poorly revealed. METHODS: GSE62254 was used as a derivation cohort, and GSE15459 was used as a validation cohort. The difference between diffuse and intestinal GC on the gene expression level was measured. Gene ontology (GO) enrichment analysis was performed for both subgroups. Hierarchical clustering and heatmap exhibition were also performed. Kaplan-Meier plot and Cox proportional hazards model were used to evaluate survival grouped by the given genes or hierarchical clusters. RESULTS: A total of 4598 genes were found differentially expressed between diffuse and intestinal GC. Immunity- and cell adhesion–related GOs were enriched for diffuse GC, whereas DNA repair– and cell cycle–related GOs were enriched for intestinal GC. We proposed a 40-gene signature (χ2 = 30.71, P < .001) that exhibits better discrimination for prognosis than Lauren classification (χ2 = 12.11, P = .002). FRZB [RR (95% CI) = 1.824 (1.115-2.986), P = .017] and EFEMP1 [RR (95% CI) = 1.537 (0.969-2.437), P = .067] were identified as independent prognostic factors only in diffuse GC but not in intestinal GC patients. KRT23 [RR (95% CI) = 1.616 (0.938-2.785), P = .083] was identified as an independent prognostic factor only in intestinal GC patients but not in diffuse GC patients. Similar results were achieved in the validation cohort. CONCLUSION: We found that GCs with different Lauren classifications had different molecular characteristics and identified FRZB, EFEMP1, and KRT23 as subtype-specific prognostic factors for GC patients.  相似文献   

7.
Genetic differences both between individuals and populations are studied for their evolutionary relevance and for their potential medical applications. Most of the genetic differentiation among populations are caused by random drift that should affect all loci across the genome in a similar manner. When a locus shows extraordinary high or low levels of population differentiation, this may be interpreted as evidence for natural selection. The most used measure of population differentiation was devised by Wright and is known as fixation index, or FST. We performed a genome-wide estimation of FST on about 4 millions of SNPs from HapMap project data. We demonstrated a heterogeneous distribution of FST values between autosomes and heterochromosomes. When we compared the FST values obtained in this study with another evolutionary measure obtained by comparative interspecific approach, we found that genes under positive selection appeared to show low levels of population differentiation. We applied a gene set approach, widely used for microarray data analysis, to detect functional pathways under selection. We found that one pathway related to antigen processing and presentation showed low levels of FST, while several pathways related to cell signalling, growth and morphogenesis showed high FST values. Finally, we detected a signature of selection within genes associated with human complex diseases. These results can help to identify which process occurred during human evolution and adaptation to different environments. They also support the hypothesis that common diseases could have a genetic background shaped by human evolution.  相似文献   

8.
Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.  相似文献   

9.
Bacterial sensing of environmental signals plays a key role in regulating virulence and mediating bacterium-host interactions. The sensing of the neuroendocrine stress hormones epinephrine (adrenaline) and norepinephrine (noradrenaline) plays an important role in modulating bacterial virulence. We used MudJ transposon mutagenesis to globally screen for genes regulated by neuroendocrine stress hormones in Salmonella enterica serovar Typhimurium. We identified eight hormone-regulated genes, including yhaK, iroC, nrdF, accC, yedP, STM3081, and the virulence-related genes virK and mig14. The mammalian α-adrenergic receptor antagonist phentolamine reversed the hormone-mediated effects on yhaK, virK, and mig14 but did not affect the other genes. The β-adrenergic receptor antagonist propranolol had no activity in these assays. The virK and mig14 genes are involved in antimicrobial peptide resistance, and phenotypic screens revealed that exposure to neuroendocrine hormones increased the sensitivity of S. Typhimurium to the antimicrobial peptide LL-37. A virK mutant and a virK mig14 double mutant also displayed increased sensitivity to LL-37. In contrast to enterohemorrhagic Escherichia coli (EHEC), we have found no role for the two-component systems QseBC and QseEF in the adrenergic regulation of any of the identified genes. Furthermore, hormone-regulated gene expression could not be blocked by the QseC inhibitor LED209, suggesting that sensing of hormones is mediated through alternative signaling pathways in S. Typhimurium. This study has identified a role for host-derived neuroendocrine stress hormones in downregulating S. Typhimurium virulence gene expression to the benefit of the host, thus providing further insights into the field of host-pathogen communication.Bacterial sensing of environmental signals plays a key role in regulating virulence gene expression and bacterium-host interactions. It is increasingly recognized that detection of host-derived molecules, such as the neuroendocrine stress hormones (catecholamines) epinephrine (adrenaline) and norepinephrine (noradrenaline), plays an important role in modulating bacterial virulence (29, 42).Physical and psychological stress has been linked to increased severity and susceptibility to infection in humans and other animals (23, 42), and epinephrine/norepinephrine levels are an important factor in this. Stress triggers an increase in plasma epinephrine levels (31), and plasma levels of epinephrine and norepinephrine have been reported to increase with patients suffering from postoperative sepsis compared to patients with no complications (32). Administration of norepinephrine and epinephrine to otherwise healthy subjects increases the severity of bacterial infections, including Clostridium perfringens in humans and enterohemorrhagic Escherichia coli (EHEC) in calves (42, 63, 65). Treatment with norepinephrine also increases the virulence of Salmonella enterica serovar Enteritidis in chicks and Salmonella enterica serovar Typhimurium in mice, with a substantial increase in bacterial numbers recovered from the cecum and liver in both cases (47, 65).Norepinephrine is found in large concentrations in the gut due to release by gastrointestinal neurones; indeed up to half the norepinephrine in the body may be produced in the enteric nervous system (ENS) (3). Epinephrine, while not normally found in the gut, is present in the bloodstream and is also produced by macrophages in response to bacteria-derived lipopolysaccharide (LPS) (12, 26). S. Typhimurium is an enteropathogen, can also cross the epithelial barrier to cause systemic infection, and will therefore encounter both these molecules in the normal infection cycle.Phenotypes induced by stress hormones in bacteria include increased adherence of EHEC to bovine intestinal mucosa (63), upregulation of type III secretion and Shiga toxin production in EHEC (22, 60), upregulation of type III secretion in Vibrio parahaemolyticus (51), increase in invasion of epithelial cells and breakdown of epithelial tight junctions by Campylobacter jejuni (15), affected motility and expression of iron uptake genes in S. Typhimurium (8, 9, 36), and modulated virulence in Borrelia burgdorferi (59). Epinephrine and norepinephrine can overcome the growth inhibition of many bacteria, including Salmonella, in serum-containing media (13, 43), due to the ability to act as a siderophore to facilitate iron uptake (13, 28, 47).Norepinephrine and epinephrine also interact with bacterial quorum-sensing (QS) systems. QS is a process of bacterial cell-cell communication in which each cell produces small signal molecules termed “autoinducers” (AIs), which regulate gene expression when a critical threshold concentration and therefore population density have been reached. QS affects diverse processes, including motility, virulence, biofilm formation, type III secretion, and luminescence (6, 64).The EHEC AI-3 QS system is important for motility and expression of the type III secretion system encoded by the locus of enterocyte effacement (LEE) (60). AI-3 sensing and signal transduction are mediated via the QseBC and QseEF two-component systems, respectively. Epinephrine and norepinephrine can substitute for AI-3, causing cross talk between the two signaling systems and induction of type III secretion and motility (57, 60). The sensor kinase QseC is autophosphorylated upon binding either epinephrine or norepinephrine (14), demonstrating the presence of adrenergic receptors in bacteria. These adrenergic phenotypes can also be blocked by the mammalian α- and β-adrenergic antagonists phentolamine and propranolol, although it should be noted that QseC is blocked only by the former (14, 60). This suggests the occurrence of cross talk between bacterial and mammalian cell signaling systems and the existence of multiple bacterial adrenergic sensors.To elucidate the role of host-derived stress hormones in the physiology and pathogenicity of S. Typhimurium, we used MudJ transposon mutagenesis to screen globally for epinephrine- and norepinephrine-regulated genes in S. Typhimurium.  相似文献   

10.
Approximately 4% of English-speaking children are affected by specific language impairment (SLI), a disorder in the development of language skills despite adequate opportunity and normal intelligence. Several studies have indicated the importance of genetic factors in SLI; a positive family history confers an increased risk of development, and concordance in monozygotic twins consistently exceeds that in dizygotic twins. However, like many behavioral traits, SLI is assumed to be genetically complex, with several loci contributing to the overall risk. We have compiled 98 families drawn from epidemiological and clinical populations, all with probands whose standard language scores fall 1.5 SD below the mean for their age. Systematic genomewide quantitative-trait–locus analysis of three language-related measures (i.e., the Clinical Evaluation of Language Fundamentals–Revised [CELF-R] receptive and expressive scales and the nonword repetition [NWR] test) yielded two regions, one on chromosome 16 and one on 19, that both had maximum LOD scores of 3.55. Simulations suggest that, of these two multipoint results, the NWR linkage to chromosome 16q is the most significant, with empirical P values reaching 10−5, under both Haseman-Elston (HE) analysis (LOD score 3.55; P=.00003) and variance-components (VC) analysis (LOD score 2.57; P=.00008). Single-point analyses provided further support for involvement of this locus, with three markers, under the peak of linkage, yielding LOD scores >1.9. The 19q locus was linked to the CELF-R expressive-language score and exceeds the threshold for suggestive linkage under all types of analysis performed—multipoint HE analysis (LOD score 3.55; empirical P=.00004) and VC (LOD score 2.84; empirical P=.00027) and single-point HE analysis (LOD score 2.49) and VC (LOD score 2.22). Furthermore, both the clinical and epidemiological samples showed independent evidence of linkage on both chromosome 16q and chromosome 19q, indicating that these may represent universally important loci in SLI and, thus, general risk factors for language impairment.  相似文献   

11.
Oral infections have a strong ethnic predilection; suggesting that ethnicity is a critical determinant of oral microbial colonization. Dental plaque and saliva samples from 192 subjects belonging to four major ethnicities in the United States were analyzed using terminal restriction fragment length polymorphism (t-RFLP) and 16S pyrosequencing. Ethnicity-specific clustering of microbial communities was apparent in saliva and subgingival biofilms, and a machine-learning classifier was capable of identifying an individual’s ethnicity from subgingival microbial signatures. The classifier identified African Americans with a 100% sensitivity and 74% specificity and Caucasians with a 50% sensitivity and 91% specificity. The data demonstrates a significant association between ethnic affiliation and the composition of the oral microbiome; to the extent that these microbial signatures appear to be capable of discriminating between ethnicities.  相似文献   

12.
13.
14.
A genome-wide association study was performed to identify genetic factors involved in susceptibility to psoriasis (PS) and psoriatic arthritis (PSA), inflammatory diseases of the skin and joints in humans. 223 PS cases (including 91 with PSA) were genotyped with 311,398 single nucleotide polymorphisms (SNPs), and results were compared with those from 519 Northern European controls. Replications were performed with an independent cohort of 577 PS cases and 737 controls from the U.S., and 576 PSA patients and 480 controls from the U.K.. Strongest associations were with the class I region of the major histocompatibility complex (MHC). The most highly associated SNP was rs10484554, which lies 34.7 kb upstream from HLA-C (P = 7.8×10−11, GWA scan; P = 1.8×10−30, replication; P = 1.8×10−39, combined; U.K. PSA: P = 6.9×10−11). However, rs2395029 encoding the G2V polymorphism within the class I gene HCP5 (combined P = 2.13×10−26 in U.S. cases) yielded the highest ORs with both PS and PSA (4.1 and 3.2 respectively). This variant is associated with low viral set point following HIV infection and its effect is independent of rs10484554. We replicated the previously reported association with interleukin 23 receptor and interleukin 12B (IL12B) polymorphisms in PS and PSA cohorts (IL23R: rs11209026, U.S. PS, P = 1.4×10−4; U.K. PSA: P = 8.0×10−4; IL12B:rs6887695, U.S. PS, P = 5×10−5 and U.K. PSA, P = 1.3×10−3) and detected an independent association in the IL23R region with a SNP 4 kb upstream from IL12RB2 (P = 0.001). Novel associations replicated in the U.S. PS cohort included the region harboring lipoma HMGIC fusion partner (LHFP) and conserved oligomeric golgi complex component 6 (COG6) genes on chromosome 13q13 (combined P = 2×10−6 for rs7993214; OR = 0.71), the late cornified envelope gene cluster (LCE) from the Epidermal Differentiation Complex (PSORS4) (combined P = 6.2×10−5 for rs6701216; OR 1.45) and a region of LD at 15q21 (combined P = 2.9×10−5 for rs3803369; OR = 1.43). This region is of interest because it harbors ubiquitin-specific protease-8 whose processed pseudogene lies upstream from HLA-C. This region of 15q21 also harbors the gene for SPPL2A (signal peptide peptidase like 2a) which activates tumor necrosis factor alpha by cleavage, triggering the expression of IL12 in human dendritic cells. We also identified a novel PSA (and potentially PS) locus on chromosome 4q27. This region harbors the interleukin 2 (IL2) and interleukin 21 (IL21) genes and was recently shown to be associated with four autoimmune diseases (Celiac disease, Type 1 diabetes, Grave''s disease and Rheumatoid Arthritis).  相似文献   

15.
Genome-wide association studies (GWAS) simultaneously investigating hundreds of thousands of single nucleotide polymorphisms (SNP) have become a powerful tool in the investigation of new disease susceptibility loci. Haplotypes are sometimes thought to be superior to SNPs and are promising in genetic association analyses. The application of genome-wide haplotype analysis, however, is hindered by the complexity of haplotypes themselves and sophistication in computation. We systematically analyzed the haplotype effects for breast cancer risk among 5,761 African American women (3,016 cases and 2,745 controls) using a sliding window approach on the genome-wide scale. Three regions on chromosomes 1, 4 and 18 exhibited moderate haplotype effects. Furthermore, among 21 breast cancer susceptibility loci previously established in European populations, 10p15 and 14q24 are likely to harbor novel haplotype effects. We also proposed a heuristic of determining the significance level and the effective number of independent tests by the permutation analysis on chromosome 22 data. It suggests that the effective number was approximately half of the total (7,794 out of 15,645), thus the half number could serve as a quick reference to evaluating genome-wide significance if a similar sliding window approach of haplotype analysis is adopted in similar populations using similar genotype density.  相似文献   

16.
Chicken is recognized as an excellent model for studies of genetic mechanism of phenotypic and genomic evolution, with large effective population size and strong human-driven selection. In the present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions employing 600K SNP Chicken chip in an F2 population of 1,534 hens, which was derived from reciprocal crosses between White Leghorn and Dongxiang chicken. Results indicated that a total of 49,151 core regions with an average length of 9.79 Kb were identified, which occupied approximately 52.15% of genome across all autosomes, and 806 significant core regions attracted us mostly. Genes in candidate regions may experience positive selection and were considered to have possible influence on beneficial economic traits. A panel of genes including AASDHPPT, GDPD5, PAR3, SOX6, GPC1 and a signal pathway of AKT1 were detected with the most extreme P-values. Further enrichment analyses indicated that these genes were associated with immune function, sensory organ development and neurogenesis, and may have experienced positive selection in chicken. Moreover, some of core regions exactly overlapped with genes excavated in our previous GWAS, suggesting that these genes have undergone positive selection may affect egg production. Findings in our study could draw a comparatively integrate genome-wide map of selection signature in the chicken genome, and would be worthy for explicating the genetic mechanisms of phenotypic diversity in poultry breeding.  相似文献   

17.
18.
Equine recurrent uveitis (ERU) is a common eye disease affecting up to 3–15% of the horse population. A genome-wide association study (GWAS) using the Illumina equine SNP50 bead chip was performed to identify loci conferring risk to ERU. The sample included a total of 144 German warmblood horses. A GWAS showed a significant single nucleotide polymorphism (SNP) on horse chromosome (ECA) 20 at 49.3 Mb, with IL-17A and IL-17F being the closest genes. This locus explained a fraction of 23% of the phenotypic variance for ERU. A GWAS taking into account the severity of ERU, revealed a SNP on ECA18 nearby to the crystalline gene cluster CRYGA-CRYGF. For both genomic regions on ECA18 and 20, significantly associated haplotypes containing the genome-wide significant SNPs could be demonstrated. In conclusion, our results are indicative for a genetic component regulating the possible critical role of IL-17A and IL-17F in the pathogenesis of ERU. The associated SNP on ECA18 may be indicative for cataract formation in the course of ERU.  相似文献   

19.
水稻株高基因及其分子生物学研究进展   总被引:6,自引:0,他引:6  
株高是水稻品种的重要农艺性状之一,植株过高容易引起倒伏而减产。近年来,随着研究手段的改进和水平的提高,在水稻株高基因的发现、定位、克隆及作用机理等方面的研究取得较大进展。本文介绍了国内外在该领域的研究进展概况。  相似文献   

20.
Development of testes or ovaries is critical to chicken breeders. Understanding the genetic mechanisms influencing the development of the testes and ovaries could enhance selection efforts which target reproductive traits. The linkage analysis was conducted within an F2 population derived from Beijing-You chickens and a commercial broiler line. The results have identified one quantitative trait loci (QTL, designated T1) for bilateral testicular weight (TW) and the percentage of TW to carcass weight, and five QTLs (designated O1–O5) for ovary weight (follicle-free, OW) and the percentage of OW to carcass weight. For the testes traits, QTL T1 is located between 6.55 and 8.56 Mb on GGA13. Especially, the gene gamma-amino butyric acid A receptor, alpha 1 (GABRA1) located near the T1 peak. For ovarian traits, QTL O2 was located at 29.31 Mb on GGA7. G protein-coupled receptor 39 (GPR39) present at the O2 peak was expressed at higher levels within the reproductive tract. It is also involved in the regulation of several reproductive functions. Other QTL peaks and the genes’ function in the ovary and testes need to be evaluated. The QTLs and the genes identified in this study could provide valuable information for establishing reproductive traits in chickens, and need further investigation.  相似文献   

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

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