Literature DB >> 25793969

Modes of interaction between individuals dominate the topologies of real world networks.

Insuk Lee1, Eiru Kim1, Edward M Marcotte2.   

Abstract

We find that the topologies of real world networks, such as those formed within human societies, by the Internet, or among cellular proteins, are dominated by the mode of the interactions considered among the individuals. Specifically, a major dichotomy in previously studied networks arises from modeling networks in terms of pairwise versus group tasks. The former often intrinsically give rise to scale-free, disassortative, hierarchical networks, whereas the latter often give rise to single- or broad-scale, assortative, nonhierarchical networks. These dependencies explain contrasting observations among previous topological analyses of real world complex systems. We also observe this trend in systems with natural hierarchies, in which alternate representations of the same networks, but which capture different levels of the hierarchy, manifest these signature topological differences. For example, in both the Internet and cellular proteomes, networks of lower-level system components (routers within domains or proteins within biological processes) are assortative and nonhierarchical, whereas networks of upper-level system components (internet domains or biological processes) are disassortative and hierarchical. Our results demonstrate that network topologies of complex systems must be interpreted in light of their hierarchical natures and interaction types.

Entities:  

Mesh:

Year:  2015        PMID: 25793969      PMCID: PMC4368763          DOI: 10.1371/journal.pone.0121248

Source DB:  PubMed          Journal:  PLoS One        ISSN: 1932-6203            Impact factor:   3.240


Introduction

Networks occur in many contexts in the real world, and the topologies of these real-world networks have potentially large practical impact in areas of biological, social, and technical importance. Topological implications of human social networks, for example, influence public policies (such as for the development of effective vaccination schemes [1]) and business strategies (such as for allocating marketing resources[2]). Likewise, the topology of the Internet affects routing protocols for robust and cost-effective communications [3], and in biology, the topologies of protein-protein interaction (PPI) networks have informed our understanding of cells and organisms [4]. As a consequence, network topologies have been extensively characterized with respect to their global topological properties, such as node degree distributions [5], node hierarchical organization [6], and assortativity (the degree correlation between connected nodes) [7]. Differences in such properties have been noted for many real world networks [8]. In the course of studying networks, we realized that many of these historical observations of contrasting network topologies could be explained by a simplifying model: that most real world networks can be categorized as one of two major classes of networks – those capturing intrinsically pairwise activities (e.g., dating or pairwise physical interactions between proteins) and those capturing intrinsically group activities (e.g., boards of directors of companies or membership in the same protein complexes). In this paper, we demonstrate that this distinction explains many of the major topology differences amongst social networks, the Internet, and biological networks, and that networks generated by the same class of activities – regardless of the precise nature of those activities – often have similar topological properties.

Materials and Methods

Global topological analyses of networks were performed as previously described for node degree distribution [5], assortativity [7,9], graphlet frequency distribution [10], and node hierarchical organization [6]. Null-model random networks for correlation profiles of assortativity test were generated by local rewiring algorithm that randomizes a network yet conserves degrees of each node [9,11]. Biological processes were defined by hierarchical clustering of YeastNet described as in [12] or by MCL clustering [13] with the granularity parameter selected so as to balance modularity and proteome coverage. For the GO biological processes network, we connected pairs of GO terms sharing at least one annotated yeast protein to generate a network of 5,587 edges among 1,066 GO biological process terms.

Results and Discussion

An intrinsic dichotomy between contact- and task-centric networks

We illustrate this key distinction among the two network classes in by introducing toy examples of two types of human social interactions similarly composed of 11 people (nodes) organized into three groups (indicated by node colors). Interpersonal relationships (edges) might be based on direct personal contact—the contact-centric model ()—such as for the cases of online dating [14] or sexual contacts [15], or alternatively based on sharing roles to perform a common task—the task-centric model ()—such as for sharing membership on company boards [16] or co-authorship of scientific papers [17]. In the contact-centric network, a few attractive individuals (represented as squares) have a large number of partners (). In contrast, a task-centric network is characterized by group activities in which the pairwise interaction reflects the tendency for individuals to participate in the same groups. Note that individuals may participate in the same task but never actually directly contact each other (represented as dotted lines)—e.g., many coauthors for scientific papers do not have a personal relationship (). The networks’ topologies can be characterized by various global topological properties [8]. Here, we will consider the three most widely studied topological properties: node degree distribution, node hierarchical organization, and assortativity.

Illustrative examples of network dichotomy.

Two broad classes of networks naturally occur in real world settings based on whether interactions between entities are intrinsically pairwise (left panel) or based on groups (right panel). Examples in human social networks might include (a) dating or sexual relationships, and (b) board memberships or co-authorships of scientific papers. (c and d) indicate small toy examples for these two types of human social networks, each composed of 11 individuals (nodes) organized into three groups (indicated by node colors). Edges in (c) indicate direct contacts between persons; edges in (d) indicate participation in shared tasks. Numbers in (c) denote the degree connectivity (count of associated edges) and clustering coefficient [18] for each node. Individuals with high node degree are marked as squares. In spite of their simplicity, the two toy networks show distinct topological features. Their cumulative probability distributions of nodes with >k degree (Pc(k)) differ; (e) is scale-free, (f) is single- or broad-scale. They differ in being (g) disassortative or (h) assortative networks, as seen by heat-map representations of the enrichment for connections between nodes of varying degrees or measured by Pearson correlation coefficient (r) of the degrees at either ends of an edge. Finally, they exhibit (i) hierarchical or (j) nonhierarchical network topologies, as judged by their relationships between node degree connectivity (k) and node clustering coefficients (C(k)). An examination of the topologies of these toy networks reveals striking differences. In general, the distribution of node degree of scale-free networks approximately follows a power-law [19], Pc(k) = k, where Pc(k) represents the cumulative probability of having nodes with >k degree. The contact-centric network is expected to be more scale-free in nature, because there is no constraint limiting the number of links for individual nodes by this modeling perspective, permitting few nodes with many partners, while most nodes have only few partners (). In contrast, the task-centric network is expected to be of single- or broad-scale [5] (). Single-scale networks are characterized by a fast-decaying tail of degree distribution, such as exponential decay. Broad-scale networks are often characterized by a power-law regime followed by exponential trimming of degree, represented by the function Pc(k) = k), where k is the cutoff degree for exponential decay. In general, networks following exponential decay or truncated power-law function exhibit an increased proportion of nodes with medium-connectivity, probably due to the constraint on the number of links by sizes of task-groups of the networks, resulting in reduced numbers of both hubs and nodes with low-connectivity. Because of the connections among entities performing the same tasks, the task-centric network shows a degree preference corresponding to the preferred sizes of the groupings. Thus, these two network classes intrinsically give rise to distinct node degree distributions. We see an equally strong dichotomy for a second well-characterized topological property that of the correlation in degree between connected nodes. A given network is termed assortative if its hub-hub connections are enriched; if depleted, the network is termed disassortative [7]. The assortativity of networks can be measured simply by Pearson correlation coefficient (r) of the degrees at either ends of an edge [7], where 09], and can be visualized as a heat map. Notably, the contact-centric network example is disassortative (r = −0.66, ), while the task-centric network example is assortative (r = 0.32, ). The third major topological property that we consider is the hierarchical organization [6]. In a hierarchically organized network, hub components bridge many disconnected regions of the network to efficiently coordinate all system components. The number of these far-reaching connections decreases as the node degree decreases, entailing increasing proportions of connections toward local neighbors. This indicates a hierarchical contribution of system components with hub components generally sharing high betweenness (a measure of network centrality). As a consequence, the clustering coefficients C(k)[18] of nodes decrease as their node degree k increases in hierarchically organized networks. It is noticeable that this definition of hierarchy assumes that real world networks are in general scale-free and modular. Thus, disassortative networks are also expected to be hierarchical networks (i.e., assortativity and hierarchy are not independent topological features). There are other hierarchy models which are independent from network modularity [20,21]. However, the hierarchy by C(k) decreasing with k has been widely used for previous network topological studies we discussed in this paper. Many real world networks exhibit this signature of hierarchical organization [6], as does the contact-centric toy network example in . Nonhierarchical networks, however, exhibit roughly equal clustering coefficients, regardless of node degree. This trend suggests enriched connections between components with similar numbers of network neighbors, and the task-centric network example is accordingly nonhierarchical (). Thus, these two major classes of networks, characterized by the nature of the relationships among the interacting entities, show distinct topological properties: a contact-centric network is a scale-free, hierarchical, disassortative network; a task-centric one is a single- or broad-scale, nonhierarchical, assortative network. The analyses of these toy examples suggest that the contrasting topological properties of many real world networks might also stem from this intrinsic dichotomy of network type. We thus next examined real world networks to test this hypothesis (listed in ).

Real-world networks exhibit the same dichotomy

We first compared the three global topological properties for two human social networks, analyzing a contact-centric online dating network [14] and a task-centric network of boards of directors of US companies [16] (). The node degree distribution of the dating network showed more scale-free character, while that of the board of directors network showed single-scale (). We observed contrasting assortativity between the two networks: the dating network is disassortative (r = −0.05), but the board of directors network is assortative (r = 0.27) (). Similarly, we found a clearly hierarchical organization for the dating network but a nonhierarchical organization for the director board network (). Thus, as for the toy network examples in , real world human social networks also show contrasting topologies according to the distinction between intrinsically pairwise vs. group activities. In the online dating network, most people online date only a few partners, while a few individuals make large numbers of online dating contacts. In contrast, membership on boards of directors carries significant obligations. Thus, directors typically participate in only a limited number of boards, and boards are limited in size; as a consequence, clustering coefficients are often similar between members belonging to large boards and those who belong to small ones.

Dichotomy of real-world networks.

Topological analysis of three real world networks (). We analyzed contact-centric networks (indicated by the letter C) and task-centric networks (indicated by the letter T) for their (a, d, g) degree distribution, (b, e, h) assortativity, and (c, f, i) hierarchical structure. We tested two different human social networks, a contact-centric online dating network (Dating) [14] and a task-centric network of shared membership on US company boards of directors (Director board) [16]. Similarly, we analyzed the Internet at two different levels, the router-level (Router) and domain level (Domain). The router-level Internet is task-centric and the domain-level Internet is contact-centric. Finally, we analyzed two alternate protein networks, testing pairwise protein interactions (CCSB-YI1) [22] and functional protein interactions (YeastNet core) [23], as examples of contact-centric and task-centric protein networks, respectively. We measured the cumulative probability of nodes with degree >k (Pc(k)) for the full range of node degree [19] to test scale-freeness of networks. We measured hierarchical connectivity by testing for decreasing clustering coefficients (C(k)) as a function of increasing node degree (k). Network assortativities were visualized as heat maps of the enrichment of connections between various ranges of degrees compared to permuted networks [9]; red indicates enrichment and blue indicates suppression of connectivity. In every case, real-world networks showed topological properties consistent with being either contact- or task-centric, as appropriate. The Internet represents another widely studied real world complex system comprising computers and other devices with IP addresses that communicate via a network of routers, each implementing routing protocols. Internet scientists have considered two major networks, those composed of routers and those of domains (or autonomous systems) (). A domain is an entity containing multiple routers under control of a network operator(s) (e.g., Internet Service Provider; ISP) using a common routing policy to the Internet. Thus, each domain is considered to be an operational unit component of the Internet. For each domain, within-domain routing protocols and between-domain protocols such as the Border Gateway Protocol (BGP) can be implemented within one or more separate routers. Therefore, the Internet, when considered at the domain-level, maps only peer-to-peer connections between BGP routers, while the router-level maps all router connections including both within-domain and between-domain connections. Thus, the domain-level Internet model can be considered to be a contact-centric network, while the router-level Internet model resembles the task-centric model. Accordingly, analysis of published router network and domain network [24] revealed the topological properties consistent with the contact-centric and task-centric models, just as for the human social networks. The contact-centric domain network is scale-free, disassortative (r = −0.17), and hierarchical, but the task-centric router network is broad-scale, relatively assortative (r = −0.01), and nonhierarchical (). Disassortativity of a domain-level Internet has been also observed in a previous study [11]. Routers with more than 300 degree connections do show exceptionally low clustering coefficients (, highest degree nodes). However, this pattern only affects a small number of routers at the highest level domain, and the overall two levels of organization share topological characteristics consistent with the trends expected for task-centric models. This trend arises because the Internet has top-level domains for mediating communication between lower-level domains. Between-domain routers of the top-level domains generally have connections to a large number of lower-level domains, resulting in much lower clustering coefficients for these few hub routers than for the rest of routers in the Internet. Thus, the Internet, considered at the level of organization for the entire router network, is nonhierarchical. Similarly, biologists have measured interactions among cellular proteins, describing complex biological networks. These networks naturally fall into contact or task centric models. For example, measurements employing the yeast two-hybrid (Y2H) methodology map the direct physical pairwise contacts between proteins. In contrast, proteins participate in groups—in pathways and complexes—in order to fulfill their functional roles within cells, and physical interactions are but one indication of functional association. We can accordingly model PPI networks considering either direct physical contacts (contact-centric) or functional associations (task-centric) between proteins (). We analyzed a high-confidence genome-wide Y2H map of yeast proteins, CCSB-YI1 [22] as a representative physical interaction network, and a published functional protein network (YeastNet core) [23] as a representative functional network. Consistent with the contact/task dichotomy, the physical network is scale-free [25], disassortative [9] (r = −0.12), and hierarchical [25], while the functional network is broad-scale, assortative (r = 0.45), and nonhierarchical (). Thus, in protein networks, as in social and technology networks, the contact/task centric models explain the dominant topological properties. We also observed a dichotomy among the real world networks for additional topological properties. One useful characteristic of network topology is the frequency distribution of the complete set of graph topologies with three, four, and five nodes, dubbed 3,4,5-node graphlets. The graphlet frequency distributions (GFD) [10] of the real world networks examined here show that task-centric networks have significantly higher frequencies for 5-node graphlets than contact-centric networks (). We also examined the distribution of shortest path lengths of the real world networks, and observed significantly higher ranges for the shortest path lengths and network diameters of task-centric real-world networks (). Taken together, we conclude that the task-centric modeling perspective generated more community structures, which increased frequencies of 5-node graphlets and shortest path lengths of the networks. These models can account for contrasting observations from previous topology studies for a variety of real world networks. Typically, many PPI networks have been implicitly modeled from the contact-centric view, while the task-centric view has been implicitly adopted for modeling many human social networks. Thus, PPI networks are generally thought to be scale-free disassortative hierarchical networks, and the human social networks as broad-scale assortative nonhierarchical networks. However, there have been reports of noncanonical topological properties for protein functional networks [23,26] and online dating networks [14]. These inconsistencies are easily explained as resulting from the contact/task centric dichotomy. Mixtures of these two models, for example by combining Y2H protein interactions and protein complexes into a single network model, may have led to inconsistent observations regarding the suppression of hub-hub connections in different PPI networks [9,27]. In fact, as we next show, the choice of how the same entities are represented in a network model markedly affects the network topology.

Alternate network representations of the same entities can also exhibit this dichotomy

A study of the toy examples in immediately suggests that different representations of the associations among the same entities could give rise to networks with either topology. We see that this is indeed the case in real-world networks and can be easily illustrated for both the internet and protein networks. In both cases, the topological properties of the resulting alternate network models match those expected for contact and task centric networks. For example, for the case of the internet, the router and domain networks differ in not only edge representation but also node representation—the task-centric router-level Internet has nodes of individual routers, whereas the contact-centric domain-level Internet has nodes of domains composed of multiple routers. As seen above (; reprinted in for clarity), the domain network is scale-free, disassortative, and hierarchical, but the router network is broad-scale, assortative, and nonhierarchical. This observation of contrasting Internet topology between router- and domain-level also has been reported by [24]. Likewise, cellular proteins can be modeled by grouping according to pairwise protein interactions or, alternatively, by clustering the proteins into functional modules or biological processes. We therefore defined biological processes by grouping proteins that were highly connected to one another in a functional network. From the resulting 333 biological processes defined by hierarchical clustering, we generated a network of processes as described in [12]. The resultant process network was revealed to be disassortative (r = −0.15) () and hierarchical (). Regardless of node degree, most proteins and routers showed clustering coefficients near the average clustering coefficient for the entire protein or router network. Notably, our findings were consistent when biological process networks were created using alternate methods (e.g. using an alternate clustering algorithm, MCL [13], or using pre-existing process annotations from the Gene Ontology [28]) (). Thus, these observations thus suggest common rules-of-organization between the Internet and the yeast proteome at equivalent system levels—assortative and nonhierarchical organization among lower-level components, and disassortative and hierarchical organization among upper-level components ().

Dichotomy of the same entities by alternative network representations.

Alternate representations of the same network can lead to different topologies, especially for networks with natural hierarchical organization. We illustrate this tendency for (a) the Internet and (b) the yeast cell proteome. Each can be modeled by networks at two different granularities, representing nodes either as upper level components (internet domains or protein processes) or lower level components (internet routers or individual proteins). For the internet, previous Internet mapping studies provide both a router-level network and a domain-level network [24]; each domain is composed on multiple routers, and domains are connected via between-domain routers. For the protein network, we defined protein processes by hierarchically clustering proteins based on their pairwise interactions as in [12]. A total of 333 biological processes were defined and connections between processes were defined based on pairwise interactions between proteins within each process. The networks’ hierarchical structure was analyzed and plotted as in , marking the mean clustering coefficient for each entire network as a horizontal solid line in the plot. The non-hierarchical router and protein networks generally exhibited clustering coefficients near this average regardless of node degree, although for the Internet router-level network, routers with >300 connections showed exceptionally low clustering coefficients primarily due to a small number of between-domain routers located at a few top-level domains of the Internet. This intrinsic dichotomy between pairwise and group-based interactions helps to clarify many previous network analyses. It has been argued that the enrichment of hub-hub connections in the board of directors network may suggest the existence of a super-group of decision-makers that influences the business of the entire nation [16]. However, given that this network represents a task-centric model, we would expect this property to be an intrinsic feature of the model. In contrast to the board of directors network, hub-hub connections are suppressed in the contact-centric online dating network (and similarly in sexual contact networks [29]). This trend has been argued as an example of the robustness of social networks, since enhanced contacts between hubs would dramatically increase the risk to society, as for example, might result from hub-hub transmission of sexually transmitted diseases [15]. Similarly, physical protein interaction networks have been argued to exhibit robustness on the basis of their disassortative topology by localizing the effects of deleterious protein perturbations [9]. However, in both of these cases, a consideration of contact vs. task centric networks reveals that the rule of dissortativeness-robustness applies only to contact-centric networks, as task-centric networks are intrinsically assortative. Finally, we note that although our focus here is on real world networks, researchers often generate random networks for study or to serve as background models for comparing to real world networks. Such randomized graphs can of course be generated by many different algorithms. We would expect the method employed to dictate the resulting topological properties. For example, one common approach is to construct graphs de novo by arbitrarily connecting randomly chosen node pairs; this approach essentially implements a contact-centric generative model, and we would expect the resulting random network’s topological properties to reflect this mechanism of generation.

Conclusions

In summary, we describe two broad classes of network models whose intrinsic topological properties explain observations on many real world networks, as we have shown with examples from human societies, the Internet, and cellular proteins. Contact-centric networks are characterized by intrinsically pairwise interactions, generating scale-free, disassortative, hierarchical networks. Task-centric networks are characterized by interactions within and between groups, and give rise to single- or broad-scale, assortative, nonhierarchical networks. Alternative representations of the same entities, such as by grouping individuals together hierarchically, can also give rise to such networks, with lower-level representations generating assortative, nonhierarchical networks, and upper-level, clustered representations generating disassortative, hierarchical networks. The intrinsic topologies of these two classes of networks may account for previous paradoxical observations among real world networks, and imply similar organizational principles among human societies, the Internet, and cellular proteomes.

Cumulative degree distribution plots for real world networks.

Three contact-centric (indicated by the letter “C”) networks and three task-centric (indicated by the letter “T”) networks were analyzed, highlighting the degree distribution tails. (TIF) Click here for additional data file.

Graphlet frequency distributions (GFD) for real world networks.

Four contact-centric (indicated by the letter “C”) networks and three task-centric (indicated by the letter “T”) networks and their randomized networks were analyzed. (TIF) Click here for additional data file.

Shortest path length distributions for real world networks.

Four contact-centric (indicated by the letter “C”) networks and three task-centric (indicated by the letter “T”) networks and their randomized networks were analyzed. (TIF) Click here for additional data file.

Heat map plotting results of assortativity analysis for the biological process network.

(TIF) Click here for additional data file.

Hierarchy analysis plots for alternate protein biological process networks.

Both a network of previously defined Gene Ontology biological processes (in which edges are defined between 7th level GO-BP annotations based on sharing of annotated proteins; labeled “Process level – GO” in the plot) and the network of biological processes inferred by MCL clustering of a functional protein network (“Process level – MCL”) show hierarchical organization, while the functional protein network itself (“YeastNet core”) shows a nonhierarchical organization, in which clustering coefficients for a given degree connectivity (C(k)) trend near the average clustering coefficient for the entire network (). (TIF) Click here for additional data file.

Summary of the networks analyzed in this study.

(PDF) Click here for additional data file.
  22 in total

1.  Specificity and stability in topology of protein networks.

Authors:  Sergei Maslov; Kim Sneppen
Journal:  Science       Date:  2002-05-03       Impact factor: 47.728

Review 2.  Sexual networks: implications for the transmission of sexually transmitted infections.

Authors:  Fredrik Liljeros; Christofer R Edling; Luis A Nunes Amaral
Journal:  Microbes Infect       Date:  2003-02       Impact factor: 2.700

Review 3.  Network biology: understanding the cell's functional organization.

Authors:  Albert-László Barabási; Zoltán N Oltvai
Journal:  Nat Rev Genet       Date:  2004-02       Impact factor: 53.242

4.  Assortative mixing in networks.

Authors:  M E J Newman
Journal:  Phys Rev Lett       Date:  2002-10-28       Impact factor: 9.161

5.  Information dynamics shape the sexual networks of Internet-mediated prostitution.

Authors:  Luis E C Rocha; Fredrik Liljeros; Petter Holme
Journal:  Proc Natl Acad Sci U S A       Date:  2010-03-15       Impact factor: 11.205

6.  A probabilistic functional network of yeast genes.

Authors:  Insuk Lee; Shailesh V Date; Alex T Adai; Edward M Marcotte
Journal:  Science       Date:  2004-11-26       Impact factor: 47.728

7.  High-quality binary protein interaction map of the yeast interactome network.

Authors:  Haiyuan Yu; Pascal Braun; Muhammed A Yildirim; Irma Lemmens; Kavitha Venkatesan; Julie Sahalie; Tomoko Hirozane-Kishikawa; Fana Gebreab; Na Li; Nicolas Simonis; Tong Hao; Jean-François Rual; Amélie Dricot; Alexei Vazquez; Ryan R Murray; Christophe Simon; Leah Tardivo; Stanley Tam; Nenad Svrzikapa; Changyu Fan; Anne-Sophie de Smet; Adriana Motyl; Michael E Hudson; Juyong Park; Xiaofeng Xin; Michael E Cusick; Troy Moore; Charlie Boone; Michael Snyder; Frederick P Roth; Albert-László Barabási; Jan Tavernier; David E Hill; Marc Vidal
Journal:  Science       Date:  2008-08-21       Impact factor: 47.728

8.  Hierarchy measure for complex networks.

Authors:  Enys Mones; Lilla Vicsek; Tamás Vicsek
Journal:  PLoS One       Date:  2012-03-28       Impact factor: 3.240

9.  Stratus not altocumulus: a new view of the yeast protein interaction network.

Authors:  Nizar N Batada; Teresa Reguly; Ashton Breitkreutz; Lorrie Boucher; Bobby-Joe Breitkreutz; Laurence D Hurst; Mike Tyers
Journal:  PLoS Biol       Date:  2006-10       Impact factor: 8.029

10.  The importance of contact network topology for the success of vaccination strategies.

Authors:  Junling Ma; P van den Driessche; Frederick H Willeboordse
Journal:  J Theor Biol       Date:  2013-01-29       Impact factor: 2.691

View more
  2 in total

1.  Network-assisted investigation of virulence and antibiotic-resistance systems in Pseudomonas aeruginosa.

Authors:  Sohyun Hwang; Chan Yeong Kim; Sun-Gou Ji; Junhyeok Go; Hanhae Kim; Sunmo Yang; Hye Jin Kim; Ara Cho; Sang Sun Yoon; Insuk Lee
Journal:  Sci Rep       Date:  2016-05-19       Impact factor: 4.379

2.  Massive migration promotes the early spread of COVID-19 in China: a study based on a scale-free network.

Authors:  Wen-Yu Song; Pan Zang; Zhong-Xing Ding; Xin-Yu Fang; Li-Guo Zhu; Ya Zhu; Chang-Jun Bao; Feng Chen; Ming Wu; Zhi-Hang Peng
Journal:  Infect Dis Poverty       Date:  2020-08-10       Impact factor: 4.520

  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.