Search

Phylogenetic and discrete-trait evolutionary inference depend heavily on an appropriate characterization of the underlying character substitution process. In this paper, we present random-effects substitution models that extend common continuous-time Markov chain models into a richer class of processes capable of capturing a wider variety of substitution dynamics. As these random-effects substitution models often require many more parameters than their usual counterparts, inference can be both statistically and computationally challenging. Thus, we also propose an efficient approach to compute an approximation to the gradient of the data likelihood with respect to all unknown substitution model parameters. We demonstrate that this approximate gradient enables scaling of sampling-based inference, namely Bayesian inference via Hamiltonian Monte Carlo, under random-effects substitution models across large trees and state-spaces. Applied to a dataset of 583 SARS-CoV-2 sequences, an HKY model with random-effects shows strong signals of nonreversibility in the substitution process, and posterior predictive model checks clearly show that it is a more adequate model than a reversible model. When analyzing the pattern of phylogeographic spread of 1441 influenza A virus (H3N2) sequences between 14 regions, a random-effects phylogeographic substitution model infers that air travel volume adequately predicts almost all dispersal rates. A random-effects state-dependent substitution model reveals no evidence for an effect of arboreality on the swimming mode in the tree frog subfamily Hylinae. Simulations reveal that random-effects substitution models can accommodate both negligible and radical departures from the underlying base substitution model. We show that our gradient-based inference approach is over an order of magnitude more time efficient than conventional approaches.

BackgroundInfection prevention (IP) measures are designed to mitigate the transmission of pathogens in healthcare. Using large-scale viral genomic and social network analyses, we determined if IP measures used during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic were adequate in protecting healthcare workers (HCWs) and patients from acquiring SARS-CoV-2.MethodsWe performed retrospective cross-sectional analyses of viral genomics from all available SARS-CoV-2 viral samples collected at UC San Diego Health and social network analysis using the electronic medical record to derive temporospatial overlap of infections among related viromes and supplemented with contact tracing data. The outcome measure was any instance of healthcare transmission, defined as cases with closely related viral genomes and epidemiological connection within the healthcare setting during the infection window. Between November 2020 through January 2022, 12 933 viral genomes were obtained from 35 666 patients and HCWs.ResultsAmong 5112 SARS-CoV-2 viral samples sequenced from the second and third waves of SARS-CoV-2 (pre-Omicron), 291 pairs were derived from persons with a plausible healthcare overlap. Of these, 34 pairs (12%) were phylogenetically linked: 19 attributable to household and 14 to healthcare transmission. During the Omicron wave, 2106 contact pairs among 7821 sequences resulted in 120 (6%) related pairs among 32 clusters, of which 10 were consistent with healthcare transmission. Transmission was more likely to occur in shared spaces in the older hospital compared with the newer hospital (2.54 vs 0.63 transmission events per 1000 admissions, P < .001).ConclusionsIP strategies were effective at identifying and preventing healthcare SARS-CoV-2 transmission.

ObjectivesThis study investigates primary peer-referral engagement (PRE) strategies to assess which strategy results in engaging higher numbers of people with HIV (PWH) who are virally unsuppressed.DesignWe develop a modeling study that simulates an HIV epidemic (transmission, disease progression, and viral evolution) over 6 years using an agent-based model followed by simulating PRE strategies. We investigate two PRE strategies where referrals are based on social network strategies (SNS) or sexual partner contact tracing (SPCT).MethodsWe parameterize, calibrate, and validate our study using data from Chicago on Black sexual minority men to assess these strategies for a population with high incidence and prevalence of HIV. For each strategy, we calculate the number of PWH recruited who are undiagnosed or out-of-care (OoC) and the number of direct or indirect transmissions.ResultsSNS and SPCT identified 256.5 [95% confidence interval (CI) 234-279] and 15 (95% CI 7-27) PWH, respectively. Of these, SNS identified 159 (95% CI 142-177) PWH OoC and 32 (95% CI 21-43) PWH undiagnosed compared with 9 (95% CI 3-18) and 2 (95% CI 0-5) for SPCT. SNS identified 15.5 (95% CI 6-25) and 7.5 (95% CI 2-11) indirect and direct transmission pairs, whereas SPCT identified 6 (95% CI 0-8) and 5 (95% CI 0-8), respectively.ConclusionWith no testing constraints, SNS is the more effective strategy to identify undiagnosed and OoC PWH. Neither strategy is successful at identifying sufficient indirect or direct transmission pairs to investigate transmission networks.

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained heterosexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

Large datasets along with sampling bias represent a challenge for phylodynamic reconstructions, particularly when the study data are obtained from various heterogeneous sources and/or through convenience sampling. In this study, we evaluate the presence of unbalanced sampled distribution by collection date, location, and risk group of human immunodeficiency virus Type 1 Subtype C using a comprehensive subsampling strategy and assess their impact on the reconstruction of the viral spatial and risk group dynamics using phylogenetic comparative methods. Our study shows that a most suitable dataset for ancestral trait reconstruction can be obtained through subsampling by all available traits, particularly using multigene datasets. We also demonstrate that sampling bias is inflated when considerable information for a given trait is unavailable or of poor quality, as we observed for the trait risk group. In conclusion, we suggest that, even if traits are not well recorded, including them deliberately optimizes the representativeness of the original dataset rather than completely excluding them. Therefore, we advise the inclusion of as many traits as possible with the aid of subsampling approaches in order to optimize the dataset for phylodynamic analysis while reducing the computational burden. This will benefit research communities investigating the evolutionary and spatio-temporal patterns of infectious diseases.

Public health surveillance data used in HIV molecular cluster analyses lack contextual information that is available from partner services (PS) data. Integrating these data sources in retrospective analyses can enrich understanding of the risk profile of people in clusters. In this study, HIV molecular clusters were identified and matched to information on partners and other information gleaned at the time of diagnosis, including coinfection with syphilis. We aimed to produce a more complete understanding of molecular cluster membership in Houston, Texas, a city ranking ninth nationally in rate of new HIV diagnoses that may benefit from retrospective matched analyses between molecular and PS data to inform future intervention. Data from PS were matched to molecular HIV records of people newly diagnosed from 2012 to 2018. By conducting analyses in HIV-TRACE (TRAnsmission Cluster Engine) using viral genetic sequences, molecular clusters were detected. Multivariable logistic regression models were used to estimate the association between molecular cluster membership and completion of a PS interview, number of named partners, and syphilis coinfection. Using data from 4,035 people who had a viral genetic sequence and matched PS records, molecular cluster membership was not significantly associated with completion of a PS interview. Among those with sequences who completed a PS interview (n = 3,869), 45.3% (n = 1,753) clustered. Molecular cluster membership was significantly associated with naming 1 or 3+ partners compared with not naming any partners [adjusted odds ratio, aOR: 1.27 (95% confidence interval, CI: 1.08-1.50), p = .003 and aOR: 1.38 (95% CI: 1.06-1.81), p = .02]. Alone, coinfection with syphilis was not significantly associated with molecular cluster membership. Syphilis coinfection was associated with molecular cluster membership when coupled with incarceration [aOR: 1.91 (95% CI: 1.08-3.38), p = .03], a risk for treatment interruption. Enhanced intervention among those with similar profiles, such as people coinfected with other risks, may be warranted.

Internally displaced persons are often excluded from HIV molecular epidemiology surveillance due to structural, behavioral, and social barriers in access to treatment. We test a field-based molecular epidemiology framework to study HIV transmission dynamics in a hard-to-reach and highly stigmatized group, internally displaced people who inject drugs (IDPWIDs). We inform the framework by Nanopore generated HIV pol sequences and IDPWID migration history. In June-September 2020, we recruited 164 IDPWID in Odesa, Ukraine, and obtained 34 HIV sequences from HIV-infected participants. We aligned them to publicly available sequences (N = 359) from Odesa and IDPWID regions of origin and identified 7 phylogenetic clusters with at least 1 IDPWID. Using times to the most recent common ancestors of the identified clusters and times of IDPWID relocation to Odesa, we infer potential post-displacement transmission window when infections likely to happen to be between 10 and 21 months, not exceeding 4 years. Phylogeographic analysis of the sequence data shows that local people in Odesa disproportionally transmit HIV to the IDPWID community. Rapid transmissions post-displacement in the IDPWID community might be associated with slow progression along the HIV continuum of care: only 63% of IDPWID were aware of their status, 40% of those were in antiviral treatment, and 43% of those were virally suppressed. Such HIV molecular epidemiology investigations are feasible in transient and hard-to-reach communities and can help indicate best times for HIV preventive interventions. Our findings highlight the need to rapidly integrate Ukrainian IDPWID into prevention and treatment services following the dramatic escalation of the war in 2022.

• Human immunodeficiency virus (HIV) drug resistance has implications for antiretroviral treatment strategies and for containing the HIV pandemic because the development of HIV drug resistance leads to the requirement for antiretroviral drugs that may be less effective, less well-tolerated, and more expensive than those used in first-line regimens.  • HIV drug resistance studies are designed to determine which HIV mutations are selected by antiretroviral drugs and, in turn, how these mutations affect antiretroviral drug susceptibility and response to future antiretroviral treatment regimens.  • Such studies collectively form a vital knowledge base essential for monitoring global HIV drug resistance trends, interpreting HIV genotypic tests, and updating HIV treatment guidelines.  • Although HIV drug resistance data are collected in many studies, such data are often not publicly shared, prompting the need to recommend best practices to encourage and standardize HIV drug resistance data sharing.  • In contrast to other viruses, sharing HIV sequences from phylogenetic studies of transmission dynamics requires additional precautions as HIV transmission is criminalized in many countries and regions.  • Our recommendations are designed to ensure that the data that contribute to HIV drug resistance knowledge will be available without undue hardship to those publishing HIV drug resistance studies and without risk to people living with HIV.

BackgroundDue to practical challenges associated with genetic sequencing in low-resource environments, the burden of hepatitis C virus (HCV) in forcibly displaced people is understudied. We examined the use of field applicable HCV sequencing methods and phylogenetic analysis to determine HCV transmission dynamics in internally displaced people who inject drugs (IDPWID) in Ukraine.MethodsIn this cross-sectional study, we used modified respondent-driven sampling to recruit IDPWID who were displaced to Odesa, Ukraine, before 2020. We generated partial and near full length genome (NFLG) HCV sequences using Oxford Nanopore Technology (ONT) MinION in a simulated field environment. Maximum likelihood and Bayesian methods were used to establish phylodynamic relationships.ResultsBetween June and September 2020, we collected epidemiological data and whole blood samples from 164 IDPWID (PNAS Nexus.2023;2(3):pgad008). Rapid testing (Wondfo® One Step HCV; Wondfo® One Step HIV1/2) identified an anti-HCV seroprevalence of 67.7%, and 31.1% of participants tested positive for both anti-HCV and HIV. We generated 57 partial or NFLG HCV sequences and identified eight transmission clusters, of which at least two originated within a year and a half post-displacement.ConclusionsLocally generated genomic data and phylogenetic analysis in rapidly changing low-resource environments, such as those faced by forcibly displaced people, can help inform effective public health strategies. For example, evidence of HCV transmission clusters originating soon after displacement highlights the importance of implementing urgent preventive interventions in ongoing situations of forced displacement.

BackgroundMonitoring the emergence and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants is an important public health objective. We investigated how the Gamma variant was established in New York City (NYC) in early 2021 in the presence of travel restrictions that aimed to prevent viral spread from Brazil, the country where the variant was first identified.MethodsWe performed phylogeographic analysis on 15 967 Gamma sequences sampled between 10 March and 1 May 2021, to identify geographic sources of Gamma lineages introduced into NYC. We identified locally circulating Gamma transmission clusters and inferred the timing of their establishment in NYC.ResultsWe identified 16 phylogenetically distinct Gamma clusters established in NYC (cluster sizes ranged 2-108 genomes); most of them were introduced from Florida and Illinois and only 1 directly from Brazil. By the time the first Gamma case was reported by genomic surveillance in NYC on 10 March, the majority (57%) of circulating Gamma lineages had already been established in the city for at least 2 weeks.ConclusionsAlthough travel from Brazil to the United States was restricted from May 2020 through the end of the study period, this restriction did not prevent Gamma from becoming established in NYC as most introductions occurred from domestic locations.

Understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in 2019 is critical to preventing future zoonotic outbreaks before they become the next pandemic. The Huanan Seafood Wholesale Market in Wuhan, China, was identified as a likely source of cases in early reports, but later this conclusion became controversial. We show here that the earliest known COVID-19 cases from December 2019, including those without reported direct links, were geographically centered on this market. We report that live SARS-CoV-2-susceptible mammals were sold at the market in late 2019 and that within the market, SARS-CoV-2-positive environmental samples were spatially associated with vendors selling live mammals. Although there is insufficient evidence to define upstream events, and exact circumstances remain obscure, our analyses indicate that the emergence of SARS-CoV-2 occurred through the live wildlife trade in China and show that the Huanan market was the epicenter of the COVID-19 pandemic.

Understanding the circumstances that lead to pandemics is important for their prevention. We analyzed the genomic diversity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) early in the coronavirus disease 2019 (COVID-19) pandemic. We show that SARS-CoV-2 genomic diversity before February 2020 likely comprised only two distinct viral lineages, denoted "A" and "B." Phylodynamic rooting methods, coupled with epidemic simulations, reveal that these lineages were the result of at least two separate cross-species transmission events into humans. The first zoonotic transmission likely involved lineage B viruses around 18 November 2019 (23 October to 8 December), and the separate introduction of lineage A likely occurred within weeks of this event. These findings indicate that it is unlikely that SARS-CoV-2 circulated widely in humans before November 2019 and define the narrow window between when SARS-CoV-2 first jumped into humans and when the first cases of COVID-19 were reported. As with other coronaviruses, SARS-CoV-2 emergence likely resulted from multiple zoonotic events.

ObjectivesWe examined whether molecular cluster membership was associated with public health identification of HIV transmission potential among named partners in Chicago.DesignHistorical cohort study.MethodsWe matched and analyzed HIV surveillance and partner services data from HIV diagnoses (2012-2016) prior to implementation of cluster detection and response interventions. We constructed molecular clusters using HIV-TRACE at a pairwise genetic distance threshold of 0.5% and identified clusters exhibiting recent and rapid growth according to the Centers for Disease Control and Prevention definition (three new cases diagnosed in past year). Factors associated with identification of partners with HIV transmission potential were examined using multivariable Poisson regression.ResultsThere were 5208 newly diagnosed index clients over this time period. Average age of index clients in clusters was 28; 47% were Black, 29% Latinx/Hispanic, 6% female and 89% MSM. Of the 537 named partners, 191 (35.6%) were linked to index cases in a cluster and of those 16% were either new diagnoses or viremic. There was no statistically significant difference in the probability of identifying partners with HIV transmission potential among index clients in a rapidly growing cluster versus those not in a cluster [adjusted relative risk 1.82, (0.81-4.06)].ConclusionPartner services that were initiated from index clients in a molecular cluster yielded similar new HIV case finding or identification of those with viremia as did interviews with index clients not in clusters. It remains unclear whether these findings are due to temporal disconnects between diagnoses and cluster identification, unobserved cluster members, or challenges with partner services implementation.

ObjectiveUnintentional drug poisoning and overdose deaths in New York City (NYC) increased 175% between 2010 and 2017, partly due to the transition from noninjectable opioids to heroin injection. This transition has led to concern of a resurgent HIV epidemic among persons who inject drugs (PWID) in NYC. Thus, we sought to characterize HIV transmission dynamics in PWID.DesignGenetic network analysis of HIV-1 public health surveillance data.MethodsWe analyzed HIV diagnoses reported to public health surveillance to determine the trajectory of the HIV epidemic among PWID in NYC, from 1985 through 2019. Genetic distance-based clustering was performed using HIV-TRACE to reconstruct transmission patterns among PWID.ResultsThe majority of the genetic links to PWID diagnosed in the 1980s and 1990s are to other PWID. However, since 2011, there has been a continued decline in new diagnoses among PWID, and genetic links between PWID have become increasingly rare, although links to noninjecting MSM and other people reporting sexual transmission risk have become increasingly common. However, we also find evidence suggestive of a resurgence of genetic links among PWID in 2018-2019. PWID who reported male-male sexual contact were not preferentially genetically linked to PWID over the surveillance period, emphasizing their distinct risk profile from other PWID.ConclusionThese trends suggest a transition from parenteral to sexual transmission among PWID in NYC, suggesting that harm reduction, syringe exchange programs, and legalization of over-the-counter syringe sales in pharmacies have mitigated HIV risk by facilitating well tolerated injection among new PWID.

During the COVID-19 pandemic, genomics and bioinformatics have emerged as essential public health tools. The genomic data acquired using these methods have supported the global health response, facilitated development of testing methods, and allowed timely tracking of novel SARS-CoV-2 variants. Yet the virtually unlimited potential for rapid generation and analysis of genomic data is also coupled with unique technical, scientific, and organizational challenges. Here, we discuss the application of genomic and computational methods for the efficient data driven COVID-19 response, advantages of democratization of viral sequencing around the world, and challenges associated with viral genome data collection and processing.

Detection of viral transmission clusters using molecular epidemiology is critical to the response pillar of the Ending the HIV Epidemic initiative. Here, we studied whether inference with an incomplete dataset would influence the accuracy of the reconstructed molecular transmission network. We analyzed viral sequence data available from ~ 13,000 individuals with diagnosed HIV (2012-2019) from Houston Health Department surveillance data with 53% completeness (n = 6852 individuals with sequences). We extracted random subsamples and compared the resulting reconstructed networks versus the full-size network. Increasing simulated completeness was associated with an increase in the number of detected clusters. We also subsampled based on the network node influence in the transmission of the virus where we measured Expected Force (ExF) for each node in the network. We simulated the removal of nodes with the highest and then lowest ExF from the full dataset and discovered that 4.7% and 60% of priority clusters were detected respectively. These results highlight the non-uniform impact of capturing high influence nodes in identifying transmission clusters. Although increasing sequence reporting completeness is the way to fully detect HIV transmission patterns, reaching high completeness has remained challenging in the real world. Hence, we suggest taking a network science approach to enhance performance of molecular cluster detection, augmented by node influence information.

Human herpes simplex virus 2 (HSV-2) is a ubiquitous, slowly evolving DNA virus. HSV-2 has two primary lineages, one found in West and Central Africa and the other found worldwide. Competing hypotheses have been proposed to explain how HSV-2 migrated out-of-Africa (i)HSV-2 followed human migration out-of-Africa 50-100 thousand years ago, or (ii)HSV-2 migrated via the trans-Atlantic slave trade 150-500 years ago. Limited geographic sampling and lack of molecular clock signal has precluded robust comparison. Here, we analyze newly sequenced HSV-2 genomes from Africa to resolve geography and timing of divergence events within HSV-2. Phylogeographic analysis consistently places the ancestor of worldwide dispersal in East Africa, though molecular clock is too slow to be detected using available data. Rates 4.2 × 10-8-5.6 × 10-8 substitutions/site/year, consistent with previous age estimates, suggest a worldwide dispersal 22-29 thousand years ago. Thus, HSV-2 likely migrated with humans from East Africa and dispersed after the Last Glacial Maximum.

BackgroundPrior to direct-acting antivirals (DAAs), HCV incidence rose among men who have sex with men (MSM) living with HIV infection in Germany despite high hepatitis C virus (HCV) treatment rates. We establish a HCV elimination modeling framework to evaluate whether existing treatment rates can achieve the World Health Organization (WHO) incidence target among MSM living with HIV in Germany.MethodsTo evaluate progress towards HCV elimination in Germany, we adapted a previously published HCV transmission model among MSM living with diagnosed HIV. We modelled HCV incidence and prevalence until 2030 (relative to 2015) under existing treatment and DAA scale-up and explored potential impacts of disruptions in treatment and behavioral risk reduction due to the COVID-19 pandemic.ResultsContinuing current treatment rates will result in stable HCV incidence among MSM living with HIV in Germany between 2015-2030. The WHO HCV incidence target is achievable under DAA scale-up to 100% treatment combined with treatment of those previously diagnosed and untreated (at a rate of 15%/year) and would result in greater reductions with early treatment (3 vs 6 months) reducing incidence from 4.0/100person-years to 0.8/100person-years by 2030. A 12-month disruption to HCV treatment (20% reduction) and risk behaviors (25%,50%,75% reduction) during the COVID-19 pandemic would result in a 15% relative increase in total HCV incidence in 2030 compared to that expected under the status quo.ConclusionsHCV elimination among MSM living with HIV in Germany requires further DAA scale-up among those newly diagnosed combined with efforts to treat those previously diagnosed but untreated. Prospective monitoring will establish whether Germany is on track for HCV microelimination.

BackgroundHIV intervention activities directed toward both those most likely to transmit and their HIV-negative partners have the potential to substantially disrupt HIV transmission. Using HIV sequence data to construct molecular transmission clusters can reveal individuals whose viruses are connected. The utility of various cluster prioritization schemes measuring cluster growth have been demonstrated using surveillance data in New York City and across the United States, by the Centers for Disease Control and Prevention (CDC).MethodsWe examined clustering and cluster growth prioritization schemes using Illinois HIV sequence data that include cases from Chicago, a large urban center with high HIV prevalence, to compare their ability to predict future cluster growth.ResultsWe found that past cluster growth was a far better predictor of future cluster growth than cluster membership alone but found no substantive difference between the schemes used by CDC and the relative cluster growth scheme previously used in New York City (NYC). Focusing on individuals selected simultaneously by both the CDC and the NYC schemes did not provide additional improvements.ConclusionGrowth-based prioritization schemes can easily be automated in HIV surveillance tools and can be used by health departments to identify and respond to clusters where HIV transmission may be actively occurring.

The use of real-time genomic epidemiology has enabled the tracking of the global spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), informing evidence-based public health decision making. Ukraine has experienced four waves of the Coronavirus Disease 2019 (COVID-19) between spring 2020 and spring 2022. However, insufficient capacity for local genetic sequencing limited the potential application of SARS-CoV-2 genomic surveillance for public health response in the country. Herein, we report local sequencing of 103 SARS-CoV-2 genomes from patient samples collected in Kyiv in July-December 2021 using Oxford Nanopore technology. Together with other published Ukrainian SARS-CoV-2 genomes, our data suggest that the third wave of the epidemic in Ukraine (June-December 2021) was dominated by the Delta Variant of Concern (VOC). Our phylogeographic analysis revealed that in summer 2021 Delta VOC was introduced into Ukraine from multiple locations worldwide, with most introductions coming from Central and Eastern European countries. The wide geographic range of Delta introductions coincides with increased volume of travel to Ukraine particularly from locations outside of Europe in summer 2021. This study highlights the need to urgently integrate affordable and easily scaled pathogen sequencing technologies in locations with less developed genomic infrastructure, in order to support local public health decision making.

Recombination is an evolutionary process by which many pathogens generate diversity and acquire novel functions. Although a common occurrence during coronavirus replication, detection of recombination is only feasible when genetically distinct viruses contemporaneously infect the same host. Here, we identify an instance of SARS-CoV-2 superinfection, whereby an individual was infected with two distinct viral variants: Alpha (B.1.1.7) and Epsilon (B.1.429). This superinfection was first noted when an Alpha genome sequence failed to exhibit the classic S gene target failure behavior used to track this variant. Full genome sequencing from four independent extracts reveals that Alpha variant alleles comprise around 75% of the genomes, whereas the Epsilon variant alleles comprise around 20% of the sample. Further investigation reveals the presence of numerous recombinant haplotypes spanning the genome, specifically in the spike, nucleocapsid, and ORF 8 coding regions. These findings support the potential for recombination to reshape SARS-CoV-2 genetic diversity.

Regional connectivity and land travel have been identified as important drivers of SARS-CoV-2 transmission. However, the generalizability of this finding is understudied outside of well-sampled, highly connected regions. In this study, we investigated the relative contributions of regional and intercontinental connectivity to the source-sink dynamics of SARS-CoV-2 for Jordan and the Middle East. By integrating genomic, epidemiological and travel data we show that the source of introductions into Jordan was dynamic across 2020, shifting from intercontinental seeding in the early pandemic to more regional seeding for the travel restrictions period. We show that land travel, particularly freight transport, drove introduction risk during the travel restrictions period. High regional connectivity and land travel also drove Jordan's export risk. Our findings emphasize regional connectedness and land travel as drivers of transmission in the Middle East.

IntroductionMolecular surveillance systems could provide public health benefits to focus strategies to improve the HIV care continuum. Here, we infer the HIV genetic network of Mexico City in 2020, and identify actively growing clusters that could represent relevant targets for intervention.MethodsAll new diagnoses, referrals from other institutions, as well as persons returning to care, enrolling at the largest HIV clinic in Mexico City were invited to participate in the study. The network was inferred from HIV pol sequences, using pairwise genetic distance methods, with a locally hosted, secure version of the HIV-TRACE tool: Seguro HIV-TRACE. Socio-demographic, clinical and behavioural metadata were overlaid across the network to design focused prevention interventions.ResultsA total of 3168 HIV sequences from unique individuals were included. One thousand and one-hundred and fifty (36%) sequences formed 1361 links within 386 transmission clusters in the network. Cluster size varied from 2 to 14 (63% were dyads). After adjustment for covariates, lower age (adjusted odds ratio [aOR]: 0.37, p34 vs.

An important component underlying the disparity in HIV risk between race/ethnic groups is the preferential transmission between individuals in the same group. We sought to quantify transmission between different race/ethnicity groups and measure racial assortativity in HIV transmission networks in major metropolitan areas in the United States. We reconstructed HIV molecular transmission networks from viral sequences collected as part of HIV surveillance in New York City, Los Angeles County, and Cook County, Illinois. We calculated assortativity (the tendency for individuals to link to others with similar characteristics) across the network for three candidate characteristics: transmission risk, age at diagnosis, and race/ethnicity. We then compared assortativity between race/ethnicity groups. Finally, for each race/ethnicity pair, we performed network permutations to test whether the number of links observed differed from that expected if individuals were sorting at random. Transmission networks in all three jurisdictions were more assortative by race/ethnicity than by transmission risk or age at diagnosis. Despite the different race/ethnicity proportions in each metropolitan area and lower proportions of clustering among African Americans than other race/ethnicities, African Americans were the group most likely to have transmission partners of the same race/ethnicity. This high level of assortativity should be considered in the design of HIV intervention and prevention strategies.

Molecular epidemiology (ME) is one tool used to end the HIV epidemic in the United States. We combined clinical and behavioral data with HIV sequence data to identify any overlap in clusters generated from different sequence datasets; to characterize HIV transmission clusters; and to identify correlates of clustering among people living with HIV (PLWH) in Washington, District of Columbia (DC). First, Sanger sequences from DC Cohort participants, a longitudinal HIV study, were combined with next-generation sequences (NGS) from participants in a ME substudy to identify clusters. Next, demographic and self-reported behavioral data from ME substudy participants were used to identify risks of secondary transmission. Finally, we combined NGS from ME substudy participants with Sanger sequences in the DC Molecular HIV Surveillance database to identify clusters. Cluster analyses used HIV-Transmission Cluster Engine to identify linked pairs of sequences (defined as distance ≤1.5%). Twenty-eight clusters of ≥3 sequences (size range: 3-12) representing 108 (3%) participants were identified. None of the five largest clusters (size range: 5-12) included newly diagnosed PLWH. Thirty-four percent of ME substudy participants (n = 213) reported condomless sex during their last sexual encounter and 14% reported a Syphilis diagnosis in the past year. Seven transmission clusters (size range: 2-19) were identified in the final analysis, each containing at least one ME substudy participant. Substudy participants in clusters from the third analysis were present in clusters from the first analysis. Combining HIV sequence, clinical and behavioral data provided insights into HIV transmission that may not be identified using traditional epidemiological methods alone. Specifically, the sexual risk behaviors and STI diagnoses reported in the substudy survey may not have been disclosed during Partner Services activities and the survey data complemented clinical data to fully characterize transmission clusters. These findings can be used to enhance local efforts to interrupt transmission and avert new infections.

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

Since the first reports of a novel severe acute respiratory syndrome (SARS)-like coronavirus in December 2019 in Wuhan, China, there has been intense interest in understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in the human population. Recent debate has coalesced around two competing ideas: a "laboratory escape" scenario and zoonotic emergence. Here, we critically review the current scientific evidence that may help clarify the origin of SARS-CoV-2.

HIV molecular epidemiology estimates the transmission patterns from clustering genetically similar viruses. The process involves connecting genetically similar genotyped viral sequences in the network implying epidemiological transmissions. This technique relies on genotype data which is collected only from HIV diagnosed and in-care populations and leaves many persons with HIV (PWH) who have no access to consistent care out of the tracking process. We use machine learning algorithms to learn the non-linear correlation patterns between patient metadata and transmissions between HIV-positive cases. This enables us to expand the transmission network reconstruction beyond the molecular network. We employed multiple commonly used supervised classification algorithms to analyze the San Diego Primary Infection Resource Consortium (PIRC) cohort dataset, consisting of genotypes and nearly 80 additional non-genetic features. First, we trained classification models to determine genetically unrelated individuals from related ones. Our results show that random forest and decision tree achieved over 80% in accuracy, precision, recall, and F1-score by only using a subset of meta-features including age, birth sex, sexual orientation, race, transmission category, estimated date of infection, and first viral load date besides genetic data. Additionally, both algorithms achieved approximately 80% sensitivity and specificity. The Area Under Curve (AUC) is reported 97% and 94% for random forest and decision tree classifiers respectively. Next, we extended the models to identify clusters of similar viral sequences. Support vector machine demonstrated one order of magnitude improvement in accuracy of assigning the sequences to the correct cluster compared to dummy uniform random classifier. These results confirm that metadata carries important information about the dynamics of HIV transmission as embedded in transmission clusters. Hence, novel computational approaches are needed to apply the non-trivial knowledge collected from inter-individual genetic information to metadata from PWH in order to expand the estimated transmissions. We note that feature extraction alone will not be effective in identifying patterns of transmission and will result in random clustering of the data, but its utilization in conjunction with genetic data and the right algorithm can contribute to the expansion of the reconstructed network beyond individuals with genetic data.

Viruses closely related to human pathogens can reveal the origins of human infectious diseases. Human herpes simplexvirus type 1 (HSV-1) and type 2 (HSV-2) are hypothesized to have arisen via host-virus codivergence and cross-species transmission. We report the discovery of novel herpes simplexviruses during a large-scale screening of fecal samples from wild gorillas, bonobos, and chimpanzees. Phylogenetic analysis indicates that, contrary to expectation, simplexviruses from these African apes are all more closely related to HSV-2 than to HSV-1. Molecular clock-based hypothesis testing suggests the divergence between HSV-1 and the African great ape simplexviruses likely represents a codivergence event between humans and gorillas. The simplexviruses infecting African great apes subsequently experienced multiple cross-species transmission events over the past 3 My, the most recent of which occurred between humans and bonobos around 1 Ma. These findings revise our understanding of the origins of human herpes simplexviruses and suggest that HSV-2 is one of the earliest zoonotic pathogens.

BackgroundClusters of HIV diagnoses in time and space and clusters of genetically linked cases can both serve as alerts for directing prevention and treatment activities. We assessed the interplay between geography and transmission across the Los Angeles County (LAC) HIV genetic transmission network.MethodsDeidentified surveillance data reported for 8186 people with HIV residing in LAC from 2010 through 2016 were used to construct a transmission network using HIV-TRACE. We explored geographic assortativity, the tendency for people to link within the same geographic region; concordant time-space pairs, the proportion of genetically linked pairs from the same geographic region and diagnosis year; and Jaccard coefficient, the overlap between geographical and genetic clusters.ResultsGeography was assortative in the genetic transmission network but less so than either race/ethnicity or transmission risk. Only 18% of individuals were diagnosed in the same year and location as a genetically linked partner. Jaccard analysis revealed that cis-men and younger age at diagnosis had more overlap between genetic clusters and geography; the inverse association was observed for trans-women and Blacks/African Americans.ConclusionsWithin an urban setting with endemic HIV, genetic clustering may serve as a better indicator than time-space clustering to understand HIV transmission patterns and guide public health action.

Understanding when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged is critical to evaluating our current approach to monitoring novel zoonotic pathogens and understanding the failure of early containment and mitigation efforts for COVID-19. We used a coalescent framework to combine retrospective molecular clock inference with forward epidemiological simulations to determine how long SARS-CoV-2 could have circulated before the time of the most recent common ancestor of all sequenced SARS-CoV-2 genomes. Our results define the period between mid-October and mid-November 2019 as the plausible interval when the first case of SARS-CoV-2 emerged in Hubei province, China. By characterizing the likely dynamics of the virus before it was discovered, we show that more than two-thirds of SARS-CoV-2-like zoonotic events would be self-limited, dying out without igniting a pandemic. Our findings highlight the shortcomings of zoonosis surveillance approaches for detecting highly contagious pathogens with moderate mortality rates.

Wide-scale SARS-CoV-2 genome sequencing is critical to tracking viral evolution during the ongoing pandemic. Variants first detected in the United Kingdom, South Africa, and Brazil have spread to multiple countries. We developed the software tool, Variant Database (VDB), for quickly examining the changing landscape of spike mutations. Using VDB, we detected an emerging lineage of SARS-CoV-2 in the New York region that shares mutations with previously reported variants. The most common sets of spike mutations in this lineage (now designated as B.1.526) are L5F, T95I, D253G, E484K or S477N, D614G, and A701V. This lineage was first sequenced in late November 2020 when it represented

Molecular cluster detection analyzes HIV sequences to identify rapid HIV transmission and inform public health responses. We describe changes in the capability to detect molecular clusters and in geographic variation in transmission dynamics. We examined the reporting completeness of HIV-1 polymerase sequences in quarterly National HIV Surveillance System datasets from December 2015 to December 2019. Priority clusters were identified quarterly. To understand populations recently affected by rapid transmission, we described the transmission risk and race/ethnicity of people in clusters first detected in 2018-2019. During December 2015 to December 2019, national sequence completeness increased from 26% to 45%. Of the 1212 people in the 136 clusters first detected in 2018-2019, 69% were men who have sex with men (MSM) and 11% were people who inject drugs (PWID). State-by-state analysis showed substantial variation in transmission risk and racial/ethnic groups in clusters of rapid transmission. HIV sequence reporting has increased nationwide. Molecular cluster analysis identifies rapid transmission in varied populations and identifies emerging patterns of rapid transmission in specific population groups, such as PWID, who, in 2015-2016, comprised only 1% of people in such molecular clusters. These data can guide efforts to focus, tailor, and scale up prevention and care services for these populations.

A disproportionate number of Greenland's Inuit population are chronically infected with Hepatitis B virus (HBV; 5-10%). HBV genotypes B and D are most prevalent in the circumpolar Arctic. Here, we report 39 novel HBV/D sequences from individuals residing in southwestern Greenland. We performed phylodynamic analyses with ancient HBV DNA calibrators to investigate the origin and relationship of these taxa to other HBV sequences. We inferred a substitution rate of 1.4 × 10-5 [95% HPD 8.8 × 10-6, 2.0 × 10-5] and a time to the most recent common ancestor of 629 CE [95% HPD 37-1138 CE]. The Greenland taxa form a sister clade to HBV/D2 sequences, specifically New Caledonian and Indigenous Taiwanese sequences. The Greenland sequences share amino acid signatures with subgenotypes D1 and D2 and ~97% sequence identity. Our results suggest the classification of these novel sequences does not fit within the current nomenclature. Thus, we propose these taxa be considered a novel quasi-subgenotype.

This study introduces an innovative methodological approach to identify potential drivers of structuring HIV-1 transmission clustering patterns between different subpopulations in the culturally and racially/ethnically diverse context of Houston, TX, the largest city in the Southern United States. Using 6332 HIV-1 pol sequences from persons newly diagnosed with HIV during the period 2010-2018, we reconstructed HIV-1 transmission clusters, using the HIV-TRAnsmission Cluster Engine (HIV-TRACE); inferred demographic and risk parameters on HIV-1 transmission dynamics by jointly estimating viral transmission rates across racial/ethnic, age, and transmission risk groups; and modeled the degree of network connectivity by using generalized estimating equations (GEE). Our results indicate that Hispanics/Latinos are most vulnerable to the structure of transmission clusters and serve as a bridge population, acting as recipients of transmissions from Whites (3.0 state changes/year) and from Blacks (2.6 state changes/year) as well as sources of transmissions to Whites (1.8 state changes/year) and to Blacks (1.2 state changes/year). There were high rates of transmission and high network connectivity between younger and older Hispanics/Latinos as well as between younger and older Blacks. Prevention and intervention efforts are needed for transmission clusters that involve younger racial/ethnic minorities, in particular Hispanic/Latino youth, to reduce onward transmission of HIV in Houston.

BackgroundPublic health action combating human immunodeficiency virus (HIV) includes facilitating navigation through the HIV continuum of care: timely diagnosis followed by linkage to care and initiation of antiretroviral therapy to suppress viral replication. Molecular epidemiology can identify rapidly growing HIV genetic transmission clusters. How progression through the care continuum relates to transmission clusters has not been previously characterized.MethodsWe performed a retrospective study on HIV surveillance data from 5226 adult cases in Los Angeles County diagnosed from 2010 through 2014. Genetic transmission clusters were constructed using HIV-TRACE. Cox proportional hazard models were used to estimate the impact of transmission cluster growth on the time intervals between care continuum events. Gamma frailty models incorporated the effect of heterogeneity associated with genetic transmission clusters.ResultsIn contrast to our expectations, there were no differences in time to the care continuum events among individuals in clusters with different growth dynamics. However, upon achieving viral suppression, individuals in high growth clusters were slower to experience viral rebound (hazard ratio 0.83, P = .011) compared with individuals in low growth clusters. Heterogeneity associated with cluster membership in the timing to each event in the care continuum was highly significant (P

In recent years, phylogenetic analysis of HIV sequence data has been used in research studies to investigate transmission patterns between individuals and groups, including analysis of data from HIV prevention clinical trials, in molecular epidemiology, and in public health surveillance programs. Phylogenetic analysis can provide valuable information to inform HIV prevention efforts, but it also has risks, including stigma and marginalization of groups, or potential identification of HIV transmission between individuals. In response to these concerns, an interdisciplinary working group was assembled to address ethical challenges in US-based HIV phylogenetic research. The working group developed recommendations regarding (1) study design; (2) data security, access, and sharing; (3) legal issues; (4) community engagement; and (5) communication and dissemination. The working group also identified areas for future research and scholarship to promote ethical conduct of HIV phylogenetic research.

Accurate understanding of the global spread of emerging viruses is critical for public health responses and for anticipating and preventing future outbreaks. Here we elucidate when, where, and how the earliest sustained severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission networks became established in Europe and North America. Our results suggest that rapid early interventions successfully prevented early introductions of the virus from taking hold in Germany and the United States. Other, later introductions of the virus from China to both Italy and Washington state, United States, founded the earliest sustained European and North America transmission networks. Our analyses demonstrate the effectiveness of public health measures in preventing onward transmission and show that intensive testing and contact tracing could have prevented SARS-CoV-2 outbreaks from becoming established in these regions.

ObjectiveTo identify correlates of incident HIV infection in rapidly growing HIV molecular clusters.DesignPhylogenetic analysis of HIV public health surveillance data.MethodsHigh-priority HIV genetic transmission clusters with evidence of rapid growth in 2012 (i.e. clusters with a pairwise genetic distance ≤0.005 substitutions/site and at least three cases diagnosed in 2012) were identified using HIV-TRACE. Then, we investigated cluster growth, defined as HIV cases diagnosed in the following 5 years that were genetically linked to these clusters. For clusters that grew during the follow-up period, Bayesian molecular clock phylogenetic inference was performed to identify clusters with evidence of incident HIV infection (as opposed to diagnosis of previously infected cases) during this follow-up period.ResultsOf the 116 rapidly growing clusters identified, 73 (63%) had phylogenetic evidence for an incident HIV case during the 5-year follow-up period. Correlates of an incident HIV case arising in clusters included a greater number of diagnosed but virally unsuppressed cases in 2012, a greater number of inferred undiagnosed cases in the cluster in 2012, and a younger time of most recent common ancestor for the cluster.ConclusionThese findings suggest that incident infections in rapidly growing clusters originate equally from diagnosed but unsuppressed cases and undiagnosed infections. These results highlight the importance of promoting retention in care and viral suppression as well as partner notification and other case-finding activities when investigating and intervening on high-priority molecular transmission clusters.

BackgroundTo develop a predictive model to prioritize persons with a transmissible HIV viral load for transmission-reduction interventions.MethodsNew York City (NYC) HIV molecular surveillance data from 2010 to 2013 were used to build a model to predict the probability that the partial pol gene of the virus of a person with a transmissible HIV viral load (>1500 copies/ml) would be genetically similar to that of a person with a new HIV infection (diagnosis at stage 0 or 1 according to the revised Centers for Disease Control and Prevention classification system). Data from 2013 to 2016 were then used to validate the model and compare it with five other selection strategies that can be used to prioritize persons for transmission-reduction interventions.ResultsA total of 10 609 persons living with HIV (PLWH) were included in the development dataset, and 8257 were included in the validation dataset. Among the six selection strategies, the predictive model had the highest area under the receiver operating characteristic curve (AUC) [0.86, 95% confidence interval (CI) 0.84--0.88], followed by the 'Young MSM' (0.79, 95% CI 0.77--0.82), 'MSM with high viral loads' (0.74, 95% CI 0.72--0.76), 'Random sample of MSM' (0.73, 95% CI 0.71--0.76), 'Persons with high viral loads' (0.56, 95% CI 0.54--0.59), and 'Random sample' (0.50, 95% CI 0.48--0.53) strategies.ConclusionsJurisdictions should consider applying predictive modeling to prioritize persons with a transmissible viral load for transmission-reduction interventions and to evaluate its feasibility and effectiveness.

SummaryIn exploring the epidemiology of infectious diseases, networks have been used to reconstruct contacts among individuals and/or populations. Summarizing networks using pathogen metadata (e.g. host species and place of isolation) and a phylogenetic tree is a nascent, alternative approach. In this paper, we introduce a tool for reconstructing transmission networks in arbitrary space from phylogenetic information and metadata. Our goals are to provide a means of deriving new insights and infection control strategies based on the dynamics of the pathogen lineages derived from networks and centrality metrics. We created a web-based application, called StrainHub, in which a user can input a phylogenetic tree based on genetic or other data along with characters derived from metadata using their preferred tree search method. StrainHub generates a transmission network based on character state changes in metadata, such as place or source of isolation, mapped on the phylogenetic tree. The user has the option to calculate centrality metrics on the nodes including betweenness, closeness, degree and a new metric, the source/hub ratio. The outputs include the network with values for metrics on its nodes and the tree with characters reconstructed. All of these results can be exported for further analysis.Availability and implementationstrainhub.io and https://github.com/abschneider/StrainHub.

HIV molecular epidemiology can identify clusters of individuals with elevated rates of HIV transmission. These variable transmission rates are primarily driven by host risk behavior; however, the effect of viral traits on variable transmission rates is poorly understood. Viral load, the concentration of HIV in blood, is a heritable viral trait that influences HIV infectiousness and disease progression. Here, we reconstruct HIV genetic transmission clusters using data from the United States National HIV Surveillance System and report that viruses in clusters, inferred to be frequently transmitted, have higher viral loads at diagnosis. Further, viral load is higher in people in larger clusters and with increased network connectivity, suggesting that HIV in the United States is experiencing natural selection to be more infectious and virulent. We also observe a concurrent increase in viral load at diagnosis over the last decade. This evolutionary trajectory may be slowed by prevention strategies prioritized toward rapidly growing transmission clusters.

Transgender people continue to be at high-risk for HIV acquisition, but little is known about the characteristics of their sexual partners. To address this gap, we examined sociodemographic and sexual characteristics of cisgender men who have sex with men (MSM) on pre-exposure prophylaxis (PrEP) reporting transgender sexual partners.A cohort of 392 MSM in southern California in a randomized clinical trial for PrEP adherence were followed from 2013 to 2016. Multivariable generalized estimating equation and logistic models identified characteristics of MSM reporting transgender sexual partners and PrEP adherence.Only 14 (4%) MSM reported having transgender sexual partners. MSM were more likely to report transgender partners if they were African American, had incident chlamydia, reported injection drug-using sexual partners, or received items for sex. Most associations remained significant in the multivariable model: African American (adjusted odds ratio [AOR] 11.20, P = .01), incident chlamydia (AOR 3.71, P = .04), and receiving items for sex (AOR 5.29, P = .04). There were no significant differences in PrEP adherence between MSM reporting transgender partners and their counterpart.MSM who report transgender sexual partners share characteristics associated with individuals with high HIV prevalence. Identifying this group distinct from larger cohorts of MSM could offer new HIV prevention opportunities for this group of MSM and the transgender community.

Molecular epidemiological analysis of viral pathogens can identify factors associated with increased transmission risk. We investigated the frequency of genetic clustering in a large data set of NS34A, NS5A, and NS5B viral sequences from patients with chronic hepatitis C virus (HCV). Within a subset of patients with longitudinal samples, Receiver Operator Characteristic (ROC) analysis was applied which identified a threshold of 0.02 substitutions/site as most appropriate for clustering. From the 7457 patients with chronic HCV infection included in this analysis, we inferred 256 clusters comprising 541 patients (7.3%). We found that HCV/HIV co-infection, young age, and high HCV viral load were all associated with increased clustering frequency, an indicator of increased transmission risk. In light of previous work on HCV/HIV co-infection in acute HCV cohorts, our results suggest that patients with HCV/HIV co-infection may disproportionately be the source of new HCV infections and treatment efforts should be geared towards viral elimination in this vulnerable population.