首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A simple theoretical model of a Darwinian system (a periodic system with a multiplication phase and a selection phase) of entities (initial form of polymer strand, primary mutant and satellite mutants) is given. First case: one mutant is considered. One individual of the mutant appears in the multiplication phase of the first generation. The probabilities to find N mutants WnM(N) after the multiplication phase M of the n-th generation (with probability δ of an error in the replication, where all possible errors are fatal errors) and WnS(N) after the following selection phase S (with probability β that one individual survives) are given iteratively. The evolutionary tree is evaluated. Averages from the distributions and the probability of extinction WS(0) are obtained. Second case: two mutants are considered (primary mutant and new form). One individual of the primary mutant appears in the multiplication phase of the first generation. The probabilities to find Np primary mutants and Nm of the new form WnM(Np, Nm) after the multiplication phase M of the n-th generation (probability ε of an error in the replication of the primary mutant giving the new form) and WnS(Np, Nm) after the following selection phase S (probabilities βp and βm that one individual each of the primary mutant and of the new form survives) are given iteratively. Again the evolutionary tree is evaluated. Averages from the distributions are obtained.  相似文献   

2.
The popular neighbor-joining (NJ) algorithm used in phylogenetics is a greedy algorithm for finding the balanced minimum evolution (BME) tree associated to a dissimilarity map. From this point of view, NJ is "optimal" when the algorithm outputs the tree which minimizes the balanced minimum evolution criterion. We use the fact that the NJ tree topology and the BME tree topology are determined by polyhedral subdivisions of the spaces of dissimilarity maps to study the optimality of the neighbor-joining algorithm. In particular, we investigate and compare the polyhedral subdivisions for n ≤ 8. This requires the measurement of volumes of spherical polytopes in high dimension, which we obtain using a combination of Monte Carlo methods and polyhedral algorithms. Our results include a demonstration that highly unrelated trees can be co-optimal in BME reconstruction, and that NJ regions are not convex. We obtain the l 2 radius for neighbor-joining for n = 5 and we conjecture that the ability of the neighbor-joining algorithm to recover the BME tree depends on the diameter of the BME tree.  相似文献   

3.
The accurate reconstruction of phylogenies from short molecular sequences is an important problem in computational biology. Recent work has highlighted deep connections between sequence-length requirements for high-probability phylogeny reconstruction and the related problem of the estimation of ancestral sequences. In Daskalakis et al. (in Probab. Theory Relat. Fields 2010), building on the work of Mossel (Trans. Am. Math. Soc. 356(6):2379–2404, 2004), a tight sequence-length requirement was obtained for the simple CFN model of substitution, that is, the case of a two-state symmetric rate matrix Q. In particular the required sequence length for high-probability reconstruction was shown to undergo a sharp transition (from O(log n) to poly(n), where n is the number of leaves) at the “critical” branch length g ML(Q) (if it exists) of the ancestral reconstruction problem defined roughly as follows: below g ML(Q) the sequence at the root can be accurately estimated from sequences at the leaves on deep trees, whereas above g ML(Q) information decays exponentially quickly down the tree.  相似文献   

4.
LetS denote the intensity of a stimulus, ν the frequency of discharge of a fiber, and ν T the integral frequency of discharge of a nerve pathway. The relation between ν T andS can be derived when the relation between ν andS, as well as the distribution function of the fiber thresholds are given. It is investigated under what conditions the functions ν T (S) and ν(s) will be both of the same exponential form, and it is shown that this can happen only ifS exceeds the highest threshold of the pathway.  相似文献   

5.
The carbon balances of whole, 21-d old French bean plants (Phaseolus vulgaris L.) grown in standard nutrient solution (1K) and its modifications without (OK) or surplus (2K) potassium were calculated from the daily photosynthetic carbon inputs of individual leaves, and the daily respiratory carbon losses by individual leaves, stalks and petioles, and roots. Under the three K concentrations, maximum net photosynthetic rates (Pn) were found in the 2nd or in the 3rd trifoliate leaves, maximum respiratory rates (Rd) in the youngest, 4th trifoliate leaves; the Pn/Rd ratio decreased with leaf age. In all leaves of 2K plants, leaf dry masses and thicknesses, Pn, Pn/Pd ratios, and stomatal and intracellular conductances were lower than in OK and IK plants. Daily whole-plant net carbon gain was highest in IK plants, whereas in OK and 2K plants it was 98.0 and 81.3 % of IK, respectively. Similar values were found in the parameters of growth analysis, namely in net assimilation rates and relative growth rates. No differences were found in water potential (Ψ w ) or water saturation deficit (Wsat) in the OK, 1K and 2K plants sufficiently supplied with water or during wilting and resaturation. The decrease in Ψw to −0.97 MPa was associated with a 19.9 %, 31.4 % and 23.4 % decrease in Pn of OK, 1K and 2K plants, respectively, but no effect on Rd was found. In the three variants, the short-time effect of mild water stress was fully reversible.  相似文献   

6.
The analysis of kinetic and thermodynamic parameters of binding of peptide and nonpeptide dimerization inhibitors of HIV protease (HIVp) to the enzyme monomers immobilized on an optical chip has been studied by surface plasmon resonance. The molecular interactions were investigated at different inhibitor concentrations (0–80 μM) and temperatures (15–35°C). Determination of kinetic (k on, k off), equilibrium (K d), and thermodynamic (ΔG, ΔH, and -TΔS) has shown that both inhibitors are characterized by similar interaction parameters and the entropic term (-TΔS) of about −20 kcal/mol is the main driving force for the HIVp complex formation with the inhibitors, while the positive value (14 kcal/mol) of the enthalpic term (ΔH) counteracted the complex formation.  相似文献   

7.
Twenty-eight specimens of the large notothenioid Dissostichus mawsoni were dissected after capture on a set line near the southern limit of its range in McMurdo Sound, Antarctica. Total length (L T) averaged 127.3 cm (range 90–162 cm) and weight (W ) was 26.7 kg (range 10.4–60.3 kg). The length-weight relationship was W=3.44×10−5 L T 2.85 (n=28, r 2=0.96). Subcutaneous lipid thickness averaged 2.6 mm and showed no difference due to sex, but a significant weak relationship to W and L T. Hepatosomatic index (I H) was 1.6% for females and 1.7% for males; gonadosomatic index (I G) was 0.9% for females and 0.2% for males. Although specimens were large and sexually mature, the histology indicated that the gonads of this November sample were in a resting stage. Testes lacked spermatids and spermatozoa. Oocytes were in the previtellogenic stage (91.2%) or in the first stage of vitellogenesis (8.8%). A few atretic oocytes and empty follicles indicated that some females in this sample had spawned previously. A summary of the life-cycle is also presented. Accepted: 27 September 1999  相似文献   

8.
Dependency of specific leaf area (SLA) on shoot diameter (x) was studied forCryptomeria japonica foliage shoot segments about 5 cm in length taken from 18 branches at different height levels on 3 trees in August 1987. The leaf area of a shoot segment (S) was defined as half the sum of the needle area (s) estimated by the allometric relationship betweens and needle length. The dry weight of woody tissue (W w) and needles (W n) of a segment and itsS were divided by the segment length (L) to give linear densities asW w/L,W n/L andS/L, respectively. The densities were related tox by power-form equations.S/L values tended to be constant around 1.2 [cm2 cm−1] within the discussed range ofx, whileW w/L andW n/L values clearly increased withx. The approximately reciprocal relationship between average needle area (s) and linear density of the number of needles supported the fact thatS/L values were roughly constant regardless ofx.SLA andSNA were defined asS/(W w+W n) andS/W n, respectively. TheSLA-x relationship expected from the average value ofS/L and theW w/L-x and theW n/L-x relationships was well fitted to the observed decrease inSLA with increasingx. SNA also decreased asx increased. Variations inSLA andSNA among the shoot segments with similarx were not systematically related to their height levels. An empirical equation with a maximum value ofx (Xmax) was also proposed in order to formulate theSLA-x relationship.  相似文献   

9.
Two approaches to determine the fraction (μ) of mitochondrial respiration sustained during illumination by measuring CO2 gas exchange are compared. In single leaves, the respiration rate in the light (`day respiration' rate Rd) is determined as the ordinate of the intersection point of A–ci curves at various photon flux densities and compared with the CO2 evolution rate in darkness (`night respiration' rate Rn). Alternatively, using leaves with varying values of CO2 compensation concentration (Γ), intracellular resistance (ri) and Rn, an average number for μ can be derived from the linear regression between Γ and the product riċRn. Both methods also result in a number c* for that intercellular CO2 concentration at which net CO2 uptake rate is equal to –Rd. c* is an approximate value of the photocompensation point Γ* (Γ in the absence of mitochondrial respiration), which is related to the CO2/O2 specificity factor of Rubisco Sc/o. The presuppositions and limitations for application of both approaches are discussed. In leaves of Nicotiana tabacum, at 22 °C, single leaf measurements resulted in mean values of μ = 0.71 and c* = 34 μmol mol−1. At the photosynthetically active photon flux density of 960 μmol quanta m−2 s−1, nearly the same numbers were derived from the linear relationship between Γ and riċRn. c* and Rd determined by single leaf measurements varied between 31 and 41 μmol mol−1 and between 0.37 and 1.22 μmol m−2 s−1, respectively. A highly significant negative correlation between c* and Rd was found. From the regression equation we obtained estimates for Γ* (39 μmol mol−1), Sc/o (96.5 mol mol−1) and the mesophyll CO2 transfer resistance (7.0 mol−1 m2 s). This revised version was published online in June 2006 with corrections to the Cover Date.  相似文献   

10.
R. K. Misra  R. Sands 《Plant and Soil》1992,140(2):269-278
Diurnal variation in sap flux (S) through stems of six trees, two each of Ulmus procera SALISB., Melaleuca styphelioides SM. and Prunus cerasifera EHRH. ‘Nigra’ (referred to hereafter by their generic names), were estimated from measurements of heat pulse velocities. Leaf water potential (ψ), stomatal conductance (g s ) and transpiration from leaves (T) of all replicate trees were measured at 1300–1500h, once during the summer. On two separate occasions measurements were made of S, ψ, (g s ) and T for one each of Ulmus and Melaleuca trees to study diurnal variations in these parameters. A 12×12 m2 area around each tree was kept covered to simulate the condition of trees growing on pavements adjacent to residential properties. Sap flux for these tree species was in the order Melaleuca>Ulmus>Prunus. It is suggested that the smaller canopy and sapwood area in Prunus compared to the other two species is responsible for lower water potential and lower transpiration rate than the other species. Detailed analysis of the diurnal variation in sap flux and water relation of leaves of Melaleuca and Ulmus indicated sap flux of Melaleuca to be greater than that of Ulmus at the same transpiration rate per unit leaf area although the sapwood area of the two species was marginally different. This may have been due either to the difference in canopy conductance or in leaf area between the two species. With the assumption that sap flux closely resembles the rate of soil water extraction for both species, results indicate that Melaleuca is likely to extract soil water at a higher rate than Ulmus and hence is capable of causing greater shrinkage and soil movement than Ulmus.  相似文献   

11.
3 S4]+, S=1/2, composed of three, antiferromagnetically coupled high-spin ferric ions) by continuous wave (CW) and pulsed EPR techniques: Azotobacter vinelandii ferredoxin I, Desulfovibrio gigas ferredoxin II, and the 3Fe forms of Pyrococcus furiosus ferredoxin and aconitase. The 35 GHz (Q-band) CW EPR signals are simulated to yield experimental g tensors, which either had not been reported, or had been reported only at X-band microwave frequency. Pulsed X- and Q-band EPR techniques are used to determine electron spin-lattice (T 1, longitudinal) relaxation times at several positions on the samples' EPR envelope over the temperature range 2–4.2 K. The T 1 values vary sharply across the EPR envelope, a reflection of the fact that the envelope results from a distribution in cluster properties, as seen earlier as a distribution in g 3 values and in 57 Fe hyperfine interactions, as detected by electron nuclear double resonance spectroscopy. The temperature dependence of 1/T 1 is analyzed in terms of the Orbach mechanism, with relaxation dominated by resonant two-phonon transitions to a doublet excited state at ∼20 cm−1 above the doublet ground state for all four of these 3Fe proteins. The experimental EPR data are combined with previously reported 57Fe hyperfine data to determine electronic spin exchange-coupling within the clusters, following the model of Kent et al. Their model defines the coupling parameters as follows: J 13=J, J 12=J(1+ε′), J 23=J(1+ε), where J ij is the isotropic exchange coupling between ferric ions i and j, and ε and ε′ are measures of coupling inequivalence. We have extended their theory to include the effects of ε′≠0 and thus derived an exact expression for the energy of the doublet excited state for any ε, ε′. This excited state energy corresponds roughly to ε J and is in the range 5–10 cm−1 for each of these four 3Fe proteins. This magnitude of the product ε J, determined by our time-domain relaxation studies in the temperature range 2–4 K, is the same as that obtained from three other distinct types of study: CW EPR studies of spin relaxation in the range 5.5–50 K, NMR studies in the range 293–303 K, and static susceptibility measurements in the range 1.8–200 K. We suggest that an apparent disagreement as to the individual values of J and ε be resolved in favor of the values obtained by susceptibility and NMR (J≳200 cm−1 and ε≳0.02 cm−1 ), as opposed to a smaller J and larger ε as suggested in CW EPR studies. However, we note that this resolution casts doubt on the accepted theoretical model for describing the distribution in magnetic properties of 3Fe clusters. Received: 23 December 1999 / Accepted: 8 March 2000  相似文献   

12.
Synopsis The routine swimming speed (S) of three groups of 4, 9 and 32 cm total length (LT) juvenile cod (Gadus morhua) was quantified in the laboratory at 6 – 10 different temperatures (T) between 3.2 and 16.7°C. At temperatures between 5 and 15°C, mean group S increased exponentially with increasing T (S=a ebT) and the effect of temperature (b = 0.082, Q10 = 2.27) was not significantly different among the groups (over the 8-fold difference in fish sizes of early- and post-settlement juveniles). Differences in mean S among individuals within each group were quite large (coefficient of variation = 40 – 80%). Swimming data for juveniles and those collected for groups of 0.4, 0.7 and 0.9 cm standard length (LS) larvae were combined to assess the effect of body size on S. At 8°C, S (mm s−1) increased with LS (mm) according to: S = 0.26LSΦ−5.28LS−1, where Φ = 1.55LS−0.08. Relative S (body lengths s−1) was related to LS by a dome-shaped relationship having a maximum value (0.49 body lengths s−1) at 18.5 – 19 mm LS corresponding to the sizes of fish at the end of larval-juvenile metamorphosis. Previous larval cod IBM’s using a cruise-predator mode likely overestimated rates of foraging (prey searching and encounters) by a factor of ~2, whereas foraging rates in pause-travel models are closer to estimates of swimming velocities obtained in this and other laboratory studies.  相似文献   

13.
Six calcareous and alluvial soil profiles differing in their texture, CaCO3 and salinity were chosen from west and middle Nile Delta for the present study. The 1st and 2nd profiles from Borg El-Arab area were sandy loam in texture and > 30% CaCO3, while the 3rd and 4th profiles (from Nubaria area) were sandy clay loam and < 30% CaCO3. The 2nd and 4th profiles were taken from cultivated area with maize. The 5th profile from Epshan area was non-saline clay alluvial soil and the 6th from El-Khamsen was saline clay alluvial soil. The relation between soil moisture content (W%) and water vapour pressure (P/P o) was established for the mentioned soils. Data showed that the specific surface area (S) values were 34–53 and 44–60 m2/g for calcareous soils of Borg El-Arab and Nubaria areas, 206–219 and 206–249 m2/g for non-saline and saline clay alluvial soils of Epshan and El-Khamsen areas, respectively. The corresponding values of the external specific surface area (S e) were 16–21, 14–22, 72–86 and 92–112 m2/g. Submitting W m+W me as an adsorption boundary of moisture films (W c) (where W m is mono-adsorbed layer of water vapour on soil particles and W me is the external mono-adsorbed layer), the maximum water adsorption capacity (W a) was found to be W c + W me or W m + 2W me. It was ranged from 1.88 to 2.70%, 1.97 to 2.95%, 9.70–10.70% and 10.80 to 13.12% while the maximum hygroscopic water (M H) values were 2.43–3.78%, 2.91–4.65%, 16–17% and 18.30–21.9% for the studied soil profiles respectively. The residual moisture content (θ r) at pF 7 and P/P o = 0 was ranged from 0.0005–0.0010%, 0.0007–0.0019% and 0.0043–0.0048% in Borg El-Arab, Nubaria and Epshan soil profiles, respectively. The inter-relations between the surface area and the hygroscopic moisture parameters of the soils under investigation were as follows Calcareous soils; W m = 0.40 M H, W c = 0.55 M H, W a = 0.70 M H, S = 14 M H Non-saline soil; W m = 0.35 M H, W c = 0.49 M H, W a = 0.63 M H, S = 13 M H Saline soil; W m = 031 M H, W c = 0.45 M H, W a = 0.59 M H, S = 12 M H These relations give possibility to deduce the soil moisture adsorption capacities and specific surface area via maximum hygroscopic water, which can be obtained through the experimental determination of water vapor adsorption isotherms.  相似文献   

14.
Fish morphometric relationships are key tools for fisheries science and studies of food web dynamics and predator foraging behaviour, but parameterisations are limited for Southern Ocean myctophids (Family Myctophidae). New standard length (LS) to total mass (MT) relationships are therefore described for the 12 biomass-dominant myctophid fish species living in the Scotia Sea, Southern Ocean, using the most comprehensive data collected in the region to date. New linear regressions for otolith size (length; OL and width; OW) and LS are also described. Significant (p < .01) LSMT relationships were established for all species using simple non-linear regression. Significant (p < .01) relationships between LS and both OL and OW were also determined for all species, with OW being the best predictor of LS in all but one species. Our study provides a comprehensive tool for reconstructing the myctophid component of marine predator diets that will improve future food web, predator behaviour and ecosystem studies in the Scotia Sea.  相似文献   

15.
The relationship between pulmonary artery pressure (P LA) and oxygen saturation of mixed venous blood (S V ) has been studied in subjects (1750 men and 1026 women) subdivided into 12 groups. Functional relationships have been found between P LA and S V , P LA = f(S V ), and S V = f(P LA), which were estimated using direct measurement of P LA and S V for each group. These factors have been found to obey the following dependences: P LA = f(S V ) and P LA = a(S V )b , where b = −0.2284a + 0.6564 in men and b = −0.285a + 1.2947 in women; S V = f(P LA) and S V = c(P LA)d , where d = −0.25131Ln(c) + 1.0212; R 2 = 0.8993 in men and d = −1.9645Ln(c) + 2.852; and R 2 = 0.9674 in women. Each group occupies a position on the curves specified by the equations. Subjects with the diagnosis of a functional heart murmur and patients with congenital aortic valve stenosis were grouped together to form the so-called normal group characterized by specific P LA = f(S V ) and S V = f(P LA) dependences. Male patients with coronary heart disease were also included in the normal group. An equation was derived to relate P LA caused by different reasons with the corresponding saturation of mixed venous blood. In the case of the changing saturation of mixed venous blood, this equation gives the corresponding value of P LA. Equilibrium between systemic and pulmonary circulations is established through interdependent changes in the physiological indices of blood circulation and gas exchange in humans.  相似文献   

16.
Chemotactic collapse for the Keller-Segel model   总被引:1,自引:0,他引:1  
 This work is concerned with the system (S) {u t u − χ∇ (uv) for x∈Ω, t>0Γ v t v+(u−1) for x∈Ω, t>0 where Γ, χ are positive constants and Ω is a bounded and smooth open set in ℝ2. On the boundary ∂Ω, we impose no-flux conditions: (N) ∂u∂n =∂v∂n =0 for x∈∂ Ω, t>0 Problem (S), (N) is a classical model to describe chemotaxis corresponding to a species of concentration u(x, t) which tends to aggregate towards high concentrations of a chemical that the species releases. When completed with suitable initial values at t=0 for u(x, t), v(x, t), the problem under consideration is known to be well posed, locally in time. By means of matched asymptotic expansions techniques, we show here that there exist radial solutions exhibiting chemotactic collapse. By this we mean that u(r, t) →Aδ(y) as t→T for some T<∞, where A is the total concentration of the species. Received 9 March 1995; received in revised form 25 December 1995  相似文献   

17.
The length to weight ratio (LWR) was used to estimate populations of four fish species sampled on the Amazon coast of the state of Maranhão, Brazil. The fishes were sampled in 2018 and 2019 by artisanal coastal fishing, using gillnets (0.20–0.60 mm). The parameters a and b of the equation WT = a LTb were estimated. The obtained LWRs were WT = 0.06LT2.93, (r2) = .989 for Oligoplites saliens; WT = 0.061LT2.57, (r2) = .982 for Peprilus crenulatus; WT = 0.014LT2.80, (r2) = .985 for Bagre bagre; and WT = 0.035LT3.21, (r2) = .981 for Nebris microps.  相似文献   

18.
 A population with birth rate function B(N) N and linear death rate for the adult stage is assumed to have a maturation delay T>0. Thus the growth equation N′(t)=B(N(tT)) N(tT) e d 1 TdN(t) governs the adult population, with the death rate in previous life stages d 1≧0. Standard assumptions are made on B(N) so that a unique equilibrium N e exists. When B(N) N is not monotone, the delay T can qualitatively change the dynamics. For some fixed values of the parameters with d 1>0, as T increases the equilibrium N e can switch from being stable to unstable (with numerically observed periodic solutions) and then back to stable. When disease that does not cause death is introduced into the population, a threshold parameter R 0 is identified. When R 0<1, the disease dies out; when R 0>1, the disease remains endemic, either tending to an equilibrium value or oscillating about this value. Numerical simulations indicate that oscillations can also be induced by disease related death in a model with maturation delay. Received: 2 November 1998 / Revised version: 26 February 1999  相似文献   

19.
Internal conductances to CO2 transfer from the stomatal cavity to sites of carboxylation (gi) in hypostomatous sun-and shade-grown leaves of citrus, peach and Macadamia trees (Lloyd et al. 1992) were related to anatomical characteristics of mesophyll tissues. There was a consistent relationship between absorptance of photosynthetically active radiation and chlorophyll concentration (mmol m?2) for all leaves, including sclerophyllous Macadamia, whose transmittance was high despite its relatively thick leaves. In thin peach leaves, which had high gi, the chloro-plast volume and mesophyll surface area exposed to intercellular air spaces (ias) per unit leaf area were similar to those in the thicker leaves of the evergreen species. Peach leaves, however, had the lowest leaf dry weight per area (D/a), the lowest tissue density (Td) and the highest chloro-plast surface area (Sc) exposed to ias. There were negative correlations between gi and leaf thickness or D/a, but positive correlations between gi and Sc or Sc/Td. We developed a one-dimensional diffusion model which partitioned gi into a gaseous diffusion conductance through the ias (gias) plus a liquid-phase conductance through mesophyll cell walls (gcw). The model accounted for a significant amount of variation (r2=0.80) in measured gi by incorporating both components. The gias component was related to the one-dimensional path-length for diffusion across the mesophyll and so was greater in thinner peach leaves than in leaves of evergreen species. The gcw component was related to tissue density and to the degree of chloroplast exposure to the ias. Thus the negative correlations between gi and leaf thickness or D/a related to gias whereas positive correlations between gi and Sc or Sc/Td, related to gcw. The gcw was consistently lower than gias, and thus represented a greater constraint on CO2 diffusion in the mesophylls of these hypostomatous species.  相似文献   

20.
Cohen SD 《Microbial ecology》2006,52(3):463-469
Discula umbrinella, a fungal endophyte of oak species, colonizes and reproduces on leaves of Quercus alba and Q. rubra in forest ecosystems. Twenty-nine isolates collected from leaves of both oak species (16 from Q. alba and 13 from Q. rubra) were assayed for oak species preference and genetic variation based on primer-specific polymerase chain reactions for the intergenic spacer region (IGS) of ribosomal DNA. DNA sequencing of the polymerase chain reaction products revealed a 10-bp insertion (237–247 bp) at the 3′ end of the IGS region present in nine isolates and absent in 20 of the isolates. Phylogenetic analysis of the IGS region using the neighbor-joining method identified IGS groups (groups I–V) based on single nucleotide sequence differences. Host selectivity and geographic origin of isolates were correlated in some instances with the IGS groups. Isolates within each IGS group were further analyzed for nucleotide polymorphisms to confirm genotype identity and genotype diversity. Ten different genotypes (Va–Vj) were identified among the isolates analyzed. Genotype diversity was greatest in IGS groups I, IV, and V. Seventy percent of the genotypes (Vc, Vd, Ve, Vf, Vg,Vi, and Vj) contained isolates with single tree species preferences.  相似文献   

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

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