首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The Icelandic horse is a breed known mainly for its ability to perform the ambling four‐beat gait ‘tölt’ and the lateral two‐beat gait pace. The natural ability of the breed to perform these alternative gaits is highly desired by breeders. Therefore, the discovery that a nonsense mutation (C>A) in the DMRT3 gene was the main genetic factor for horses' ability to perform gaits in addition to walk, trot and canter was of great interest. Although several studies have demonstrated that homozygosity for the DMRT3 mutation is important for the ability to pace, only about 70% of the homozygous mutant (AA) Icelandic horses are reported to pace. The aim of the study was to genetically compare four‐ and five‐gaited (i.e. horses with and without the ability to pace) AA Icelandic horses by performing a genome‐wide association (GWA) analysis. All horses (= 55) were genotyped on the 670K Axiom Equine Genotyping Array, and a GWA analysis was performed using the genabel package in r . No SNP demonstrated genome‐wide significance, implying that the ability to pace goes beyond the presence of a single gene variant. Despite its limitations, the current study provides additional information regarding the genetic complexity of pacing ability in horses. However, to fully understand the genetic differences between four‐ and five‐gaited AA horses, additional studies with larger sample materials and consistent phenotyping are needed.  相似文献   

2.
A previous study revealed a strong association between the DMRT3:Ser301STOP mutation in horses and alternate gaits as well as performance in harness racing. Several follow‐up studies have confirmed a high frequency of the mutation in gaited horse breeds and an effect on gait quality. The aim of this study was to determine when and where the mutation arose, to identify additional potential causal mutations and to determine the coalescence time for contemporary haplotypes carrying the stop mutation. We utilized sequences from 89 horses representing 26 breeds to identify 102 SNPs encompassing the DMRT3 gene that are in strong linkage disequilibrium with the stop mutation. These 102 SNPs were genotyped in an additional 382 horses representing 72 breeds, and we identified 14 unique haplotypes. The results provided conclusive evidence that DMRT3:Ser301STOP is causal, as no other sequence polymorphisms showed an equally strong association to locomotion traits. The low sequence diversity among mutant chromosomes demonstrated that they must have diverged from a common ancestral sequence within the last 10 000 years. Thus, the mutation occurred either just before domestication or more likely some time after domestication and then spread across the world as a result of selection on locomotion traits.  相似文献   

3.
Horse owners choose whether or not to train American Saddlebred horses (ASHs) to perform the 4‐beat gaits called rack and slow gait. The rack and slow gait are similar to ambling gaits shown to be associated with variation in the DMRT3 gene in other breeds but are trained rather than naturally occurring gaits. A premature stop codon in the DMRT3 gene (DMRT3_Ser301STOP) caused by the Ch23:g.22999655C>A SNP has an effect on the pattern of locomotion in horses and allows for the pacing gait and strong association with performance of ambling gaits in diverse breeds. We used horse show records to identify ASHs trained to perform as 5‐gaited horses and compared them to other Saddlebreds for the prevalence of the A allele of the Ch23:g.22999655C>A SNP of DMRT3. Genomic DNA was typed using a PCR‐RFLP technique for 5‐gaited ASHs (= 55), 3‐gaited or harness ASHs (= 64) and ASHs not shown (= 128). We did not observe differences in the frequencies of the Ch23:g.22999655C>A SNP of DMRT3 when comparing 5‐gaited ASHs to other ASHs (> 0.22). These results suggest that DMRT3 alleles do not affect the choice of breeders to train and show ASHs as 5‐gaited horses.  相似文献   

4.
The Mangalarga Marchador (MM) is a Brazilian horse breed known for a uniquely smooth gait. A recent publication described a mutation in the DMRT3 gene that the authors claim controls the ability to perform lateral patterned gaits (Andersson et al. 2012). We tested 81 MM samples for the DMRT3 mutation using extracted DNA from hair bulbs using a novel RFLP. Horses were phenotypically categorized by their gait type (batida or picada), as recorded by the Brazilian Mangalarga Marchador Breeders Association (ABCCMM). Statistical analysis using the plink toolset (Purcell, 2007) revealed significant association between gait type and the DMRT3 mutation (= 2.3e‐22). Deviation from Hardy–Weinberg equilibrium suggests that selective pressure for gait type is altering allele frequencies in this breed (= 1.00e‐5). These results indicate that this polymorphism may be useful for genotype‐assisted selection for gait type within this breed. As both batida and picada MM horses can perform lateral gaits, the DMRT3 mutation is not the only locus responsible for the lateral gait pattern.  相似文献   

5.
For centuries, domestic horses have represented an important means of transport and served as working and companion animals. Although their role in transportation is less important today, many horse breeds are still subject to intense selection based on their pattern of locomotion. A striking example of such a selected trait is the ability of a horse to perform additional gaits other than the common walk, trot and gallop. Those could be four‐beat ambling gaits, which are particularly smooth and comfortable for the rider, or pace, used mainly in racing. Gaited horse breeds occur around the globe, suggesting that gaitedness is an old trait, selected for in many breeds. A recent study discovered that a nonsense mutation in DMRT3 has a major impact on gaitedness in horses and is present at a high frequency in gaited breeds and in horses bred for harness racing. Here, we report a study of the worldwide distribution of this mutation. We genotyped 4396 horses representing 141 horse breeds for the DMRT3 stop mutation. More than half (2749) of these horses also were genotyped for a SNP situated 32 kb upstream of the DMRT3 nonsense mutation because these two SNPs are in very strong linkage disequilibrium. We show that the DMRT3 mutation is present in 68 of the 141 genotyped horse breeds at a frequency ranging from 1% to 100%. We also show that the mutation is not limited to a geographical area, but is found worldwide. The breeds with a high frequency of the stop mutation (>50%) are either classified as gaited or bred for harness racing.  相似文献   

6.
Previously, we performed a genome scan for type 2 diabetes (T2DM) using 638 African-American (AA) affected sibling pairs from 247 families; non-parametric linkage analysis suggested evidence of linkage at 6q24–27 (LOD 2.26). To comprehensively evaluate this region, we performed a two-stage association study by first constructing a SNP map of 754 SNPs selected from HapMap on the basis of linkage disequilibrium (LD) in 300 AAT2DM end-stage renal disease (ESRD) subjects, 311 AA controls, 43 European American controls and 45 Yoruba Nigerian samples (Set 1). Replication analyses were conducted in an independent population of 283 AA T2DM-ESRD subjects and 282 AA controls (Set 2). In addition, we adjusted for the impact of admixture on association results by using ancestry informative markers (AIMs). In Stage 1, 137 (18.2%) SNPs showed nominal evidence of association (P < 0.05) in one or more of tests of association: allelic (n = 33), dominant (n = 36), additive (n = 29), or recessive (n = 34) genotypic models, and 2- (n = 47) and 3-SNP (n = 43) haplotypic analyses. These SNPs were selected for follow-up genotyping. Stage 2 analyses confirmed association with a predicted 2-SNP “risk” haplotype in the PARK2 gene. Also, two intergenic SNPs showed consistent genotypic association with T2DM-ESRD: rs12197043 and rs4897081. Combined analysis of all subjects from both stages revealed nominal associations with 17 SNPs within genes, including suggestive associations in ESR1 and PARK2. This study confirms known diabetic nephropathy loci and identifies potentially novel susceptibility variants located within 6q24–27 in AA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

7.
With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (= 129) and TC <2 points (= 140); controls were those with PB >0.90 (= 58) and TC >4 points (= 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (< 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.  相似文献   

8.
We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3rd and 4th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3rd and 4th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.  相似文献   

9.
All species of the genus Populus (poplar, aspen) are dioecious, suggesting an ancient origin of this trait. Despite some empirical counter examples, theory suggests that nonrecombining sex‐linked regions should quickly spread, eventually becoming heteromorphic chromosomes. In contrast, we show using whole‐genome scans that the sex‐associated region in Populus trichocarpa is small and much younger than the age of the genus. This indicates that sex determination is highly labile in poplar, consistent with recent evidence of ‘turnover’ of sex‐determination regions in animals. We performed whole‐genome resequencing of 52 P. trichocarpa (black cottonwood) and 34 Populus balsamifera (balsam poplar) individuals of known sex. Genomewide association studies in these unstructured populations identified 650 SNPs significantly associated with sex. We estimate the size of the sex‐linked region to be ~100 kbp. All SNPs significantly associated with sex were in strong linkage disequilibrium despite the fact that they were mapped to six different chromosomes (plus 3 unmapped scaffolds) in version 2.2 of the reference genome. We show that this is likely due to genome misassembly. The segregation pattern of sex‐associated SNPs revealed this to be an XY sex‐determining system. Estimated divergence times of X and Y haplotype sequences (6–7 Ma) are much more recent than the divergence of P. trichocarpa (poplar) and Populus tremuloides (aspen). Consistent with this, in P. tremuloides, we found no XY haplotype divergence within the P. trichocarpa sex‐determining region. These two species therefore have a different genomic architecture of sex, suggestive of at least one turnover event in the recent past.  相似文献   

10.
Four single nucleotide polymorphisms (SNPs, rs2237892, rs2237895, rs2237897, rs2283228) in KCNQ1 are associated with type 2 diabetes mellitus in different ancestral groups. We investigated whether these 4 genetic markers are determinants of type 2 diabetes and premature coronary artery disease (CAD) in a Chinese population. We studied 398 consecutive patients, including 180 with coronary stenosis ≥50% or previous myocardial infarction (male <55 years, female <65 years) and 218 controls without documented CAD. CAD cases and controls were genotyped for 4 SNPs by using the ligase detection reaction method. The 3 genotypes AA, AC, and CC were present in rs2283228 and rs2237895, and the 3 genotypes CC, CT, TT were present in rs2237897 and rs2237892. No differences were found in genotype distribution and allele frequencies of these 4 SNPs between subjects with and without type 2 diabetes. Logistic regression showed that the risk of premature CAD in subjects carrying the CC genotype at rs2237892 was reduced by 90% in relation to individuals carrying the TT genotype (OR = 0.100, 95% CI: 0.018–0.564, P = 0.009). The association of other 3 SNPs with premature CAD could not be detected, nor did there exist any association of these 4 SNPs among groups of patients with 0, 1, 2, and 3-vessel disease (all P > 0.05). Our data implicate rs2237892 in KCNQ1 as a protective gene variant against premature CAD and we couldn’t replicate any association of these 4 SNPs with T2DM or extent of coronary lesions in a Chinese population.  相似文献   

11.
In this study, we performed a new genome‐wide association study using SLAF‐seq technology. A total of 19 single nucleotide polymorphism effects involving nine different SNP markers reached 5% Bonferroni‐corrected genome‐wide significance. In addition, a 5‐Mb region spanning 72.9–77.9 Mb on GGA4, exhibiting many significant SNP effects, was identified. The LDB2 gene in this region had a very strong association with body weight. Another SNP on GGA1, located in the INTS6 gene, had the strongest association with late body weight (weeks 10–16). Some of the SNPs that reached suggestive significance level overlapped with previously reported quantitative trait locus regions.  相似文献   

12.
2′-5′-oligoadenylate synthetases are a family of interferon-induced enzymes playing an important role in antiviral defense in mammals. In the human genome, three genes encoding functional synthetases (OAS1, OAS2 and OAS3) form a cluster. Previously, we found that particular genotypes and/or alleles of five single nucleotide polymorphisms (SNPs) of OAS2 and OAS3 are associated with predisposition to severe forms of tick-borne encephalitis (TBE) in Russians. In the current study, we investigated the distribution of three of the above SNPs, OAS3 rs2285932 (C/T, Ile438Ile), OAS3 rs2072136 (G/A, Ser567Ser), and OAS2 rs15895 (G/A, Trp720Ter relative to p71 isoform), in seven populations of North Eurasia: Caucasians (Russians, Germans from Altai region), Central Asian Mongoloids (Altaians, Khakass, Tuvinians, and Shorians), and Arctic Mongoloids (Chukchi). Interpopulational differences in genotype, allele and haplotype frequencies and in linkage disequilibrium structure for these SNPs were detected. These frequencies correlated with the ethnicity of the populations and with their supposed differential exposure to the TBE virus. In particular, the lowest frequencies of G/G genotype for OAS3 SNP rs2072136 (which, according to our earlier results, is associated with predisposition to severe forms of TBE) were found in Altaians, Khakass, Tuvinians, and Shorians, who commonly contact with the TBE virus in their habitation regions. Thus, the data obtained suggest that the TBE virus might act as a selection factor for particular OAS variants in Central Asian Mongoloids.  相似文献   

13.
L. Shi  L. Liu  Z. Ma  X. Lv  C. Li  L. Xu  B. Han  Y. Li  F. Zhao  Y. Yang  D. Sun 《Animal genetics》2019,50(5):430-438
Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (= 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (= <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (= <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.  相似文献   

14.
Neuroblastoma ranks as the most commonly seen and deadly solid tumour in infancy. The aberrant activity of m6A‐RNA methyltransferase METTL3 is involved in human cancers. Therefore, functional genetic variants in the METTL3 gene may contribute to neuroblastoma risk. In the current nine‐centre case‐control study, we aimed to analyse the association between the METTL3 gene single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility. We genotyped four METTL3 gene SNPs (rs1061026 T>G, rs1061027 C>A, rs1139130 A>G, and rs1263801 G>C) in 968 neuroblastoma patients and 1814 controls in China. We found significant associations between these SNPs and neuroblastoma risk in neither single‐locus nor combined analyses. Interestingly, in the stratified analysis, we observed a significant risk association with rs1061027 AA in subgroups of children ≤ 18 months of age (adjusted OR = 1.87, 95% CI = 1.03‐3.41, P = .040) and females (adjusted OR = 1.86, 95% CI = 1.07‐3.24, P = .028). Overall, we identified a significant association between METTL3 gene rs1061027 C>A polymorphism and neuroblastoma risk in children ≤18 months of age and females. Our findings provide novel insights into the genetic determinants of neuroblastoma.  相似文献   

15.
A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (< 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.  相似文献   

16.
Genetic studies in Turkish, Native American, European American, and African American (AA) families have linked chromosome 18q21.1–23 to susceptibility for diabetes-associated nephropathy. In this study, we have carried out fine linkage mapping in the 18q region previously linked to diabetic nephropathy in AAs by genotyping both microsatellite and single nucleotide polymorphisms (SNPs) for linkage analysis in an expanded set of 223 AA families multiplexed for type 2 diabetes associated ESRD (T2DM-ESRD). Several approaches were used to evaluate evidence of linkage with the strongest evidence for linkage in ordered subset analysis with an earlier age of T2DM diagnosis compared to the remaining pedigrees (LOD 3.9 at 90.1 cM, ∆P = 0.0161, NPL P value = 0.00002). Overall, the maximum LODs and LOD-1 intervals vary in magnitude and location depending upon analysis. The linkage mapping was followed up by performing a dense SNP map, genotyping 2,814 SNPs in the refined LOD-1 region in 1,029 AA T2DM-ESRD cases and 1,027 AA controls. Of the top 25 most associated SNPs, 10 resided within genic regions. Two candidate genes stood out: NEDD4L and SERPINB7. SNP rs512099, located in intron 1 of NEDD4L, was associated under a dominant model of inheritance [P value = 0.0006; Odds ratio (95% Confidence Interval) OR (95% CI) = 0.70 (0.57–0.86)]. SNP rs1720843, located in intron 2 of SERPINB7, was associated under a recessive model of inheritance [P value = 0.0017; OR (95% CI) = 0.65 (0.50–0.85)]. Collectively, these results suggest that multiple genes in this region may influence diabetic nephropathy susceptibility in AAs.  相似文献   

17.
Neuroblastoma ranks the most common seen solid tumour in childhood. Overexpression of LIN28A gene has been linked to the development of multiple human malignancies, but the relationship between LIN28A single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility is still under debate. Herein, we evaluated the correlation of four potentially functional LIN28A SNPs (rs3811464 G>A, rs3811463 T>C, rs34787247 G>A, and rs11247957 G>A) and neuroblastoma susceptibility in 505 neuroblastoma patients and 1070 controls from four independent hospitals in China. The correlation strengths were determined by using odds ratios (ORs) and corresponding 95% confidence intervals (CIs). Among these SNPs, rs34787247 G>A exhibited a significant association with increased susceptibility in neuroblastoma (GA vs GG: adjusted OR = 1.30, 95% CI = 1.03‐1.64; AA vs GG: adjusted OR = 2.51, 95% CI = 1.36‐4.64, AA/GA vs GG: adjusted OR = 1.42, 95% CI = 1.12‐1.80, AA vs GG/GA: adjusted OR = 2.39, 95% CI = 1.29‐4.42). Furthermore, the combined analysis of risk genotypes revealed that subjects carrying three risk genotypes (adjusted OR = 1.64, 95% CI = 1.02‐2.63) are more inclined to develop neuroblastoma than those without risk genotype, and so do carriers of 1‐4 risk genotypes (adjusted OR = 1.26, 95% CI = 1.01‐1.56). Stratification analysis further revealed risk effect of rs3811464 G>A, rs34787247 G>A and 1‐4 risk genotypes in some subgroups. Haplotype analysis of these four SNPs yields two haplotypes significantly correlated with increased neuroblastoma susceptibility. Overall, our finding indicated that LIN28A SNPs, especially rs34787247 G>A, may increase neuroblastoma risk.  相似文献   

18.
American Standardbred horses are divided into two groups based upon gait: the trot and the pace. The tendency to trot (diagonally opposite legs moving forward together) or pace (the two legs on the same side of the body moving forward together) appears to be genetically determined, although no formal genetic analysis has been undertaken. There is nearly complete assortative mating for gait; however, about 20% of the offspring sired by trotters are registered as pacers, while fewer than 1% of those sired by pacers are registered as trotters. Electrophoretically detectable genic variation at 13 protein loci has been analysed for 371 trotters and 856 pacers, and 10 blood group loci have been examined for 600 trotters and 1227 pacers. Trotters and pacers shared common alleles at 20 of the 23 loci; however, there were significant differences in allele frequencies at 21 of the 23 loci. Highly significant fixation indices (FSTS) were observed for 17 of the loci. The extent of genetic difference between Standardbred trotters and pacers was as great as or greater than that seen between some distinct horse breeds.  相似文献   

19.
Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.  相似文献   

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

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