Introduction
Humankind is currently facing the effects of an unparalleled pandemic. COVID-19, a disease caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARSCoV-2), has gravely impacted health and socio-economic systems in practically every country in the world (WHO, 2020). While the critical effects of this pandemic were unexpected, the SARS-Cov-2 virus is just one of a long list of infectious agents, primarily viral, that are responsible for high impact human diseases, including acquired immune deficiency syndrome (AIDS), severe acute respiratory syndrome (SARS-CoV-1), ebola, and influenza, among others (Taylor et al., 2001; Christou, 2011; Plowright et al., 2017). These are all examples of recently emerged zoonotic diseases that have shaken healthcare systems at a global scale, which lead to reexamine how the entire health field should be approached, and present the need to include disciplines beyond medicine and public health, such as ecology, evolutionary biology, and anthropology, to better control and understand disease. Given that epidemics and pandemics are the result of inter-species interactions (usually among wildlife, domestic animal species, and human populations) which are influenced by ecosystem transformation and multiple sociocultural factors, there is an urgent need for transdisciplinary collaboration that allows for a better understanding and a stronger approach to each problem (Lewis, 2005; Wilcox and Colwell, 2005).
To understand the complexity of the current coronavirus emergence, it is necessary to contextualize the factors that are consistently associated with these types of events (Johnson et al., 2015; Webster et al., 2016). Zoonosis emergence includes anthropic transformation of natural ecosystems, wildlife community modification, and the alteration of the biogeographic and evolutionary history of the microorganisms that have been established in their natural hosts for millions of years (Taylor et al., 2001; Hoberg and Brooks, 2015; White and Razgour, 2020). Thus, disease ecology provides precise information to contextualize the causes of emergence and to appraise the role of interactions among pathogens, domestic animals, wildlife, and humans, in an integrated manner.
In particular, the emergence of diseases caused by coronaviruses (CoVs) resulted from increased direct and indirect interactions between humans and other animal species. Human activities have modified practically all ecosystems on the planet, which has led to changes in their structure and function, as well as changes in geographic distribution of wild species, including disease vectors and hosts, and modifications in the dynamics of infectious agents (Cunningham et al., 2017; Daszak et al., 2001). Ecosystem alterations modify patterns of occurrence and distribution of emerging infectious diseases (EIDs). This phenomenon has been exemplified by the recent emergence of all three of the novel coronaviruses (SARSCoV-1, Middle East Respiratory Syndrome (MERS), and COVID-19).
The One Health/EcoHealth perspective integrates public health, animal health, and ecosystem health, simultaneously emphasizing the necessity for international connectivity and interinstitutional collaboration, for a unifying management approach (Barrett et al., 2011; Errecaborde et al., 2019). The SARS-CoV-1 outbreak, that originated in China and that lasted eight months, infected 8,098 people and killed 774 (9.5%) worldwide (WHO, 2003), becoming the first pandemic of the 21st century. This, pandemic revealed interactions between different bat species, civets, raccoon dogs, and humans, that led to changes in infection dynamics (Luk et al., 2019). Later, in 2012, during the emergence of MERS, a disease associated with another coronavirus variant, 2,494 infected people and 858 deaths in 27 countries were recorded (WHO, 2019). These new CoVs, as the recently described SARS-CoV-2, were noteworthy because of their high pathogenicity in humans. In contrast, other previously reported CoVs infections in human populations, such as HCoV-229E, HCoV-OC43, HCoV-NL63 were considered to present low pathogenicity (Walsh et al., 2013), emphasizing particular variations.
Disease ecology is a discipline that studies the ecological and evolutionary factors that determine the occurrence, dynamics, and distribution of infectious diseases (Johnson et al., 2015). Over the past 25 years, disease ecology has allowed for the identification of the factors that determine the occurrence of emerging diseases in space and time (Cook et al., 2004; Gibbs and Gibbs, 2013; Kelly et al., 2017). This transdisciplinary approach contextualizes information related to the causes of disease emergence and appraises the role of interactions among pathogens, domestic animals, wildlife, and humans, in an integrated manner. This review aims to present a disease ecology perspective by including complementary and transdisciplinary insight and methodologies while illustrating some of their input in results related to coronavirus investigations conducted in Mexico. We draw from various disciplines, including ecology, molecular epidemiology, biogeography and evolutionary biology, to address the emergence of coronaviruses, including SARSCoV-2 (Figure 1).
Coronavirus taxonomy and evolution
The CoVs infect a wide diversity of animals, including humans. Their clinical presentation varies, affecting respiratory, enteric, renal, hepatic, and neurological tissues with fluctuating degrees of severity (Masters, 2006; Quinn, et al. 2011; Perlman, 2013). CoVs belong to the order Nidovirales, which contain the Coronaviridae and the Torovirinae families. CoVs are genetically diverse but share a large genomic size with a similar genetic organization, that encompass comparable replication strategies. However, they differ in morphology and genome length. These viruses are globally distributed, are highly infectious, and are extremely difficult to control due to their genetic diversity, short generation time, and high mutation rates. Bats and other wild species are known to have a large reservoir of previously unidentified CoVs (Anthony et al., 2017; Cui et al., 2019; Luk et al., 2019).
CoVs were initially classified into three distinct genera based on genotypic and serological characterization (CoV 1, 2, and 3; (Masters, 2006). Later, the Coronaviridae Study Group of the International Committee on Taxonomy of Viruses proposed three new genera-Alfacoronavirus (α-CoV), Betacoronavirus (ß-CoV) and Gammacoronavirus (ϒ-CoV)-to replace the traditional CoV 1, 2 and 3 groups (de Groot RJ et al., 2011). Recently, a new genus has been identified, the Deltacoronavirus (Δ-CoV), which includes the smallest genome viruses among all CoVs (25 to 26 kb); (de Groot RJ, 2011; Woo et al., 2012a).
The α-CoV and ß-CoV infect mammals, while the ϒ-CoV and Δ-CoV mainly infect birds, though some can also infect mammals. The α-CoV and ß-CoV generally cause respiratory diseases in humans and gastroenteritis in animals. Both groups of viruses are highly pathogenic, with SARS-CoV-1, MERS-CoV and SARS-CoV-2 causing severe respiratory syndrome in humans, while the human-endemic CoVs (HCoV-NL63, HCoV-229E, HCoV-OC43 and HKU1) mainly induce mild respiratory conditions in immunosuppressed individuals, being the origin of close to 15% of seasonal colds. Certain human-endemic CoVs can also cause severe infections in infants, young children, and the elderly (Masters, 2006; Perlman, 2013).
In general, CoVs found in humans originated in animals, with SARS-CoV-1, MERS, SARS-CoV-2, HCoV-NL63 and HCoV-229E coming from bats, while HCoVOC43 and HKU1 probably emerging from rodents (Anthony et al., 2017; Cui et al., 2019). Based on molecular clock analyses, it has been estimated that the most recent common ancestor of all CoVs originated approximately 8,100 years ago (YoE, years of evolution), and those of the α-CoV, ß-CoV, ϒ-CoV and Δ-CoV encompass approximately 2400 YoE, 3300 YoE, 2800 YoE, and 3000 YoE, respectively. It has been proposed that bats and birds are the ideal host for CoVs, with bats giving rise to the α-CoV and ß-CoVs, and birds originating the ϒ-CoV and Δ-CoV groups (Woo et al., 2012b). CoVs that affect humans are found within the α-CoV and ß-CoV genera (Luk et al., 2019). There are four distinct sub-genera in the ß-CoV genus: ßa, ßb, ßc, and ßd. In particular, the two SARS CoVs that affect humans are found in the subgenus ßb, while MERS-CoV is classified within the ßc subgenus. CoVs described in bats belong to the ßb, ßc and ßd subgenera, as well as to a specific clade of the α-CoV genus (Figure 2).
It has been shown both experimentally and in nature, that CoVs have a genetic divergence that favors host and niche switching and adaptation (Jackwood, 2006; Woo et al., 2009, 2006). In particular, CoVs exhibit continuous genetic recombination through the genomic mechanism of template switching (Woo et al., 2006). These changes generate point mutations at a similar speed to other RNA viruses, including type A influenza viruses (Brian and Baric, 2005). This tendency for recombination and high mutation rates increase likelihood for regular host switching between animals and humans and multi-host capacity, leading to potential epidemics and epizootics.
Coronavirus structure and molecular biology
CoVs are enveloped viruses with a positive-sense single-stranded RNA genome and a helical nucleocapsid. The CoVs genome size ranges from 26 to 32 kilobases, the largest among all RNA viruses. The shape of the virion is pleomorphic, varying in diameter from 80 to 100 nm. These viruses have characteristic projections on their surface (spikes), which are reminiscent of the solar corona when observed by electron microscope imaging, hence acquiring their name (Jackwood, 2006; Masters, 2006).
The genomes of all CoVs have a 5′ methylated cap, a 5′- UTR leader region, and a UTR region with a polyadenylated tail at the 3′ end. The main proteins encoded by the CoV genome from the 5’ end are: the viral polymerase, the polyprotein replicase/transcriptase (which is coded by Open Reading Frames, ORF 1a and 1b, that occupy the first two thirds of the genome, and are cut to form non-structural proteins), the hemagglutinin-esterase (HE, which is present in some viruses), the spike glycoprotein (S, located on the surface of the virus), a small envelope protein (E), a membrane integrated glycoprotein (M), and the nucleocapsid protein (N, which is closely associated with the viral RNA). Interspersed between these reading frames are additional reading frames coding for accessory proteins that determine a number and functions that are unique to each specific coronavirus (Jackwood, 2006; Masters, 2006) (Figure 3).
At a molecular level, CoVs employ various unusual strategies to complete a complex gene expression program. Their replication involves ribosomal frameshifting during translation of the genome, synthesis of multiple genomic and subgenomic RNA species, and the assembly of progeny virions via a pathway that is unique among enveloped RNA viruses (Brian and Baric, 2005; Masters, 2006). In infected cells, CoVs use a nested transcription strategy in which gene expression occurs through the synthesis of several 3’-coterminal subgenomic mRNA. This unique strategy was recognized as distinctive for the establishment of the order Nidovirales (from the Latin, nidus, nest) (de Vries et al., 1997; Jackwood, 2006; Masters, 2006).
The α-CoVs do not have the HE protein; the M protein is N-glycosylated, and the S glycoprotein is not cut. Most of the viruses from the ß-CoV do have the HE protein; the M protein is O-glycosylated, and the S glycoprotein is cut into two subunits. The ϒ-CoV lack HE protein; the M protein is N-glycosylated, and the S glycoprotein undergoes proteolytic cuts (Cui et al., 2019; Jackwood, 2006; Woo et al., 2012b). The Δ-CoVs have a single protease domain, similar to papain at gene nsp3; the accessory ORF NS6 gene (ORF) situated between genes M and N), and a variable number of additional accessory genes (up to four) (Woo et al., 2012b).
Factors that favor infectious disease emergence: the coronavirus epidemics
In addition to recognizing the particular features of coronaviruses, it is necessary to describe additional factors that favor viral emergence, as well as the associated ecological and anthropogenic processes, that are critical to understand disease development.
The zoonotic viral emergence process
The diseases caused by novel CoVs during the 21st century are representative of zoonotic transmission from wild animal species. While the source of SARS-CoV-2 has not been 100% confirmed, it is highly probable that like SARS-CoV-1 and MERS, the origin of SARS- CoV-2 is associated with bats (H. Zhou et al., 2020; P. Zhou et al., 2020) and probably pangolins (Lam et al., 2020; Xiao et al., 2020).
The ecological process of a zoonotic virus emergence is complex, since it involves intrinsic factors of both a particular virus and its host species, as well as extrinsic factors related to the environment in which the virus and hosts interact (Boots et al., 2009; Day and Proulx, 2004). To understand this process, it is necessary to acquire information related to viral evolution mechanisms, permanence of the virus in natural hosts, and transmission among taxonomically distinct species, like humans. It is also necessary to understand how the disease process arises and how, through repeated infections in the new species, the virus becomes capable to generate sufficient morbidity or mortality rates to be detected and characterized as a new health problem of local, regional, or global importance.
The ecological process of viral emergence must go through four fundamental stages, the first two being a prerequisite for the emergence of a new disease. First, contact must occur between the natural reservoir and a new susceptible species. Second, there must be an interspecies transmission, in which the virus crosses a species barrier and completes its cycle of infection within a new host. In some diseases, such as the Hantavirus pulmonary syndrome for example (Ermonval et al., 2016), the new host may become infected but cannot transmit infection, becoming an accidental host, and limiting infection to this single individual. Some infections may require the participation of several species as intermediate hosts, as is the case of camels and alpacas in MERS, or of civets in SARS-CoV-1 (Gong and Bao, 2018). In other instances the reservoir can transmit the disease directly to the new species, without the need for an intermediate host. The role of intermediate hosts in the emergence of COVID-19 is not yet clear, though the Malayan Pangolin (Manis javanica) has been purported as a possibility (Lam et al., 2020; Liu et al., 2019; Xiao et al., 2020). It is also important to recognize that intermediate hosts can maintain diseases caused by CoVs, while being reservoirs that favor viral recombination, hence having the capacity to infect other species, such as humans (Mohd et al., 2016; Rabaan et al., 2017).
The third stage of the viral emergence process is sustained transmission among individuals of the new species. In the SARS-CoV-2 situation for instance, the individual or individuals that were in contact with the reservoirs and became infected, were also able to transmit the agent to their conspecifics. The fourth and final stage of disease emergence is genetic adaptation following transmission. The evolutionary adaptation of the virus and the new host (and vice versa) results in a new virus that is qualitatively and quantitatively distinct from the original virus (Agudelo-Romero et al., 2008; Cressler et al., 2016).
Ecosystem transformation
Anthropogenic ecosystem transformation includes land use change, expansion of the agricultural and livestock boundary, urbanization, changes to environmental and climate conditions, legal and illegal wildlife trafficking, as well as water, soil, and air pollution (Ellis et al., 2010). All of these factors modify the interactions between organisms within the ecosystem and broaden the contact interface between wildlife, domestic animals, and humans (Mackenstedt et al., 2015; Patz et al., 2004). Increasing the rates of inter-species contact also alters the behavioral and population ecology of the hosts, as well as that of their infectious agents, and the characteristics of the host-virus interaction (Daszak et al., 2000; Jones et al., 2008; Morens et al., 2004; Morse, 1995). Therefore, anthropogenic activities can have a profound impact on the assembly of host species communities (i.e., diversity, composition and abundance), which can subsequently impact pathogen prevalence (Johnson et al., 2015).
To study how diseases are maintained within ecosystems, disease ecology incorporates ecological and phylogenetic indexes, such as taxonomic and phylogenetic beta diversity metrics, that are used in community ecology to understand the interactions within multi-host and multi-pathogen systems (Behdenna et al., 2019; Woolhouse et al., 2001). These indexes are also used to measure changes in the viral composition within the ecosystem community, including CoVs, considering hosts at different spatial scales, over several environmental gradients (Anthony et al., 2013; Rico-Chávez et al., 2015; Nieto-Rabiela et al., 2019, 2018). A Mexican study described a strong relationship between bat diversity and richness of four viral families of medical importance (Hantavirus, Paramyxovirus, Pegivirus and Coronavirus) in three regions of the country. Results of this work suggest that generalist bat species, namely species that are not affected by habitat modifications, are related with more viral families than bat species affected by habitat modifications (Rico-Chávez et al., 2015). Thus habitat modification can bias the composition of the potential host communities, at different taxonomic and phylogenetic levels, toward species that tend to carry more viruses, generating zoonotic disease emergence risk areas, often near human populations.
Focusing on spatial scales has allowed for the identification of space and timeframe levels where infections operate in multi-host systems. Most of the processes that favor a viral taxonomic jump occur at different spatial scales and are contextspecific. Nonertheless, identifying the scale at which an infectious disease operates is fundamental for its prevention and modeling (Milholland et al., 2018). The landscape scale is where most of the interaction processes occur, such as inter-species contact or community composition modifications, and also where the human-wildlife-domestic animal interface takes place. This scale can therefore be considered as the scale of operation of emerging zoonotic diseases, yet it is seldomly evaluated in infection modeling and simulation (Suzán et al., 2015).
Nowadays it has become increasingly crucial to comprehensively understand how food production methods, land-use planning, and other human activities can affect spatial patterns and ecological processes in ways that impact infection emergence. Landscape ecology studies of interactions between flora and fauna, as well as social components, including spatial and temporal variation, are imperative to adequately assess disease emergence scenarios. Addressing biotic and abiotic interactions from an earth science perspective and characterizing the geometric conformation, functionality, and heterogeneity of the landscape becomes essential (Wiens et al., 1993; Urban, 2006). Landscape ecology is a recent discipline, whose development has been largely driven by technological advances in remote sensing. In the context of infectious disease, there is a growing number of studies that use the characterization of landscapes to relate landscape transformation to infection dynamics (Cobb et al., 2012; Gottdenker et al., 2014). The analysis of landscape heterogeneity and connectivity show how dynamic transmission scenarios can be and how infections can persist.
Anthropized environments are characterized by alternating environmental heterogeneity and homogeneity which often include combinations of grazing pasture, crop fields, fragments of natural and secondary vegetation of variable size and shape, in rural, semi-rural, suburban and urban human settlements. Each of these landscape variables has been shown to impact disease transmission. The shape, size, and connectivity of environmental patches affect infection dynamics. Empirical studies in Costa Rica (Daily et al., 2003), Venezuela (Utrera and Duno, 2007), the United States (Diffendorfer et al., 1995), Chile (Torres-Pérez et al., 2004), Panama (Suzán and Ceballos, 2005) and Paraguay (Yahnke, 2006), as well as studies at continental and global scales (Rubio et al., 2014; López-Pérez et al., 2019) have shown that habitat loss and landscape modification favor pathogen reservoir species that tolerate fragmented environments. For infectious agents associated with bats, different land use types and landscape variations, including linear fragmentation and edges, favor viral reservoir species (Rico-Chávez et al., 2015; Lanzagorta-Valencia et al., 2019) and a differential distribution of viral diversity, including that of CoVs that we have reported in southern Mexico (Anthony et al., 2013; Rico-Chávez et al., 2015).
The occurrence of emerging disease outbreaks with high human health impact, has been linked to proximity of people and animals. The emergence of MERS for example is attributed to contact between people and domesticated camels, which cohabitate with other species of domestic animals as well as bats that live in nearby refuges (Goldstein and Weiss, 2017; Smith and Wang, 2013). Similarly, appearance of , the Hendra virus in people in Australia in 1994 was associated with bat (genus Pteropus) transmission to horses, and then from horses to humans (Halpin et al., 2000). Also, , the Nipah virus in several Asian countries emerged from interactions between bats, pigs, palm plantations, and humans. This process was favored by ecosystem fragmentation and niche invasion (Kulkarni et al., 2013). The case of SARS-CoV-2 is similar. Phylogenetic analysis of viral variants suggest that bats are the primary source of this CoV (H. Zhou et al., 2020; P. Zhou et al., 2020). While the precise location and mechanisms of transmission are not yet known, it has been suggested that Malayan Pangolins (Manis javanica) acted as intermediate hosts, contracting the virus from bats, then transmitting the infection to humans (Lam et al., 2020). The interaction among these three species is likely linked to a confluence of conditions that favored more frequent encounters in the human-animal interface, ecosystem fragmentation and deterioration, illegal markets, niche invasion, poaching, and consumption of bushmeat without sanitary regulations (VanWormer et al., 2013; Hassell et al., 2017; Paige et al., 2017). Landscape ecology assesses landscape transformation and infection dynamics, including sociocultural consumption habits and other conducts, as well as aspects of utilization of natural resources. These scenarios, where biological and cultural processes occur, signify the ideal landscape scale for modeling multi-factor scenarios that resemble the current pandemic emergence events.
Wildlife consumption and legal and illegal trafficking
Every year, hundreds of thousands of plants, animals, as well as their derivatives are legally traded to satisfy food demand of different societies (Rosen and Smith, 2010). Meanwhile, the illicit wildlife trade is considered the second most lucrative illegal market after drug trafficking (Wyler and Sheikh, 2016). The unlawful nature of this trade prohibits a full knowledge of its magnitude, and hinders the understanding of its impact on ecosystems and health. Many pathogens have been introduced into new ecosystems and communities when host species are transported, legally or illegally, to sites out of their original distribution boundaries. Introduced pathogens can potentially lead to local species extinctions, including endemic and/ or endangered species, and can trigger ecosystem collapse. Pathogens that have been repeatedly mobilized by this mechanism include rabies, Bacillus anthracis, Mycobacterium tuberculosis, Echinococcus spp. and Leptospira spp., among others (Pavlin et al., 2009).
Systematic monitoring of species trafficking in markets in Guangzhou, China for instance have revealed constant trading of masked palm civets (Paguma larvata), ferret-badgers (Melogale spp.), wild boar (Sus scrofa), and several endangered species (Lau et al., 2010). After the emergence of SARS-CoV-1 in 2003, 838,500 wildlife animals were confiscated from markets in this city (British Broadcasting Corporation, 2003). It has been recognized that this type of trafficking favors interactions among species, such that wild-life trade markets serve as points of contact that facilitate taxonomic jumps of different infectious agents to new hosts. The origin of the Human Immunodeficiency Virus is also attributed to human consumption of non-human primates (Gao et al., 1999). Similarly, the recent Ebola outbreaks are associated with human consumption of great apes (Leroy et al., 2004), and the SARS-CoV-1 has been related to international trade of small carnivores and bats (May et al., 2004; Lau et al., 2005).
Consumption of bushmeat is known to increase food availability for low-income human communities in many countries around the world, particularly in Africa (Friant et al., 2020). Unfortunately, this practice is not regulated, and endangered species are often overexploited. Certainly, a report to to the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), informed of the trade of pangolins (mostly from Asia), mainly for their skin. The report estimates that from 1977 to 2014, approximately 50,000 pangolins were trafficked, most being of the Manis javanica species (Heinrich, 2017). This is the same species that has been associated with the transmission of SARS-CoV-2. Current data estimate that nearly 900,000 pangolins have been traded across Southeast Asia in the past two decades (WHO, 2020).
Wildlife and emerging coronaviruses
Coronaviruses in wildlife
The human epidemic outbreaks of SARS-CoV-1 (and MERS-CoV) in China (2002-2003) and Saudi Arabia (2012), (Drosten et al., 2003; Zaki et al., 2012), have led to an increased interest in this viral family and its relationship to wildlife on a global scale. At the beginning of the SARS-CoV-1 outbreak, the viral reservoir in wild animal species was unknown, and in initial studies, the palm civet (Paguma larvata), the raccoon dog (Nyctereutes procyonoides) and the Chinese ferret-badger (Melogale moschata) were proposed as the virus natural hosts (Song et al., 2005; Xu et al., 2004). Ensuing field studies and extensive scale phylogenetic analyses revealed that horseshoe bats (Rhinolophidae) acted as the natural reservoir of SARS-CoV-1, and that the civet, raccoon and ferret were intermediate hosts that facilitated amplification of the virus (Lau et al., 2005; Li et al., 2005). The first CoVs report in bats came from China in 2005, as part of a monitoring program established to determine viral circulation in wildlife (Poon et al., 2005). Since then, a growing number of wildlife studies, mainly in bats, have increased knowledge of CoVs diversity, allowing for new viral species discoveries in a wide range of mammal and avian hosts.
Bats appear to be particularly important reservoirs for CoVs and have been identified as natural hosts for several CoVs genoypes after the SARS epidemic (Drexler et al., 2014). Moreover, bat CoVs are more diverse than those identified in other wild mammal species; bats being natural hosts to the largest diversity of CoVs described to date (Drexler et al., 2014; Anthony et al., 2017). Consequently, search efforts have increasingly focused on bats, as they continue to be established as reservoirs for numerous viral agents. The abundant descriptions of new viruses in bats and other wild animals have drastically changed our perception of the relevance of reservoir species in emerging zoonoses (Ang et al., 2018; Bengis et al., 2004; Hayman, 2016). However, knowledge of the ecological and evolutionary factors that affect viral diversity is often overlooked. For example, there is a geographic bias in CoVs studies in bats, with a notable lack of information in high biodiversity regions including Africa, Asia, and Latin America (Drexler et al., 2014). Given that bats are central CoV reservoirs, an entire line of research in this mammal species has become critical to understand factors that affect viral prevalence in wildlife populations. Certainly, it has been suggested that host ecology could be determinant in defining the degree of susceptibility and viral prevalence, which could limit or favor viral transmission within or between different species (Luis et al., 2013).
Coronavirus evolutionary patterns
Before the SARS-CoV-1 epidemic in 2002, most of the scientific knowledge regarding CoVs was the result of animal health research, mainly in domestic animals, leaving ecological and evolutionary aspects practically unexplored. It has however become increasingly evident that phylogenetic and coevolutionary studies of emerging pathogens are central in identifying the origin of the virus, its current distribution and associated host species, as well as for providing the bases for understanding viral ecology and zoonotic potential in different risk areas (Davies and Pedersen, 2008; Lei and Olival, 2014; Chu et al., 2011).
After the SARS-CoV-1 epidemic, nine new species of CoVs were identified in hosts where they had not been previously described. Given that the interactions between infectious agents and species that act as their natural reservoirs occur over evolutionary time scales (millions of years), development of a large diversity of well-adapted pathogens can take place (Carroll et al., 2007). For CoVs in particular, most genotypes naturally infect a single host species or a small number of closely related species. High genotype diversity of CoVs is mainly due to their elevated mutation rate (0.80 - 2.38 × 103 mutations per site per year) and their evolutionary dynamics, including recombination and other viral strategies (Drexler et al., 2014; Durães-Carvalho et al., 2015). These characteristics could explain the proposed coevolutionary relationship between CoVs and different species of bats (Anthony et al., 2013; Cui et al., 2019).
Bat CoVs have a large genetic divergence and are the oldest identified CoVs considering the evolutionary scale. Interestingly, it has been shown that CoVs present stable viral loads in populations of bats, while an epidemic population growth can be seen in other host species (Vijaykrishna et al., 2007; Gloza-Rausch et al., 2008; Pfefferle et al., 2009). Given that there is a close bat-CoV relationship between specific species of both types of organisms, it has been important to consider how the diversity and distribution of bat species affects the richness and diversity of CoVs. In the Americas, for example, the relationship between similar CoV genotypes and closely related bat species has supported the hypothesis of coevolution (Ojeda, 2015). This association has been observed when comparing CoV genotypes present in populations from geographically distant areas (Gloza-Rausch et al., 2008; Tang et al., 2006). Nonetheless, reports also show that within the same refuge, two host species from different families of bats can harbor the same CoV genotype, contrary to the theory of host specificity (Osborne et al., 2011). The coevolutionary relationship found between bats and CoV populations in Mexico, indicates that hosts drive CoV selection and evolution, highlighting the importance of studying these evolutionary associations.
Among the different approaches pertinent for the study of the emergence of zoonotic diseases, focusing on evolutionary associations between host and viruses could thus be important for viral reservoir identification and related inter-species jumps. This methodology considers that emerging viruses circulating in new species share sufficient similarity with the viruses circulating in reservoir communities.
In Mexico, monitoring bat CoVs has revealed a coevolutionary relationship in which host species act as drivers for CoVs selection and evolution, such that the host geographic distribution pattern could determine the distribution and epidemiology of the viruses (Anthony et al., 2013). For example, a bat related CoV monophyletic group has been identified in Mexico City (located at the center of the country), which has different genome sequences than those identified for a group in the southeast of the country, suggesting that the difference in phylogenetic composition of CoV communities could reflect the variations in the composition of the host communities distributed in the regions of the Nearctic, Neotropic, and transition zone located in Mexico (Ojeda, 2015).
Mexico has an extraordinary bat diversity and abundance that are unusually greater than those found for other groups of mammals. In 2013, eight genotypes of α-CoVs and four of β-CoVs detected in Mexican bats showed phylogenetic distinctiveness, suggesting that the host species could act as key drivers of the CoV evolution (Anthony et al., 2013). Interestingly, a 6.9% CoVs prevalence was detected in bats sampled in the transition zone located in Central Mexico in 2015. Since the most abundant captured bat species in the sample coincided with the highest prevalence, this parameter could be partly explained by sampling bias. An alternate explanation could however be that gregarious species (the most abundant in this region) with higher rates of contact among individuals could have higher CoV prevalences, especially when considering species from large, dense colonies (Ortega and Arita, 1998; Schountz, 2014; Plowright et al., 2015; Martínez-Duque, 2017). CoVs positive individuals from species with lower abundances within the sampled refuge were also detected, indicating that high host population density is not required for CoVs to circulate within a bat population. This suggests that CoV genotypes may persist if adapted to a host population with consistent and frequent bat reproductive rates (Papkou et al., 2016). Hence, an integrated approach is required to understand aspects of ecology and evolution in local bat-CoVs dynamics.
Coronaviruses biogeographic patterns
Due to the variation in ecological processes and dynamics, the drivers of disease emergence also present spatial variations, such that the environmental and biogeographic characteristics of a region are significant determining factors (Peterson et al., 2004; Scheiner, 2009; Reperant, 2010). In this sense, the interaction between the virus, the host, and the environment, occur on an ecological time scale (hundreds of years), which leads to the confrontation to constant changes in short time periods. This causes alterations at different temporal and spatial scales with implications for both ecological and evolutionary processes (Carroll et al., 2007).
Based on our previous results (Anthony et al., 2013), and given that bat studies in the Americas have been linked to CoV presence in relatively few species (less than 15%), we studied the CoV-bat relationship in terms of its biogeographic affinity, comparing the CoV-bat local associations between the Neotropical region and the Nearctic - Neotropic transition zone. We found that the CoV geographic distribution seems to be influenced by host distribution, meaning that the viral genotype allocation could be mostly restricted to the host species distribution area (Martínez-Duque 2017). The identification of a high genotype divergence of CoVs in restricted geographic areas highlights CoVs adaptability towards different bat species. Since both host and habitat attributes are strong selective drivers in CoV evolution, even in geographically separated populations, the combination of evolutionary and biogeographic approaches is useful to broaden the CoVs study framework. While it is important to investigate coronavirus diversity in bats, it is also important to remember that the majority of viruses they carry will not pose any clinical risk, and bats should not be stigmatized ubiquitously as significant threats to public health.
Research perspectives
Community-level viral diversity
Measuring the genetic diversity of viruses is essential to understand their evolutionary dynamics, but reconstructing the genetic characteristics of a population, such as the composition and relative frequencies of viral genomes, represents a problem at a different scale from the experimental design to the analysis of information (Iwasa et al., 2004). The community level approach has been employed for the study of infectious diseases at a macroecological scale, requiring the identification of geographic distribution patterns of organisms (host, vectors and pathogens) involved in the ecology of a disease (Murray et al., 2015). Currently, the use of distribution models as a tool for identifying disease emergence in risk zones continues to increase, and advantages and disadvantages of the use of these models have been described (Anderson, 2012; Beck et al., 2014). The accuracy of distribution models depends strongly on data availability and quality, and lack of information can seriously bias their performance (Anderson, 2012; Beck et al., 2014). Nonetheless, distribution model outcomes help to determine priority areas to investigate CoVs, increasing efficient use of available resources and directing efforts towards regions where information is most needed.
Network theory is a useful community level approach that has been recently incorporated to facilitate the study and description of complex systems, as the associations between viruses and hosts, considering their respective diversities. Efforts have been undertaken to describe CoVs distribution patterns at a macroecological scale to detect communities within the network, represented as a set of closely related nodes. A global study of CoVs distribution revealed that these viruses show a high level of continental modularity. That means that CoVs genotypes are not shared among continents, with the exception of Asia and Africa, where only two are common to both continents (Anthony et al., 2017). In a recent study using network theory to estimate viral diversity and to describe the relationship of CoVs with other viral families and their hosts, we found that the CoVs have high centrality values, which are higher than those found for other commonly studied viruses such as Lyssavirus and Paramyxovirus. This is probably due to the large diversity of host species with which CoVs have been associated, as well as to a high host specificity (Nieto-Rabiela et al., 2019). In terms of disease ecology, the community level approach permits to explore host-host, virus-host and virus-virus interactions, by measuring network attributes, such as connectivity. Network analyses thus offer the opportunity to recognize highly diverse viruses and hosts based on a high degree of connectedness. This approach can be used to better understand of how interactions take place within pathogen communities, how hosts are connected with pathogens, their preferred associations and patterns of pathogen transmission. It is thus a useful tool for identifying hosts and viruses of potential importance in the zoonotic disease transmission dynamics.
Molecular analyses
Depending on the scale of the study, number of species (phylogenetic), or involved population (phylogeographic), molecular tools and phylogenetic analyses allow for the development and testing of hypotheses related to the dynamics of viruses and their hosts at macro- and micro-ecological scales (Denef et al., 2010; Hunt et al., 2008). These techniques are necessary to broaden our understanding of the mechanisms involved in the emergence of CoVs. Through the molecular identification and description of the evolutionary history of microorganisms, it is possible to establish the origins and patterns of dispersal and evolution involved in the emergence of disease in different populations. In relation to SARS-CoV-2 specifically, this approach is being intensively explored, mostly in relation to the evolution of the virus in the human population throughout different regions, and as a function of the time of local, regional, and worldwide spread (da Silva Candido et al., 2020; Deng et al., 2020; Lemey et al., 2020). Approaches based on molecular epidemiology have allowed for genetic and functional data analyses on SARS-CoV-2 with the aim of characterizing its evolutionary mechanisms, pathogenicity, and response to vaccines, as well as possible antiviral medications that are under development (Liu et al., 2019; Amanat and Krammer, 2020; Mavian et al., 2020). The molecular identification and characterization of CoVs in wildlife can also be useful for reservoir host identification and possible intermediate host SARS-CoV-2 association. An integrated evaluation of CoVs that are closely related to SARS-CoV-2 and that circulate in wildlife is essential to understand specific viral functions. While relevance of establishing the exact mechanisms by which SARS-CoV-2 arose and evolved in human populations through natural selection by molecular analyses is important, continuous monitoring of COVID-19 in humans and other animals is also of the utmost importance.
Conclusions
The current COVID-19 pandemic highlights an imperative need to improve the healthcare and environmental sectors, as well as the relationship of humankind with nature. The emerging One Health/Ecohealth paradigm must be transdisciplinary, encouraging intersectoral strategies at a national level, as well as international collaboration. Rethinking conservation, food production, and economic models that relate human and other species to natural resources has now high priority.
Disease ecology provides an integrated framework to develop novel research for infectious disease dynamics, leading to a better understanding of pathogen, wildlife hosts, domestic animals, and human interactions, that can only be accomplished by identifying the causes and processes associated with disease emergence. Transdisciplinary work that integrates biological, physical, and social sciences is essential for developing innovative and adaptive prevention and control strategies. A constant research and monitoring of potential emerging viruses at the interface among wildlife, domestic animals, and humans is also required. Simultaneously, evaluation of environmental and landscape characteristics that favor emerging infectious diseases, including land use and widening of the agricultural border, environmental degradation, and socioeconomic conditions, is fundamental to the development of adequate responses to future crises. Multi-host-multi-pathogen model integration in the study of CoVs is hence critical to identify high risk areas, potential hosts, and adequate prevention strategies.