首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
BackgroundSpike protein is the surface glycoprotein of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) necessary for the entry of the virus via the transmembrane receptors of the human respiratory cells causing COVID-19 disease.AimHere, we aimed to predict the three-dimensional monomer structure of spike protein of SARS-CoV-2 from 20 Jordanian nasopharyngeal samples and to determine the percentage of single amino acid variants (SAV) in the spike protein of SARS-CoV-2.MethodsThe output of the Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) found four single amino acid variants in the spike gene.ResultsThe first variant represented by 5% of samples that showed tyrosine deletion at Y144 located in the N terminal domain. The second and the dominant variant, represented by 62%, showed aspartate a coil amino acid substitution to glycine an extracellular amino acid at D614G located in the spike recognition binding site. The third variant, represented by 5%, showed aspartate substitution to tyrosine at D1139Y, and the fourth variant, represented by 5% glycine substitution to serine at G1167S.ConclusionOur results have shown low mutational sensitivity in all variants except to D614G the one with the most likely neutral mutational sensitivity that all variants might not explicitly affect the function of spike glycoprotein. However, D614G might change the viral conformational plasticity and hence a potential viral fitness gain but one must be cautious about drawing any concrete conclusions about the severity of symptoms and viral transmission from genomic data only.General significanceStudying mutations such as D614G in deep is essential to control the pandemic in terms of immune systems, antibodies, or even vaccines.  相似文献   

2.
Mutations in the spike protein of SARS-CoV-2 are the major causes for the modulation of ongoing COVID-19 infection. Currently, the D614G substitution in the spike protein has become dominant worldwide. It is associated with higher infectivity than the ancestral (D614)variant. We demonstrate using Gaussian network model-based normal mode analysis that the D614G substitution occurs at the hinge region that facilitates domain-domain motions between receptor binding domain and S2 region of the spike protein. Computer-aided mutagenesis and inter-residue energy calculations reveal that contacts involving D614 are energetically frustrated. However, contacts involving G614 are energetically favourable, implying the substitution strengthens residue contacts that are formed within as well as between protomers. We also find that the free energy difference (ΔΔG) between two variants is -2.6 kcal/mol for closed and -2.0 kcal/mol for 1-RBD up conformation. Thus, the thermodynamic stability has increased upon D614G substitution. Whereas the reverse mutation in spike protein structures having G614 substitution has resulted in the free energy differences of 6.6 kcal/mol and 6.3 kcal/mol for closed and 1-RBD up conformations, respectively, indicating that the overall thermodynamic stability has decreased. These results suggest that the D614G substitution modulates the flexibility of spike protein and confers enhanced thermodynamic stability irrespective of conformational states. This data concurs with the known information demonstrating increased availability of the functional form of spikeprotein trimer upon D614G substitution.  相似文献   

3.
In-silico studies on SARS-CoV-2 genome are considered important to identify the significant pattern of variations and its possible effects on the structural and functional characteristics of the virus. The current study determined such genetic variations and their possible impact among SARS-CoV-2 variants isolated in India. A total of 546 SARS-CoV-2 genomic sequences (India) were retrieved from the gene bank (NCBI) and subjected to alignment against the Wuhan variant (NC_045512.2), the corresponding amino acid changes were analyzed using NCBI Protein-BLAST. These 546 variants revealed 841 mutations; most of these were non-synonymous 464/841 (55.1%), there was no identical variant compared to the original strain. All genes; coding and non-coding showed nucleotide changes, most of the structural genes showed frequent nonsynonymous mutations. The most affected genes were ORF1a/b followed by the S gene which showed 515/841 (61.2%) and 120/841 (14.3%) mutations, respectively. The most frequent non-synonymous mutation 486/546 (89.01%) occurred in the S gene (structural gene) at position 23,403 where A changed to G leading to the replacement of aspartic acid by glycine in position (D614G). Interestingly, four variants also showed deletion. The variants MT800923 and MT800925 showed 12 consecutive nucleotide deletion in position 21982–21993 resulting in 4 consecutive amino acid deletions that were leucine, glycine, valine, and tyrosine in positions 141, 142, 143, and 144 respectively. The present study exhibited a higher mutations rate per variant compared to other studies carried out in India.  相似文献   

4.
Severe COVID‐19 is characterized by lung abnormalities, including the presence of syncytial pneumocytes. Syncytia form when SARS‐CoV‐2 spike protein expressed on the surface of infected cells interacts with the ACE2 receptor on neighboring cells. The syncytia forming potential of spike variant proteins remain poorly characterized. Here, we first assessed Alpha (B.1.1.7) and Beta (B.1.351) spread and fusion in cell cultures, compared with the ancestral D614G strain. Alpha and Beta replicated similarly to D614G strain in Vero, Caco‐2, Calu‐3, and primary airway cells. However, Alpha and Beta formed larger and more numerous syncytia. Variant spike proteins displayed higher ACE2 affinity compared with D614G. Alpha, Beta, and D614G fusion was similarly inhibited by interferon‐induced transmembrane proteins (IFITMs). Individual mutations present in Alpha and Beta spikes modified fusogenicity, binding to ACE2 or recognition by monoclonal antibodies. We further show that Delta spike also triggers faster fusion relative to D614G. Thus, SARS‐CoV‐2 emerging variants display enhanced syncytia formation.  相似文献   

5.
The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct “S-R/x3” to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics.  相似文献   

6.
As novel SARS-CoV-2 variants continue to emerge, it is critical that their potential to cause severe disease and evade vaccine-induced immunity is rapidly assessed in humans and studied in animal models. In early January 2021, a novel SARS-CoV-2 variant designated B.1.429 comprising 2 lineages, B.1.427 and B.1.429, was originally detected in California (CA) and it was shown to have enhanced infectivity in vitro and decreased antibody neutralization by plasma from convalescent patients and vaccine recipients. Here we examine the virulence, transmissibility, and susceptibility to pre-existing immunity for B 1.427 and B 1.429 in the Syrian hamster model. We find that both variants exhibit enhanced virulence as measured by increased body weight loss compared to hamsters infected with ancestral B.1 (614G), with B.1.429 causing the most marked body weight loss among the 3 variants. Faster dissemination from airways to parenchyma and more severe lung pathology at both early and late stages were also observed with B.1.429 infections relative to B.1. (614G) and B.1.427 infections. In addition, subgenomic viral RNA (sgRNA) levels were highest in oral swabs of hamsters infected with B.1.429, however sgRNA levels in lungs were similar in all three variants. This demonstrates that B.1.429 replicates to higher levels than ancestral B.1 (614G) or B.1.427 in the oropharynx but not in the lungs. In multi-virus in-vivo competition experiments, we found that B.1. (614G), epsilon (B.1.427/B.1.429) and gamma (P.1) dramatically outcompete alpha (B.1.1.7), beta (B.1.351) and zeta (P.2) in the lungs. In the nasal cavity, B.1. (614G), gamma, and epsilon dominate, but the highly infectious alpha variant also maintains a moderate size niche. We did not observe significant differences in airborne transmission efficiency among the B.1.427, B.1.429 and ancestral B.1 (614G) and WA-1 variants in hamsters. These results demonstrate enhanced virulence and high relative oropharyngeal replication of the epsilon (B.1.427/B.1.429) variant in Syrian hamsters compared to an ancestral B.1 (614G) variant.  相似文献   

7.
The high transmissibility of SARS-CoV-2 is related to abundant replication in the upper airways, which is not observed for the other highly pathogenic coronaviruses SARS-CoV and MERS-CoV. We here reveal features of the coronavirus spike (S) protein, which optimize the virus towards the human respiratory tract. First, the S proteins exhibit an intrinsic temperature preference, corresponding with the temperature of the upper or lower airways. Pseudoviruses bearing the SARS-CoV-2 spike (SARS-2-S) were more infectious when produced at 33°C instead of 37°C, a property shared with the S protein of HCoV-229E, a common cold coronavirus. In contrast, the S proteins of SARS-CoV and MERS-CoV favored 37°C, in accordance with virus preference for the lower airways. Next, SARS-2-S-driven entry was efficiently activated by not only TMPRSS2, but also the TMPRSS13 protease, thus broadening the cell tropism of SARS-CoV-2. Both proteases proved relevant in the context of authentic virus replication. TMPRSS13 appeared an effective spike activator for the virulent coronaviruses but not the low pathogenic HCoV-229E virus. Activation of SARS-2-S by these surface proteases requires processing of the S1/S2 cleavage loop, in which both the furin recognition motif and extended loop length proved critical. Conversely, entry of loop deletion mutants is significantly increased in cathepsin-rich cells. Finally, we demonstrate that the D614G mutation increases SARS-CoV-2 stability, particularly at 37°C, and, enhances its use of the cathepsin L pathway. This indicates a link between S protein stability and usage of this alternative route for virus entry. Since these spike properties may promote virus spread, they potentially explain why the spike-G614 variant has replaced the early D614 variant to become globally predominant. Collectively, our findings reveal adaptive mechanisms whereby the coronavirus spike protein is adjusted to match the temperature and protease conditions of the airways, to enhance virus transmission and pathology.  相似文献   

8.
The non-synonymous mutations of SARS-CoV-2 isolated from across the world have been identified during the last few months. The surface glycoprotein spike of SARS-CoV-2 forms the most important hotspot for amino acid alterations followed by the ORF1a/ORF1ab poly-proteins. It is evident that the D614G mutation in spike glycoprotein and P4715L in RdRp is the important determinant of SARS-CoV-2 evolution since its emergence. P4715L in RdRp, G251V in ORF3a and S1498F of Nsp3 is associated with the epitope loss that may influence pathogenesis caused by antibody escape variants. The phylogenomics distinguished the ancestral viral samples from China and most part of Asia, isolated since the initial outbreak and the later evolved variants isolated from Europe and Americas. The evolved variants have been found to predominant globally with the loss of epitopes from its proteins. These have implications for SARS-CoV-2 transmission, pathogenesis and immune interventions.  相似文献   

9.
The pandemic of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is far from being controlled despite the great effort that have been taken throughout the world. Herd immunity through vaccination is our major expectation to rein the virus. However, the emergence of widespread genetic variants could potentially undermine the vaccines. The evidence that some variants could evade immune responses elicited by vaccines and previous infection is growing. In this review, we summarized the current understanding on five notable genetic variants, i.e., D614G, Cluster 5, VOC 202012/01, 501Y.V2 and P.1, and discussed the potential impact of these variants on the virus transmission, pathogenesis and vaccine efficacy. We also highlight that mutations in the N-terminal domain of spike protein should be considered when evaluating the antibody neutralization abilities. Among these genetic variants, a concern of genetic variant 501Y.V2 to escape the protection by vaccines was raised. We therefore call for new vaccines targeting this variant to be developed.  相似文献   

10.
The D614G mutation in the spike protein of SARS-CoV-2 alters the fitness of the virus, leading to the dominant form observed in the COVID-19 pandemic. However, the molecular basis of the mechanism by which this mutation enhances fitness is not clear. Here we demonstrated by cryo-electron microscopy that the D614G mutation resulted in increased propensity of multiple receptor-binding domains (RBDs) in an upward conformation poised for host receptor binding. Multiple substates within the one RBD-up or two RBD-up conformational space were determined. According to negative staining electron microscopy, differential scanning calorimetry, and differential scanning fluorimetry, the most significant impact of the mutation lies in its ability to eliminate the unusual cold-induced unfolding characteristics and to significantly increase the thermal stability under physiological pH. The D614G spike variant also exhibited exceptional long-term stability when stored at 37 °C for up to 2 months. Our findings shed light on how the D614G mutation enhances the infectivity of SARS-CoV-2 through a stabilizing mutation and suggest an approach for better design of spike protein-based conjugates for vaccine development.  相似文献   

11.
High-resolution melting (HRM) analysis was conducted to discriminate between SARS-CoV-2 Omicron variant BA.1 (B.1.1.529.1) and subvariant BA.2 (B.1.1.529.2). We performed two-step PCR consisting of the first PCR and the second nested PCR to prepare the amplicon for HRM analysis, which detected G339D, N440K, G446S and D796Y variations in the SARS-CoV-2 spike protein. The melting temperatures (Tms) of the amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 receptor binding domain (RBD) in spike protein were the same: 75.2 °C (G339D variation) and 73.4 °C (D796Y variation). These Tms were distinct from those of SARS-CoV-2 isolate Wuhan-Hu-1, and were specific to the Omicron variant. In HRM analyses that detected the N440K and G446S variations, the Tms of amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 RBDs were 73.0 °C (N440K and G446S variations) and 73.5 °C (G446S variation). This difference indicates that the SARS-CoV-2 Omicron variants BA.1 and BA.2 can be clearly discriminated. Our study demonstrates the usefulness of HRM analysis after two-step PCR for the discrimination of SARS-CoV-2 variants.  相似文献   

12.
Detection of pathogens with single-nucleotide variations is indispensable for the disease tracing, but remains technically challenging. The D614G mutation in the SARS-CoV-2 spike protein is known to markedly enhance viral infectivity but is difficult to detect. Here, we report an effective approach called “synthetic mismatch integrated crRNA guided Cas12a detection” (symRNA-Cas12a) to detect the D614 and G614 variants effectively. Using this method, we systemically screened a pool of crRNAs that contain all the possible nucleotide substitutions covering the -2 to +2 positions around the mutation and identify one crRNA that can efficiently increase the detection specificity by 13-fold over the ancestral crRNA. With this selected crRNA, the symRNA-Cas12a assay can detect as low as 10 copies of synthetic mutant RNA and the results are confirmed to be accurate by Sanger sequencing. Overall, we have developed the symRNA-Cas12a method to specifically, sensitively and rapidly detect the SARS-CoV-2 D614G mutation.  相似文献   

13.
The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.  相似文献   

14.
Severe acute respiratory syndrome coronavirus 2(SARS-CoV-2), the etiologic agent of the current coronavirus disease 2019(COVID-19) pandemic, has evolved to adapt to human host and transmission over the past 12 months. One prominent adaptive mutation is the asparagine-to-glycine substitution at amino acid position 614 in the viral spike protein(D614G), which has become dominant in the currently circulating virus strains. Since spike protein determines host ranges, tissue tropism, and pathogenesis through binding to the cellular receptor of angiotensin converting enzyme 2(ACE2), the D614G mutation is hypothesized to enhance viral fitness in human host, leading to increased transmission during the global pandemic. Here we summarize the recent progress on the role of the D614G mutation in viral replication, pathogenesis, transmission, and vaccine and therapeutic antibody development. These findings underscore the importance in closely monitoring viral evolution and defining their functions to ensure countermeasure efficacy against newly emerging variants.  相似文献   

15.
In mammalian cells, the repair of DNA bases that have been damaged by reactive oxygen species is primarily initiated by a series of DNA glycosylases that include OGG1, NTH1, NEIL1, and NEIL2. To explore the functional significance of NEIL1, we recently reported that neil1 knock-out and heterozygotic mice develop the majority of symptoms of metabolic syndrome (Vartanian, V., Lowell, B., Minko, I. G., Wood, T. G., Ceci, J. D., George, S., Ballinger, S. W., Corless, C. L., McCullough, A. K., and Lloyd, R. S. (2006) Proc. Natl. Acad. Sci. U. S. A. 103, 1864-1869). To determine whether this phenotype could be causally related to human disease susceptibility, we have characterized four polymorphic variants of human NEIL1. Although three of the variants (S82C, G83D, and D252N) retained near wild type levels of nicking activity on abasic (AP) site-containing DNA, G83D did not catalyze the wild type beta,delta-elimination reaction but primarily yielded the beta-elimination product. The AP nicking activity of the C136R variant was significantly reduced. Glycosylase nicking activities were measured on both thymine glycol-containing oligonucleotides and gamma-irradiated genomic DNA using gas chromatography/mass spectrometry. Two of the polymorphic variants (S82C and D252N) showed near wild type enzyme specificity and kinetics, whereas G83D was devoid of glycosylase activity. Although insufficient quantities of C136R could be obtained to carry out gas chromatography/mass spectrometry analyses, this variant was also devoid of the ability to incise thymine glycol-containing oligonucleotide, suggesting that it may also be glycosylase-deficient. Extrapolation of these data suggests that individuals who are heterozygous for these inactive variant neil1 alleles may be at increased risk for metabolic syndrome.  相似文献   

16.
The outbreak of SARS-CoV-2 (SARS2) has caused a global COVID-19 pandemic. The spike protein of SARS2 (SARS2-S) recognizes host receptors, including ACE2, to initiate viral entry in a complex biomechanical environment. Here, we reveal that tensile force, generated by bending of the host cell membrane, strengthens spike recognition of ACE2 and accelerates the detachment of spike’s S1 subunit from the S2 subunit to rapidly prime the viral fusion machinery. Mechanistically, such mechano-activation is fulfilled by force-induced opening and rotation of spike’s receptor-binding domain to prolong the bond lifetime of spike/ACE2 binding, up to 4 times longer than that of SARS-S binding with ACE2 under 10 pN force application, and subsequently by force-accelerated S1/S2 detachment which is up to ~103 times faster than that in the no-force condition. Interestingly, the SARS2-S D614G mutant, a more infectious variant, shows 3-time stronger force-dependent ACE2 binding and 35-time faster force-induced S1/S2 detachment. We also reveal that an anti-S1/S2 non-RBD-blocking antibody that was derived from convalescent COVID-19 patients with potent neutralizing capability can reduce S1/S2 detachment by 3 × 106 times under force. Our study sheds light on the mechano-chemistry of spike activation and on developing a non-RBD-blocking but S1/S2-locking therapeutic strategy to prevent SARS2 invasion.Subject terms: Molecular biology, Structural biology  相似文献   

17.
Exposure histories to SARS-CoV-2 variants and vaccinations will shape the specificity of antibody responses. To understand the specificity of Delta-elicited antibody immunity, we characterize the polyclonal antibody response elicited by primary or mRNA vaccine-breakthrough Delta infections. Both types of infection elicit a neutralizing antibody response focused heavily on the receptor-binding domain (RBD). We use deep mutational scanning to show that mutations to the RBD’s class 1 and class 2 epitopes, including sites 417, 478, and 484–486 often reduce binding of these Delta-elicited antibodies. The anti-Delta antibody response is more similar to that elicited by early 2020 viruses than the Beta variant, with mutations to the class 1 and 2, but not class 3 epitopes, having the largest effects on polyclonal antibody binding. In addition, mutations to the class 1 epitope (e.g., K417N) tend to have larger effects on antibody binding and neutralization in the Delta spike than in the D614G spike, both for vaccine- and Delta-infection-elicited antibodies. These results help elucidate how the antigenic impacts of SARS-CoV-2 mutations depend on exposure history.  相似文献   

18.
In our previous attempt at in vitro selection of a trans - acting human hepatitis delta virus (HDV) ribozyme, we found that one of the variants, G10-68-725G, cleaved a 13 nt substrate, HDVS1, at two sites [Nishikawa,F., Kawakami,J., Chiba,A., Shirai,M., Kumar,P.K.R. and Nishikawa,S. (1996) Eur. J. Biochem., 237, 712-718]. One site was the normal cleavage site and the other site was shifted 1 nt toward the 3'-end. To clarify the interactions between nucleotides around the cleavage site of the trans -acting HDV ribozyme, we analyzed the efficiency of the reaction for every possible base pair between the substrate and the ribozyme at positions -1 (-1N:726N) and +1 (+1N:725N) relative to the cleavage site using the genomic HDV ribozyme, TdS4(Xho), and derivatives of the most active variant, G10-68. These mutagenesis analyses revealed that the +1 base of the substrate affects the structure of the catalytic core in the complex with G10-68-725G, substrate and divalent metal ions, and it shifts the cleavage site. In a comparison with other variants of the trans -acting HDV ribozyme, we found that this cleavage site shift occurred only with G10-68-725G.  相似文献   

19.
The repair of free-radical oxidative DNA damage is carried out by lesion-specific DNA glycosylases as the first step of the highly conserved base excision repair (BER) pathway. In humans, three orthologs of the prototypical endonuclease VIII (Nei), the Nei-like NEIL1-3 enzymes are involved in the repair of oxidized DNA lesions. In recent years, several genome and cancer single-nucleotide polymorphic variants of the NEIL1 glycosylase have been identified. In this study we characterized four variants of human NEIL1: S82C, G83D, P208S, and ΔE28, and tested their ability to excise pyrimidine-derived lesions such as thymine glycol (Tg), 5-hydroxyuracil (5-OHU), and dihydrouracil (DHU) and the purine-derived guanidinohydantoin (Gh), spiroiminodihydantoin 1 (Sp1), and methylated 2,6-diamino-4-hydroxy-5-formamidopyrimidine (MeFapyG). The P208S variant has near wild-type activity on all substrates tested. The S82C and ΔE28 variants exhibit decreased Tg excision compared to wild-type. G83D displays little to no activity with any of the substrates tested, with the exception of Gh and Sp1. Human NEIL1 is known to undergo editing whereby the lysine at position 242 is recoded into an arginine. The non-edited form of NEIL1 is more efficient at cleaving Tg than the R242 form, but the G83D variant does not cleave Tg regardless of the edited status of NEIL1. The corresponding G86D variant in Mimivirus Nei1 similarly lacks glycosylase activity. A structure of a G86D–DNA complex reveals a rearrangement in the β4/5 loop comprising Leu84, the highly-conserved void-filling residue, thereby providing a structural rationale for the decreased glycosylase activity of the glycine to aspartate variant.  相似文献   

20.
Human alpha1-antitrypsin-deficient variants may aggregate in the liver, with subsequent deficiency in the plasma, which can lead to emphysema. The structural and functional characteristics of 10 dysfunctional alpha1-antitrypsin variants (R39C, S53F, V55P, I92N, G115S, N158K, E264V, A336T, P369S, and P369L) were analyzed in detail. Most of them were unstable, as compared to the wild-type molecule, and many of the variants folded into an intermediate form. When five thermostable mutations (T68A, A70G, M374I, S381A, and K387R) were introduced into dysfunctional alpha1-antitrypsin variants, the stabilities and inhibitory activities of most of the variants were restored to levels comparable to those of the wild-type molecule. However, the extremely unstable S53F variant was not stabilized sufficiently by these mutations so as to exhibit function. N158K variant, which carries a mutation in the region critical for the reactive site loop insertion into beta-sheet A, exhibited a reduced level of inhibitory activity, despite conformational stabilization. These results show that aberrant folding caused by conformational destabilization due to mutations can be compensated for by increasing the overall stability of the alpha1-antitrypsin molecule, with exception of a mutation in the highly localized region critical for functional execution.  相似文献   

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

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