ACE2 and TMPRSS2 variation in savanna monkeys (Chlorocebus spp.): Potential risk for zoonotic/anthroponotic transmission of SARS-CoV-2 and a potential model for functional studies

doi:10.21203/rs.3.rs-27287/v1

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has devastated health infrastructure around the world. Both ACE2 (an entry receptor) and TMPRSS2 (used by the virus for spike protein priming) are key proteins to SARS-CoV-2 cell entry, enabling progression to COVID-19 in humans. Comparative genomic research into critical ACE2 binding sites, associated with the spike receptor binding domain, has suggested that Old World primates may also be susceptible to disease from SARS-CoV-2 infection. Savanna monkeys (Chlorocebus spp.) are a widespread non-human primate with well-established potential as a bi-directional zoonotic/anthroponotic agent due to high levels of human interaction throughout their range in sub-Saharan Africa and the Caribbean. To characterize potential functional variation in savanna monkey ACE2 and TMPRSS2, we inspected recently published genomic data from 245 savanna monkeys, including 163 wild monkeys from Africa and the Caribbean and 82 captive monkeys from the Vervet Research Colony (VRC). We found several missense variants. One missense variant in ACE2 (X:14,077,550; Asp30Gly), common in Ch. sabaeus, causes a change in amino acid residue that has been inferred to reduce binding efficiency of SARS-CoV-2, suggesting potentially reduced susceptibility. The remaining populations appear as susceptible as humans, based on these criteria for receptor usage. All missense variants observed in wild Ch. sabaeus populations are also present in the VRC, along with two splice acceptor variants (at X:14,065,076) not observed in the wild sample that are potentially disruptive to ACE2 function. The presence of these variants in the VRC suggests a promising model for SARS-CoV-2 infection and vaccine and therapy development. In keeping with a One Health approach, characterizing actual susceptibility and potential for bi-directional zoonotic/anthroponotic transfer in savanna monkey populations may be an important consideration for controlling COVID-19 epidemics in communities with frequent human/non-human primate interactions that, in many cases, may have limited health infrastructure.

Computational Biology

Zoonoses

anthroponoses

zoonoses

COVID-19

SARS-CoV-2

ACE2

TMPRSS2

Chlorocebus

vervet

The emergence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV–2) in China in December 2019 has since led to a devastating human pandemic of COVID–19[1], the potentially fatal respiratory disease resulting from infection with the novel coronavirus [2]. Due to the ongoing severity of the pandemic in humans, little attention has yet focused on human/animal transmission beyond the initial presumed zoonotic event with a bat and/or pangolin host [3]. Although potential anthroponotic transmissions to domesticated animals has been detected, these transmissions have been relatively rare (e.g., ferrets, cats, dogs) [4] and human-to-human transmission is the primary risk [5] [6]. Aspects of the biology of SARS-CoV–2 infection call for stricter scrutiny of such transmissions among non-human primates kept as pets, among wild non-human primates, and other populations that may have frequent exposure to potentially infected humans.

The angiotensin-converting enzyme 2 (ACE2), previously identified as the key receptor for the SARS coronavirus [7], has since also been identified as one of the key binding receptors by which SARS-CoV–2 gains entry to cells [8] [9] [10]. The spike protein of SARS-CoV–2 contains a variable receptor binding domain (RBD) which binds to the ACE2 receptor facilitating entry of the virus into cells. Mutations at critical sites in the RBD appear to alter the binding affinity to ACE2 receptors, facilitating spill-over events into other species. SARS-CoV–2 likely originated in bats with zoonotic spill-over, via an intermediate host, to humans [3] [11]. In the same way that mutations in the spike protein have altered binding efficacy to the ACE2 receptor, evidence suggests that mutations or natural variants in the ACE2 receptor may influence susceptibility or resistance of a species to infection. Identifying key receptors associated with infectivity and transmission will be important for understanding host range and potentially susceptible hosts. It has been shown that changes in key residues of ACE2 influence the receptor efficiency and consequently entry of severe acute respiratory syndrome coronavirus 1 (SARS-CoV) in to cells [12]. To gain insight into host range for SARS-CoV–2 based on receptor usage, alignment of predicted amino acid sequences for human, non-human primates and selected domestic animals showed shared identity at binding regions and critical residues [13].

Recent research has supported the idea that, due to ACE2 receptor congruence–particularly at critical binding sites for the SARS-CoV–2 spike receptor binding domain–catarrhine primates may be as susceptible to infection with COVID–19-like disease as humans [13] [14] [15], presenting the possibility of bi-directional human-non-human primate infections in areas with human-wildlife contact. Although rhesus macaques (Macaca mulatta) have been experimentally infected and do show COVID–19-like symptoms upon infection [16] [17] [18] zoonotic/anthroponotic potential in wild populations has not yet been investigated. This is cause for great concern, not just from a global and community health perspective, but also from a conservation perspective, as many critically endangered non-human primates may also be vulnerable to COVID–19 [14].

Cell entry is facilitated by binding of SARS-CoV–2 to ACE2 receptors and on spike protein priming by host cell proteases. Transmembrane serine protease 2 (TMPRSS2) is a cellular protease that had previously been shown to be a proteolytic activator of both SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV), responsible for cleaving the spike protein of each into a conformation that facilitates cell membrane fusion and increases infectivity [19] [20]. TMPRSS2 also cleaves the ACE2 binding site itself, which further increases coronavirus infectivity [21]. It has since been shown to similarly prime SARS-CoV–2 for infection [22]. As in SARS-CoV [20], a TMPRSS2 inhibitor has been shown to block SARS-CoV–2 from cellular entry, and may serve as a treatment for COVID–19 [10]. This suggests that key polymorphisms in TMPRSS2 protein structure might also yield differential susceptibility to infection with SARS-CoV–2. This gene region shows pleiotropy in primates, and has previously shown signs of positive selection across the order while investigating its putative role in influencing the properties of seminal fluid [23]. Androgens regulate expression of TMPRSS2, and its functionality is autocatalyzed by cleavage of the protein at the R255 residue, yielding the functional form [24]. Previous research has identified the proteolytic activation function of TMPRSS2 to be dependent on the coordination of a catalytic triad comprised of the H296, D345, and S441 residues [20].

Savanna monkeys (Chlorocebus spp.) are a catarrhine primate that is ubiquitous throughout sub-Saharan African savanna environments [25]. They are often characterized as problem animals within this range due to their easy co-habitation with human settlements ranging from rural farmland to densely packed urban areas [26] [27] [28]. Green monkeys (Ch. sabaeus) are also endemic to the Caribbean islands of St. Kitts, Nevis, and Barbados, having been brought to the islands over 300 years ago [29], where they are densely populated and have frequent contact with humans [27]. The established Caribbean populations have since spread to the islands of St. Martin/St. Maarten and Tortola as escaped pets [30]. There is also a small feral population of green monkeys established near Ft. Lauderdale, Florida, in the United States [31].

Savanna monkeys have been shown to be more susceptible than rhesus macaques to respiratory infection and viral replication when exposed to coronavirus SARS-CoV, but when infected they show no evident external symptoms of disease [32]. Green monkeys (Ch. sabaeus) exposed to SARS-CoV (Urbani strain) via intratracheal swab have high titers of the virus in tracheal and lung tissue accompanied by pneumonitis on days 2 through 4 of infection, followed by eventual resolution [32]. Given the increased severity of response to SARS-CoV in savanna monkeys compared to rhesus, recent evidence that rhesus macaques show COVID–19-like symptoms upon SARS-CoV–2 infection [16] [17] [18] strongly suggests that savannah monkeys will also respond to infection with COVID–19-like disease. This not only raises the potential of the savanna monkey as a potential model for developing vaccines and therapeutics against COVID–19, but also for anthroponotic and zoonotic infection where savanna monkeys and humans meet.

Contact between savanna monkeys and humans with transmissible COVID–19 infections seem an inevitability, even with social distancing and stay-at-home orders in place in many African and Caribbean nations showing high infection rates. This is primarily due to their frequent raiding of crops and interactions with human refuse and other potential fomites in the wild [26] [27] [33]. Use of human resources as ‘hotspots’ of interest could potentially exacerbate animal disease propagation in human-frequented areas [34]. Urban areas present a special risk, as urban savanna monkeys are more habituated to human proximity and less likely to flee when humans are near [35]. Indeed, several zoonoses and anthroponoses have already been documented to pass between humans and wild savanna monkeys throughout their range (e.g.,simian foamy virus [36], Staphylococcus aureus [37], gastrointestinal parasites [38]). The abundance of sanctuary and rehabilitation centers that focus on caring for those savanna monkeys that have been adversely affected by human impacts serves as another key mediator of potential transmission [39] [40] [41].

Given the high zoonotic/anthroponotic transmission potential between savanna monkeys and humans, we used recently published genome sequence data from 163 wild savanna monkeys [42] and 82 captive green monkeys from the Vervet Research Colony (VRC) at Wake Forest School of Medicine [43] to document variation in the ACE2 and TMPRSS2 gene regions in this widespread genus. Our goals were 1) to assess the likelihood of anthroponotic infection of wild savanna monkeys upon exposure to infected humans, in order to evaluate the risk of savanna monkey populations becoming 2a) as highly impacted as human populations by potential epidemics, 2b) vectors of infection to humans via zoonotic exposure and 2c) potential reservoirs of infection with potential for the post-pandemic initiation of novel COVID–19 outbreaks among humans; 3) identify natural genetic models for functional studies of the host proteins critical to SARS-CoV–2 entry in the context of transmission of the virus and development of treatments for COVID–19; and 4) assess if captive savanna monkey colonies might be appropriate for modeling COVID–19 disease infection and therapeutics studies, as with SARS-CoV.

There is substantial variation in the ACE2 sequence of wild savanna monkeys. None of the detected variants are predicted to cause frameshifts or premature stops in the transcription process. Exonic variants include 6 missense variants, 14 synonymous variants, and two potential canonical splice region variants (Table 1), which are the only examples of putative protein-truncating variants in the vervet ACE2 discovered so far. Other detected variants include 43 in the 3’ UTR region, and 680 intronic variants including 3 potential splice region variants. Missense variants in savanna monkeys appear to occur only in the Ethiopian Ch. aethiops and West African/Caribbean Ch. sabaeus populations sampled. Only one of the missense-mediated changes in amino acid residues, also in Ch. sabaeus (X:14,077,550; Asp30Gly), involves a conversion at one of the critical binding sites for SARS-CoV–2 identified by Melin et al. [14]. This substitution is predicted to destabilize the binding of SARS-CoV–2 to ACE2 receptors (ΔΔG_bind = 0.958 kcal/mol), and of all the point mutations analyzed in a recent survey of primates [14], Asp30Gly reduces the predicted binding affinity more than all but the Tyr21His variant, which is found in most strepsirrhines (lemurs and lorisoids) and all platyrrhines (monkeys from the Americas) considered.

Savanna monkeys in the extended pedigree of Caribbean-origin Ch. sabaeus in the VRC showed marked variation in ACE2. A total of eight ACE2 variants are segregating in the pedigree: 2 splice acceptor site variants, 4 missense variants, and 2 synonymous variants. The two splice acceptor site variants are located at X:14,065,076 (at the end of a polyT tract) and were not observed in the wild monkeys. Both variants are common in the VRC (AAF_VRC = 0.08 and AAF_VRC = 0.34), and both are expected to be protein truncating variants with potential deleterious effects on the protein. These novel variants might have emerged in the VRC as de novo mutations, or may be present in wild populations at very low frequencies not captured with our moderate sample sizes. That the VRC pedigree is founded from a relatively small number of individuals, themselves originating from the founder population on St. Kitts, followed by the expansion of the pedigree might have led to an increased frequency of such rare variants. All four missense variants observed in wild Ch. sabaeus (X:14,035,353, X:14,035,374, X:14,045,023, and X:14,077,550causingAsp30Gly) are also present the VRC.

Table 1: Potential functional variants in ACE2 gene region sequence among wild savanna monkeys.

Position	Variant	Consequence	AAF	AA	Pos_AA	Notes
X:14035311	T/C	Missense	0.09	I/V	753	Alt. allele prevalent in Ethiopia (AAF=0.94); absent in all others
X:14035353*	C/T	Missense	0.33	V/I	739	Alt. allele only in Ch. sabaeus (AAF=0.75-1; AAF_VRC=0.84).
X:14035354	G/A	Synonymous	0.09	P	738	Alt. allele prevalent in Ethiopia (AAF=0.94); absent in all others
X:14035357	G/A	Synonymous	0.02	S	737	Alt. allele only found in Ethiopia (AAF=0.16)
X:14035374*	T/C	Missense	0.97	I/V	732	Ref. allele only found in St. Kitts (AAF=0.75; AAF_VRC=0.85).
X:14041898	T/C	Synonymous	0.13	G	629	Alt. allele only found in southern Africa (AAF=0.36-1.00)
X:14043283	G/C	Synonymous	0.04	L	585	Alt. allele predominantly found in Central African Republic (AAF=0.36)
X:14043289*	T/G	Synonymous	0.06	P	583	AAF_VRC =0.01
X:14043304	G/A	Synonymous	0.03	N	578	Alt. allele only found in South Africa (AAF=0.10)
X:14043773	T/A	Intronic SRV	0.45	-	-	Fixation of minor allele in Ch. sabaeus
X:14043797	T/A	Missense	0.02	E/V	549	Alt. allele only found in Ethiopia (AAF=0.19)
X:14045023*	C/G	Missense	0.62	E/D	483	Ref. allele near fixation in Ch. sabeaus (AAF=0.00-0.50; AAF_VRC =0.01).
X:14049928	G/A	Synonymous	0.01	I	358
X:14052949	A/G	Synonymous	0.18	F	315	Alt. allele only found in southern Africa (AAF=0.22-0.52)
X:14061092	C/T	Intronic SRV	0.02	-	-	Alt. allele only present in The Gambia (AAF=0.14)
X:14063390	A/G	Synonymous	0.02	L	266	Alt. allele only found in Ethiopia (AAF=0.22)
X:14064963	C/T	Synonymous; SRV	0.02	E	232	Alt. allele only present in The Gambia (AAF=0.14)
X:14065002	G/A	Synonymous	0.01	R	219	Alt. allele only present in East Africa (AAF=0.13-0.25)
X:10465076	C/CT CT/C	Splice Acceptor	0.08 0.34	-	-	Only present in the VRC, expected to be deleterious.
X:14067819	T/C	Intronic SRV	0.04	-	-	Alt. allele only found in southern Africa (AAF=0.03-0.13)
X:14077504	A/G	Synonymous	0.02	L	45	Alt. allele only present in Ethiopia (AAF=0.25)
X:14077524*	A/G	Synonymous	0.61	L	39	Fixation of reference allele in Ch. sabaeus (AAF_VRC=0.01)
X:14077531	G/A	Synonymous	0.01	A	36	Alt. allele only present in Zambia (AAF=0.06)
X:14077550*	T/C	Missense	0.13	D/G	30	*Alt. allele only found in Ch. sabaeus* (AAF=0.01-1.00; AAF_VRC=0.36)**

Emboldened, shaded text indicates coding regions or residues critical to SARS-CoV-2 binding. AAF = alternative allele frequency for the full sample. AA = change in amino acid residue predicted to accompany sequence variation. Pos_AA = amino acid position in the protein. SRV = splice region variant. Asterisks (*) = variant also present in the Vervet Research Colony at Wake Forest School of Medicine. AAF_VRC = alternative allele frequency in the VRC.

We also examined the genetic variation in TMPRSS2 both in the wild populations and VRC (Supplementary Table 1). A total of 15 missense and 21 synonymous variants were observed in the wild populations. In the VRC, one missense, one splice region and 9 synonymous variants were identified, including several variants shared with the East, Central, and South African populations, and all Ch. sabaeus variants (all of which are synonymous) observed in the wild populations, except one rare variant found in the Gambian population. An additional missense variant and splice variant were observed only in the VRC. Although there is substantially more variation in TMPRSS2 compared to ACE2 among savanna monkeys, all retain the functional residues for the cleavage site (R255) and for the catalytic triad (H296, D345, and S441).

The results of this first investigation of SARS-CoV–2 susceptibility in savanna monkeys in both the wild populations and in the VRC pedigree suggest that there are population differences in ACE2 sequence variation among savanna monkeys, some of which are estimated to yield changes in amino acid residues or exert deleterious effects on the transcript and protein function. The majority of the variants observed in this study in the wild monkeys are unlikely to dramatically alter susceptibility to disease upon exposure to SARS-CoV–2. The lack of frameshift mutations is in agreement with the paucity of protein-truncating variants in the human ortholog observed in gnomAD population data, apparently reflecting the intolerance of this gene against this type of variation in humans [44]. One of these changes (X:14077550; Asp30Gly) occurs at a previously identified critical binding site, part of the SARS-CoV–2 spike receptor binding domain [14]. A similar variant at this site (Asp30Glu, E30) is seen in the Microcebus murinus reference genome, and is predicted to confer reduced binding affinity of SARS-CoV–2 to ACE2 receptors (ΔΔG_bind = 0.692 kcal/mol) [14]. The Asp30Gly mutation, found in West African savanna monkeys and their Caribbean descendant populations, is predicted to destabilize the bond between the virus and the receptor to a greater degree (ΔΔG_bind = 0.958 kcal/mol). This suggests that a sizable proportion of the populations of Ch. sabaeus in West Africa and the Caribbean, those possessing the G30 residue, may have reduced susceptibility to SARS-CoV–2 infection and developing COVID–19-like illness.

Given the lack of clear impact on known functional regions in TMPRSS2, it remains to be seen what effect, if any, the population-specific missense variants observed might have on SARS-CoV–2 susceptibility. That expression of TMPRSS2 is dependent on androgen (AR) expression [24], however, may explain why human males appear to show significantly higher mortality with COVID–19 [45]. We might expect to see a similar, sex-specific pattern of severity in non-human primates. That male savanna monkeys are often bolder than females [46] suggests that they may both be more likely to interact with humans and more vulnerable to severe outcomes from SARS-CoV–2 infection.

Overall, our analysis suggests that savanna monkeys are likely susceptible to SARS-CoV–2 infection, though some populations may possess additional resistance to the disease, and that bi-directional infections between humans and savanna monkeys are possible. This high probability for susceptibility suggests that, as in other infected non-human primates (e.g., rhesus macaques) [16] [17] [18], savanna monkeys will likely become ill upon contracting the disease, and, like macaques, may suffer similar symptoms while being infectious to humans. Although this means there is high potential for savanna monkeys to experience concurrent outbreaks of COVID–19-like illness with sympatric human populations, they are unlikely to serve as a long-term reservoir of the disease, and may be less susceptible in the case of re-infection [32]. Additionally, strict quarantines in response to epidemics in nations with high contact between savanna monkeys and humans, like those implemented in March and April of 2020 in South Africa [47], will likely prevent such bi-directional transmissions from occurring as they will also isolate humans from non-human primate populations. Although personal protective equipment is the norm in captive primate facilities and many rescue centers and sanctuaries, similar isolation measures should be considered among those populations when human infection is likely.

We acknowledge that identifying susceptibility or resistance to infection based on receptor usage has limitations, especially with regards to describing severity of disease which has multiple contributing factors, for example the role of inflammatory cytokine responses in immunopathogenesis [48]. However, receptor usage provides an indication of species that could potentially be susceptible [11] [14]. The close proximity of humans and wild non-human primates provides potential for cross-species transmission of pathogens, and for some endangered species this could have devastating effects [14]. Similarly, identifying if non-human primates have the potential to act as intermediates hosts, as for Ebola virus [49], would be important for understanding zoonotic transmission and use of animal models.

Clear steps that must be taken to better understand this potential threat include characterizing the clinical signs of SARS-CoV–2 infection in savanna monkeys, and assessing whether they can develop a COVID–19-like disease state as seen in rhesus macaques. Established captive research colonies, such as those at the Vervet Research Colony of the Wake Forest School of Medicine in the United States (Ch. sabaeus),, the St. Kitts Biomedical Research Foundation in St. Kitts & Nevis (Ch. sabaeus),, and the Medical Research Council in South Africa (Ch. pygerythrus pygerythrus),, would be ideal for such assessments [25] [50]. These colonies may also be critical to the establishment of vaccines and therapies for COVID–19 should the green monkeys serve as a better model of SARS-CoV–2 infection and symptoms than macaque models, as is the case with SARS-CoV [32].

To facilitate such assessments, we characterized ACE2 and TMPRSS2 variation in the Vervet Research Colony. All four missense variants present in the wild Ch. sabaeus populations, including the potential modifying variant Asp30Gly (AAF_VRC = 0.36) are segregating in the VRC. The monkeys in the VRC may be useful models for better understanding variable susceptibility to infection with SARS-CoV–2 and developing appropriate interventions in these nonhuman primate populations. Human data collected so far lacks common natural mutations conferring differing efficiencies in coronavirus S-protein binding [51]. From this perspective, a model species with natural genetic variation with potential modulatory effects on the ACE2 role in the SARS-CoV–2 infection is of particular interest. In the VRC, we observed two splice acceptor variants potentially resulting in the truncation of the ACE2 protein. While they may provide a potential natural model for studies of the pathomechanism of the SARS-CoV–2 infection, they should be treated with caution as these variants are located in the repetitive region and require further validation through genotyping animals of interest.

Efforts should also be made to sample widely from wild savanna monkeys–with a focus on those in close contact with humans–in order to assess whether the virus has infected wild primate populations, and what characteristics of the local environment are important to human-monkey transmission. Researchers affiliated with the International Vervet Research Consortium already have contacts and protocols in place to potentially conduct such research in St. Kitts & Nevis, South Africa, and The Gambia [25] [50]. Such assessments could best characterize whether anthroponotic/zoonotic transfer of SARS-CoV–2 is a genuine public health and conservation concern, and predict areas where risk is high for such transmission. This work could better inform local policy with respect to human wildlife interactions. These steps will be essential to assess the potential risk of COVID–19 outbreaks mediated by human-monkey contact or conflict with savanna monkeys in Africa and the Caribbean.

Such assessments are critical, and would be the first steps towards a One Health approach to future monitoring of SARS-CoV–2 [52] [53]. Savanna monkeys predominantly live in nations with generally moderate health infrastructure, which have rapidly prepared but are also predicted to be among the hardest-hit by the current pandemic [54] [55]. Previous cycles of zoonotic disease transmission have already been well-established for other fatal outbreaks in such settings, in non-human primates with significantly lower frequencies of human interaction (e.g., wild apes and Ebola [49]). A clear understanding of potential coincident outbreaks and zoonotic cycles, especially in highly human impacted non-human primate populations, may be important to the long-term management of SARS-CoV–2 in both human health and conservation contexts.

Whole genome sequence data (median 4.4X coverage, with at least one member per taxon sequenced to at least 10X coverage) was previously generated for 163 individual wild savanna monkeys as part of research conducted by the International Vervet Research Consortium (for details, see [42] [50]. Raw sequence data were aligned to the Vervet-AGM reference assembly (ChlSab1.1 [56]), and variant calls are archived and publicly available online via Ensembl (Release 99 [57]). These data include resequencing from populations across sub-Saharan Africa and the Caribbean, including: 16 grivets (Ch. aethiops) from Ethiopia; 14 vervets (Ch. pygerythrus hilgerti) from Kenya and Tanzania; 16 malbroucks (Ch. cynosuros) from Zambia; 51 vervets (Ch. pygerythrus pygerythrus) from South Africa and Botswana; 11 tantalus monkeys (Ch. tantalus) from the Central African Republic; and 63 African green monkeys (Ch. sabaeus),, 22 from The Gambia, 2 from Ghana, 34 from St. Kitts & Nevis, and 5 from Barbados [42]. Although the taxonomy of these populations is contentious and currently unresolved [25], for simplicity we have here chosen to follow the taxonomic designations of Groves [58]. Variant calls in in the ACE2 region in wild savanna monkeys (X:14,030,233–14,077,785) were inspected using the Variant Table function, and filtered by annotations produced by Ensembl genebuild [59].

Allele frequencies in the VRC monkeys are reported here for 82 monkeys with moderate to high sequencing coverage depth (10–30X). Variant discovery and genotype calling steps are described in detail by Ramensky et al. [43]. For variant annotation, we used Ensembl Variant Effect Predictor [60]. Sequenced monkeys in the VRC are highly related, so in this case MAF should be interpreted as alternative allele occurrence rather than frequency in the strict statistical sense, which implies independent observations.

Finally, we used SSIPe [61] to assess the change in protein binding affinity caused by the discovered Asp30Gly mutation at the critical residue identified via alanine scanning mutagenesis analyses performed by Melin et al. [14]. The change in free energy due to the mutation was inferred by comparing against a reference interaction of human ACE2 receptor bound to the SARS-CoV–2 spike receptor binding domain (PBD 6M0J).

Acknowledgements

For the generation of this data, the authors gratefully acknowledge assistance in the field by Jennifer Danzy Cramer, Oliver “Pess” Morton, Yoon Jung, Stephanie Groman, David Jentsch, Eugene Redmond and the St. Kitts Biomedical Research Foundation in St. Kitts; veterinarians Lizanne Meiring, Murray Stokoe and numerous field assistants in South Africa including especially Gabriel Coetzer, Magali Jacquier, Helene DeNys, JD Pampush, Elzet Answegen, Tegan Gaetano, Dewald DuPlessis, Micah Beller, David Beller and Pess Morton; the University of Limpopo and the University of the Free State; and Mamkumba Sanneh, Katie Camfield, Michel Dione, Martin Antonio and the Medical Research Council in the Gambia. Collections of these data were conducted under permits and permissions issued by the following: Ethiopian Wildlife Conservation Authority; Gambia Department of Parks and Wildlife Management; Ministry of Tourism and Wildlife, Kenya; in South Africa the Department of Economic Development and Environmental Affairs, Eastern Cape; Department of Tourism, Environmental and Economic Affairs, Free State Province; Ezemvelo KZN Wildlife; Department of Economic Development, Environment and Tourism, Limpopo; the Department of Agriculture, Conservation and Environment, Mpumalanga and the Department of Environment and Nature Conservation, Northern Cape, as well as the South African National Department of Environmental Affairs.

Data Availability

All wild sequence data is publicly available in the European Variation Archive (EVA: PRJEB22988), and on Ensembl. Sequence data for the VRC are available in the NCBI Sequence Read Archive under BioProject Accession number PRJNA240242 and in the European Nucleotide Archive under the same accession number. Individual genotypes are available via the vervet genome browser at https://coppolalab.ucla.edu/vgb/home.

Author Contributions

CAS conceptualized the research and wrote the original manuscript, with additional writing by CMB, AJJ, VR, FB, and editorial input from HS, JPG, TRT. Analyses conducted by CAS, AJJ, VR, and protein binding analysis by CMB. HS, AJJ, NBF, MJJ are responsible for data curation. These data were generated with funding acquired by AJJ, MJJ, NBF, TRT.

Competing interests: The authors declare no competing interests.

Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020. doi: 1016/S1473-3099(20)30120-1

Xu Z, Shi L, Wang Y, Zhang J, Huang L, Zhang C, et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir Med. 2020;8(4): P420-422.

Zhang T, Wu Q, Zhang Z. Probably pangolin origin of SARS-CoV-2 associated with the COVID-19 outbreak. Curr Biol. 2020;30(7): 1346-1351.e2

Shi J, Wen Z, Zhong G, Yang H, Wang C, Liu R, et al. (2020). Susceptibility of ferrets, cats, dogs, and different domestic animals to SARS-coronavirus-2. bioRxiv [Preprint]. 2020. [cited 2020 April 17]. doi: 10.1101/2020.03.30.015347

American Veterinary Medical Association. COVID-19: FAQs for Pet Owners. Available from: https://www.avma.org/resources-tools/animal-health-and-welfare/covid-19/covid-19-faqs-pet-owners [Accessed 20 April 2020].

Centers for Disease Control and Prevention. If You Have Animals. Available from: https://www.cdc.gov/coronavirus/2019-ncov/daily-life-coping/animals.html [Accessed 20 April 2020].

Jia HP, Look DC, Hickey M, Shi L, Pewe L, Netland J, et al. Infection of human airway epithelia by SARS coronavirus is associated with ACE2 expression and localization. In: The Nidoviruses. Boston, MA: Springer; 2006. pp. 479-484.

Yan R, Zhang Y, Li Y, Xia L, Guo Y, Zhou Q. Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2. Science. 2020;367(6485): 1444-1448.

Zhang H, Penninger JM, Li Y, Zhong N, Slutsky AS. Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: molecular mechanisms and potential therapeutic target. Intensive Care Med. 2020; 1-5.

Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell. 2020;181(2): 271-280.e8. doi: 1016/j.cell.2020.02.052

Liu Z, Xiao X, Wei X, Li J, Yang J, Tan H, et al. Composition and divergence of coronavirus spike proteins and host ACE2 receptors predict potential intermediate hosts of SARS-CoV-2. J Med Virol. 2020;92(6): 595-601.

Li F, Li W, Farzan M, Harrison SC. Structure of SARS coronavirus spike receptor binding domain complexed with receptor. Science. 2005;309(5742): 1864-1868.

Li R, Qiao S, Zhang G. Analysis of angiotensin-converting enzyme 2 (ACE2) from different species sheds some light on cross-species receptor usage of a novel coronavirus 2019-nCoV. J Infect. 2020;80(4): 469-496.

Melin AD, Janiak MC, Marrone F, Arora PS, Higham JP. Comparative ACE2 variation and primate COVID-19 risk. bioRxiv [Preprint]. 2020. [cited 2020 April 24]. doi: 10.1101/2020.04.09.034967

Wu, C., Zheng, M., Yang, Y., Li, M., Liu, Y., Zhang, P. et al. In silico analysis of intermediate hosts and susceptible animals of SARS-CoV-2. ChemRxiv [Preprint]. 2020. [cited 2020 April 20]. doi: 10.26434/chemrxiv.12057996.v1

Shan C, Yao YF, Yang XL, Zhou YW, Wu J, et al. Infection with novel coronavirus (SARS-CoV-2) causes pneumonia in the rhesus macaques. ResearchSquare [Preprint]. 2020. [cited 2020 April 17]. doi: 21203/rs.2.25200/v1

Bao L, Deng W, Gao H, Xiao C, Liu J, Xue J, et al. Reinfection could not occur in SARS-CoV-2 infected rhesus macaques. bioRxiv [Preprint]. 2020 [cited 2020 April 17]. doi: 10.1101/2020.03.13.990226

Yu P, Qi F, Xu Y, Li F, Liu P, Liu J, et al. Age-related rhesus macaque models of COVID-19. Anim Model Exp Med. 2020;3(1): 93-97.

Matsuyama S, Nagata N, Shirato K, Kawase M, Takeda M, Taguchi F. Efficient activation of the Severe Acute Respiratory Syndrome coronavirus spike protein by the transmembrane protease TMPRSS2. J Virol 2010;84: 12658-12664.

Simmons G, Zmora P, Gierer S, Heurich A, Pöhlmann. Proteolytic activation of the SARS-coronavirus spike protein: Cutting enzymes at the cutting edge of antiviral research. Antiviral Res. 2013;100(3): 605-614.

Shulla A, Heald-Sargent T, Subramanya G, Zhao J, Perlman S, Gallagher T. A transmembrane serine protease is linked to the severe acute respiratory syndrome coronavirus receptor and activates virus entry. J Virol. 2011;85: 873-882.

Matsuyama S, Nao N, Shirato K, Kawase M, Saito S, Takayama I, et al. Enhanced isolation of SARS-CoV-2 by TMPRSS2-expressing cells. PNAS. 2020:117(13);7001-7003.

Clark N, Swanson WJ. Pervasive adaptive evolution in primate seminal proteins. PLoS Genet. 2005;1(3): e35.

Afar DE, Vivanco I, Hubert RS, Kuo J, Chen E, Saffran DC, et al. Catalytic cleavage of the androgen-regulated TMPRSS2 protease results in its secretion by prostate and prostate cancer epithelia. Cancer Res. 2001;61(4): 1686-1692.

Turner TR, Schmitt CA, Cramer JD. Savanna Monkeys: The Genus Chlorocebus. Cambridge, MA: Cambridge University Press; 2019.

Saj TL, Sicotte P, Paterson JD. The conflict between vervet monkeys and farmers at the forest edge in Entebbe, Uganda. Afr J Ecol. 2001;39(2): 195-199.

Dore KM. Identity construction and symbolic association in farmer-vervet monkey (Chlorocebus aethiops sabaeus) interconnections in St. Kitts. Folia Primatol. 2018;89: 63-80.

Turner TR, Schmitt CA, Cramer JD, Lorenz J, Grobler JP, Jolly CJ, Freimer NB. Morphological variation in the genus Chlorocebus: Ecogeographic and anthropogenically mediated variation in body mass, postcranial morphology, and growth. Am J Phys Anthropol. 2018;166: 682-707.

Denham WW. History of green monkeys in the West Indies. I. Migration from Africa. J Barbados Museum Historical Soc. 1982;36: 211-228.

Dore KM. Vervets in the Caribbean. International Encyclopedia of Primatology. 2017: 1-3. doi: 1002/9781119179313.wbprim0358

Anderson CJ, Hostetler ME, Johnson SA. History and status of introduced non-human primate populations in Florida. Southeastern Nat. 2017;16(1): 19-36.

McAuliffe J, Vogel L, Roberts A, Fahle G, Fischer S, Shieh WJ, et al. Replication of SARS coronavirus administered into the respiratory tract of African Green, rhesus and cynomolgus monkeys. Virology. 2004;330(1): 8-15.

Schmitt CA, Rich AM, Parke SAR, Blaszczyk MB, Cramer JD, Freimer NB, et al. Anthropogenic food enhancement alters the timing of maturational landmarks among wild savanna monkeys (Chlorocebus pygerythrus). bioRxiv [Preprint]. 2020. [cited 2020 April 17]. doi: 10.1101/2020.04.18.048314

Benavides J, Walsh PD, Meyers LA, Raymond M, Caillaud D. Transmission of infectious diseases en route to habitat hotspots. PLoS One. 2012;7(2): e31290.

Mikula P, Saffa G, Nelson E, Tryjanowski P. Risk perception of vervet monkeys Chlorocebus pygerythrusto humans in urban and rural environments. Behav Processes. 2018;147: 21-27.

Wolfe ND, Switzer WM, Carr JK, Bhullar VB, Shanmugam V, Tamoufe U, et al. Naturally acquired simian retrovirus infections in central African hunters. Lancet. 2004;363(9413): 932-937.

Senghore M, Bayliss S, Kwambana-Adams B, Foster-Nyarko E, Manneh J, Dione M, et al. Whole-genome sequencing reveals transmission of Staphylococcus aureusfrom humans to green monkeys in The Gambia. Appl Environ Microbiol. 2016;82(19): 5910-5917.

Gaetano TJ, Danzy J, Mtshali MS, Theron N, Schmitt CA, Grobler JP, et al. Mapping correlates of parasitism in wild South African vervet monkeys (Chlorocebus aethiops). S African J Wildlife Res 2014;44(1): 56-70.
Healy A, Nijman V. Pets and pests: vervet monkey intake at a specialist South African rehabilitation centre. Anim Welf. 2014;23: 353-360.

Wimberger K, Downs CT, Boyes RS. A survey of wildlife rehabilitation in South Africa: is there a need for improved management? Anim Welf. 2010;19: 481-499.

Grobler P, Jacquier M, deNys H, Blair M, Whitten PL, Turner TR. Primate sanctuaries, taxonomy and survival: A case study from South Africa. Ecol Environ Anthropol. 2006;2: 12–16.

Svardal H, Jasinska AJ, Apetrei C, Coppola G, Huang Y, Schmitt CA, et al. Ancient hybridization and strong adaptation to viruses across African vervet monkey populations. Nat Genet.2017;49: 1705–1713.

Ramensky V, Jasinska, AJ, Deverasetty, S, Svardal, H, Zelaya, I, Jorgensen, et al. The burden of deleterious variants in a non-human primate biomedical model. bioRxiv [Preprint]. 2019. [cited 2020 April 22]. doi: 10.1101/784132

Karczewski KJ, Francoli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. bioRxiv [Preprint]. 2019. [cited 2020 April 25]. doi: 10.1101/531210

Nasiri MJ, Haddadi S, Tahvildari A, Farsi Y, Arbabi M, Hasanzadeh S et al. COVID-19 clinical characteristics, and sex-specific risk of mortality: Systematic Review and Meta-analysis. medRxiv [Preprint]. 2020. [Accessed 21 April 2020]. doi: 10.1101/2020.03.24.20042903

Blaszczyk MB. Boldness towards novel objects predicts predator inspection in wild vervet monkeys. Anim Behav. 2017;123:91-100.

South African Government. Regulations and Guidelines – Coronavirus Covid-19. Available from: https://www.gov.za/coronavirus/guidelines [Accessed 24 April 2020].

Mehta P, McAuley DF, Brown M, Sanchez E, Tattersall RS, Manson JJ, et al. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet. 2020;395(10229): 1033-1034.

Leendertz SA, Wich S, Ancrenaz M, Bergl RA, Gonder MK, Humle T, Leendertz FH. Ebola in great apes – current knowledge, possibilities for vaccination, and implications for conservation and human health. Mammal Rev. 2017;47(2): 98-111.

Jasinska AJ, Schmitt CA, Service SK, Cantor RM, Dewar K, Jentsch DJ, et al. Systems biology of the vervet monkey. ILAR / Natl Res Council. 2013;54: 122-143.

Cao Y, Li L, Feng Z, Wan S, Huang P, Sun X, et al. Comparative genetic analysis of the novel coronavirus (2019-nCoV/SARS-CoV-2) receptor ACE2 in different populations. Cell Discov. 2020;6: 11. doi: 10.1038/s41421-020-0147-1

Kelly TR, Machalaba C, Karesh WB, Cook PZ, Gilardi K, Nziza J, et al. Implementing One Health approaches to confront emerging and re-emerging zoonotic disease threats: lessons from PREDICT. One Health Outlook. 2020;2(1). doi: 1186/s42522-019-0007-9

Gillespie TR, Nunn CL, Leendertz FH. Integrative approaches to the study of primate infectious disease: implications for biodiversity conservation and global health. Yrbk Phys Anthropol. 2008;51: 53-69.

Hopman J, Allegranzi B, Mehtar S. Viewpoint: Managing COVID-19 in low- and middle-income countries. JAMA. 2020. doi: 1001/jama.2020.4169

Gilbert M, Pullano G, Pinotti F, Valdano E, Poletto C, Böelle PY, et al. Preparedness and vulnerability of African countries against importations of COVID-19: a modelling study. Lancet. 2020;395(10227): 871-877.

Warren WC, Jasinska AJ, Garcia-Perez R, Svardal H, Tomlinson C, Rocchi M, et al. The genome of the vervet (Chlorocebus aethiops sabaeus). Genome Res.2015;25: 1921-1933.

Cunningham F, Achuthan P, Akanni W, Allen J, Amode MR, Armean IM, et al. Ensembl 2019. Nucleic Acids Res. 2019;47(D1): D745-D751.

Groves CP. Primate Taxonomy. Washington, DC: Smithsonian Institution Press; 2001.

Curwen V, Eyras E, Andrews TD, Clarke L, Mongin E, Searle SMJ, Clamp M. The Ensembl automatic gene annotation system. Genome Res. 2004;14(5): 942-950.

McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, et al.The Ensembl Variant Effect Predictor. Genome Biol. 2016;17(1): 122.

Huang X, Zheng W, Pearce R, Zhang Y. SSIPe: accurately estimating protein-protein binding affinity change upon mutations using evolutionary profiles in combination with an optimized physical energy function. Bioinformatics. 2019. doi: 10.1093/bioinformatics/btz926.36

SupplementaryTable1.docx

ACE2 and TMPRSS2 variation in savanna monkeys (Chlorocebus spp.): Potential risk for zoonotic/anthroponotic transmission of SARS-CoV-2 and a potential model for functional studies

Status:

Journal Publication

Version 1

Abstract

Introduction

Results

Discussion

Materials And Methods

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1