On-Campus Dormitories as Viral Transmission Sinks: Phylodynamic Insights into Student Housing Networks During the COVID-19 Pandemic =================================================================================================================================== * Juan Bolanos * Alex Dornburg * April Harris * Samuel Kunkleman * Jannatul Ferdous * William Taylor * Jessica Schlueter * Cynthia Gibas ## Abstract University student housing environments are often viewed as hotspots for infectious disease transmission due to their high-density living conditions and high frequency of interpersonal interactions. During the COVID-19 pandemic, concerns arose that on-campus dormitories could serve as amplifiers of viral spread, seeding outbreaks into surrounding off-campus student residences. However, whether on-campus housing acts as a primary driver of transmission or as a recipient of infections introduced from the broader off-campus community remains unresolved. Here, we analyzed 1,431 SARS-CoV-2 genomes collected from students residing on and off campus at the University of North Carolina at Charlotte (UNCC) between September 2020 and May 2022. Using Bayesian phylodynamic and ancestral state reconstruction approaches, we traced viral transmission pathways to determine the directionality of spread between residential settings. Our results indicate that transmission from off-campus housing consistently seeded on-campus dormitory outbreaks. In contrast, viral movement from on-campus to off-campus housing was minimal. These patterns persisted across all major pandemic waves, regardless of shifting mitigation strategies, and suggest that on-campus residences acted as transmission sinks rather than sources of broader student outbreaks. These findings raise the possibility that on-campus residences may be more vulnerable than often considered, functioning as epidemiological ‘islands’ that primarily receive infections from off-campus sources. **Author Summary** University dormitories are often considered high-risk environments for disease spread due to their crowded living conditions and the inter-connectivity of individuals. This risk has led to the perception that university student housing could act as a major driver of outbreaks, raising concerns that on-campus dormitories fuel transmission into surrounding off-campus residences and the broader community. However, the actual transmission dynamics between on and off campus communities are often unclear. To address this, we analyzed over 1,400 SARS-CoV-2 genomes from students living on and off campus at the University of North Carolina at Charlotte between September 2020 and May 2022. Using Bayesian phylodynamic and ancestral state reconstruction approaches, we mapped viral transmission pathways and found that infections consistently flowed from off-campus residences into on-campus dormitories—not the other way around. This pattern persisted across multiple waves of the pandemic, regardless of changing mitigation strategies such as masking, testing, and vaccination policies. Our findings suggest on-campus housing functioned more like transmission “islands”—receiving infections from off-campus networks but rarely exporting them. These results suggest that universities should shift their focus toward managing transmission in off-campus student communities rather than concentrating mitigation efforts solely on dormitory settings. By reframing on-campus housing as a transmission *sink* rather than a *source*, our study provides critical insights for improving future outbreak and pandemic response strategies in university environments. ## Introduction The COVID-19 pandemic has underscored the critical importance of understanding viral transmission dynamics in densely populated environments. University campuses, which serve as microcosms of broader societal interactions, provide a unique lens for studying the spread and evolution of infectious diseases [1]. These environments are characterized by high mobility [2], tightly connected contact networks [3], and varied mitigation strategies in the face of disease outbreaks [4], making them ideal models for understanding how interventions shape transmission within and between residential settings. However, the relationship between on-campus and off-campus transmission is complex and context-dependent. Some studies suggest that reopening campuses can significantly elevate transmission rates in surrounding communities [5,6] with estimates indicating potential increases of over 50% [7]. Conversely, robust on-campus interventions such as testing, contact tracing, and wastewater surveillance can successfully limit spillover between residential groups [8–10]. These heterogeneous outcomes highlight the need for nuanced, data-driven analyses of campus-community transmission patterns to inform context-specific policies to optimize public health outcomes. Mathematical models and empirical datasets analyzing COVID-19 transmission have identified University campuses as possible superspreaders [11]. These findings are supported by additional studies demonstrating that university-related gatherings and close-contact environments played a significant role in local outbreaks during peak transmission periods [12]. However, the extent to which on-campus dormitories contribute to spread across dispersed off-campus university residential settings is not clear. The dynamics of dormitory populations may differ from the dynamics of the university population as a whole. For example, some universities implemented intervention strategies that significantly lowered transmission within campus housing communities [13]. These outcomes align with mathematical simulations that suggest robust contact tracing and demonstrate that intervention programs can successfully mitigate outbreaks in university settings [14]. College campuses also vary dramatically in their demographics and geographical settings, both of which are highly correlated to incidence patterns [15–17]. University residential settings span the full urban to rural gradient, from campuses located within major population centers to those located within small isolated communities. Most studies have focused on isolated extremes of this spectrum of settings [12], raising the question of how universities in urban environments with off-campus residences scattered throughout high density metropolitan housing areas differ in transmission risks and dynamics. Urban campuses are characterized by high population densities, extensive student clustering, and substantial daily movement between on campus dormitories, off-campus apartments, and other communal spaces [18,19]. The high frequency of interpersonal interactions, shared spaces, and proximity to common gathering areas could position on-campus housing as viral incidence amplification sites that facilitate streams of transmission to off-campus residences. An alternative hypothesis, though, is that on-campus residences may act as epidemiological ‘islands’, where infections are introduced by contact with individuals residing off-campus who have had greater exposure to contacts in the broader community. In this scenario,viral strains are introduced to campus, and may circulate within the campus population, but contribute little to overall community spread beyond the campus [8,20]. Supporting this latter scenario, Purdue University detected at least 10 independent introductions of the gamma variant onto its campus in 2021 [21], emphasizing the possibility of multiple entry points into on campus residences. Disentangling these contrasting hypotheses–akin to classic source-sink dynamics in biogeography [22–25]– is critical for understanding the degree to which student housing acts as a self-contained transmission network or as a gateway for broader community spillover. To address this question, we apply a methodological approach capable of tracing viral movement with high resolution. Traditional epidemiological metrics such as incidence data have been widely used to reveal overall trends in infection rates between campus and other environments [5]. However, these metrics cannot resolve the history of transmission between individuals in the absence of detailed contact tracing. Recent advances in phylodynamic approaches facilitate the integration of viral genomic data with epidemiological models permit reconstruction of transmission pathways with robust resolution [26,27]. Such approaches have been widely used to track the spread of infectious diseases at global and national scales [28–30]. Instances of their application to university campuses is less common [12], creating a knowledge gap in our understanding of general trends in the transmission dynamics of viruses between campus populations. Here we integrated extensive SARS-CoV-2 genomic sequence data (n=1431 genomes) with public health and intervention data from the University of North Carolina at Charlotte (UNCC), a mid-sized urban institution serving over 31,000 students within the densely populated Charlotte metropolitan region [31]. We analyzed these genomic and epidemiological data in a Bayesian phylodynamic framework. We apply time-calibrated phylogenies and discrete ancestral state reconstructions [32] to test three key hypotheses regarding the role of student housing settings in transmission dynamics. First, we assessed the degree to which transmission dynamics between on- and off-campus populations reflect alternate possible source-sink scenarios. Specifically, we determined whether the university primarily drives spillover into surrounding off-campus housing or functions as a recipient of infections introduced from off-campus residences. Second, we evaluated whether shifts in on-campus mitigation strategies, such as masking mandates, asymptomatic and symptomatic testing, and dormitory occupancy limits, corresponded to measurable changes in transmission patterns within and between campus and off-campus residences. Lastly, we examined whether viral lineage dynamics in residences represented early indicators of broader outbreaks or if infections lagged behind those in the greater Charlotte metropolitan area. Collectively, our results challenge the prevailing popular media notion that on campus residences act as primary superspreader hubs or reservoirs amplifying off-campus transmission. Instead, they show a higher rate of transmission into the campus under all mitigation conditions examined. ## Results & Discussion ### Campuses dormitories can act as sinks of viral transmission We sequenced viral genomes from 1444 SARS-CoV-2 positive clinical samples collected at the University of North Carolina at Charlotte (UNCC) between September, 2020 and May, 2022. These samples were from 470 and 961 patients living on-campus or off-campus, respectively. Using this data, we estimated a time calibrated Bayesian phylogenetic framework to model the transmission pathways of SARS-CoV-2 between on-campus and off-campus student populations via ancestral state reconstruction (ASR).(**Figure 1A**). Given the location of the school and the large number of commuter students, our initial hypothesis was that on campus dormitories could act as a primary source of infections, fueling outbreaks among off-campus students and facilitating spread into the broader community. However, our results strongly support that this is not the case. (**Figure 1**). Assessing the distribution of off and on campus incidence across the phylogenetic history of transmission revealed numerous shifts between on and off campus residences (**Figure 1A**). In most cases, transmission pathways reflect bursts of infections within either the on or off campus community, reflecting a higher rate of transmission within either community than between communities (**Figure 1B**). In cases of transmission between on and off campus residences, transmission from off-campus residences into campus residences dramatically outnumber the reverse (**Figure 1C).** This asymmetry suggests that viral introductions into the university environment were primarily seeded from off-campus sources rather than spreading outward into the community, with as many as one fifth of the total number of cases on-campus being a direct consequence of transmission from off-campus residents. These results were robust to bootstrapping of the data (**Supplemental Figure 2**), and they challenge assumptions that university dormitory populations act as amplifiers of community transmission. Rather, our findings suggest the on-campus student population is a transmission sink within the larger context of the student housing network.  [Figure 1.](http://medrxiv.org/content/early/2025/03/10/2025.03.07.25323493/F1) Figure 1. Transition history and transition dynamics between on and off campus communities during the study period. **(A)** Maximum likelihood tree topology of SARS-CoV-2 sequences collected during this study. Tip labels indicate whether samples were derived from on campus (cool shading) or off campus (warm shading) locations. * symbols in the delta strain portion of the phylogeny demarcate lineages identified as alpha by diagnostic mutations, but are phylogenetically resolved within delta (see supplemental materials). **(B)** Box plots depict the estimated quantiles of the log-transformed number of transitions occurring between communities. Top row reflects the frequency of transmissions occurring between the campus and off campus community (left) or between the off campus community and the campus (right). Bottom row reflects the frequency of transmissions within the campus (left) or within the off campus community (right). Cartoons above each box plot indicate the transmission mode and correspond to the legend in **(A)**. **(C)** The tempo of transmission between communities. Violin plots depict the estimated distribution of evolutionary transition rates estimated in BEAST, truncated to the 3rd quartile (75%). **(D)** Relative frequency of transmissions between and within communities. Chord diagram depicts the relative frequency of transmissions, shaded by community as in **(A)**. Outer bands correspond to total relative cases for each community, and inner chords illustrate transmission mode scaled to their frequency. Arrows and cartoons illustrate the direction of transmission. Universities are often assumed to be amplifiers of viral transmission due to high-density living and frequent interpersonal interactions [11]. However, our results suggest that the off-campus student housing was the primary driver of infection trends at UNC Charlotte, with the on-campus population largely reflecting, rather than shaping, viral transmission dynamics among off-campus students. The clustering of on-campus infections following introductions from the off-campus population indicates that once introduced, transmission among those students living in the university environment was sustained at lower overall rates than among those living in the surrounding community. This pattern suggests that on-campus student residences are vulnerable to external seeding events and may not act as primary drivers of viral transmission in the community. As viral genomic data from COVID-19 pandemic continues to be analyzed, evidence for multiple introductions to urban college campuses that seeded in situ outbreaks continue to accrue, as does evidence for these environments generating only minimal spillover into their urban surroundings [21,33]. More broadly, these results suggest that on campus housing may follow similar transmission dynamics as other highly structured sub-populations such as nursing homes or long term care facilities [34]. It is becoming increasingly evident that community-wide mobility networks are predictive of infection pathways in the urban landscape [35], suggesting that the degree of connection of students to the community outside of campus is likely a major determinant of transmission into the campus. Our data lack fine-scale resolution concerning mobility between neighborhoods and individuals, precluding an analysis of whether certain areas have disproportionately seeded the campus. Regardless, our findings emphasize the importance of considering universities within the broader landscape of urban disease ecology [36,37]. Rather than viewing university dormitories as isolated entities shaping incidence patterns across a metropolitan region, integration of urban campus residence data into the broader urban community structure [38] will allow communities to harness information from universities’ healthcare platforms to understand how these settings are embedded within complex mobility networks that mediate pathogen flow across the urban landscape. ### Mitigation strategies impact transmission rates, not source-sink dynamics To assess whether shifts in mitigation strategies corresponded to changes in transmission dynamics between on-campus and off-campus student housing, we analyzed SARS-CoV-2 lineage movement across four major phases of the pandemic: Wuhan (Fall 2020), Alpha (Spring–Summer 2021), Delta (Fall 2021), and Omicron (Spring 2022). Throughout this period, mitigation strategies transitioned from stringent controls—such as reducing dormitory occupancy to one-third capacity, enforcing masking and social distancing, and relocating SARS-CoV-2-positive students to designated isolation dormitories—to increasingly relaxed control measures that culminated in the resumption of in-person learning with many previous mitigation strategies optional (**see supplemental materials**). As measures relaxed, so did compliance with self-reported symptom tracking, mirroring broader trends in pandemic fatigue and shifting public health priorities [39–41]. These changes in campus policy largely aligned with evolving state and federal guidelines for North Carolina (**supplemental materials**), and raise the possibility that viral transmission dynamics between off and on campus may experience concomitant shifts that amplify spillover between these communities. However, our analyses strongly suggest this to not be the case. Considering transmission dynamics across major variant-driven case surges of the pandemic reveals incidence trends that consistently aligned with mitigation shifts (**Figure 2**). In general, relaxation of restrictions led to higher frequencies of incidence in the on campus population, as well as a separation between the distribution of transitions between on campus cases, and transitions from off campus onto campus, with higher levels of on campus circulation (**Figure 2**). For example, during the initial Wuhan phase, when stringent mitigation measures such as mask mandates, social distancing, and routine testing were in place, on-campus transmission was equivalent to transitions from off to on campus (**Figure 2A**). In contrast, incidence levels during Alpha, when many restrictions were lifted, led to on campus transmission events outpacing the transitions from off to on campus (**Figure 2B**). This same pattern was repeated when comparing delta and omicron, the latter of which vastly exceeded transmission levels observed in other phases of the pandemic, both on and off campus (**Figure 2C&2D; Supplemental Table 1**). These incidence patterns reflect a global trend, as Omicron’s increased transmissibility and immune evasion properties contributed to surges of infection worldwide [42–44].  [Figure 2.](http://medrxiv.org/content/early/2025/03/10/2025.03.07.25323493/F2) Figure 2. State transition distributions across four major waves of the pandemic: (A) Wuhan; (B) Alpha; (C) Delta; and (D) Omicron. The frequency (Y axis) of transitions (X axis) between and within residence types are indicated by the shadings in the figure legend with areas indicating the direction of transmission for each shading (i.e., off-campus to on-campus). However, regardless of time period, movement of viruses from the campus residences into the off-campus community remained unchanged. Transmission occurred primarily from the off-campus residences into on-campus residences, rather than the reverse (**Figure 3**). Even as on-campus transmission increased during Omicron, there was no corresponding rise in outward viral movement from campus into the community (**Figure 2D**; **Figure 3**). These patterns were robust to resampling, indicating this result is not an artifact of differential incidence levels between on-campus and off-campus populations alone (**Supplemental Figure 2**). Instead, the persistent asymmetry between off-campus to on-campus and on-campus to off-campus transitions suggests that the primary drivers of viral introduction into the university were external, rather than originating within campus populations. This pattern mirrors findings from other university settings during Omicron, where campuses exhibited high in situ transmission but limited spillover into the surrounding community [9,45]. These results suggest that, regardless of mitigation measures, UNC Charlotte continued to act primarily as a sink rather than a source of community transmission.  [Figure 3.](http://medrxiv.org/content/early/2025/03/10/2025.03.07.25323493/F3) Figure 3. Transmission dynamics between on and off campus across variant waves: (A) Wuhan; (B) Alpha; (C) Delta; and (D) Omicron. The relative frequency of transmissions between and within communities is depicted in chord diagram, shaded by community as in Figure 1. Outer bands correspond to total relative cases for each community, and inner chords illustrate transmission mode scaled to their frequency. Arrows and cartoons illustrate the direction of transmission. There is no doubt that initial pandemic restrictions and mitigation strategies saved millions of lives across the globe [46]. Similarly, targeted university interventions such as masking mandates and routine testing are estimated to have reduced infections on University campuses by an order of magnitude or more [47]. However, universities, like other sectors [48], had to navigate a complex trade-off between the immediate health benefits of mitigation and the long-term academic, social, and economic consequences of restrictive policies [49,50]. Advances in genomic surveillance methods [51,52] as well as computational models have come to play a crucial role in balancing disease suppression with the need for in-person education [53,54]. Our results suggest that placing the university within the broader framework of urban disease ecology could provide vital additional context to multiple urban sectors. In particular, our finding of a marked asymmetry in viral transmission, that persisted across mitigation strategies, suggests that on campus residences such as those at UNC Charlotte may be epidemiologically insular, acting as relatively self-contained transmission units rather than perpetual community wide spreaders. Further work leveraging computational [55] and phylodynamic approaches [26,56] are critically needed to assess whether similar trends hold across other urban campuses. ### On campus infection waves lag behind off campus waves If on-campus student housing were a major source of viral transmission, we would expect infection waves among dormitory residents to precede or at least coincide with those in off-campus housing. However, our phylodynamic analyses reveal a consistent lag in on-campus infection waves relative to those in off-campus housing, with peaks in on-campus lineage counts occurring after surges in off-campus populations (**Figure 4A**). This pattern is particularly pronounced for the Delta and Omicron waves (**Figure 4**), where variant-driven outbreaks first became established in off-campus housing before spreading into dormitories. Incidence peaks largely align with the broader emergence of new variants of concern in the surrounding community (**Figure 4B**), and the delayed pattern of spread is consistent with the expected transmission lag when considering the average incubation period of SARS-CoV-2, which ranges from 5 to 14 days [57,58]. The variation between off and on campus incidence patterns over time was particularly pronounced during periods of higher mitigation on campus, which resulted in an overall muting of peak incidence patterns (**Supplemental Figure 3**), and likely reflects the efficacy of mitigation measures. Collectively, these results imply that consistent transmission from source off-campus residences seeded numerous on campus cases and outbreaks during all waves of the pandemic, with times of heightened use of on-campus mitigation strategies, controlled environments, and testing preventing further in situ outbreaks.  [Figure 4.](http://medrxiv.org/content/early/2025/03/10/2025.03.07.25323493/F4) Figure 4. Trends in lineages through time and incidence. (A) SARS-COV-2 relative frequency of lineages through time for off-campus (warm shading) and on-campus (cool shading) residences plotted in comparison to incidence in Mecklenburg county (gray). Circles in lineage through time plots represent specific time points. **(B)**. Abundance of circulating strains during the time period depicted in **A**. Strains are indicated in the legend. Throughout the COVID-19 pandemic, transmission dynamics often varied spatially, with some areas experiencing infection waves weeks ahead of others [59]. This staggered pattern of viral spread can occur both at larger geographic scales (e.g., hemisphere, country, etc) as well as at more localized scales (e.g., state, county, etc) [60,61]. Understanding variation in the spread of disease is of particular concern in densely populated urban environments where social connectivity, mobility patterns, and mitigation strategies influence the timing and intensity of outbreaks [62]. If universities were major sources of community spread, we would expect campus surges to precede infection waves in off-campus residences, that, in turn, would seed transmission through the larger metropolitan area. However, our findings strongly suggest the reverse. The on campus residences at UNC Charlotte largely acted as recipient of community viral spread, rather than an independent epicenter of transmission. This pattern of transmission raises the possibility that off campus housing acts as a gateway enabling viruses circulating in the community to enter and circulate on campus. Whether there are asymmetrical pathways of transmission such as this between the broader community and the campus remains unknown and represents an exciting area of future research. ## Conclusion Our findings demonstrate that on-campus student housing at UNC Charlotte functioned as a transmission sink rather than a source of community-wide viral spread throughout the COVID-19 pandemic. By leveraging whole-genome sequencing and Bayesian phylogenetic analyses, we reconstructed SARS-CoV-2 transmission pathways and found strong evidence that viral introductions into dormitories overwhelmingly originated from off-campus housing rather than spreading outward from campus into the broader student population. This pattern persisted across multiple pandemic waves and mitigation phases, suggesting a consistent epidemiological role for on-campus residences as recipients of community transmission rather than amplifiers of broader outbreaks. Importantly, even during periods of high incidence—such as the Omicron surge—there was no corresponding increase in outward viral movement from dormitories to off-campus residences. These results suggest on-campus housing to possibly be similar in expected viral transmission dynamics to other structured, high-density subpopulations that have limited interactions with the broader community, such as long-term care or military facilities, in which external introductions predominantly seed internal outbreaks [63–65]. Future work assessing additional universities with similar on and off campus residence structures is critically needed to determine the degree to which these transmission source-sink dynamics persist across institutional and geographic contexts. ## Materials and Methods ### Sample Collection and RNA extraction Samples were obtained from the Covid Testing Center (later absorbed into the Student Health Center) at the University of North Carolina Charlotte, and StarMed, an urgent care network serving Mecklenburg and surrounding counties. Collection involved buccal nasal swabs from symptomatic as well as asymptomatic individuals, followed by viral RNA extraction by the respective facility. StarMed utilized the Kingfisher system with the MagMAX™-96 Total RNA Isolation Kit (Catalog No: AM1830) for RNA extraction, while the UNC Charlotte testing center employed the KingFisher Flex with the MagMAX™ Viral/Pathogen II (MVP II) Nucleic Acid Isolation Kit (Catalog No: A48383). After extraction, qPCR was conducted to detect SARS-CoV-2 positive samples. UNC Charlotte used three probes (N1, N2, and S primers, designed by the CDC) with the TaqPath™ COVID-19 Combo Kit (Catalog No: A47814), whereas StarMed relied solely on the N1 primer for detection, using the same kit. Only RNA from SARS-CoV-2 positive samples was forwarded for sequencing. This sampling resulted in 1444 unique positive cases of which 1431 were successfully sequenced. This sample represents nearly 90% of all known campus associated SARS-CoV-2 cases during the study period, with only approximately 200 additional cases voluntarily reported by individuals who were tested in locations away from campus without a surveillance sequencing program. ### Library preparation and sequencing RNA extracts of clinical samples were reverse-transcribed and amplified in PCR according to a tiling amplicon protocol modified from Pater et al. [66]. A fully detailed guide to this modified protocol is available on protocols.io [67]. Briefly, first we transformed extracted RNA into complementary DNA (cDNA), with cDNA enrichment of the SARS-CoV-2 genome performed by multiplexing PCR with two pools of ARTIC primers. For each sample, the ARTIC PCR products were pooled and amplicons were isolated through magnetic Solid Phase Reversible Immobilization (SPRI) beads and quantified by Qubit fluorescence. Subsequent sequencing library preparation comprised amplicon end-prep, sample barcoding, sample pooling and cleaning, and sequencing adapter ligation. Finally, the sample library was loaded onto an R9 flow cell (Oxford Nanopore Technologies, ONT) and sequenced using the ONT PromethION instrument. Manufacturer-recommended settings for the LSK-109 genomic sequencing by ligation kit, Native Barcode Expansion Kit (Exp-NBD196), and FLO-PRO002 flow cell were applied at the start of the run. ### Sequence assembly and variant calling Basecalling was initially carried out in real time using the Oxford Nanopore guppy software with high accuracy setting, via the MinKnow user interface on the ONT PromethION instrument. During spring of 2021, updates to ONT basecalling algorithms became available and were applied as they were released. In order to avoid biases arising due to changes in basecalling algorithms, sequences were re-basecalled from raw fast5 files using the ONT 1.2.1 “super high accuracy” basecalling method. Assembly and variant calling of SARS-CoV-2 sequences was performed using a modified version of the ARTIC ‘fieldbioinformatics’ pipeline ([https://github.com/artic-network/fieldbioinformatics](https://github.com/artic-network/fieldbioinformatics)). Briefly, reads were filtered and trimmed with guppyplex to remove V3 adapter sequences and to exclude chimeric sequences smaller than 250 bases and longer than 700 bases. Alignments were constructed using minimap2, followed by medaka ([https://github.com/nanoporetech/medaka](https://github.com/nanoporetech/medaka)) for polishing, variant calling, and consensus building using the ‘r941\_prom\_high_g360’ model. Consensus sequences were classified using Pangolin ([https://github.com/cov-lineages/pangolin](https://github.com/cov-lineages/pangolin)) versions 3.1.11, 3.1.14, 3.1.16, 3.1.17, 3.1.19, 3.1.20, and 4.1.2 to ensure accurate lineage assignments over the course of the study (For example, earlier versions would not be able to accurately classify Delta or Omicron due to their unique mutation signature). Sequences were inspected for quality based on the pass/fail metrics in Pangolin that assess minimum sequence length and maximum allowable ambiguous bases (N content).. Of the sequenceable samples, 13 failed quality control at this classification step and were excluded from downstream analyses. Additional quality control was conducted using Nextclade [68], which assesses missing data, ambiguous sites, mutation clusters, stop codons, and frame shifts for SARS-CoV-2 ([https://github.com/nextstrain/ncov](https://github.com/nextstrain/ncov)). ### Phylogenetic Analyses A total of 1444 viral samples were collected from both on-campus and off-campus populations between September, 2020 and May, 2022. Of these, 1431 passed the sequence quality control assessments described above. The final set of sequences (available in GISAID with EpiSet ID: EPI\_SET\_250220my) is composed of 470 samples originating from on-campus students and 961 from off-campus individuals. We used the Wuhan strain (accession [NC\_045512](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=NC_045512&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom)) to serve as a reference genome and outgroup for analysis. Sequences were aligned using MAFFT v7.525 [69]. To estimate the evolutionary history of our collected samples, we used MAPLE (MAximum Parsimonious Likelihood Likelihood Estimation), a computational tool designed for large-scale phylogenetic inference of epidemiological datasets that has been shown to significantly reduce computational analysis time while increasing accuracy relative to other likelihood approaches [70]. We used RelTime with Dated Tips (RTDT) [71] in Mega X [72] to time calibrate the MAPLE inferred maximum likelihood topology, allowing us to estimate the transmission timing of SARS-CoV-2 across on-campus and off-campus student housing. RTDT is an algebraic method designed for estimating divergence times from temporally sampled molecular sequences, making it particularly well-suited for analyzing pandemic-scale datasets with precise sampling time information. RTDT incorporates a relative rate framework, enabling computational efficiency while maintaining a level of accuracy that is par with more computationally complex Bayesian approaches to divergence time estimation such as those implemented in BEAST [73–75], including for the analysis of SARS-CoV-2 data [76]. Divergence times were calibrated using the time each viral genome was sampled. ### Reconstructing the history of transmission We used the discrete phylogeographic diffusion model implemented in BEAST v2.6.7 [32,77,78] to reconstruct the transmission history of SARS-CoV-2 between student housing communities. This model integrates ancestral state reconstruction (ASR), a well-established method in evolutionary inference [79–82], with a diffusion model that models viral spread through distinct transmission events between predefined geographic units [83]. This approach allowed us to quantify the directionality and frequency of transitions between on-campus and off-campus housing while integrating branch-specific rate variation to accommodate heterogeneity in transmission dynamics [83]. We fixed the time calibrated phylogeny estimated above, and ran two independent markov chains for up to 250,000,000 generations to estimate transmission histories. Runs were stopped at 232,979,000 generations as inspection of the likelihood and other parameter values in Tracer v1.7.2 [84] indicated convergence between runs, with all effective sample size (ESS) values above 200, indicating effective sampling of the target distributions. To verify that posterior estimates were informed by the empirical data, we conducted an additional BEAST analysis in which parameters were sampled exclusively from the prior distributions.. We used the Babel package ([https://github.com/rbouckaert/Babel](https://github.com/rbouckaert/Babel)) to quantify the distribution of on or off campus lineages through time with the LineagesThroughTimeCounter function, and also quantified transitions between on and off campus using the transitionStateCounter function. Resulting transition histories were imported into R visualized using Circos plots with functions from circlize v0.4.16 [85,86] treedataverse v0.0.1 [87,88], and tidyverse v0.4.6 [89] libraries for additional visualizations. Transmission histories were visually compared to incidence frequencies in Mecklenburg County using public incidence data available from the North Carolina Respiratory Virus Summary Dashboard ([https://covid19.ncdhhs.gov/dashboard](https://covid19.ncdhhs.gov/dashboard)). Our dataset captured approximately 90% of all cases known within these populations during the study period. However, it is possible that the uneven levels of incidence between off and on campus might bias estimation of transmission dynamics. To assess the degree to which possible sampling biases of campus populations impact estimated transmission dynamics, the above analyses were repeated on randomly subsampled data sets with predefined on-campus:off-campus ratios (1:5, 2:3, 1:1, 3:2, 5:1). For each set of datasets at a given ratio (e.g. 5 on campus sampler per very 1 off campus sample), total transmissions between and within population (off-campus to off-campus, on-campus to on-campus, off-campus to on-campus, and on-campus to off-campus) were quantified to ensure a consistent signal of source-sink dynamics relative to that based on the observed incidences of COVID-19. ## Supporting information Supplemental Materials [[supplements/323493_file02.pdf]](pending:yes) ## Data Availability All data produced in the study are available as a GISAID EpiSet or as part of a Zenodo data package. DOI/URL information for supplementary materials shared via EpiCoV, Zenodo and GitHub are provided. [https://doi.org/10.55876/gis8.250220my](https://doi.org/10.55876/gis8.250220my) [https://doi.org/10.5281/zenodo.14946855](https://doi.org/10.5281/zenodo.14946855) [https://github.com/jbolanos93/Covid-Campus-Phylodynamics.git](https://github.com/jbolanos93/Covid-Campus-Phylodynamics.git) ## Funding This research was supported by COVID-19 contract funding from The NC Collaboratory CORVASEQ project (CG and JS), Mecklenburg County Public Health (CG, JS, WT, AH, JF), and the State of North Carolina (CG, JS and AD) and by National Science Foundation of the United States of America RAPID 2031204 (AD). JB was supported by a graduate teaching assistantship from the UNC Charlotte Department of Bioinformatics and Genomics. ## Conflict of Interest The authors declare no conflict of interest. ## Declaration of generative AI and AI-assisted technologies in the writing process During the preparation of this work the author(s) used ChatGPT in order to obtain feedback and enhance the manuscript. After using this tool/service, the author(s) reviewed and edited the content as needed and take(s) full responsibility for the content of the publication. ## CRediT author roles Bolanos: Conceptualization, Methodology, Investigation, Formal Analysis, Visualization, Software, Writing – Original Draft, Writing – Review and Edit Dornburg: Conceptualization, Methodology, Visualization, Funding Acquisition, Writing – Original Draft, Writing – Review and Edit Harris: Resources, Data Curation, Investigation Kunkleman: Data Curation, Software, Methodology Ferdous: Methodology, Investigation Taylor: Project Administration, Methodology, Investigation Schlueter: Funding Acquisition, Project Administration, Supervision, Writing – Review and Edit Gibas: Conceptualization, Funding Acquisition, Project Administration, Supervision, Visualization, Writing – Original Draft, Writing – Review and Edit ## Data Availability SARS-CoV-2 sequence has been deposited to GISAID (doi: 10.55876/gis8.250220my). All code used for analyses has been archived on Zenodo (DOI:10.5281/zenodo.14946855; URL: [https://doi.org/10.5281/zenodo.14946855](https://doi.org/10.5281/zenodo.14946855)). ## Acknowledgements Dr. Robert Jones at the UNC Charlotte Student Health Center and Dr. Angelica Martins in the UNC Charlotte Division of Research provided samples and COVID-19 diagnostic test metadata as specified under UNC Charlotte protocol IRB 21-0113. We thank them, their teams, and the other staff of the UNC Charlotte Division of Research for their collaboration and tireless public health efforts during the COVID-19 pandemic emergency. * Received March 7, 2025. * Revision received March 7, 2025. * Accepted March 9, 2025. * © 2025, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Stubbs CW, Springer M, Thomas TS. The impacts of testing Cadence, mode of instruction, and student density on Fall 2020 COVID-19 rates on campus. medRxiv. 2020. p. 2020.12.08.20244574. doi:10.1101/2020.12.08.20244574 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMi4wOC4yMDI0NDU3NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjUvMDMvMTAvMjAyNS4wMy4wNy4yNTMyMzQ5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 2. 2.Bhandari R, Blumenthal P. Global student mobility and the twenty-first century silk road: National trends and new directions. International Students and Global Mobility in Higher Education. New York: Palgrave Macmillan US; 2011. pp. 1–23. 3. 3.Israel U, Koester BP, McKay TA. Campus connections: Student and course networks in higher education. Innov High Educ. 2020;45: 135–151. 4. 4.Klein B, Generous N, Chinazzi M, Bhadricha Z, Gunashekar R, Kori P, et al. Higher education responses to COVID-19 in the United States: Evidence for the impacts of university policy. PLOS Digit Health. 2022;1: e0000065. 5. 5.Andersen MS, Bento AI, Basu A, Marsicano CR, Simon KI. College openings in the United States increase mobility and COVID-19 incidence. PLoS One. 2022;17: e0272820. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36037207&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 6. 6.Walke HT, Honein MA, Redfield RR. Preventing and Responding to COVID-19 on College Campuses. JAMA. 2020;324: 1727–1728. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 7. 7.Leidner AJ, Barry V, Bowen VB, Silver R, Musial T, Kang GJ, et al. Opening of large institutions of higher education and county-level COVID-19 incidence - United States, July 6-September 17, 2020. MMWR Morb Mortal Wkly Rep. 2021;70: 14–19. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33411699&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 8. 8.Valesano AL, Fitzsimmons WJ, Blair CN, Woods RJ, Gilbert J, Rudnik D, et al. SARS-CoV-2 genomic surveillance reveals little spread from a large University campus to the surrounding community. Open Forum Infect Dis. 2021;8: ofab518. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34805437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 9. 9.Park SW, Daskalaki I, Izzo RM, Aranovich I, Te Velthuis AJW, Notterman DA, et al. Relative role of community transmission and campus contagion in driving the spread of SARS-CoV-2: Lessons from Princeton University. PNAS Nexus. 2023;2: gad201. 10. 10.Clancey E, Mietchen MS, McMichael C, Lofgren ET. Unexpected Transmission Dynamics in a University Town: Lessons from COVID-19. medRxiv. 2024. doi:10.1101/2024.01.10.24301116 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyNC4wMS4xMC4yNDMwMTExNnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjUvMDMvMTAvMjAyNS4wMy4wNy4yNTMyMzQ5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. 11.Lu H, Weintz C, Pace J, Indana D, Linka K, Kuhl E. Are college campuses superspreaders? A data-driven modeling study. Comput Methods Biomech Biomed Engin. 2021;24: 1136–1145. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33439055&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 12. 12.Andrews KR, New DD, Gour DS, Francetich K, Minnich SA, Robison BD, et al. Genomic surveillance identifies potential risk factors for SARS-CoV-2 transmission at a mid-sized university in a small rural town. Sci Rep. 2023;13: 7902. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37193760&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 13. 13.Lee J, Acosta N, Waddell BJ, Du K, Xiang K, Van Doorn J, et al. Campus node-based wastewater surveillance enables COVID-19 case localization and confirms lower SARS-CoV-2 burden relative to the surrounding community. Water Res. 2023;244: 120469. 14. 14.Brown RA. A simple model for control of COVID-19 infections on an urban campus. Proc Natl Acad Sci U S A. 2021;118. doi:10.1073/pnas.2105292118 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxODoiMTE4LzM2L2UyMTA1MjkyMTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjUvMDMvMTAvMjAyNS4wMy4wNy4yNTMyMzQ5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 15. 15.Karmakar M, Lantz PM, Tipirneni R. Association of Social and Demographic Factors With COVID-19 Incidence and Death Rates in the US. JAMA Netw Open. 2021;4: e2036462. 16. 16.Pijls BG, Jolani S, Atherley A, Derckx RT, Dijkstra JIR, Franssen GHL, et al. Demographic risk factors for COVID-19 infection, severity, ICU admission and death: a meta-analysis of 59 studies. BMJ Open. 2021;11: e044640. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMjoiMTEvMS9lMDQ0NjQwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjUvMDMvMTAvMjAyNS4wMy4wNy4yNTMyMzQ5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 17. 17.Zhang CH, Schwartz GG. Spatial Disparities in Coronavirus Incidence and Mortality in the United States: An Ecological Analysis as of May 2020. J Rural Health. 2020;36: 433–445. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 18. 18.Doyle-Baker PK, Ladle A, Rout A, Galpern P. Smartphone GPS Locations of Students’ Movements to and from Campus. ISPRS International Journal of Geo-Information. 2021;10: 517. 19. 19.Kirk CM, Lewis RK. Sense of community on an urban, commuter campus. International Journal of Adolescence and Youth. 2015 [cited 30 Jan 2025]. doi:10.1080/02673843.2013.763833 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/02673843.2013.763833&link_type=DOI) 20. 20.Bharti N, Lambert B, Exten C, Faust C, Ferrari M, Robinson A. Large university with high COVID-19 incidence is not associated with excess cases in non-student population. Sci Rep. 2022;12: 3313. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-022-07155-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35228585&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 21. 21.Ciubotariu II, Dorman J, Perry NM, Gorenstein L, Kattoor JJ, Fola AA, et al. Genomic surveillance of SARS-CoV-2 in a university community: Insights into tracking variants, transmission, and spread of Gamma (P.1) variant. Open Forum Infect Dis. 2022;9: ofac268. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35818365&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 22. 22.Dornburg A, Federman S, Lamb AD, Jones CD, Near TJ. Cradles and museums of Antarctic teleost biodiversity. Nat Ecol Evol. 2017;1: 1379–1384. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29046532&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 23. 23.Gilroy JJ, Edwards DP. Source-sink dynamics: A neglected problem for landscape-scale biodiversity conservation in the tropics. Curr Landsc Ecol Rep. 2017;2: 51–60. 24. 24.Goldberg EE, Roy K, Lande R, Jablonski D. Diversity, endemism, and age distributions in macroevolutionary sources and sinks. Am Nat. 2005;165: 623–633. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/430012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15937743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229047500005&link_type=ISI) 25. 25.Schurr FM, Pagel J, Cabral JS, Groeneveld J, Bykova O, O’Hara RB, et al. How to understand species’ niches and range dynamics: a demographic research agenda for biogeography. J Biogeogr. 2012;39: 2146–2162. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2699.2012.02737.x&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000311384100008&link_type=ISI) 26. 26.Volz EM, Koelle K, Bedford T. Viral phylodynamics. PLoS Comput Biol. 2013;9: e1002947. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1002947&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23555203&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 27. 27.Baele G, Dellicour S, Suchard MA, Lemey P, Vrancken B. Recent advances in computational phylodynamics. Curr Opin Virol. 2018;31: 24–32. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.COVIRO.2018.08.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30248578&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 28. 28.Attwood SW, Hill SC, Aanensen DM, Connor TR, Pybus OG. Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic. Nat Rev Genet. 2022;23: 547–562. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41576-022-00483-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35459859&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 29. 29.Lin SYC, Toyoda H, Kumada T, Liu H-F. Molecular Evolution and Phylodynamics of Acute Hepatitis B Virus in Japan. PLoS One. 2016;11: e0157103. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0157103&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27280441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 30. 30.Truong Nguyen P, Kant R, Van den Broeck F, Suvanto MT, Alburkat H, Virtanen J, et al. The phylodynamics of SARS-CoV-2 during 2020 in Finland. Commun Med (Lond). 2022;2: 65. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35698660&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 31. 31.Brown P. UNC Charlotte sets historic enrollment record. In: Inside UNC Charlotte [Internet]. 5 Sep 2024 [cited 28 Jan 2025]. Available: [https://inside.charlotte.edu/2024/09/05/unc-charlotte-sets-historic-enrollment-record/](https://inside.charlotte.edu/2024/09/05/unc-charlotte-sets-historic-enrollment-record/) 32. 32.Bouckaert R, Vaughan TG, Barido-Sottani J, Duchêne S, Fourment M, Gavryushkina A, et al. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2019;15: e1006650. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1006650&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30958812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 33. 33.Soda KJ, Chen X, Feinn R, Hill DR. Monitoring and responding to emerging infectious diseases in a university setting: A case study using COVID-19. PLoS One. 2023;18: e0280979. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37196023&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 34. 34.Wilmink G, Summer I, Marsyla D, Sukhu S, Grote J, Zobel G, et al. Real-Time Digital Contact Tracing: Development of a System to Control COVID-19 Outbreaks in Nursing Homes and Long-Term Care Facilities. JMIR Public Health Surveill. 2020;6: e20828. 35. 35.Wen T-H, Hsu C-S, Hu M-C. Evaluating neighborhood structures for modeling intercity diffusion of large-scale dengue epidemics. Int J Health Geogr. 2018;17: 9. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29724243&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 36. 36.Kache PA, Santos-Vega M, Stewart-Ibarra AM, Cook EM, Seto KC, Diuk-Wasser MA. Bridging landscape ecology and urban science to respond to the rising threat of mosquito-borne diseases. Nat Ecol Evol. 2022;6: 1601–1616. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36303000&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 37. 37.LaDeau SL, Allan BF, Leisnham PT, Levy MZ. The ecological foundations of transmission potential and vector-borne disease in urban landscapes. Funct Ecol. 2015;29: 889–901. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26549921&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 38. 38.Hoen AG, Hladish TJ, Eggo RM, Lenczner M, Brownstein JS, Meyers LA. Epidemic Wave Dynamics Attributable to Urban Community Structure: A Theoretical Characterization of Disease Transmission in a Large Network. J Med Internet Res. 2015;17: e169. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26156032&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 39. 39.Du Z, Wang L, Shan S, Lam D, Tsang TK, Xiao J, et al. Pandemic fatigue impedes mitigation of COVID-19 in Hong Kong. Proc Natl Acad Sci U S A. 2022;119: e2213313119. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36417445&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 40. 40.MacIntyre CR, Nguyen P-Y, Chughtai AA, Trent M, Gerber B, Steinhofel K, et al. Mask use, risk-mitigation behaviours and pandemic fatigue during the COVID-19 pandemic in five cities in Australia, the UK and USA: A cross-sectional survey. Int J Infect Dis. 2021;106: 199–207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2021.03.056&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33771668&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 41. 41.Haktanir A, Can N, Seki T, Kurnaz MF, Dilmaç B. Do we experience pandemic fatigue? current state, predictors, and prevention. Curr Psychol. 2022;41: 7314–7325. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34690475&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 42. 42.Jalali N, Brustad HK, Frigessi A, MacDonald EA, Meijerink H, Feruglio SL, et al. Increased household transmission and immune escape of the SARS-CoV-2 Omicron compared to Delta variants. Nat Commun. 2022;13: 5706. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-33233-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36175424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 43. 43.Shrestha LB, Foster C, Rawlinson W, Tedla N, Bull RA. Evolution of the SARS-CoV-2 omicron variants BA.1 to BA.5: Implications for immune escape and transmission. Rev Med Virol. 2022;32: e2381. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35856385&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 44. 44.McCallum M, Czudnochowski N, Rosen LE, Zepeda SK, Bowen JE, Walls AC, et al. Structural basis of SARS-CoV-2 Omicron immune evasion and receptor engagement. Science. 2022;375: 864–868. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abn8652&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35076256&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 45. 45.Arnold CRK, Srinivasan S, Rodriguez S, Rydzak N, Herzog CM, Gontu A, et al. A longitudinal study of the impact of university student return to campus on the SARS-CoV-2 seroprevalence among the community members. Sci Rep. 2022;12: 8586. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-022-12499-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35597780&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 46. 46.Yakusheva O, van den Broek-Altenburg E, Brekke G, Atherly A. Lives saved and lost in the first six month of the US COVID-19 pandemic: A retrospective cost-benefit analysis. PLoS One. 2022;17: e0261759. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35061722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 47. 47.Losina E, Leifer V, Millham L, Panella C, Hyle EP, Mohareb AM, et al. College Campuses and COVID-19 Mitigation: Clinical and Economic Value. Ann Intern Med. 2021;174: 472–483. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33347322&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 48. 48.Miles D, Stedman M, Heald A. LIVING WITH COVID-19: BALANCING COSTS AGAINST BENEFITS IN THE FACE OF THE VIRUS. National Institute Economic Review. 2020;253: R60–R76. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/nie.2020.30&link_type=DOI) 49. 49.Doidge S, Doyle J. Australian universities in the age of Covid. Educational Philosophy and Theory. 2022 [cited 2 Feb 2025]. doi:10.1080/00131857.2020.1804343 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00131857.2020.1804343&link_type=DOI) 50. 50.Ahlburg DA. Covid-19 and UK Universities. Polit Q. 2020;91: 649–654. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32836411&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 51. 51.Ciubotariu II, Wilkes RP, Kattoor JJ, Christian EN, Carpi G, Kitchen A. Investigating the rise of Omicron variant through genomic surveillance of SARS-CoV-2 infections in a highly vaccinated university population. Microb Genom. 2024;10: 001194. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=38334271&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 52. 52.Gibas C, Lambirth K, Mittal N, Juel MAI, Barua VB, Roppolo Brazell L, et al. Implementing building-level SARS-CoV-2 wastewater surveillance on a university campus. Sci Total Environ. 2021;782: 146749. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2021.146749&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33838367&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 53. 53.Johnson KE, Pasco R, Woody S, Lachmann M, Johnson-Leon M, Bhavnani D, et al. Optimizing COVID-19 testing strategies on college campuses: Evaluation of the health and economic costs. PLoS Comput Biol. 2023;19: e1011715. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=38134223&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 54. 54.Hekmati A, Luhar M, Krishnamachari B, Matarić M. Simulating COVID-19 classroom transmission on a university campus. Proc Natl Acad Sci U S A. 2022;119: e2116165119. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35609196&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 55. 55.Heesterbeek H, Anderson RM, Andreasen V, Bansal S, De Angelis D, Dye C, et al. Modeling infectious disease dynamics in the complex landscape of global health. Science. 2015;347: aaa4339. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE2OiIzNDcvNjIyNy9hYWE0MzM5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjUvMDMvMTAvMjAyNS4wMy4wNy4yNTMyMzQ5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 56. 56.Pybus OG, Rambaut A. Evolutionary analysis of the dynamics of viral infectious disease. Nat Rev Genet. 2009;10: 540–550. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2583&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19564871&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268038200012&link_type=ISI) 57. 57.Wu Y, Kang L, Guo Z, Liu J, Liu M, Liang W. Incubation Period of COVID-19 Caused by Unique SARS-CoV-2 Strains: A Systematic Review and Meta-analysis. JAMA Netw Open. 2022;5: e2228008. 58. 58.Jiang X, Rayner S, Luo M-H. Does SARS-CoV-2 has a longer incubation period than SARS and MERS? J Med Virol. 2020;92: 476–478. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.25708&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 59. 59.Bergeri I, Whelan MG, Ware H, Subissi L, Nardone A, Lewis HC, et al. Global SARS-CoV-2 seroprevalence from January 2020 to April 2022: A systematic review and meta-analysis of standardized population-based studies. PLoS Med. 2022;19: e1004107. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1004107&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36355774&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 60. 60.Townsend JP, Hassler HB, Lamb AD, Sah P, Alvarez Nishio A, Nguyen C, et al. Seasonality of endemic COVID-19. mBio. 2023;14: e0142623. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37937979&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 61. 61.Townsend JP, Hassler HB, Dornburg A. Optimal Annual COVID-19 Vaccine Boosting Dates Following Previous Booster Vaccination or Breakthrough Infection. Clin Infect Dis. 2024. doi:10.1093/cid/ciae559 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciae559&link_type=DOI) 62. 62.Bansal S, Chowell G, Simonsen L, Vespignani A, Viboud C. Big Data for Infectious Disease Surveillance and Modeling. J Infect Dis. 2016;214: S375–S379. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiw400&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28830113&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 63. 63.Russell KL, Broderick MP, Franklin SE, Blyn LB, Freed NE, Moradi E, et al. Transmission dynamics and prospective environmental sampling of adenovirus in a military recruit setting. J Infect Dis. 2006;194: 877–885. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/507426&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16960774&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240318300003&link_type=ISI) 64. 64.Gardner W, States D, Bagley N. The Coronavirus and the Risks to the Elderly in Long-Term Care. J Aging Soc Policy. 2020;32: 310–315. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/08959420.2020.1750543&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32245346&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 65. 65.Dykgraaf SH, Matenge S, Desborough J, Sturgiss E, Dut G, Roberts L, et al. Protecting Nursing Homes and Long-Term Care Facilities From COVID-19: A Rapid Review of International Evidence. J Am Med Dir Assoc. 2021;22: 1969–1988. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jamda.2021.07.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34428466&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 66. 66.Pater AA, Bosmeny MS, White AA, Sylvain RJ, Eddington SB, Parasrampuria M, et al. High throughput nanopore sequencing of SARS-CoV-2 viral genomes from patient samples. J Biol Methods. 2021;8: e155. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.14440/jbm.2021.360&link_type=DOI) 67. 67.Ferdous J, Weathers T, Bharati Barua V, Stiers E, France A, C Lambirth K, et al. A SARS-CoV-2 Surveillance Sequencing Protocol Optimized for Oxford Nanopore PromethION v1. 2021. doi:10.17504/protocols.io.butbnwin [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.17504/protocols.io.butbnwin&link_type=DOI) 68. 68.Aksamentov I, Roemer C, Hodcroft E, Neher R. Nextclade: clade assignment, mutation calling and quality control for viral genomes. J Open Source Softw. 2021;6: 3773. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21105/joss.03773&link_type=DOI) 69. 69.Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30: 772–780. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/mst010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23329690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317002300004&link_type=ISI) 70. 70.De Maio N, Kalaghatgi P, Turakhia Y, Corbett-Detig R, Minh BQ, Goldman N. Maximum likelihood pandemic-scale phylogenetics. Nat Genet. 2023;55: 746–752. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-023-01368-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37038003&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 71. 71.Miura S, Tamura K, Tao Q, Huuki LA, Kosakovsky Pond SL, Priest J, et al. A new method for inferring timetrees from temporally sampled molecular sequences. PLoS Comput Biol. 2020;16: e1007046. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31951607&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 72. 72.Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35: 1547–1549. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msy096&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29722887&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 73. 73.Tamura K, Tao Q, Kumar S. Theoretical Foundation of the RelTime Method for Estimating Divergence Times from Variable Evolutionary Rates. Mol Biol Evol. 2018;35: 1770–1782. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msy044&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29893954&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 74. 74.Tamura K, Battistuzzi FU, Billing-Ross P, Murillo O, Filipski A, Kumar S. Estimating divergence times in large molecular phylogenies. Proc Natl Acad Sci U S A. 2012;109: 19333–19338. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzQ3LzE5MzMzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjUvMDMvMTAvMjAyNS4wMy4wNy4yNTMyMzQ5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 75. 75.Mello B, Tao Q, Tamura K, Kumar S. Fast and Accurate Estimates of Divergence Times from Big Data. Mol Biol Evol. 2017;34: 45–50. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msw247&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27836983&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 76. 76.Townsend JP, Hassler HB, Wang Z, Miura S, Singh J, Kumar S, et al. The durability of immunity against reinfection by SARS-CoV-2: a comparative evolutionary study. Lancet Microbe. 2021;2: e666–e675. 77. 77.Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4: e88. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.0040088&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16683862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 78. 78.Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009;5: e1000520. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1000520&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19779555&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 79. 79.Townsend JP, Hassler HB, Sah P, Galvani AP, Dornburg A. The durability of natural infection and vaccine-induced immunity against future infection by SARS-CoV-2. Proc Natl Acad Sci U S A. 2022;119: e2204336119. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1073/pnas.2204336119&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35858382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 80. 80.Pagel M, Meade A, Barker D. Bayesian estimation of ancestral character states on phylogenies. Syst Biol. 2004;53: 673–684. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/10635150490522232&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15545248&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224960100001&link_type=ISI) 81. 81.Chazot N, Condamine FL, Dudas G, Peña C, Kodandaramaiah U, Matos-Maraví P, et al. Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies. Nat Commun. 2021;12: 5717. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-25906-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34588433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 82. 82.Cincotta A, Nicolaï M, Campos HBN, McNamara M, D’Alba L, Shawkey MD, et al. Pterosaur melanosomes support signalling functions for early feathers. Nature. 2022;604: 684–688. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-04622-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35444275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 83. 83.Lemey P, Rambaut A, Welch JJ, Suchard MA. Phylogeography takes a relaxed random walk in continuous space and time. Mol Biol Evol. 2010;27: 1877–1885. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msq067&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20203288&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280296700015&link_type=ISI) 84. 84.Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst Biol. 2018;67: 901–904. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/sysbio/syy032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29718447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 85. 85.Gu Z, Gu L, Eils R, Schlesner M, Brors B. *circlize*implements and enhances circular visualization in R. Bioinformatics. 2014;30: 2811–2812. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu393&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24930139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000343082900017&link_type=ISI) 86. 86.Wickham H. Ggplot2. Wiley Interdiscip Rev Comput Stat. 2011;3: 180–185. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/wics.147&link_type=DOI) 87. 87.Yu G. Data integration, manipulation and visualization of phylogenetic trees. Boca Raton: Chapman and Hall/CRC; 2022. 88. 88.Wang L-G, Lam TT-Y, Xu S, Dai Z, Zhou L, Feng T, et al. Treeio: An R Package for Phylogenetic Tree Input and Output with Richly Annotated and Associated Data. Mol Biol Evol. 2020;37: 599–603. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msz240&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31633786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2025%2F03%2F10%2F2025.03.07.25323493.atom) 89. 89.Wickham H, Averick M, Bryan J, Chang W, McGowan L, François R, et al. Welcome to the tidyverse. J Open Source Softw. 2019;4: 1686. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21105/joss.01686&link_type=DOI)