Measuring Spatial Dependence for Infectious Disease Epidemiology.

Global spatial clustering is the tendency of points, here cases of infectious disease, to occur closer together than expected by chance. The extent of global clustering can provide a window into the spatial scale of disease transmission, thereby providing insights into the mechanism of spread, and informing optimal surveillance and control. Here the authors present an interpretable measure of spatial clustering, τ, which can be understood as a measure of relative risk. When biological or temporal information can be used to identify sets of potentially linked and likely unlinked cases, this measure can be estimated without knowledge of the underlying population distribution. The greater our ability to distinguish closely related (i.e., separated by few generations of transmission) from more distantly related cases, the more closely τ will track the true scale of transmission. The authors illustrate this approach using examples from the analyses of HIV, dengue and measles, and provide an R package implementing the methods described. The statistic presented, and measures of global clustering in general, can be powerful tools for analysis of spatially resolved data on infectious diseases

Understanding the drivers of preferential migration of people living with HIV to fishing communities of Lake Victoria in Uganda.

Fishing communities around Lake Victoria have among the highest burdens of HIV globally. Growing evidence suggests that high HIV prevalence is partially due to selective migration of people living with HIV to fishing communities. However, the reasons for this preferential migration are unclear. We recruited 60 men and women for qualitative in-depth interviews (30% living with HIV; 70% recent migrants of unknown HIV status) from seven Ugandan fishing communities. Interviews discussed mobility histories and the social context surrounding migration. Interviews were audio-recorded, transcribed, and translated. A version of the 'Push-Pull' theory of migration helped structure a conceptual thematic framework for data analysis. Unfavourable conditions related primarily to stigma, social discrimination, humiliation, rejection or HIV labelling, and violence, induced individuals to leave their home communities. Factors which eventually resulted in migration to fishing communities included anticipating less HIV-related stigma and a safe, friendly environment that accommodates all people. Access to healthy food (fish) and the perceived availability of community-based HIV care services were also attractions. We found that stigma is the major social phenomenon shaping preferential migration to fishing communities in Uganda

Ethical considerations in global HIV phylogenetic research.

Phylogenetic analysis of pathogens is an increasingly powerful way to reduce the spread of epidemics, including HIV. As a result, phylogenetic approaches are becoming embedded in public health and research programmes, as well as outbreak responses, presenting unique ethical, legal, and social issues that are not adequately addressed by existing bioethics literature. We formed a multidisciplinary working group to explore the ethical issues arising from the design of, conduct in, and use of results from HIV phylogenetic studies, and to propose recommendations to minimise the associated risks to both individuals and groups. We identified eight key ethical domains, within which we highlighted factors that make HIV phylogenetic research unique. In this Review, we endeavoured to provide a framework to assist researchers, public health practitioners, and funding institutions to ensure that HIV phylogenetic studies are designed, done, and disseminated in an ethical manner. Our conclusions also have broader relevance for pathogen phylogenetics

Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in Europe

HIV-1 set-point viral load-the approximately stable value of viraemia in the first years of chronic infection-is a strong predictor of clinical outcome and is highly variable across infected individuals. To better understand HIV-1 pathogenesis and the evolution of the viral population, we must quantify the heritability of set-point viral load, which is the fraction of variation in this phenotype attributable to viral genetic variation. However, current estimates of heritability vary widely, from 6% to 59%. Here we used a dataset of 2,028 seroconverters infected between 1985 and 2013 from 5 European countries (Belgium, Switzerland, France, the Netherlands and the United Kingdom) and estimated the heritability of set-point viral load at 31% (CI 15%-43%). Specifically, heritability was measured using models of character evolution describing how viral load evolves on the phylogeny of whole-genome viral sequences. In contrast to previous studies, (i) we measured viral loads using standardized assays on a sample collected in a strict time window of 6 to 24 months after infection, from which the viral genome was also sequenced; (ii) we compared 2 models of character evolution, the classical "Brownian motion" model and another model ("Ornstein-Uhlenbeck") that includes stabilising selection on viral load; (iii) we controlled for covariates, including age and sex, which may inflate estimates of heritability; and (iv) we developed a goodness of fit test based on the correlation of viral loads in cherries of the phylogenetic tree, showing that both models of character evolution fit the data well. An overall heritability of 31% (CI 15%-43%) is consistent with other studies based on regression of viral load in donor-recipient pairs. Thus, about a third of variation in HIV-1 virulence is attributable to viral genetic variation.Peer reviewe

Phylogenetic estimation of the viral fitness landscape of HIV-1 set-point viral load

Set-point viral load (SPVL), a common measure of human immunodeficiency virus (HIV)-1 virulence, is partially determined by viral genotype. Epidemiological evidence suggests that this viral property has been under stabilising selection, with a typical optimum for the virus between 10$^{4}$ and 10$^{5}$ copies of viral RNA per ml. Here we aimed to detect transmission fitness differences between viruses from individuals with different SPVLs directly from phylogenetic trees inferred from whole-genome sequences. We used the local branching index (LBI) as a proxy for transmission fitness. We found that LBI is more sensitive to differences in infectiousness than to differences in the duration of the infectious state. By analysing subtype-B samples from the Bridging the Evolution and Epidemiology of HIV in Europe project, we inferred a significant positive relationship between SPVL and LBI up to approximately 10$^{5}$ copies/ml, with some evidence for a peak around this value of SPVL. This is evidence of selection against low values of SPVL in HIV-1 subtype-B strains, likely related to lower infectiousness, and perhaps a peak in the transmission fitness in the expected range of SPVL. The less prominent signatures of selection against higher SPVL could be explained by an inherent limit of the method or the deployment of antiretroviral therapy

Phylogenetic estimation of the viral fitness landscape of HIV-1 set-point viral load

Set-point viral load (SPVL), a common measure of human immunodeficiency virus (HIV)-1 virulence, is partially determined by viral genotype. Epidemiological evidence suggests that this viral property has been under stabilising selection, with a typical optimum for the virus between 104 and 105 copies of viral RNA per ml. Here we aimed to detect transmission fitness differences between viruses from individuals with different SPVLs directly from phylogenetic trees inferred from whole-genome sequences. We used the local branching index (LBI) as a proxy for transmission fitness. We found that LBI is more sensitive to differences in infectiousness than to differences in the duration of the infectious state. By analysing subtype-B samples from the Bridging the Evolution and Epidemiology of HIV in Europe project, we inferred a significant positive relationship between SPVL and LBI up to approximately 105 copies/ml, with some evidence for a peak around this value of SPVL. This is evidence of selection against low values of SPVL in HIV-1 subtype-B strains, likely related to lower infectiousness, and perhaps a peak in the transmission fitness in the expected range of SPVL. The less prominent signatures of selection against higher SPVL could be explained by an inherent limit of the method or the deployment of antiretroviral therapy.Peer Reviewe

Quantifying HIV transmission flow between high-prevalence hotspots and surrounding communities: a population-based study in Rakai, Uganda.

BACKGROUND: International and global organisations advocate targeting interventions to areas of high HIV prevalence (ie, hotspots). To better understand the potential benefits of geo-targeted control, we assessed the extent to which HIV hotspots along Lake Victoria sustain transmission in neighbouring populations in south-central Uganda. METHODS: We did a population-based survey in Rakai, Uganda, using data from the Rakai Community Cohort Study. The study surveyed all individuals aged 15-49 years in four high-prevalence Lake Victoria fishing communities and 36 neighbouring inland communities. Viral RNA was deep sequenced from participants infected with HIV who were antiretroviral therapy-naive during the observation period. Phylogenetic analysis was used to infer partial HIV transmission networks, including direction of transmission. Reconstructed networks were interpreted through data for current residence and migration history. HIV transmission flows within and between high-prevalence and low-prevalence areas were quantified adjusting for incomplete sampling of the population. FINDINGS: Between Aug 10, 2011, and Jan 30, 2015, data were collected for the Rakai Community Cohort Study. 25 882 individuals participated, including an estimated 75·7% of the lakeside population and 16·2% of the inland population in the Rakai region of Uganda. 5142 participants were HIV-positive (2703 [13·7%] in inland and 2439 [40·1%] in fishing communities). 3878 (75·4%) people who were HIV-positive did not report antiretroviral therapy use, of whom 2652 (68·4%) had virus deep-sequenced at sufficient quality for phylogenetic analysis. 446 transmission networks were reconstructed, including 293 linked pairs with inferred direction of transmission. Adjusting for incomplete sampling, an estimated 5·7% (95% credibility interval 4·4-7·3) of transmissions occurred within lakeside areas, 89·2% (86·0-91·8) within inland areas, 1·3% (0·6-2·6) from lakeside to inland areas, and 3·7% (2·3-5·8) from inland to lakeside areas. INTERPRETATION: Cross-community HIV transmissions between Lake Victoria hotspots and surrounding inland populations are infrequent and when they occur, virus more commonly flows into rather than out of hotspots. This result suggests that targeted interventions to these hotspots will not alone control the epidemic in inland populations, where most transmissions occur. Thus, geographical targeting of high prevalence areas might not be effective for broader epidemic control depending on underlying epidemic dynamics. FUNDING: The Bill & Melinda Gates Foundation, the National Institute of Allergy and Infectious Diseases, the National Institute of Mental Health, the National Institute of Child Health and Development, the Division of Intramural Research of the National Institute for Allergy and Infectious Diseases, the World Bank, the Doris Duke Charitable Foundation, the Johns Hopkins University Center for AIDS Research, and the President's Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention

The Transient Receptor Potential Ion Channel TRPV6 Is Expressed at Low Levels in Osteoblasts and Has Little Role in Osteoblast Calcium Uptake

Background: TRPV6 ion channels are key mediators of regulated transepithelial absorption of Ca2+ within the small intestine. Trpv6-/- mice were reported to have lower bone density than wild-type littermates and significant disturbances in calcium homeostasis that suggested a role for TRPV6 in osteoblasts during bone formation and mineralization. TRPV6 and molecules related to transepithelial Ca2+ transport have been reported to be expressed at high levels in human and mouse osteoblasts. Results: Transmembrane ion currents in whole cell patch clamped SaOS-2 osteoblasts did not show sensitivity to ruthenium red, an inhibitor of TRPV5/6 ion channels, and 45Ca uptake was not significantly affected by ruthenium red in either SaOS-2 (P = 0.77) or TE-85 (P = 0.69) osteoblastic cells. In contrast, ion currents and 45Ca uptake were both significantly affected in a human bronchial epithelial cell line known to express TRPV6. TRPV6 was expressed at lower levels in osteoblastic cells than has been reported in some literature. In SaOS-2 TRPV6 mRNA was below the assay detection limit; in TE-85 TRPV6 mRNA was detected at 6.90±1.9 × 10−5 relative to B2M. In contrast, TRPV6 was detected at 7.7±3.0 × 10−2 and 2.38±0.28 × 10−4 the level of B2M in human carcinoma-derived cell lines LNCaP and CaCO-2 respectively. In murine primary calvarial osteoblasts TRPV6 was detected at 3.80±0.24 × 10−5 relative to GAPDH, in contrast with 4.3±1.5 × 10−2 relative to GAPDH in murine duodenum. By immunohistochemistry, TRPV6 was expressed mainly in myleocytic cells of the murine bone marrow and was observed only at low levels in murine osteoblasts, osteocytes or growth plate cartilage. Conclusions: TRPV6 is expressed only at low levels in osteoblasts and plays little functional role in osteoblastic calcium uptake

Easy and accurate reconstruction of whole HIV genomes from short-read sequence data with shiver

Studying the evolution of viruses and their molecular epidemiology relies on accurate viral sequence data, so that small differences between similar viruses can be meaningfully interpreted. Despite its higher throughput and more detailed minority variant data, next-generation sequencing has yet to be widely adopted for HIV. The difficulty of accurately reconstructing the consensus sequence of a quasispecies from reads (short fragments of DNA) in the presence of large between-and within-host diversity, including frequent indels, may have presented a barrier. In particular, mapping (aligning) reads to a reference sequence leads to biased loss of information; this bias can distort epidemiological and evolutionary conclusions. De novo assembly avoids this bias by aligning the reads to themselves, producing a set of sequences called contigs. However contigs provide only a partial summary of the reads, misassembly may result in their having an incorrect structure, and no information is available at parts of the genome where contigs could not be assembled. To address these problems we developed the tool shiver to pre-process reads for quality and contamination, then map them to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences. Run with two commands per sample, it can easily be used for large heterogeneous data sets. We used shiver to reconstruct the consensus sequence and minority variant information from paired-end short-read whole-genome data produced with the Illumina platform, for sixty-five existing publicly available samples and fifty new samples. We show the systematic superiority of mapping to shiver's constructed reference compared with mapping the same reads to the closest of 3,249 real references: median values of 13 bases called differently and more accurately, 0 bases called differently and less accurately, and 205 bases of missing sequence recovered. We also successfully applied shiver to whole-genome samples of Hepatitis C Virus and Respiratory Syncytial Virus. shiver is publicly available from https://github.com/ChrisHIV/shiver.Peer reviewe