首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Life history patterns are usually identified by comparisons of extant species. Because of inferences regarding phylogenetic constraints, comparative data are often not statistically independent. In order to remove phylogenetic patterns embedded in life history data completely, we adopted a phylogenetic autoregressive method to reanalyse a data set of the ovipositional and developmental rates of 45 Phytoseiidae species. We first calculated the phylogenetic correlation in relation to different taxonomic levels using Moran's I statistics. Significant and positive phylogenetic correlations were found at the subgenus and subfamily levels. This indicates that some variation in both of these life-history traits could be accounted for by phylogeny. Phylogenetic associations, therefore, were removed by a phylogenetic autoregressive method. Using corrected data from this method, the specific components of the ovipositional rate are positively correlated with the specific components of th e developmental rate. The method that we have used obtains the same conclusion as others but differs from the phylogenetic effect in the way that it influences the relationship between comparative data. Because of no data reduction in the phylogenetic autoregressive method, the specific components are more useful than the mean values derived from the higher taxonomic nodes for testing ecological and evolutionary hypotheses about life history patterns. © Rapid Science Ltd. 1998  相似文献   

2.
Sahas Barve  Nicholas A. Mason 《Ibis》2015,157(2):299-311
The ecology of cavity nesting in passerine birds has been studied extensively, yet there are no phylogenetic comparative studies that quantify differences in life history traits between cavity‐ and open‐nesting birds within a passerine family. We test existing hypotheses regarding the evolutionary significance of cavity nesting in the Old World flycatchers (Muscicapidae). We used a multi‐locus phylogeny of 252 species to reconstruct the evolutionary history of cavity nesting and to quantify correlations between nest types and life history traits. Within a phylogenetic generalized linear model framework, we found that cavity‐nesting species are larger than open‐nesting species and that maximum clutch sizes are larger in cavity‐nesting lineages. In addition to differences in life history traits between nest types, species that breed at higher latitudes have larger average and maximum clutch sizes and begin to breed later in the year. Gains and losses of migratory behaviour have occurred far more often in cavity‐nesting lineages than in open‐nesting taxa, suggesting that cavity nesting may have played a crucial role in the evolution of migratory behaviour. These findings identify important macro‐evolutionary links between the evolution of cavity nesting, clutch size, interspecific competition and migratory behaviour in a large clade of Old World songbirds.  相似文献   

3.
Evaluating trait correlations across species within a lineage via phylogenetic regression is fundamental to comparative evolutionary biology, but when traits of interest are derived from two sets of lineages that coevolve with one another, methods for evaluating such patterns in a dual‐phylogenetic context remain underdeveloped. Here, we extend multivariate permutation‐based phylogenetic regression to evaluate trait correlations in two sets of interacting species while accounting for their respective phylogenies. This extension is appropriate for both univariate and multivariate response data, and may use one or more independent variables, including environmental covariates. Imperfect correspondence between species in the interacting lineages can also be accommodated, such as when species in one lineage associate with multiple species in the other, or when there are unmatched taxa in one or both lineages. For both univariate and multivariate data, the method displays appropriate type I error, and statistical power increases with the strength of the trait covariation and the number of species in the phylogeny. These properties are retained even when there is not a 1:1 correspondence between lineages. Finally, we demonstrate the approach by evaluating the evolutionary correlation between traits in fig species and traits in their agaonid wasp pollinators. R computer code is provided.  相似文献   

4.
Morphological divergence among species may be constrained by the pattern of genetic variances and covariances among traits within species. Assessing the existence of such a relationship in nature requires information on the stability of intraspecific correlation and covariance structure and the correspondence of this structure to the pattern of evolutionary divergence within a lineage. Here, we investigate these issues for nine morphological traits and 15 species of stalk-eyed flies in the genus Diasemopsis. Within-species matrices for these traits were generated from phenotypic data for all the Diasemopsis species and from genetic data for a single Diasemopsis species, D. dubia. The among-species pattern of divergence was assessed by calculating the evolutionary correlations for all pairwise combinations of the morphological traits along the phylogeny of these species. Comparisons of intraspecific matrices reveal significant similarity among all species in the phenotypic correlations matrices but not the covariance matrices. In addition, the differences in correlation structure that do exist among species are not related to their phylogenetic placement or change in the means of the traits. Comparisons of the phenotypic and phylogenetic matrices suggest a strong relationship between the pattern of evolutionary change among species and both the intraspecific correlation structure and the stability of this structure among species. The phenotypic and the phylogenetic matrices are significantly similar, and pairs of traits whose intraspecific correlations are more stable across taxa exhibit stronger coevolution on the phylogeny. These results suggest either the existence of strong constraints on the pattern of evolutionary change or a consistent pattern of correlated selection shaping both the phenotypic and phylogenetic matrices. The genetic correlation structure for D. dubia, however, does not correspond with patterns found in the phenotypic and phylogenetic data. Possible reasons for this disagreement are discussed.  相似文献   

5.
The question is often raised whether it is statistically necessary to control for phylogenetic associations in comparative studies. To investigate this question, we explore the use of a measure of phylogenetic correlation, lambda, introduced by Pagel (1999), that normally varies between 0 (phylogenetic independence) and 1 (species' traits covary in direct proportion to their shared evolutionary history). Simulations show lambda to be a statistically powerful index for measuring whether data exhibit phylogenetic dependence or not and whether it has low rates of Type I error. Moreover, lambda is robust to incomplete phylogenetic information, which demonstrates that even partial information on phylogeny will improve the accuracy of phylogenetic analyses. To assess whether traits generally show phylogenetic associations, we present a quantitative review of 26 published phylogenetic comparative data sets. The data sets include 103 traits and were chosen from the ecological literature in which debate about the need for phylogenetic correction has been most acute. Eighty-eight percent of data sets contained at least one character that displayed significant phylogenetic dependence, and 60% of characters overall (pooled across studies) showed significant evidence of phylogenetic association. In 16% of tests, phylogenetic correlation could be neither supported nor rejected. However, most of these equivocal results were found in small phylogenies and probably reflect a lack of power. We suggest that the parameter lambda be routinely estimated when analyzing comparative data, since it can also be used simultaneously to adjust the phylogenetic correction in a manner that is optimal for the data set, and we present an example of how this may be done.  相似文献   

6.
The quantitative characterization of the ecology of individual phytoplankton taxa is essential for model resolution of many aspects of aquatic ecosystems. Existing literature cannot directly parameterize all phytoplankton taxa of interest, as many traits and taxa have not been sampled. However, valuable clues on the value of traits are found in the evolutionary history of species and in common correlations between traits. These two resources were exploited with an existing, statistically consistent method built upon evolutionary concepts. From a new data set with >700 observations on freshwater phytoplankton traits and a qualitative phytoplankton phylogeny, estimates were derived for the size, growth rate, phosphate affinity, and susceptibility to predation of 277 phytoplankton types, from evolutionary ancestors to present‐day species. These estimates account simultaneously for phylogenetic relationships between types, as imposed by the phylogeny, and approximate power‐law relationships (e.g., allometric scaling laws) between traits, as reconstructed from the data set. Results suggest that most phytoplankton traits are to some extent conserved in evolution: cross‐validation demonstrated that the use of phylogenetic information significantly improves trait value estimates. By providing trait value estimates as well as uncertainties, these results could benefit most quantitative studies involving phytoplankton.  相似文献   

7.
The aim of this study was to evaluate the levels of phylogenetic heritability of the geographical range size, shape and position for 88 species of fiddler crabs of the world, using phylogenetic comparative methods and simulation procedures to evaluate their fit to the neutral model of Brownian motion. The geographical range maps were compiled from literature, and range size was based on the entire length of coastline occupied by each species, and the position of each range was calculated as its latitudinal and longitudinal midpoint. The range shape of each species was based in fractal dimension (box‐counting technique). The evolutionary patterns in the geographical range metrics were explored by phylogenetic correlograms using Moran’s I autocorrelation coefficients, autoregressive method (ARM) and phylogenetic eigenvector regression (PVR). The correlograms were compared with those obtained by simulations of Brownian motion processes across phylogenies. The distribution of geographical range size of fiddler crabs is right‐skewed and weak phylogenetic autocorrelation was observed. On the other hand, there was a strong phylogenetic pattern in the position of the range (mainly along longitudinal axis). Indeed, the ARM and PVR evidenced, respectively, that ca. 86% and 91% of the longitudinal midpoint could be explained by phylogenetic relationships among the species. The strong longitudinal phylogenetic pattern may be due to vicariant allopatric speciation and geographically structured cladogenesis in the group. The traits analysed (geographical range size and position) did not follow a Brownian motion process, thus suggesting that both adaptive ecological and evolutionary processes must be invoked to explain their dynamics, not following a simple neutral inheritance in the fiddler‐crab evolution.  相似文献   

8.
Species distributions are often constrained by climatic tolerances that are ultimately determined by evolutionary history and/or adaptive capacity, but these factors have rarely been partitioned. Here, we experimentally determined two key climatic niche traits (desiccation and cold resistance) for 92–95 Drosophila species and assessed their importance for geographic distributions, while controlling for acclimation, phylogeny, and spatial autocorrelation. Employing an array of phylogenetic analyses, we documented moderate‐to‐strong phylogenetic signal in both desiccation and cold resistance. Desiccation and cold resistance were clearly linked to species distributions because significant associations between traits and climatic variables persisted even after controlling for phylogeny. We used different methods to untangle whether phylogenetic signal reflected phylogenetically related species adapted to similar environments or alternatively phylogenetic inertia. For desiccation resistance, weak phylogenetic inertia was detected; ancestral trait reconstruction, however, revealed a deep divergence that could be traced back to the genus level. Despite drosophilids’ high evolutionary potential related to short generation times and high population sizes, cold resistance was found to have a moderate‐to‐high level of phylogenetic inertia, suggesting that evolutionary responses are likely to be slow. Together these findings suggest species distributions are governed by evolutionarily conservative climate responses, with limited scope for rapid adaptive responses to future climate change.  相似文献   

9.
Phylogenetic comparative methods have long considered phylogenetic signal as a source of statistical bias in the correlative analysis of biological traits. However, the main life-history strategies existing in a set of taxa are often combinations of life history traits that are inherently phylogenetically structured. In this paper, we present a method for identifying evolutionary strategies from large sets of biological traits, using phylogeny as a source of meaningful historical and ecological information. Our methodology extends a multivariate method developed for the analysis of spatial patterns, and relies on finding combinations of traits that are phylogenetically autocorrelated. Using extensive simulations, we show that our method efficiently uncovers phylogenetic structures with respect to various tree topologies, and remains powerful in cases where a large majority of traits are not phylogenetically structured. Our methodology is illustrated using empirical data, and implemented in the adephylo package for the free software R.  相似文献   

10.
The statistical estimation of phylogenies is always associated with uncertainty, and accommodating this uncertainty is an important component of modern phylogenetic comparative analysis. The birth–death polytomy resolver is a method of accounting for phylogenetic uncertainty that places missing (unsampled) taxa onto phylogenetic trees, using taxonomic information alone. Recent studies of birds and mammals have used this approach to generate pseudoposterior distributions of phylogenetic trees that are complete at the species level, even in the absence of genetic data for many species. Many researchers have used these distributions of phylogenies for downstream evolutionary analyses that involve inferences on phenotypic evolution, geography, and community assembly. I demonstrate that the use of phylogenies constructed in this fashion is inappropriate for many questions involving traits. Because species are placed on trees at random with respect to trait values, the birth–death polytomy resolver breaks down natural patterns of trait phylogenetic structure. Inferences based on these trees are predictably and often drastically biased in a direction that depends on the underlying (true) pattern of phylogenetic structure in traits. I illustrate the severity of the phenomenon for both continuous and discrete traits using examples from a global bird phylogeny.  相似文献   

11.
A method is proposed to conduct phylogenetic analyses of comparative or interspecific data when the true phylogeny is not known. Standard models of speciation and/or extinction or other methods are used to generate a sample from the set of all possible phylogenies for the measured species. The comparative data are then analyzed on each of the possible trees to obtain a distribution of possible evolutionary statistics for these data. The mean of this distribution is proposed as a reasonable estimate of the true evolutionary statistic of interest. Ways of obtaining confidence intervals and of developing hypothesis tests for this mean statistic are also proposed. The method can be used with any comparative method or phylogenetic analysis technique when phylogenetic relationships among species are not known or when branch lengths for a phylogeny in units of expected character change (as required by most methods) are not available. Computer programs to conduct the analyses are available on request.  相似文献   

12.
Aim Wood properties are related to tree physiology and mechanical stability and are influenced by both phylogeny and the environment. However, it remains unclear to what extent geographical gradients in wood traits are shaped by either phylogeny or the environment. Here we aimed to disentangle the influences of phylogeny and the environment on spatial trends in wood traits. Location China. Methods We compiled a data set of 11 wood properties for 618 tree species from 98 sampling sites in China to assess their phylogenetic and spatial patterns, and to determine how many of the spatial patterns in wood properties are attributable to the environment after correction for phylogenetic influences. Result All wood traits examined exhibited significant phylogenetic signal. The widest divergence in wood traits was observed between gymnosperms and angiosperms, Rosids and Asterids, Magnoiliids and Eudicots, and in Lamiales. For most wood traits, the majority of trait variation was observed at genus and species levels. The mechanical properties of wood showed correlated evolution with wood density. Most of the mechanical properties of wood exhibited significant latitudinal variation but limited or no altitudinal variation, and were positively correlated with mean annual precipitation based on both Pearson's correlation analysis and the phylogenetic comparative method. Correlations at family level between mean annual temperature and wood density, compression strength, cross‐section hardness, modulus of elasticity and volumetric shrinkage coefficient became significant after correction for phylogenetic influences. Main conclusions Phylogeny interacted with the environment in shaping the spatial patterns of wood traits of trees across China because most wood properties showed strong phylogenetic conservatism and thus affected environmental tolerances and distributions of tree species. Mean annual precipitation was a key environmental factor explaining the spatial patterns of wood traits. Our study provides valuable insights into the geographical patterns in productivity, distribution and ecological strategy of trees linking to wood traits.  相似文献   

13.
Recent developments in the analysis of comparative data   总被引:5,自引:0,他引:5  
Comparative methods can be used to test ideas about adaptation by identifying cases of either parallel or convergent evolutionary change across taxa. Phylogenetic relationships must be known or inferred if comparative methods are to separate the cross-taxonomic covariation among traits associated with evolutionary change from that attributable to common ancestry. Only the former can be used to test ideas linking convergent or parallel evolutionary change to some aspect of the environment. The comparative methods that are currently available differ in how they manage the effects brought about by phylogenetic relationships. One method is applicable only to discrete data, and uses cladistic techniques to identify evolutionary events that depart from phylogenetic trends. Techniques for continuous variables attempt to control for phylogenetic effects in a variety of ways. One method examines the taxonomic distribution of variance to identify the taxa within which character variation is small. The method assumes that taxa with small amounts of variation are those in which little evolutionary change has occurred, and thus variation is unlikely to be independent of ancestral trends. Analyses are then concentrated among taxa that show more variation, on the assumption that greater evolutionary change in the character has taken place. Several methods estimate directly the extent to which ancestry can predict the observed variation of a character, and subtract the ancestral effect to reveal variation of phylogeny. Yet another can remove phylogenetic effects if the true phylogeny is known. One class of comparative methods controls for phylogenetic effects by searching for comparative trends within rather than across taxa. With current knowledge of phylogenies, there is a trade-off in the choice of a comparative method: those that control phylogenetic effects with greater certainty are either less applicable to real data, or they make restrictive or untestable assumptions. Those that rely on statistical patterns to infer phylogenetic effects may not control phylogeny as efficiently but are more readily applied to existing data sets.  相似文献   

14.
There is an urgent need to reduce drastically the rate at which biodiversity is declining worldwide. Phylogenetic methods are increasingly being recognised as providing a useful framework for predicting future losses, and guiding efforts for pre-emptive conservation actions. In this study, we used a reconstructed phylogenetic tree of angiosperm species of the Eastern Arc Mountains – an important African biodiversity hotspot – and described the distribution of extinction risk across taxonomic ranks and phylogeny. We provide evidence for both taxonomic and phylogenetic selectivity in extinction risk. However, we found that selectivity varies with IUCN extinction risk category. Vulnerable species are more closely related than expected by chance, whereas endangered and critically endangered species are not significantly clustered on the phylogeny. We suggest that the general observation for taxonomic and phylogenetic selectivity (i.e. phylogenetic signal, the tendency of closely related species to share similar traits) in extinction risks is therefore largely driven by vulnerable species, and not necessarily the most highly threatened. We also used information on altitudinal distribution and climate to generate a predictive model of at-risk species richness, and found that greater threatened species richness is found at higher altitude, allowing for more informed conservation decision making. Our results indicate that evolutionary history can help predict plant susceptibility to extinction threats in the hyper-diverse but woefully-understudied Eastern Arc Mountains, and illustrate the contribution of phylogenetic approaches in conserving African floristic biodiversity where detailed ecological and evolutionary data are often lacking.  相似文献   

15.
We have presented a formal model for the quantitative analysis of phylogenetic and specific effects on the distribution of trait values among species. Total trait values are divided into phylogenetic values, inherited from an ancestral species, and specific values, the result of independent evolution. This allows a quantitative assessment of the strength of the phylogenetic inertia, or burden, displayed by a character in a lineage, so that questions concerning the relative importance of phylogenetic constraints in evolution can be answered. The separation of phylogenetic from specific effects proposed here also allows phylogenetic factors to be explicitly included in cross-species comparative analyses of adaptation. This solves a long-standing problem in evolutionary comparative studies. Only species' specific values can provide information concerning the independent evolution of characters in a set of related species. Therefore, only correlations among specific values for traits may be used as evidence for adaptation in cross-species comparative analyses. The phylogenetic autocorrelation model was applied to a comparative analysis of the determinants of sexual dimorphism in weight among 44 primate species. In addition to sexual dimorphism in weight, mating system, habitat, diet, and size (weight itself) were included in the analysis. All of the traits, except diet, were substantially influenced by phylogenetic inertia. The comparative analysis of the determinants of sexual dimorphism in weight indicates that 50% of the variation among primate species is due to phylogeny. Size, or scaling, could account for a total of 36% of the variance, making it almost as important as phylogeny in determining the level of dimorphism displayed by a species. Habitat, mating system, and diet follow, accounting for minor amounts of variation. Thus, in attempting to explain why a particular modern primate species is very dimorphic compared to other primates, we would say first because its ancestor was more dimorphic than average, second because it is a relatively large species, and third because it is terrestrial, polygynous, and folivorous.  相似文献   

16.
Avian acoustic communication has resulted from evolutionary pressures and ecological constraints. We therefore expect that auditory detectability in birds might be predictable by species traits and phylogenetic relatedness. We evaluated the relationship between phylogeny, species traits, and field‐based estimates of the two processes that determine species detectability (singing rate and detection distance) for 141 bird species breeding in boreal North America. We used phylogenetic mixed models and cross‐validation to compare the relative merits of using trait data only, phylogeny only, or the combination of both to predict detectability. We found a strong phylogenetic signal in both singing rates and detection distances; however the strength of phylogenetic effects was less than expected under Brownian motion evolution. The evolution of behavioural traits that determine singing rates was found to be more labile, leaving more room for species to evolve independently, whereas detection distance was mostly determined by anatomy (i.e. body size) and thus the laws of physics. Our findings can help in disentangling how complex ecological and evolutionary mechanisms have shaped different aspects of detectability in boreal birds. Such information can greatly inform single‐ and multi‐species models but more work is required to better understand how to best correct possible biases in phylogenetic diversity and other community metrics.  相似文献   

17.
We show that a new comparative method that sheds light on evolutionary trends among species may also illuminate trends within species. This finding comes from a phylogenetic autocorrelation analysis of morphological traits among individuals sampled from ten populations of a cooperatively breeding songbird, the Grey-crowned Babbler (Pomatostomus temporalis). Highly variable mitochondrial DNA (mtDNA) sequences from both the eastern (Pomatostomus temporalis temporalis) and western (Pomatostomus temporalis rubeculus) lineages were used to define genetic distances among 120 individuals and to estimate correlations among individuals in wing length, tarsus length, and body weight via an intraspecific weighting matrix. The autoregressive model effectively removed intraspecific correlations for all three morphological variables, and the proportion of the total phenotypic variance due to genealogical relationships varied from 0.68 (weight) to 0.23 (tarsus). The analysis revealed correlations among the specific components of traits, in which none were previously detected (type-I error) and diminished correlations that appeared strong when phylogeny was ignored. Group size was the only trait for which the autoregressive model failed to remove intraspecific correlations, a result likely due to the plasticity, convergence, and clinal variation in this trait in both the eastern and western lineages. The autocorrelation analysis weakened significant negative correlations between group size and total values for wing length and body weight, but the interpretation of this result depends on the adaptive significance ascribed to the “phylogenetic component” of trait values removed by the analysis. Comparative methods employing distance matrices are one useful way of summarizing the pattern of nonhierarchical relationships among conspecific individuals (“tokogenetic” relationships, sensu Hennig).  相似文献   

18.
The presence of phylogenetic signal is assumed to be ubiquitous. However, for microorganisms, this may not be true given that they display high physiological flexibility and have fast regeneration. This may result in fundamentally different patterns of resemblance, that is, in variable strength of phylogenetic signal. However, in microbiological inferences, trait similarities and therewith microbial interactions with its environment are mostly assumed to follow evolutionary relatedness. Here, we tested whether indeed a straightforward relationship between relatedness and physiological traits exists for aerobic methane‐oxidizing bacteria (MOB). We generated a comprehensive data set that included 30 MOB strains with quantitative physiological trait information. Phylogenetic trees were built from the 16S rRNA gene, a common phylogenetic marker, and the pmoA gene which encodes a subunit of the key enzyme involved in the first step of methane oxidation. We used a Blomberg's K from comparative biology to quantify the strength of phylogenetic signal of physiological traits. Phylogenetic signal was strongest for physiological traits associated with optimal growth pH and temperature indicating that adaptations to habitat are very strongly conserved in MOB. However, those physiological traits that are associated with kinetics of methane oxidation had only weak phylogenetic signals and were more pronounced with the pmoA than with the 16S rRNA gene phylogeny. In conclusion, our results give evidence that approaches based solely on taxonomical information will not yield further advancement on microbial eco‐evolutionary interactions with its environment. This is a novel insight on the connection between function and phylogeny within microbes and adds new understanding on the evolution of physiological traits across microbes, plants and animals.  相似文献   

19.

Background

Species number, functional traits, and phylogenetic history all contribute to characterizing the biological diversity in plant communities. The phylogenetic component of diversity has been particularly difficult to quantify in species-rich tropical tree assemblages. The compilation of previously published (and often incomplete) data on evolutionary relationships of species into a composite phylogeny of the taxa in a forest, through such programs as Phylomatic, has proven useful in building community phylogenies although often of limited resolution. Recently, DNA barcodes have been used to construct a robust community phylogeny for nearly 300 tree species in a forest dynamics plot in Panama using a supermatrix method. In that study sequence data from three barcode loci were used to generate a well-resolved species-level phylogeny.

Methodology/Principal Findings

Here we expand upon this earlier investigation and present results on the use of a phylogenetic constraint tree to generate a community phylogeny for a diverse, tropical forest dynamics plot in Puerto Rico. This enhanced method of phylogenetic reconstruction insures the congruence of the barcode phylogeny with broadly accepted hypotheses on the phylogeny of flowering plants (i.e., APG III) regardless of the number and taxonomic breadth of the taxa sampled. We also compare maximum parsimony versus maximum likelihood estimates of community phylogenetic relationships as well as evaluate the effectiveness of one- versus two- versus three-gene barcodes in resolving community evolutionary history.

Conclusions/Significance

As first demonstrated in the Panamanian forest dynamics plot, the results for the Puerto Rican plot illustrate that highly resolved phylogenies derived from DNA barcode sequence data combined with a constraint tree based on APG III are particularly useful in comparative analysis of phylogenetic diversity and will enhance research on the interface between community ecology and evolution.  相似文献   

20.
We assessed whether (1) arbuscular mycorrhizal colonization of roots (RC) and/or plant responses to arbuscular mycorrhizae (MR) vary with plant phylogeny and (2) MR and RC can be more accurately predicted with a phylogenetic predictor relative to a null model and models with plant trait and taxonomic predictors. In a previous study, MR and RC of 95 grassland species were measured. We constructed a phylogeny for these species and found it explained variation in MR and RC. Next, we used multiple regressions to identify the models that most accurately predicted plant MR. Models including either phylogenetic or phenotypic and taxonomic information similarly improved our ability to predict MR relative to a null model. Our study illustrates the complex evolutionary associations among species and constraints of using phylogenetic information, relative to plant traits, to predict how a plant species will interact with AMF.  相似文献   

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

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