首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 562 毫秒
1.
The aim of this study was to study the population structure, to characterize the LD structure and to define core regions based on low recombination rates among SNP pairs in the genome of Piétrain pigs using data from the PorcineSNP60 BeadChip. This breed is a European sire line and was strongly selected for lean meat content during the last decades. The data were used to map signatures of selection using the REHH test. In the first step, selection signatures were searched genome‐wide using only core haplotypes having a frequency above 0.25. In the second step, the results from the selection signature analysis were matched with the results from the recently conducted genome‐wide association study for economical relevant traits to investigate putative overlaps of chromosomal regions. A small subdivision of the population with regard to the geographical origin of the individuals was observed. The extent of LD was determined genome‐wide using r2 values for SNP pairs with a distance ≤5 Mb and was on average 0.34. This comparable low r2 value indicates a high genetic diversity in the Piétrain population. Six REHH values having a p‐value < 0.001 were genome‐wide detected. These were located on SSC1, 2, 6 and 17. Three positional candidate genes with potential biological roles were suggested, called LOC100626459, LOC100626014 and MIR1. The results imply that for genome‐wide analysis especially in this population, a higher marker density and higher sample sizes are required. For a number of nine SNPs, which were successfully annotated to core regions, the REHH test was applied. However, no selection signatures were found for those regions (p‐value < 0.1).  相似文献   

2.
Coat color is one of the important factors characterizing breeds for domestic animals. Melanocortin 1 receptor (MC1R) is a representative responsible gene for this phenotype. Two single‐nucleotide polymorphisms (SNPs) in bovine MC1R gene, c.296T > C and c.310G>‐, have been well characterized, but these SNPs are not enough to explain cattle coat color. As far as we know, MC1R genotypes of Kumamoto sub‐breed of Japanese Brown cattle have not been analyzed. In the current study, genotyping for c.296T > C and c.310G>‐ was performed to elucidate the role of MC1R in determining the coat color of this sub‐breed. As a result, most animals were e/e genotype, suggesting the coat color of this sub‐breed is derived from the e allele of MC1R gene. However, we found six animals with E/e genotype, which coat color would be black theoretically. Subsequently, sequence comparison was performed with these animals to identify other polymorphisms affecting coat color, elucidating that these animals possessed the A allele of c.871G > A commonly. c.871G > A was a non‐synonymous mutation in the seventh transmembrane domain, suggesting alteration of the function and/or the structure of MC1R protein. Our data indicated that the A allele of c.871G > A might be a loss‐of‐function mutation.  相似文献   

3.
  相似文献   

4.
The genetic identification of the population of origin of individuals, including animals, has several practical applications in forensics, evolution, conservation genetics, breeding and authentication of animal products. Commercial high‐density single nucleotide polymorphism (SNP) genotyping tools that have been recently developed in many species provide information from a large number of polymorphic sites that can be used to identify population‐/breed‐informative markers. In this study, starting from Illumina BovineSNP50 v1 BeadChip array genotyping data available from 3711 cattle of four breeds (2091 Italian Holstein, 738 Italian Brown, 475 Italian Simmental and 407 Marchigiana), principal component analysis (PCA) and random forests (RFs) were combined to identify informative SNP panels useful for cattle breed identification. From a PCA preselected list of 580 SNPs, RFs were computed using ranking methods (Mean Decrease in the Gini Index and Mean Accuracy Decrease) to identify the most informative 48 and 96 SNPs for breed assignment. The out‐of‐bag (OOB) error rate for both ranking methods and SNP densities ranged from 0.0 to 0.1% in the reference population. Application of this approach in a test population (10% of individuals pre‐extracted from the whole data set) achieved 100% of correct assignment with both classifiers. Linkage disequilibrium between selected SNPs was relevant (r2 > 0.6) only in few pairs of markers indicating that most of the selected SNPs captured different fractions of variance. Several informative SNPs were in genes/QTL regions that affect or are associated with phenotypes or production traits that might differentiate the investigated breeds. The combination of PCA and RF to perform SNP selection and breed assignment can be easily implemented and is able to identify subsets of informative SNPs useful for population assignment starting from a large number of markers derived by high‐throughput genotyping platforms.  相似文献   

5.
Despite the wide range of observed phenotypic diversities and adaptation to different ecological conditions, little has been studied regarding the genetics of adaptation in the genome of indigenous cattle breeds of developing countries. Here, we investigated the linkage disequilibrium (LD) and identified the subset of outlier loci that are highly differentiated among cattle populations adapted to different ecological conditions in Ethiopia. Specifically, we genotyped 47 unrelated animals sampled from high‐ versus low‐altitude environments using a Bovine 50K SNP BeadChip. Linkage disequilibrium was assessed using both D′ and r2 between adjacent SNPs. We calculated FST and heterozygosity at different significance levels as measures of genetic differentiation for each locus between high‐ and low‐altitude populations following the hierarchical island model approach. We identified 816 loci (p < 0.01) showing selection signals and are associated with genes that might have roles in local adaptation. Some of them are associated with candidate genes that are involved in metabolism (ATP2A3, CA2, MYO18B, SIK3, INPP4A, and IREB2), hypoxia response (BDNF, TFRC, and PML) and heat stress (PRKDC, CDK1, and TFDC). Average r2 and D′ values were 0.14 ± 0.21 and 0.57 ± 0.34, respectively, for a minor allele frequency (MAF) ≥ 0.05 and were found to increase with increasing MAF value. The outlier loci identified in the studied Ethiopian cattle populations indicate the presence of genetic variation produced/shaped by adaptation to different environmental conditions and provide a basis for further validation and functional analysis using a reasonable sample size and high‐density markers.  相似文献   

6.
Bamaxiang pig is from Guangxi province in China, characterized by its small body size and two‐end black coat colour. It is an important indigenous breed for local pork market and excellent animal model for biomedical research. In this study, we performed genomewide association studies (GWAS) on 43 growth and carcass traits in 315 purebred Bamaxiang pigs based on a 1.4 million SNP array. We observed considerable phenotypic variability in the growth and carcass traits in the Bamaxiang pigs. The corresponding SNP based heritability varied greatly across the 43 traits and ranged from 9.0% to 88%. Through a conditional GWAS, we identified 53 significant associations for 35 traits at p value threshold of 10?6. Among which, 26 associations on chromosome 3, 7, 14 and X passed a genomewide significance threshold of 5 × 10?8. The most remarkable loci were at around 30.6 Mb on chromosome 7, which had growth stage‐dependent effects on body lengths and cannon circumferences and showed large effects on multiple carcass traits. We discussed HMGA1 NUDT3, EIF2AK1, TMEM132C and AFF2 that near the lead SNP of significant loci as plausible candidate genes for corresponding traits. We also showed that including phenotypic covariate in GWAS can help to reveal additional significant loci for the target traits. The results provide insight into the genetic architecture of growth and carcass traits in Bamaxiang pigs.  相似文献   

7.
Chinese indigenous pigs display marked genetic and phenotypic differences compared with western commercial pigs. In this study, we tested the association between 660K SNPs and longissimus muscle fatty acid composition traits in Chinese Erhualian (n = 331) and Bamaxiang (n = 315) pigs based on a customized 1.4 million SNP array. We identified a total of 64 significant associations for 20 fatty acid composition traits at the p‐value threshold of 1 × 10?6 among which 42 associations in low linkage disequilibrium (r2 < .2) with previously reported loci were considered novel. We substantially improved the strength and precision of the associations at four previously detected loci near FADS2, ELOVL7, ELOVL6 and FASN genes, facilitating follow‐up candidate gene studies. Moreover, we also identified loci near ABCD2, ACSBG1, ELOVL5, HPGDS, DAGT2, ACAD10 and ACSL1 genes with function relevant to metabolism of fatty acids. In this study, valuable genetic variants and candidate genes associated with fatty acid composition traits were identified in Erhualian and Bamaxiang pigs. Some identified loci could be used to improve pork nutrition in pig breeding practice. Using the SNP array with higher marker density and less ascertainment bias improved QTL detection power and precision in Chinese indigenous pigs.  相似文献   

8.
Linkage disequilibrium (LD) plays an important role in genomic selection and mapping of quantitative trait loci (QTL). This study investigated the pattern of LD and effective population size (Ne) in Gir cattle selected for yearling weight. For this purpose, 173 animals with imputed genotypes (from 18 animals genotyped with the Illumina BovineHD BeadChip and 155 animals genotyped with the Bovine LDv4 panel) were analysed. The LD was evaluated at distances of 25–50 kb, 50–100 kb, 100–500 kb and 0.5–1 Mb. The Ne was estimated based on 5 past generations. The r2 values (a measure of LD) were, respectively, .35, .29, .18 and .032 for the distances evaluated. The LD estimates decreased with increasing distance of SNP pairs and LD persisted up to a distance of 100 kb (r2 = .29). The Ne was greater in generations 4 and 5 (24 and 30 animals, respectively) and declined drastically after the last generation (12 animals). The results showed high levels of LD and low Ne, which were probably due to the loss of genetic variability as a consequence of the structure of the Gir population studied.  相似文献   

9.
Boar reproductive traits are economically important for the pig industry. Here we conducted a genome‐wide association study (GWAS) for 13 reproductive traits measured on 205 F2 boars at day 300 using 60 K single nucleotide polymorphism (SNP) data imputed from a reference panel of 1200 pigs in a White Duroc × Erhualian F2 intercross population. We identified 10 significant loci for seven traits on eight pig chromosomes (SSC). Two loci surpassed the genome‐wide significance level, including one for epididymal weight around 60.25 Mb on SSC7 and one for semen temperature around 43.69 Mb on SSC4. Four of the 10 significant loci that we identified were consistent with previously reported quantitative trait loci for boar reproduction traits. We highlighted several interesting candidate genes at these loci, including APN, TEP1, PARP2, SPINK1 and PDE1C. To evaluate the imputation accuracy, we further genotyped nine GWAS top SNPs using PCR restriction fragment length polymorphism or Sanger sequencing. We found an average of 91.44% of genotype concordance, 95.36% of allelic concordance and 0.85 of r2 correlation between imputed and real genotype data. This indicates that our GWAS mapping results based on imputed SNP data are reliable, providing insights into the genetic basis of boar reproductive traits.  相似文献   

10.
Meat quality of Iberian pigs is defined by the combination of their genetic characteristics and the particular production system. To carry out a genetic analysis of the main meat quality traits, we estimated their heritabilities, genetic correlations and the association effects of 32 selected SNPs of 12 candidate genes. A total of ten traits were measured in longissimus dorsi samples from 1,199 Iberian pigs fattened in the traditional free-range system: water holding capacity (thawing, cooking and centrifuge force water losses), instrumental colour (lightness L*, redness a* and yellowness b*), myoglobin content, shear force on cooked meat, and shear force and maximum compression force on dry-cured loin. Estimated heritability values were low to moderate (0.01 to 0.43) being the lowest for L* and the highest for cooking loss. Strong genetic correlations between water holding capacity traits (0.93 to 0.96) and between myoglobin content and a* (0.94) were observed. The association analyses revealed 19 SNPs significantly associated with different traits. Consistent and strong effects were observed between PRKAG3 SNPs (rs319678464G > C and rs330427832C > T), MYH3_rs81437544T > C, CASP3_rs319658214G > T and CTSL_rs332171512A > G and water losses. Also for CAPN1_rs81358667G > A and CASP3_rs319658214G > T and shear force. The SNPs mapping on PRKAG3 showed the highest effects on Minolta colour traits. Genotyping of these SNPs could be useful for the selection of Iberian young boars with similar estimated breeding values for productive traits.  相似文献   

11.
Circulating leptin concentrations could potentially be used as a predictor of production traits in cattle. This study aimed to clarify the correlations between circulating leptin concentrations and growth performance, carcass traits, and meat quality indexes in finishing bulls fed high‐concentrate diets (concentrate‐to‐forage ratio 70:30). Fifty‐seven Simmental × Luxi F1 crossbred bulls were used for 112‐day finishing experiment. Circulating leptin concentrations and relevant indexes of growth performance, and carcass traits and meat quality were measured during or after finishing trail. The results indicated that the leptin concentrations tended to be negatively correlated with dry matter intake (DMI) (r = ?.233, p = .081), and were positively correlated with 12th‐rib fat thickness (r = .330, p = .012), marbling score (r = .336, p = .011), and intramuscular fat content (r = .368, p = .021). Moreover, the leptin concentrations were negatively correlated with cholesterol content (r = ?.339, p = .037) and were not correlated with sensory indexes including tenderness, juiciness, and like flavor (p > .05). In conclusion, circulating leptin concentrations may potentially be used as a predictor of carcass traits related to content of fat and beef quality traits related to content of cholesterol in finishing bulls fed high‐concentrate diets.  相似文献   

12.
Conservation decisions based on neutral genetic diversity have been observed to promote retention of useful quantitative variation in biological populations. An experiment was undertaken to determine the association between microsatellite marker polymorphisms and phenotypic variation in semen production and cryosurvival traits in bulls. Thirty-five ejaculates were collected from ten bulls of two breeds and evaluated before and after cryopreservation for several semen traits. The bulls were also genotyped using a set of sixteen bovine-specific microsatellite marker loci. Fixation indices (FST), heterozygosity and Nei's genetic distance measures were computed from allele frequency data for each of the bulls. Molecular and phenotypic data were used to compute tri-distance matrices for the ten bulls and correlated using Mantel's test in GenAIEx 6.5. The study revealed extensive heterogeneity in semen traits, heterozygosity and FST values among the bulls. Large pairwise phenotypic and genetic distances were also observed. Correlation between pairwise genetic distances and phenotypic distances was significant and highly positive for sperm viability (r = .61, p < .001) and moderately positive for sperm motility (r = .40–42, p < .05) variables. For sperm morphology, ejaculate volume and sperm concentration, correlation with genetic distances was positive, low and not significantly different from zero (p > .05). A tendency for a triangular-shaped relationship between genetic and phenotypic distances for post-thaw motility and viability traits was also observed. Accordingly, association with neutral genetic diversity was absent for semen production traits and moderate to highly positive for sperm cryosurvival traits. Given these findings, conservation decisions based on neutral genetic diversity may capture variation in some adaptive traits, but not others.  相似文献   

13.
The aim of the present study was to detect quantitative trait loci affecting fatty acid composition in back fat and intramuscular fat in a Duroc pig population comprising seventh‐generation pedigrees using genome‐wide association studies (GWAS). In total, 305 animals were genotyped using single nucleotide polymorphisms (SNPs) array and five selected SNPs from regions containing known candidate genes related to fatty acid synthesis or metabolism. In total, 24 genome‐wide significant SNP regions were detected in 12 traits, and 76 genome‐wide suggestive SNP regions were detected in 33 traits. The Sus scrofa chromosome (SSC) 7 at 10.3 Mb was significantly associated with C17:0 in intramuscular fat, while the SSC9 at 13.6 Mb was significantly associated with C14:0 in intramuscular fat. The SSC12 at 1.0 Mb was significantly associated with C14:0 in back fat and the SSC14 at 121.0 Mb was significantly associated with C18:0 in intramuscular fat. These regions not only replicated previously reported loci containing some candidate genes involved in fatty acid composition (fatty acid synthase and stearoyl‐CoA desaturase) but also included several additional related loci.  相似文献   

14.
Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8–59.3 Mb on BTA 4; 67.5–69.3 on BTA 5; 39.5–41.0 Mb on BTA 6; and 26.4–29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.  相似文献   

15.
The aim of this study was to separate marked additive genetic variability for three quantitative traits in chickens into components associated with classes of minor allele frequency (MAF), individual chromosomes and marker density using the genomewide complex trait analysis (GCTA) approach. Data were from 1351 chickens measured for body weight (BW), ultrasound of breast muscle (BM) and hen house egg production (HHP), each bird with 354 364 SNP genotypes. Estimates of variance components show that SNPs on commercially available genotyping chips marked a large amount of genetic variability for all three traits. The estimated proportion of total variation tagged by all autosomal SNPs was 0.30 (SE 0.04) for BW, 0.33 (SE 0.04) for BM, and 0.19 (SE 0.05) for HHP. We found that a substantial proportion of this variation was explained by low frequency variants (MAF <0.20) for BW and BM, and variants with MAF 0.10–0.30 for HHP. The marked genetic variance explained by each chromosome was linearly related to its length (R2 = 0.60) for BW and BM. However, for HHP, there was no linear relationship between estimates of variance and length of the chromosome (R2 = 0.01). Our results suggest that the contribution of SNPs to marked additive genetic variability is dependent on the allele frequency spectrum. For the sample of birds analysed, it was found that increasing marker density beyond 100K SNPs did not capture additional additive genetic variance.  相似文献   

16.
The objective of this study was to evaluate, using three different genotype density panels, the accuracy of imputation from lower‐ to higher‐density genotypes in dairy and beef cattle. High‐density genotypes consisting of 777 962 single‐nucleotide polymorphisms (SNP) were available on 3122 animals comprised of 269, 196, 710, 234, 719, 730 and 264 Angus, Belgian Blue, Charolais, Hereford, Holstein‐Friesian, Limousin and Simmental bulls, respectively. Three different genotype densities were generated: low density (LD; 6501 autosomal SNPs), medium density (50K; 47 770 autosomal SNPs) and high density (HD; 735 151 autosomal SNPs). Imputation from lower‐ to higher‐density genotype platforms was undertaken within and across breeds exploiting population‐wide linkage disequilibrium. The mean allele concordance rate per breed from LD to HD when undertaken using a single breed or multiple breed reference population varied from 0.956 to 0.974 and from 0.947 to 0.967, respectively. The mean allele concordance rate per breed from 50K to HD when undertaken using a single breed or multiple breed reference population varied from 0.987 to 0.994 and from 0.987 to 0.993, respectively. The accuracy of imputation was generally greater when the reference population was solely comprised of the breed to be imputed compared to when the reference population comprised of multiple breeds, although the impact was less when imputing from 50K to HD compared to imputing from LD.  相似文献   

17.
Previously accurate genomic predictions for Bacterial cold water disease (BCWD) resistance in rainbow trout were obtained using a medium‐density single nucleotide polymorphism (SNP) array. Here, the impact of lower‐density SNP panels on the accuracy of genomic predictions was investigated in a commercial rainbow trout breeding population. Using progeny performance data, the accuracy of genomic breeding values (GEBV) using 35K, 10K, 3K, 1K, 500, 300 and 200 SNP panels as well as a panel with 70 quantitative trait loci (QTL)‐flanking SNP was compared. The GEBVs were estimated using the Bayesian method BayesB, single‐step GBLUP (ssGBLUP) and weighted ssGBLUP (wssGBLUP). The accuracy of GEBVs remained high despite the sharp reductions in SNP density, and even with 500 SNP accuracy was higher than the pedigree‐based prediction (0.50–0.56 versus 0.36). Furthermore, the prediction accuracy with the 70 QTL‐flanking SNP (0.65–0.72) was similar to the panel with 35K SNP (0.65–0.71). Genomewide linkage disequilibrium (LD) analysis revealed strong LD (r2 ≥ 0.25) spanning on average over 1 Mb across the rainbow trout genome. This long‐range LD likely contributed to the accurate genomic predictions with the low‐density SNP panels. Population structure analysis supported the hypothesis that long‐range LD in this population may be caused by admixture. Results suggest that lower‐cost, low‐density SNP panels can be used for implementing genomic selection for BCWD resistance in rainbow trout breeding programs.  相似文献   

18.
The objective was to assess goodness of fit and predictive ability of subsets of single nucleotide polymorphism (SNP) markers constructed based on minor allele frequency (MAF), effect sizes and varying marker density. Target traits were body weight (BW), ultrasound measurement of breast muscle (BM) and hen house egg production (HHP) in broiler chickens. We used a 600 K Affymetrix platform with 1352 birds genotyped. The prediction method was genomic best linear unbiased prediction (GBLUP) with 354 564 single nucleotide polymorphisms (SNPs) used to derive a genomic relationship matrix ( G ). Predictive ability was assessed as the correlation between predicted genomic values and corrected phenotypes from a threefold cross‐validation. Predictive ability was 0.27 ± 0.002 for BW, 0.33 ± 0.001 for BM and 0.20 ± 0.002 for HHP. For the three traits studied, predictive ability decreased when SNPs with a higher MAF were used to construct G . Selection of the 20% SNPs with the largest absolute effect sizes induced a predictive ability equal to that from fitting all markers together. When density of markers increased from 5 K to 20 K, predictive ability enhanced slightly. These results provide evidence that designing a low‐density chip using low‐frequency markers with large effect sizes may be useful for commercial usage.  相似文献   

19.
This study investigated genomic predictions across Nordic Holstein and Nordic Red using various genomic relationship matrices. Different sources of information, such as consistencies of linkage disequilibrium (LD) phase and marker effects, were used to construct the genomic relationship matrices (G‐matrices) across these two breeds. Single‐trait genomic best linear unbiased prediction (GBLUP) model and two‐trait GBLUP model were used for single‐breed and two‐breed genomic predictions. The data included 5215 Nordic Holstein bulls and 4361 Nordic Red bulls, which was composed of three populations: Danish Red, Swedish Red and Finnish Ayrshire. The bulls were genotyped with 50 000 SNP chip. Using the two‐breed predictions with a joint Nordic Holstein and Nordic Red reference population, accuracies increased slightly for all traits in Nordic Red, but only for some traits in Nordic Holstein. Among the three subpopulations of Nordic Red, accuracies increased more for Danish Red than for Swedish Red and Finnish Ayrshire. This is because closer genetic relationships exist between Danish Red and Nordic Holstein. Among Danish Red, individuals with higher genomic relationship coefficients with Nordic Holstein showed more increased accuracies in the two‐breed predictions. Weighting the two‐breed G‐matrices by LD phase consistencies, marker effects or both did not further improve accuracies of the two‐breed predictions.  相似文献   

20.
The effectiveness of five selection methods for genetic improvement of net merit comprising trait 1 of low heritability (h2 = 0.1) and trait 2 of high heritability (h2 = 0.4) was examined: (i) two‐trait quantitative trait loci (QTL)‐assisted selection; (ii) partial QTL‐assisted selection based on trait 1; (iii) partial QTL‐assisted selection based on trait 2; (iv) QTL‐only selection; and (v) conventional selection index without QTL information. These selection methods were compared under 72 scenarios with different combinations of the relative economic weights, the genetic correlations between traits, the ratio of QTL variance to total genetic variance of the trait, and the ratio of genetic variances between traits. The results suggest that the detection of QTL for multiple‐trait QTL‐assisted selection is more important when the index traits are negatively correlated than when they are positively correlated. In contrast to literature reports that single‐trait marker‐assisted selection (MAS) is the most efficient for low heritability traits, this study found that the identified QTL of the low heritability trait contributed negligibly to total response in net merit. This is because multiple‐trait QTL‐assisted selection is designed to maximize total net merit rather than the genetic response of the individual index trait as in the case of single‐trait MAS. Therefore, it is not economical to identify the QTL of the low heritability traits for the improvement of total net merit. The efficient, cost‐effective selection strategy is to identify the QTL of the moderate or high heritability traits of the QTL‐assisted selection index to facilitate total economic returns. Detection of the QTL of the low h2 traits for the QTL‐assisted index selection is justified when the low h2 traits have high negative genetic correlation with the other index traits and/or when both economic weights and genetic variances of the low h2 traits are larger as compared to the other index traits of higher h2. This study deals with theoretical efficiency of QTL‐assisted selection, but the same principle applies to SNP‐based genomic selection when the proportion of the genetic variance ‘explained by the identified QTLs’ in this study is replaced by ‘explained by SNPs’.  相似文献   

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

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