Assessing biases in phylodynamic inferences in the presence of super-spreaders

Reconstructing the phylodynamic history and geographic spread of the CRF01_AE-predominant HIV-1 epidemic in the Philippines from PR/RT sequences sampled from 2008 to 2018

The Philippines has had a rapidly growing human immunodeficiency virus (HIV) epidemic with a shift in the prevalent subtype from B to CRF01_AE. However, the phylodynamic history of CRF01_AE in the Philippines has yet to be reconstructed. We conducted a descriptive retrospective study reconstructing the history of HIV-1 CRF01_AE transmissions in the Philippines through molecular epidemiology. Partial polymerase sequences (n = 1144) collected between 2008 and 2018 from three island groups were collated from the Research Institute for Tropical Medicine drug resistance genotyping database. Estimation of the time to the most recent common ancestor (tMRCA), effective reproductive number (Re), effective viral population size (Ne), relative migration rates, and geographic spread of CRF01_AE was performed with BEAST. Re and Ne were compared between CRF01_AE and B. Most CRF01_AE sequences formed a single clade with a tMRCA of June 1996 [95 per cent highest posterior density (HPD): December 1991, October 1999]. An increasing CRF01_AE Ne was observed from the tMRCA to 2013. The CRF01_AE Re reached peaks of 2.46 [95 per cent HPD: 1.76, 3.27] in 2007 and 2.52 [95 per cent HPD: 1.83, 3.34] in 2015. A decrease of CRF01_AE Re occurred in the intervening years of 2007 to 2011, reaching as low as 1.43 [95 per cent HPD: 1.06, 1.90] in 2011, followed by a rebound. The CRF01_AE epidemic most likely started in Luzon and then spread to the other island groups of the country. Both CRF01_AE and Subtype B exhibited similar patterns of Re fluctuation over time. These results characterize the subtype-specific phylodynamic history of the largest CRF01_AE cluster in the Philippines, which contextualizes and may inform past, present, and future public health measures toward controlling the HIV epidemic in the Philippines.

Using multiple sampling strategies to estimate SARS-CoV-2 epidemiological parameters from genomic sequencing data

The choice of viral sequences used in genetic and epidemiological analysis is important as it can induce biases that detract from the value of these rich datasets. This raises questions about how a set of sequences should be chosen for analysis. We provide insights on these largely understudied problems using SARS-CoV-2 genomic sequences from Hong Kong, China, and the Amazonas State, Brazil. We consider multiple sampling schemes which were used to estimate Rt and rt as well as related R0 and date of origin parameters. We find that both Rt and rt are sensitive to changes in sampling whilst R0 and the date of origin are relatively robust. Moreover, we find that analysis using unsampled datasets result in the most biased Rt and rt estimates for both our Hong Kong and Amazonas case studies. We highlight that sampling strategy choices may be an influential yet neglected component of sequencing analysis pipelines.