首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Turnera sidoides is a complex of outcrossing, perennial, rhizomatous herbs that is widely distributed in southern South America. Five subspecies are recognized taxonomically based on morphological features and geographical distribution. In certain regions, the areas of distribution of the subspecies overlap partially. In such contact zones, the extent of reproductive barriers among subspecies is still largely unknown, but morphologically intermediate individuals have been found in the field, indicating that hybridization may actually occur between subspecies. Crossability among subspecies of T. sidoides has been shown by experimental studies with cultivated plants, but the mechanisms involved in natural populations are still unknown. To investigate the mechanisms that underlie gene flow within the T. sidoides complex, in this paper we analyze the morphological and genetic variation, as well as the crossability among taxa in a contact zone between subspecies pinnatifida and sidoides, in southeastern Uruguay. Our results constitute the first evidences of ongoing natural hybridization between subspecies of T. sidoides and suggest that, although hybridization may not have been of significance in the early phase of the species differentiation, reticulate evolution is ongoing enhancing the current morphological and genetic variability of the complex.  相似文献   

2.
Irrespective of its causes, strong population genetic structure indicates a lack of gene flow. Understanding the processes that underlie such structure, and the spatial patterns it causes, is valuable for conservation efforts such as restoration. On the other hand, when a species is invasive outside its native range, such information can aid management in the non-native range. Here we explored the genetic characteristics of the Australian tree Acacia dealbata in its native range. Two subspecies of A. dealbata have previously been described based on morphology and environmental requirements, but recent phylogeographic data raised questions regarding the validity of this taxonomic subdivision. The species has been widely planted within and outside its native Australian range and is also a highly successful invasive species in many parts of the world. We employed microsatellite markers to investigate the population genetic diversity and structure among 42 A. dealbata populations from across the species’ native range. We also tested whether environmental variables purportedly relevant for the putative separation of subspecies are linked with population genetic differentiation. We found no relationship between population genetic structure of A. dealbata in Australia and these environmental features. Rather, we identified two geographically distinct genetic clusters that corresponded with populations in the northeastern part of mainland Australia, and the southern mainland and Tasmanian range of the species. Our results do not support the taxonomic subdivision of the species into two distinct subspecies based on environmental features. We therefore assume that the observed morphological differences between the putative subspecies are plastic phenotypic responses. This study provides population genetic information that will be useful for the conservation of the species within Australia as well as to better understand the invasion dynamics of A. dealbata.  相似文献   

3.
The northern Neotropical region is characterized by a heterogeneous geological and climatic history. Recent studies have shown contrasting patterns regarding the role of geographic elements as barriers that could have determined phylogeographic structure in various species. Recently, the phylogeography and biogeography of Quercus species have been studied intensively, and the patterns observed so far suggest contrasting evolutionary histories for Neotropical species in comparison with their Holarctic relatives. The goal of this study was to describe the phylogeographic structure of two Neotropical oak species (Quercus insignis and Quercus sapotifolia) in the context of the geological and palaeoclimatic history of the northern Neotropics. Populations through the distribution range of both species were collected and characterized using nine chloroplast DNA microsatellite loci. Both oak species showed high levels of genetic diversity and strong phylogeographic structure. The distribution of genetic variation in Q. insignis suggested an influence of two major barriers, the Isthmus of Tehuantepec and the Nicaraguan Depression, while Q. sapotifolia exhibited a genetic structure defined by the heterogeneity of the Chortis highlands. The haplotype networks of both species indicated complex histories, suggesting that colonization from the Sierra Madre de Chiapas to central Mexico and from the north of the Nicaraguan Depression to the Costa Rican mountains may have occurred during different stages, and apparently more than one time. In conclusion, the phylogeographic structure of Neotropical oak species seems to be defined by a combination of geological and climatic events.  相似文献   

4.
Peach belongs to the genus Prunus, which includes Prunus persica and its relative species, P. mira, P. davidiana, P. kansuensis, and P. ferganensis. Of these, P. ferganensis have been classified as a species, subspecies, or geographical population of P. persica. To explore the genetic difference between P. ferganensis and P. persica, high-throughput sequencing was used in different peach accessions belonging to different species. First, low-depth sequencing data of peach accessions belonging to four categories revealed that similarity between P. ferganensis and P. persica was similar to that between P. persica accessions from different geographical populations. Then, to further detect the genomic variation in P. ferganensis, the P. ferganensis accession “Xinjiang Pan Tao 1” and the P. persica accession “Xia Miao 1” were sequenced with high depth, and sequence reads were assembled. The results showed that the collinearity of “Xinjiang Pan Tao 1” with the reference genome “Lovell” was higher than that of “Xia Miao 1” and “Lovell” peach. Additionally, the number of genetic variants, including single nucleotide polymorphisms (SNPs), structural variations (SVs), and the specific genes annotated from unmapped sequence in “Xia Miao 1” was higher than that in “Xinjiang Pan Tao 1” peach. The data showed that there was a close distance between “Xinjiang Pan Tao 1” (P. ferganensis) and reference genome which belong to P. persica, comparing “Xia Miao 1” (P. persica) and reference ones. The results accompany with phylogenetic tree and structure analysis confirmed that P. ferganensis should be considered as a geographic population of P. persica rather than a subspecies or a distinct species. Furthermore, gene ontology analysis was performed using the gene comprising large-effect variation to understand the phenotypic difference between two accessions. The result revealed that the pathways of gene function affected by SVs but SNPs and insertion-deletions markedly differed between the two peach accessions.  相似文献   

5.
Ecological speciation has long been noted as a central topic in the field of evolutionary biology, and investigation into the relative importance of ecological and geographical factors is becoming increasingly emphasized. We surveyed genetic variation of 277 samples from 25 populations of nine Rhododendron species within Tsutsusi subgenus in Taiwan using simple sequence repeats of expressed sequence tags. Bayesian clustering revealed four genetic lineages: (1) the Rhododendron simsii, Rhododendron kanehirai, and Rhododendron nakaharae lineage (lineage 1); (2) the Rhododendron longiperulatum, Rhododendron breviperulatum, and Rhododendron noriakianum lineage (lineage 2); (3) the Rhododendron rubropilosum lineage (lineage 3); and (4) the Rhododendron oldhamii lineage (lineage 4). Asymmetric introgressions were found from lineage 3 into lineages 1 and 2 (introgressed lineages). Genetic admixture of non-R. oldhamii species was also revealed by a neighbor-joining tree. Variation partitioning showed that environment explained much larger portions of genetic variation than geography between non-introgressed lineages (i.e., between R. oldhamii and other lineages). However, the Mantel and partial Mantel tests and the multiple matrix regression with randomization found that isolation-by-distance played a more important role than isolation-by-environment (IBE) in contributing to genetic variation in most between lineage comparisons. Nevertheless, strong IBE was found when compared between non-introgressed lineages of R. oldhamii and R. rubropilosum, suggesting post-speciation ecological divergence. Several environmental variables, including annual mean temperature, aspect, isothermality, seasonal precipitation, slope, and soil pH, could be important ecological drivers involved in reproductive isolation between R. oldhamii and non-R. oldhamii species within the Tsutsusi subgenus.  相似文献   

6.
The Neotropical otter, Lontra longicaudis, is considered near threatened by the IUCN, and with declining populations, it is still one of the lesser known species of otter in the world. Like other otter species, L. longicaudis is associated with water courses in a large variety of habitats and biomes across its area of distribution. Here, we present a review of the available literature on L. longicaudis aiming to update information on its distribution, ecology and systematics. We performed online searches using combinations of the words L. longicaudis, river otter, Neotropical otter, nutria and lobito del río, limiting our search to 1972–2016. Recent information on evolutionary history and biogeographic patterns, distribution patterns, as well as on the ecology of the species across a large part of its distribution range, including trophic ecology and intraspecific communication, was retrieved. We found severe gaps in knowledge including aspects of life history, ground validation of new potential distribution areas and long-term population monitoring aiming at designing more efficient conservation measures and techniques. In this context, we conclude by presenting perspectives for future studies.  相似文献   

7.
Oceanic archipelagos are typically rich in endemic taxa, because they offer ideal conditions for diversification and speciation in isolation. One of the most remarkable evolutionary radiations on the Canary Islands comprises the 16 species included in Limonium subsection Nobiles, all of which are subject to diverse threats, and legally protected. Since many of them are single-island endemics limited to one or a few populations, there exists a risk that a loss of genetic variation might limit their long-term survival. In this study, we used eight newly developed microsatellite markers to characterize the levels of genetic variation and inbreeding in L. macrophyllum, a species endemic to the North-east of Tenerife that belongs to Limonium subsection Nobiles. We detected generally low levels of genetic variation over all populations (H T = 0.363), and substantial differentiation among populations (F ST = 0.188; R ST = 0.186) coupled with a negligible degree of inbreeding (F?=?0.042). Obligate outcrossing may have maintained L. macrophyllum relatively unaffected by inbreeding despite the species’ limited dispersal ability and the genetic bottlenecks likely caused by a prolonged history of grazing. Although several factors still constitute a risk for the conservation of L. macrophyllum, the lack of inbreeding and the recent positive demographic trends observed in the populations of this species are factors that favour its future persistence.  相似文献   

8.
Mass occurrences of scypozoan medusae have become increasingly common in recent decades in European seas, including species in the genera Aurelia and Rhizostoma. We inferred the phylogeographic patterns of metagenetic scyphozoa Aurelia spp. and Rhizostoma pulmo from mitochondrial COI and nuclear ITS regions. No genetic structure was detected in R. pulmo over the Mediterranean Sea. By contrast, the phylogeographic analyses confirmed the separation of Aurelia spp. to several proposed cryptic species. Our results do not support the null hypothesis that both genera have concordant phylogeographic patterns. The resolvable parsimony network of haplotypes was retrieved for Aurelia aurita, Aurelia sp. 5, and Aurelia sp. 8 without connectivity between them and no genetic structure were found within those groups. Even though evidence of hybridization was found between A. aurita and Aurelia sp. 5, that did not break down the phylogenetic separation among them. The lowest haplotype and nucleotide diversity were found in samples of Aurelia sp. 8 and R. pulmo from the northern Adriatic, which acts as a sink area due to strong genetic drift. These new findings will facilitate linking the phenotype of the organism and its ability to survive in a particular environment—which shapes phylogeographic patterns.  相似文献   

9.
Far Eastern daces, genus Tribolodon (Cyprinidae), are thought to have diversified and developed unique diadromous life histories under changing conditions in the Sea of Japan and the surrounding environment. To examine the relationships between life history traits, distribution, and genetic population structures, we conducted a comparative phylogeographic analysis using partial mtDNA sequence data from samples collected over almost the full ranges of all four Tribolodon species. Phylogenetic analyses revealed several intraspecific haplotype groups that differentiated in the early Pleistocene to the Pliocene with or without geographic overlaps. A time-calibrated phylogeny suggested that the relatively smaller geographic ranges of the strictly freshwater species, T. sachalinensis and T. nakamurai, were explained not by the recent origins of these species, but by their limited dispersal abilities and smaller historical population sizes. The wider-ranging diadromous species, T. brandtii and T. hakonensis, exhibited similar major phylogeographic structures in their distributions, but the chronological order and timing of formation of this structure largely differed between the two species. In addition to those differences, the overlapping patterns of the differentiated intraspecific lineages in these species suggest dynamic, but somewhat restricted dispersal during the Plio-Pleistocene. Tribolodon hakonensis, one of the most widespread species of East Asian freshwater fishes, included both common and unique phylogeographic patterns compared to other fish species; the unique patterns (i.e., its wide range across freshwater biogeographic boundaries like the sea and mountains) would reflect its ecological features as a remarkable generalist inhabiting lakes, upper and lower reaches of rivers, and even coastal areas.  相似文献   

10.
The genus Jatropha (Euphorbiaceae) contains species that are of significant economic and ornamental value. However, Jatropha breeding material is rather limited due to incomplete information regarding phylogenetic relationships among germplasm resources. Phylogenetic analyses were performed based on the internal transcribed spacer of nuclear ribosomal DNA (nrDNA ITS), two chloroplast regions (trnL-F and rbcL), and the combined (ITS+trnL-F+rbcL) dataset among twenty-five specimens representing six key Jatropha species. Phylogenetic relationships of Jatropha were well resolved between subgenus Curcas and subgenus Jatropha, and demonstrated the intermediate position of section Polymorphae among sections of both subgenera. Jatropha curcas and J. integerrima demonstrated a close phylogenetic relationship. The molecular data agreed with the morphological classification that recognized J. multifida and J. podagrica in sec. Peltatae. The distinct intraspecific divergence that occurred in J. curcas could be attributed to restricted gene flow caused by geographical isolation and different ecological conditions. Phylograms produced with trnL-F and rbcL sequence data suggested slow rates of sequence divergence among Jatropha spp., while the ITS gene tree had good resolution suggesting high genetic variation of ITS among Jatropha species.  相似文献   

11.
Sandy beach ecosystems are decreasing worldwide and organisms living there are becoming threatened. The burrowing wolf spider Lycosa ishikariana is one such example. To establish effective conservation strategies under habitat fragmentation, we examined population genetic structure of L. ishikariana from mitochondrial cytochrome oxidase I gene and 6 microsatellite loci. Mitochondrial sequence data revealed 6 population subgroups with very high fixation indices, indicating that L. ishikariana has a clear phylogeographic structure and that the level of differentiation among regions is considerable. In particular, one subgroup in the western Honshu mainland (clade G) has a highly distinct genetic structure, despite having no clear geographic barriers from its parapatric population. Moreover, the distribution ranges of the other two subgroups (clades D and E) were highly restricted, suggesting their vulnerability to local human impacts and highlighting their high conservation priorities. Microsatellite data revealed 10 subgroups that were compatible with the clades identified from the mitochondrial data. Fixation indices among these groups were very high, indicating a limited gene flow induced by male spiders. Based on these results, we proposed six conservation units of L. ishikariana and effective conservation/restoration strategies in the face of ongoing coastal armoring.  相似文献   

12.
To characterize the molecular genetic diversity of the genus Darevskia, several populations were examined by the inter-SINE-PCR method, reporting the number and sizes of the spacers between individual copies of SINE-like interspersed repeats. Examination of 17 D. raddei geographical populations and several reference species revealed unequal genetic differences, measured as Nei and Li’s genetic distances (DNL), for different groups of samples. The highest homogeneity was observed for the apparently panmictic D. raddei nairensis population from the basin of the Hrazdan River: genetic differences within each of the five samples and between them were similarly low (less than 0.1). The difference between ten samples of D. raddei raddei from Armenia and Karabakh (0.2–0.3) was somewhat higher than the interindividual difference within each sample (0.1–0.2), indicating that the samples belonged to different populations. The assumption was supported by the phylogenetic tree topology and multidimensional scaling. The differences between samples from the morphological subspecies D. raddei raddei and D. raddei nairensis ranged 0.3–0.4. The difference of two D. raddei raddei samples of Talysh (Azerbaijan) from other samples of the same subspecies corresponded to the subspecific level. The genetic distances between the good species D. raddei and D. rudis was 0.6–0.7. In terms of DNL, a questionable population from northwestern Turkey (“D. tristis”) was closer to D. rudis (DNL = 0.45), probably representing its subspecies. The phylogeography of D. raddei is discussed.  相似文献   

13.
Competitive allele-specific PCR (KASPar) assay is a user-friendly system that provides flexibility in the numbers of single nucleotide polymorphisms (SNPs) and genotypes. Based on Illumina-GA-IIx genomic data from 10 genotypes with a broad genetic background, 3183 SNPs were selected for KASPar assays development, and 568 were finally converted and selected for Brassica rapa germplasm characterization (17.8%) on the basis of reproducibility, missing data rate, and uniform genetic distribution. High levels of polymorphism of these markers across 231 B. rapa genotypes were verified, illustrating by high polymorphic information content (averaged 0.34), minor allele frequency (0.37), genetic diversity (0.45), and the low observed heterozygosity (0.10). Based on the SNP dataset, structure and principal coordinates analysis, and neighbor-joining phylogenetic methods were used to examine the population structure and gave highly consistent results. The 231 accessions were divided into the four primary subspecies, representing 99 accessions from B. rapa ssp. pekinensis, 85 from B. rapa ssp. chinensis, 30 from B. rapa ssp. rapifera, and 17 from B. rapa ssp. oleifera and were further subdivided into 12 lower-order clusters according to different morphotypes. The genetic variability and pairwise fixation index analysis revealed that the ssp. pekinensis accessions possess the most extensive genetic variation among the four subspecies. The KASPar system is highly useful for validating SNPs and will be valuable for genetics research and breeding applications in B. rapa.  相似文献   

14.
Zanthoxylum is an economically and ecologically important genus of the Rutaceae family, of which Z. bungeanum and Z. armatum have a long history of cultivation in China. However, how the natural processes such as selection and drift and agriculture practices have influenced the genetic variation of cultivated Zanthoxylum species during long-term domestication remains elusive. Herein, we determined the population genetic structure of current widely cultivated Zanthoxylum species, Z. bungeanum and Z. armatum. Microsatellite markers revealed a high level of genetic variation and significant genetic differentiation for both species despite Z. bungeanum showed higher genetic diversity than Z. armatum. AMOVA indicated that most of the genetic variation exists within individuals rather than among provenances for both species. Population structure analyses generated three distinct groups within the entire accessions. All Z. bungeanum accessions were distinguished into two major geographic groups, north and south groups, with Qinling Mountains as the main geographic barrier to gene flow while a significant genetic differentiation was observed between cultivated and wild Z. armatum accessions. Mantel test of Z. bungeanum displayed a significant correlation between genetic and geographic distances within each inferred group but no correlation between genetic and geographic distance was observed when comparing genetic and geographic distances focusing only on pairwise of north vs. south provenances, ruling out the hypothesis that gene flow between north and south provenances followed an isolation-by-distance model. Our research provided a fundamental genetic profile that will improve the conservation and responsible exploitation of the extant germplasm of Zanthoxylum.  相似文献   

15.
Elaborate and expensive endeavours are underway worldwide to understand and manage biological invasions. However, the success of such efforts can be jeopardised due to taxonomic uncertainty. We highlight how unresolved native range taxonomy can complicate inferences in invasion ecology using the invasive tree Acacia dealbata in South Africa as an example. Acacia dealbata is thought to comprise two subspecies based on morphological characteristics and environmental requirements within its native range in Australia: ssp. dealbata and spp. subalpina. Biological control is the most promising option for managing invasive A. dealbata populations in South Africa, but it remains unknown which genetic/taxonomic entities are present in the country. Resolving this question is crucial for selecting appropriate biological control agents and for identifying areas with the highest invasion risk. We used species distribution models (SDMs) and phylogeographic approaches to address this issue. The ability of subspecies-specific and overall species SDMs to predict occurrences in South Africa was also explored. Furthermore, as non-overlapping bioclimatic niches between the two taxonomic entities may translate into evolutionary distinctiveness, we also tested genetic distances between the entities using DNA sequencing data and network analysis. Both approaches were unable to differentiate the two putative subspecies of A. dealbata. However, the SDM approach revealed a potential niche shift in the non-native range, and DNA sequencing results suggested repeated introductions of different native provenances into South Africa. Our findings provide important information for ongoing biological control attempts and highlight the importance of resolving taxonomic uncertainties in invasion ecology.  相似文献   

16.
Wood mice of the genus Apodemus belong to the most frequent and epidemiologically important rodents of Europe and adjacent regions. Previous studies showed that in the Middle East region species of this genus exhibit extraordinary morphological similarity precluding their proper determination without application of molecular characters. In order to determine the species of the studied populations and to obtain an insight into their phylogeographic history, we analyzed their genetic variation. We sequenced 1139 bp fragment of the mitochondrial DNA control region and flanking tRNA genes in samples collected from six localities. Phylogenetic analyses revealed presence of distinct clades corresponding to species A. uralensis and A. witherbyi. In most localities we confirmed presence of both species which suggests their large sympatric and syntopic occurrence. We recognized an extensive genetic variability, 38 specimens of A. uralensis belong to 32 distinct haplotypes, while 19 specimens of A. witherbyi to 14 haplotypes. We confirmed presence of several distinct haplotypes that may originate from multiple wood mouse colonization waves from distinct geographic regions.  相似文献   

17.
Transposable elements account for up to 85% of the maize genome and have significant implications in crop-improvement and evolutionary analyses. The Mutator (Mu) transposon superfamily, a class of DNA transposons, comprises the most complex and active elements in the maize genome, suggesting a special role in plant evolution. Here, we designed a set of Mu-specific primers based on terminal invert repeats and used a transposon display (TD) method for genotyping. We analyzed the distribution pattern of Mu insertions in teosinte (wild relative), sorghum (distant relative), and domesticated maize accessions (dent, sweet, and waxy). The MU-TD analysis suggested the presence of high polymorphic insertions among the species and subspecies, indicating the utility of the method in studying genetic variation and species relationships. Furthermore, we analyzed 80 maize recombinant inbred line populations. Mu-TD generated an average of 60% Mu-anchored polymorphic fragments in which insertions appeared to be segregating in significantly high numbers. The amplification profile was highly reproducible, confirming the utility of Mu elements as a new set of TD markers for developing high-density genetic maps.  相似文献   

18.

Key message

Using a much higher number of SNP markers and larger sample sizes than all the previous studies, we characterized the genetic relationships among wild and cultivated plants of section Beta.

Abstract

We analyzed the genetic variation of Beta section Beta, which includes wild taxa (Beta macrocarpa, B. patula, B. vulgaris subsp. adanensis and B. vulgaris subsp. maritima) and cultivars (fodder beet, sugar beet, garden beet, leaf beet, and swiss chards), using 9724 single nucleotide polymorphism markers. The analyses conducted at the individual level without a priori groups confirmed the strong differentiation of B. macrocarpa and B. vulgaris subsp. adanensis from the other taxa. B. vulgaris subsp. maritima showed a complex genetic structure partly following a geographical pattern, which confounded the differences between this taxon and the cultivated varieties. Cultivated varieties were structured into three main groups: garden beets, fodder and sugar beets, and leaf beets and swiss chards. The genetic structure described here will be helpful to correctly estimate linkage disequilibrium and to test for statistical associations between genetic markers and environmental variables.
  相似文献   

19.
Subspecies were traditionally defined by identifying gaps between phenotypes across the geographic range of a species, and may represent important units in the development of conservation strategies focused on preserving genetic diversity. Previous taxonomic research proposed that phenotypic variation between scattered Indri indri populations warranted the naming of two distinct subspecies, I. i. indri and I. i. variegatus. We tested these subspecific designations using mitochondrial sequence data generated from the control region or D-loop (569 bp) and a large section (2362 bp) of multiple genes and tRNAs known as Pastorini’s fragment and nuclear microsatellite markers. This study used 114 samples of I. indri from 12 rainforest sites in eastern Madagascar, encompassing the entire range of the species. These genetic samples represent multiple populations from low- and high-elevation forests from both putative subspecies. Molecular analyses of the mitochondrial sequence data did not support the two proposed subspecies. Furthermore, the microsatellite analyses showed no significant differences across the range beyond population level differentiation. This study demonstrates the utility of incorporating multiple lines of evidence in addition to phenotypic traits to define species or subspecies.  相似文献   

20.

Key message

Genetic structure among M. azedarach populations was detected and two subpopulations were present among them. A significant ‘isolation by distance’ was found in M. azedarach population in North-Western Plains of India.

Abstract

Melia azedarach is an important forest tree with pharmaceutical, insecticidal, pesticidal, and commercial significance. It is a good reforestation tree because of its fast growth and drought hardy nature. Genetic variation in a species allows itself to adapt, evolve and respond to environmental stress. It provides the basis for survival of a species and critically influences its evolutionary potential. Assessment of genetic diversity is necessary for improvement and conservation of a species. For this, microsatellite markers are of particular interest given the attributes like co-dominance, reproducibility, hyper variability and abundance throughout the genome. In the present study, we analyzed the genetic diversity and population structure of M. azedarach, an ecologically imperative species growing in the North-Western Plains of India. We developed 43 microsatellite markers, of which 20 were subsequently employed for analysis of diversity and population structure among 33 populations encompassing 318 genotypes representing North-Western Plains of India. A moderate level of diversity (Na = 5.1, Ho = 0.506, He = 0.712, I = 1.386) was assessed. The highest value of ΔK estimated using STRUCTURE indicated 2 subpopulations (K = 2). AMOVA exhibited 73 % variation within populations and 12 % variation was found among regions. Significant positive correlation between geographical and genetic distance was found (Rxy = 0.365, P = 0.010). The present study lays a foundation on a better understanding of genetic dynamics of the species and reveals its diversity and population structure in North-Western Plains of India.
  相似文献   

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

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