Zika virus (ZIKV) exhibits unique transmission dynamics in that it is concurrently spread by a mosquito vector and through sexual contact. We show that this sexual component of ZIKV transmission induces novel processes on networks through the highly asymmetric durations of infectiousness between males and females -it is estimated that males are infectious for periods up to ten times longer than females -leading to an asymmetric percolation process on the network of sexual contacts. We exactly solve the properties of this asymmetric percolation on random sexual contact networks and show that this process exhibits two epidemic transitions corresponding to a core-periphery structure. This structure is not present in the underlying contact networks, which are not distinguishable from random networks, and emerges because of the asymmetric percolation. We provide an exact analytical description of this double transition and discuss the implications of our results in the context of ZIKV epidemics. Most importantly, our study suggests a bias in our current ZIKV surveillance as the community most at risk is also one of the least likely to get tested. A bstract modeling of epidemics on networks remains an active field because some of the most basic features of epidemics are still misunderstood. The classic model is quite simple (1): disease spreads stochastically, with a fixed transmission probability, T , through contacts around a given patient zero. The outbreak dies quickly if T is too small, but spreads to a macroscopic fraction S of the entire population if T is larger than a threshold Tc. At Tc, most of the typical insights from phase transition theory are valuable. For instance, the sizes of microscopic outbreaks follow a power-law distribution such that the expected size of microscopic outbreaks, s , indicates the position of a phase transition. Indeed, as T increases, s monotonically increases, diverges exactly at Tc, and then monotonically goes down; meanwhile the expected macroscopic epidemic size, S, starts increasing monotonically at Tc. However, simple modifications to this model can dramatically alter its phenomenology. The epidemic threshold can vanish in networks with a scale-free degree distribution (2) or in growing networks (3). The phase transition can be discontinuous in the case of complex contagions with threshold exposition or reinforcement (4), interacting epidemics (5, 6), or adaptive networks (7-9). Recently, a unique phenomenon of double phase transitions has also been observed numerically when networks have a very heterogeneous and clustered structure (10, 11). The current Zika virus (ZIKV) epidemic exhibits two unique properties. First, while the main transmission pathway for ZIKV is through a mosquito vector (predominantly Aedes aegypti or Aedes albopictus (13, 14) ), a feature which has its own type of well-studied model and behavior (14-16), it can also spread through sexual contacts (17, 18) . Second, the probability of sexual transmission is highly asymmetric between males and females. Although this is also true for other sexually transmitted infections such as HIV (19), it reaches an extreme level of asymmetry in the case of ZIKV. Indeed, males can be infectious for over 180 days (20) while females are infectious for less than 20 days (21). Assuming a symmetric risk of transmission per contact, males would be 10 times more likely to transmit to a partner than females. This is, however, a rather conservative estimate since male-tofemale transmissions tend to be more likely than the opposite (19, 22) . The dynamics of the ZIKV epidemic is well understood in countries where the vector-borne pathway dominates (23). However, with travelers moving to and from endemic regions, the potential of ZIKV as an emerging STI in regions without the mosquito vector remains to be fully assessed. Indeed, with only few reported cases of sexual transmission of ZIKV -including male-to-male, male-to-female, and female-to-male (18) -the scientific community still struggles to reach a consensus on the impact of sexual transmission of ZIKV (24, 25) . It is therefore imperative to investigate the extent to which canonical knowledge about emerging infectious diseases applies to the threat assessment of ZIKV as an STI. We model the ZIKV sexual transmission through asym- Zika virus (ZIKV) continues to be a threat to countries with conditions suitable for transmission, namely, adequate temperatures and the presence of competent mosquito vectors. Estimates of risk in other countries, based on the sexual transmission of ZIKV, may be underestimated due to inadequate surveillance. Here we formulate random network models of sexual transmission of ZIKV with asymmetric transmission (men being infectious for longer than women) and show that, contrary to previous work, there exists two epidemic thresholds and that certain men-who-have-sex-with-men communities could sustain transmission on their own. Our results also shed light on a new class of processes on random networks by providing a complete analysis of dynamics with multiple critical points. A.A. and L.H.-D. performed the calculations; all authors conceived of the project, interpreted the results, and produced the final manuscript. The authors declare no conflict of interest.  c , corresponds to the global epidemic threshold of the population. The second, T c , corresponds to the invasion of the large heterosexual subpopulation. The threshold T (1) c corresponds to the value of T such that the largest eigenvalue of the Jacobian matrix of Eq. (1) equals 1. The second threshold T (2) c is obtained similarly but with the probability of transmission between homosexual males set to zero. The homo-/bi-/heterosexual subpopulations represent 5%, 3% and 92%, respectively and are equally split between males and females. The degrees are distributed according to a Poisson distribution, metric percolation on random sexual contact networks and solve it exactly using a multitype (multivariate) generating function formalism (26). We then show how the asymmetric percolation leads to a double transition. Interestingly, the formulation of our model allows us to provide a first analytic framework for the aforementioned numerical results on double transitions. More importantly, this allows us to identify two different thresholds for ZIKV to be endemic as an STI in regions where the mosquito vector is absent, but where travellers to/from endemic regions can spark a sexual epidemic when they return/visit. We also find that, in the large interval of parameter space between those two thresholds, the asymmetric percolation creates a core-periphery structure in a system where there was none. Finally, we discuss the implications of this core-periphery structure for the surveillance and control of the ZIKV epidemic, and provide policy guidelines. Inspired by the sexual transmission of ZIKV, we investigate the effect of asymmetry on bond percolation on networks, and show that it yields outcomes akin to the double phase transitions observed numerically in other contexts (10, 11). To isolate the effect of asymmetry alone and thus provide a clear proof of concept, we consider a very simple model in which nodes belong to one of 6 types based on their sex and sexual orientation (i.e., female/male and homo-/bi-/heterosexual). Each node is assigned a number of contacts, k, independently of its type (i.e., all nodes have the same degree distribution {p k } k≥0 ), and links are created randomly via a simple stub-matching scheme constrained by the sexual orientations (12, 26). For instance, bisexual males choose their partners randomly in the pools of heterosexual females, bisexual males and females, and homosexual males. This implies that there is no correlation between the type of a node and its number of contacts, and consequently no core-periphery structure. In fact, this model generates well-mixed contact networks that are indistinguishable from networks generated with the configuration model and the same degree distribution (see Supporting Information). While these networks are originally undirected, asymmetric percolation implies that links can be more likely to exist (i.e., transmit) in one direction than in the other, thus inducing an effective semi-directed structure to the networks (27). In other words, Tij = Tji in general, with Tij being the probability of transmission from a node of type i to a node of type j (hereafter we denote N the set of the six possible types of nodes). In particular, we set Tij = T for every i, j ∈ N except when i corresponds to a female, in which case we set Tij = T /a to enforce asymmetric probabilities of transmission (i.e., females are a times less likely to transmit ZIKV than males). We adapt the formalism presented in Ref. (26) to compute the epidemic threshold and the expected final size of outbreaks in the limit of large networks. It is worth pointing out that since asymmetric percolation (i.e., whenever a = 1) induces an effective semi-directed structure to the networks, the probability for the existence of an extensive connected component (i.e, an epidemic) does not equal to its relative size as for symmetric, traditional bond percolation (i.e., a = 1). Here we focus on the relative size for the sake of conciseness; we refer the readers to the Supporting Information for full details of the analysis and numerical validation. To obtain the relative size of the extensive component, we define vi as the probability that a neighbor of type i is not in the extensive component, which we solve by a self-consistent argument. If the neighbor of a node is not in the extensive component, then none of its other neighbors should be in it either. The probability that the neighbor has a degree equal to k being kp k / k , with k = k kp k , this self-consistent argument can be written as [1] where α j|i is the probability that a neighbor of a node of type i is of type j (i.e., j α j|i = 1 for any i). Solving this equation for every i ∈ N , the probability that a node of type i is part of the extensive component, Si, corresponds to the probability that at least one of its neighbors is in it as well [2] The relative size of the extensive component is then S = i∈N wiSi, where wi is the fraction of the nodes that are of type i. Below the epidemic or percolation threshold, every vi is equal to 1 since there is no extensive component. The percolation threshold corresponds to the point where the largest eigenvalue of the Jacobian matrix of Eq. (1) equals 1. The distribution of the composition of the small, nonextensive components can be computed in a similar fashion (see Supporting Information for full details). Let us define the probability generating function (pgf) Hi(x) whose coefficients correspond to the probability that a neighbor of type i leads to a small component of a given composition (i.e., the number of nodes of type j is given by the exponent of xj). Invoking the same self-consistency argument as above, the pgfs are the solution of where the extra xi has been added to account for the neighbor of type i itself. Similarly, the small component that can be reached from a node of type i is therefore given by [4] The distribution of the composition of the small components is K(x) = i∈N wiKi(x). It is worth noting that whenever S > 0, the distribution generated by K(x) is no longer normalized, K(1) < 1, such that the average number of nodes of type i in the small components is An example of the general phenomenology is shown in Fig. 1 . Unlike the classic epidemic transition picture, where s diverges at the epidemic threshold where the macroscopic epidemic emerges, we now find two peaks in s . This double transition is similar to numerical results from Ref. (10), but here observed without the need for either strong clustering nor heterogeneity in degree distribution. In fact, we used the homogeneous Poisson degree to ensure that the asymmetry in the transmission is the only salient feature of the model. Interestingly, as shown in Fig. 2, T (1) c and T (2) c are virtually equal for small values of the asymmetry. As asymmetry increases, the peak separates thus yielding a double transition corresponding to an effective core-periphery organization in the network of infections. The core then corresponds to the men-having-sex-with-men (MSM) population where infections are more frequent than in the remaining population. Figure 3 shows the network of who infected whom for two values of c , the extensive component is mostly composed of one type of nodes and any spillover in the other types quickly dies out. However, at T = T (2) c , these spillovers now cause cascades into other types with truncated power-law distributed sizes (see Fig. 4 ). For T > T (2) c , the extensive component recovers the well-mixed structure of the original underlying network. Altogether, the second peak in the average size of outbreaks, s , corresponds to a transition between subcritical and supercritical spillover in a less susceptible sub-population, but not to a second phase transition in the classic sense. Indeed, the analytical nature of our results allows confirm the null critical exponent observed in Ref. (10) for the scaling of the height of the second susceptibility peak with regards to system size. Even in the infinite system considered by our calculations, the peak saturates, which is the only possible outcome for a system whose order parameter is already non-zero. Interestingly, a critical power-law-like behavior is nonetheless observed in the heterosexual population at both thresholds. Moreover, our results suggest that the asymmetry in transmission probability is reflected in the asymmetric prevalence within the male and c , the power-law tail in the component size distribution is mainly due to the critical core of homosexual males while the exponential behavior is mainly due to heterosexuals. The power-law tail in the distributions of heterosexuals is due to spillovers from the critical core. (lower right) At T (2) c , the power-law portion of the distribution is due to heterosexuals now forming a critical core while homosexual males, being already almost exclusively in the extensive component, do not contribute. All curves were obtained by solving Eq. (4) with asymmetry a = 3 and the parameters given in the caption of Fig. 1. female heterosexual populations, which is reminiscent of recent empirical results (28). Based on our results, we can summarize the phase diagram of the ZIKV epidemic in 3 possible outcomes. First, with T < T (1) c , all outbreaks are microscopic, quickly die out, and mostly infect MSM. Second, with T (1) c (2) c , we now see a macroscopic epidemic within the network of homosexual contacts between males, with microscopic spill-over into the rest of the population via bisexual males. Third, with T > T (2) c , we now find a more classic epidemic scenario in the sense that it is of macroscopic scale in most of the population. It is also worth mentioning that this phenomenology is robust to the presence of multiple infectious seeds sparking outbreaks (see Supporting Information). Our results are thus valid beyond ZIKV for any infections with asymmetry in probabilities of direct transmission, regardless of whether or not there is also vector transmission. We developed a network model of ZIKV transmission highlighting the importance of asymmetric sexual transmission between males and females. We find a double transition generated by a core group of MSM that could maintain ZIKV transmission without the presence of a viable mosquito vector, such as in regions where people may have brought back ZIKV with them after a trip to endemic regions. These results are unique as previous models showing double transitions relied on the need for strong clustering and heterogeneity in degree distribution. Our study carries important consequences for the ongoing ZIKV epidemic and stresses the large knowledge gap in the sexual transmission of ZIKV (25). The aim of our work is to present the epidemiological consequences of possible sustained sexual transmission. While there are many unknowns, recent work demonstrates a) multiple anecdotal cases of sexual transmission of ZIKV in humans (25, (29) (30) (31) , b) multiple separate animal models demonstrating sexual transmission (32-34), c) strong asymmetries between durations of ZIKV shedding in semen and vaginal secretions (20, 21), and d) differential risk between sexes for ZIKV infection in sexually active populations. Indeed, recent work has identified 90% more ZIKV infections in women between 15 and 65 than men of the same age in Rio de Janero (28) adjusted for gender-related healthseeking behavior and pregnancy status. Importantly, this risk difference was not seen in women less than 15 years of age or greater than 65, indicating the potentially large impact of sexual transmission of ZIKV in a country with known ongoing vectored transmission of ZIKV. A similar situation has also been observed in Colombia (35) and in the Dominican Republic (36). While more research on the epidemiological impacts and basic biology of sexual ZIKV transmission is needed, there is compelling need to be prepared with epidemiological studies examining transmission on a population-scale. We demonstrated that potential ZIKV persistence in MSM, even if barely critical within that sub-population, could cause subcritical but dramatic spillover into the heterosexual community. ZIKV infections in adults are largely asymptomatic (37) and, therefore, most testing occurs in the roughly 20% of cases that are symptomatic or in individuals seeking to have children (38). The vast majority of these individuals will be outside of the MSM community (38). This means that the community most at risk is also one of the least likely to get tested. To avoid underestimating the spread of ZIKV, it is therefore important for health officials and policy makers to keep its unique behavior and phenomenology in mind. Given the extent of foreign travel to locations endemic with ZIKV, public health practitioners should be aware of the potential for infectious introduction into local MSM communities. Travel history as well as sexual history should be employed when evaluating an occult fever. Cities which have a viable vector for ZIKV should be doubly aware of the potential transmission routes of ZIKV. As it stands, current estimates of the basic reproductive number, R0, of ZIKV may be too low as they fail to account for sustained sexual transmission (17, 18, 39, 40) . Important future work will be to accurately estimate R0 of ZIKV across various settings with differing sexual practices and mosquito fauna. ACKNOWLEDGMENTS. A.A. acknowledges support from the
In this section, we will briefly describe the model present in Asymmetric percolation drives a double transition in sexual contact networks REF , the basis for our work.
Asymmetric percolation drives a double transition in sexual contact networks
Proceedings of the National Academy of Sciences USA 114, 8969-8973 (2017)
Physics Biology Medicine