Next Article in Journal
Isolation of Subtype 3c, 3e and 3f-Like Hepatitis E Virus Strains Stably Replicating to High Viral Loads in an Optimized Cell Culture System
Previous Article in Journal
Attenuation of Bluetongue Virus (BTV) in an in ovo Model Is Related to the Changes of Viral Genetic Diversity of Cell-Culture Passaged BTV
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Limited Sustained Local Transmission of HIV-1 CRF01_AE in New South Wales, Australia

1
The Kirby Institute, The University of New South Wales, Sydney, New South Wales 2052, Australia
2
Department of Infectious Diseases & Microbiology, Royal Prince Alfred Hospital, Camperdown, New South Wales 2050, Australia
3
New South Wales State Reference Laboratory for HIV/AIDS, Darlinghurst, New South Wales 2010, Australia
4
New South Wales Health Pathology, Royal Prince Alfred Hospital, Camperdown, New South Wales 2050, Australia
5
Health Protection New South Wales, New South Wales Health, NSW, North Sydne, New South Wales 2060, Australia
6
Positive Life New South Wales, Surry Hills, New South Wales 2010, Australia
7
ACON Health Ltd., Surry Hills, New South Wales 2010, Australia
8
Centre for Population Health, New South Wales Ministry of Health, North Sydney, New South Wales 2059, Australia
9
Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, Institute of Medical Virology, University of Zurich, 8091 Zurich, Switzerland
10
Department of Clinical Immunology & Allergy, Royal Prince Alfred Hospital, Camperdown, New South Wales 2050, Australia
11
Sydney Medical School, University of Sydney, Sydney, New South Wales 2006, Australia
12
New South Wales Health Pathology-Institute of Clinical Pathology and Medical Research, Westmead Hospital, Westmead, New South Wales 2145, Australia
13
Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, Victoria 3000, Australia
14
Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Life and Environmental Sciences, The University of Sydney, Sydney, New South Wales 2006, Australia
*
Author to whom correspondence should be addressed.
Viruses 2019, 11(5), 482; https://doi.org/10.3390/v11050482
Submission received: 12 April 2019 / Revised: 21 May 2019 / Accepted: 24 May 2019 / Published: 27 May 2019
(This article belongs to the Section Animal Viruses)

Abstract

:
Australia’s response to the human immunodeficiency virus type 1 (HIV-1) pandemic led to effective control of HIV transmission and one of the world’s lowest HIV incidence rates—0.14%. Although there has been a recent decline in new HIV diagnoses in New South Wales (NSW), the most populous state in Australia, there has been a concomitant increase with non-B subtype infections, particularly for the HIV-1 circulating recombinant form CRF01_AE. This aforementioned CRF01_AE sampled in NSW, were combined with those sampled globally to identify NSW-specific viral clades. The population growth of these clades was assessed in two-year period intervals from 2009 to 2017. Overall, 109 NSW-specific clades were identified, most comprising pairs of sequences; however, five large clades comprising ≥10 sequences were also found. Forty-four clades grew over time with one or two sequences added to each in different two-year periods. Importantly, while 10 of these clades have seemingly discontinued, the remaining 34 were still active in 2016/2017. Seven such clades each comprised ≥10 sequences, and are representative of individual sub-epidemics in NSW. Thus, although the majority of new CRF01_AE infections were associated with small clades that rarely establish ongoing chains of local transmission, individual sub-epidemics are present and should be closely monitored.

1. Introduction

The human immunodeficiency virus type 1 (HIV-1) is still a major burden worldwide with approximately 36.9 million individuals infected with HIV globally (WHO). The Joint United Nations Programme on HIV/AIDS (UNAIDS) set the 90/90/90 targets with the aim that by 2020 90% of all people living with HIV-1 know their status, 90% of these will be receiving antiretroviral therapy, and 90% of these will have an undetectable plasma viral load [1]. The latter is of particular importance, as infected people that are virally suppressed do not develop the acquired immunodeficiency syndrome (AIDS) [2], and cannot transmit the virus further while on treatment [3]. Australia has an estimated HIV-1 prevalence of 0.14%, with approximately 1000 new infections annually [4]. New South Wales (NSW) is the most populous state in Australia, and has the highest HIV-1 burden, with approximately 10,100 people living with HIV-1 [5]. Importantly, NSW achieved the 90/90/90 targets in 2016 [5,6], and has experienced the first decline in new infections in ten years [7]. However, this decline was only observed in Australian-born men-who-have-sex-with-men (MSM), and an increase in new infections was observed in overseas-born MSM [6], indicating that there are differences in how different risk groups have been impacted by public health interventions.
HIV-1 circulates globally as different genetically defined subtypes, although the subtype prevalence varies across different countries; subtype C is most prevalent in Africa (and globally), the circulating recombinant form 01_AE (CRF01_AE) is most prevalent in Southeast Asia, and subtype B is most prevalent in the USA, Europe, and Australia [8]. This is of current importance because the HIV-1 epidemic in Australia is experiencing a shift in subtype incidence with a major increase in non-B infections, as an estimated 20% of all new infections in Australia were of non-B origin in 2006, and that increased to 38% in 2016 [9]. CRF01_AE is of particular interest for Australia due to the close proximity to Southeast Asia, with both travel and migration between these regions representing a conduit for HIV-1 importation [10]. Yet, little is known about the geographic origins and epidemiological characteristics of CRF01_AE in Australia.
To better understand the phylodynamics of CRF01_AE in NSW we (i) identified independent transmission clades, (ii) estimated clade growth over time, and (iii) assessed the impact of transmission clade growth on the NSW epidemic.

2. Materials and Methods

2.1. Ethics

This study was approved by the NSW Population and Health Services Research Ethics Committee and the ACON Research Ethics Review Committee (RERC) [AU RED Reference: HREC/15/CIHS/38, Cancer Institute NSW reference number: 2015/08/605]. The HIV/AIDS Legal Centre was consulted for legal advice on data anonymity.

2.2. HIV-1 Sequence Data

Protease and reverse transcriptase sequences sampled between 2004 and 2017 were obtained from the three laboratories performing HIV genotypic drug resistance testing for HIV-1 in NSW; St. Vincent’s HIV State Reference Laboratory, New South Wales Health Pathology-Royal Prince Alfred Hospital, and New South Wales Health Pathology Institute of Clinical Pathology and Medical Research (ICPMR).
The final data set was created by the Center for Health Record Linkage (CHeReL) by linking the individual sequence identifiers obtained from the laboratories performing resistance testing to HIV notification case identifiers from the NSW HIV registry. Patient data always remained anonymized. The earliest available sequence for each patient was used for the analyses to avoid including multiple sequences from the same patient (n = 7000). HIV-1 subtype was determined using the REGA HIV-1 Subtyping Tool Version 3.0 [11].

2.3. Phylogenetic Analysis

A maximum likelihood tree including 83 HIV-1 reference sequences from the REGA tool and the 7000 NSW sequences was estimated using FastTreeMP v2.1.10 [12] via the CIPRES Science Gateway [13], implementing a GTR+Γ model of nucleotide substitution, and accounting for invariable sites, and NNI branch-swapping was implemented for topology search. Sequences forming a monophyletic CRF01_AE group (n = 758) were used for subsequent analyses. All HIV-1 sequences were retrieved from the Los Alamos HIV database (http://www.hiv.lanl.gov/) and used for comparison. NSW protease and reverse transcriptase nucleotide sequences were blasted against the global database, and the best 50 hits for each NSW sequence were retrieved from the global database [14]. Accordingly, a total of 2807 global sequences were used as background sequence data; three HIV-1 subtype G sequences (accession numbers DQ168576 and EU786670) were added as an outgroup to root the phylogeny (Table S1). Pairwise alignment was performed in Mafft, implementing the L-INS-I algorithm [15]. The sequences from the Los Alamos HIV database, as well as the sequences obtained from the three reference laboratories in NSW differ slightly in their length and genome position. As a consequence, the alignment was reduced to the nucleotide positions 2253–3384, based on the reference sequence HXB2 (accession number K03455.1). The alignment was visually inspected in Geneious 11.1.3 (https://www.geneious.com), and sequences that were too short were removed as well as codons associated with drug resistance mutations (according to reference sequence HXB2; accession number K03455.1) [16]. The final alignment comprised 3510 sequences of 1132 nucleotides in length, of which 702 were from NSW. Next, the NSW CRF01_AE sequences were split into four groups according to sampling year, with two new years of sampling added successively: (i) 2004–2009 (n = 149), (ii) 2004–2011 (n = 262), (iii) 2004–2013 (n = 387), (iv) 2004–2015 (n = 542), and (v) 2004–2017 (n = 702), the latter representing the complete data set. A maximum likelihood (ML) tree was estimated, employing FastTreeMP for each sub-data set, and using the same background data as described above. An approximate Shimodaira-Hasegawa test was implemented in PhyML v3.1. to test for node support [17,18]. Phylogenetic trees were visualized in FigTree v1.4.3. The high sampling density in this study and in potential subsequent studies enables the identification of complete transmission networks which would compromise patient privacy. Thus, a random subset of 10% of sequence data (n = 70) from the 702 CRF01_AE samples used in this study are available via NCBI under the accession numbers MK941065-MK941134 (Table S1).

2.4. Clade Identification

Clade identification was performed using R v3.4.3 [19], and a NSW-specific HIV-1 clade was defined as a node on the CRF01_AE phylogeny in which at least 80% of all tips were NSW sequences collected in this study [14,20,21]. ClusterPicker v1.2.5. [22] was used to identify potential active transmission clusters with that where characterized by a node support of 0.9 or higher and a genetic distance of 1.5% or lower [23,24].

3. Results and Discussion

3.1. Increase in non-B HIV-1 Subtype in NSW

Virus subtype was determined and changes in subtype incidence were measured over time for NSW HIV-1 sequences sampled between 2004–2017 (n = 7000). The proportion of non-B subtype sequences increased over time: Non-B subtypes comprised an average of 15% of all samples during 2004–2008, an average of 28% of all samples between 2009–2013, and 42% between 2014–2017 (Figure 1). Notably, CRF01_AE became the most prevalent non-B subtype and second most prevalent subtype overall in NSW, i.e., approximately 17% of all samples between 2014–2017 were CRF01_AE.

3.2. Identification of Large NSW-Specific Clades

Virus sequences were available from 2004 onwards. Although sequence numbers were low in the first three years, they increased steadily between 2007 and 2017, with an average of 62 sequences available for each year after 2007 (Figure 2). These data were combined with a background data set consisting of global sequences to identify NSW-specific independent clades (Figure S1). For the complete data set including all 702 NSW sequences, a total of 109 NSW-specific clades were identified, comprising 463 sequences, representing 66% of all NSW sequences in this study. A total of 50% (n = 54) of these clades were sequence pairs, with no evidence of additional onwards transmission. However, 355 sequences fell within 55 clades characterized by three or more sequences, and were distributed as follows: 19 clades contained three sequences, nine clades contained four sequences, 13 clades containing five sequences, three clades contained six and eight sequences each, two clades contained nine and ten sequences each, and four clades contained more than 10 sequences, i.e., 14, 16, 19, and 68 sequences, respectively. This clade size pattern is typical for HIV-1 phylogenies, with the majority of clades being small (two sequences), and only a handful of larger clades representing sub-epidemics [25,26]. This pattern also strongly supports multiple independent virus entries, i.e., individual sub-epidemics, of CRF01_AE in NSW, which arrived in NSW either via acquisition outside the state due to travel, or migration to NSW.
We next mapped clade duration—that is, the time between the earliest and oldest sequences present in a clade—onto a time scale (Figure 2). Notably, while some small clades with two to four sequences were sampled within a short time frame, other small clades contain sequences sampled up to ten years apart. In fact, of the 100 clades identified, 51 spanned a time frame of three to 12 years, and 44 clades were sampled within one year. Thus, clade number and size appear to be growing over time.
Notably, the background data set contained 361 sequences sampled in Australia; however, only two of these formed a cluster with the NSW sequences, such that the contribution of inter-state transmission within Australia is seemingly limited, which is in concordance with a previous study showing that there is only limited clustering between sequences from different states [27].

3.3. Clade Growth is Slow over Time

To assess the rate of clade growth, NSW sequences were split into four different data subsets according to their year of sampling and combined with the global sequences. The first sub-data set comprised of 149 NSW sequences sampled between 2004–2009. A wider time range was chosen in this case to enable a solid starting point for clade identification. Two additional sampling years were added for each subsequent sub-data set. NSW-specific clades were identified for each sub-data set. Clades where the most recent sequence was sampled before 2016 were classified as temporarily discontinued, although these may include more recent undiagnosed cases (11% of all HIV-1 infected people living in Australia are thought to be undiagnosed [9]). Ten such discontinued clades were found (Figure 3A). Notably, 34 clades were still active in 2017 (ongoing clades) and have been steadily growing over time, with an addition of one or two sequences for each two-year period. This includes larger clades found comprising more than ten sequences, and including the largest clade with 68 sequences, all of which represent independent sub-epidemics. This is consistent with a previous study showing that clades grow slowly and irregularly [28].
Within the starting data set with sequences up to 2009, approximately 42% of the sequences represented singletons, and 58% were part of a clade (Figure 3B). Notably, this was the highest ratio of singleton sequences as the proportion of singleton sequences decreased over time in the following sub-data sets from 38% in the 2011 sub-data set to 24% in the 2017 data set.
In contrast, the proportion of sequences forming a new cluster increased initially from 15% in the 2011 sub-data set to 17% in the 2013 sub-data set, but then decreased again to 14% in the complete data set (Figure 3B). This could simply be associated with an increase in the extent of sampling (HIV-1 testing). However, the number of new HIV notifications actually decreased in NSW from 332 in 2011, to 317 in 2016, and 313 notifications in 2017 [6]. Hence, while existing clades are growing more slowly, new clades are arising, such that public health interventions to stop virus transmission may be impacting on different populations unevenly. Indeed, epidemiological differences have been observed between HIV-1 positive individuals born in Australia, and those born overseas; for example a higher proportion of late diagnoses were observed in overseas born individuals [9]. However, from the 702 NSW CRF01_AE sequences used in this study, only 51% had available meta data including transmission risk and country of birth, so that no correlation between demographic factor and clade size and duration could be performed.

3.4. Potential Active Transmission Networks

Finally, potential direct transmission events were investigated by implementing a distance threshold of 1.5% and a node support of 99%. Sequences fulfilling these parameters form a transmission network, representing either direct virus transmission between individuals (i.e., a transmitter-recipient pair), or virus transmission to two recipients from the same transmitter. Sixteen singleton sequences had a high similarity to one or more of the global sequences, representing a potential transmitter-recipient pair between a sequence from overseas and a sequence sampled in NSW, although the directionality of transmission could not be established. In addition, four NSW-specific clades had a low genetic distance with one or more sequences from the global data set, indicating that transmission between sequences from different countries is not uncommon (Figure S2).
Forty NSW-specific sequence pairs had a low genetic distance and high node support, indicative of a transmission network. Similarly, 13 networks with three sequences, four networks with five sequences, and one network with four, six, and eight sequences, respectively, were also identified. One of these has been active for at least six years, which is indicative of active clade growth, and which could be prevented by focused public health interventions. Notably, in one of the networks all three sequences share the drug resistance mutations K101E, Y181C and G190S, which can lead to reduced susceptibility to the non-nucleoside analog reverse transcriptase inhibitor drug class [16].

4. Conclusions

Our study has revealed the depth of active HIV-1 transmission for the CRF01_AE subtype in NSW, Australia. Importantly, the marked increase in CRF01_AE infections might be associated with an increase of small local transmission networks forming local sub-epidemics within the state. The slow growth in clade size and increase in clade number over time identified here highlights the likely effectiveness of interventions facilitating early diagnosis and contact-tracing.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/5/482/s1, Figure S1: Global phylogeny of CRF 01_AE; Figure S2: Potential active transmission networks; Table S1: Accession numbers.

Author Contributions

Conceptualization, F.D.G., A.N.P., C.C., K.P., J.H., E.C.H. and A.D.K.; Data curation, A.N.P., A.S., A.C., H.S., B.T., C.S. and A.D.K.; Funding acquisition, F.D.G., A.E.G. and A.D.K.; Methodology, F.D.G., A.N.P., N.B., and S.D.; Supervision, A.N.P. and A.D.K.; Visualization, F.D.G.; Writing—original draft, F.D.G., A.N.P., P.K., F.J.L., D.E.D., E.C.H. and A.D.K.

Funding

This work was supported by the Swiss National Science Foundation (P300PA_174462 to F.D.G.), the Practitioner Fellowship, the NHMRC Postgraduate Scholarship (APP1074467 to A.N.P.), and the Australian Laureate Fellowship (FL170100022 to E.C.H.). This analysis was further supported by funding from NHMRC Partnership (APP1092852) and Practitioner Fellowship (APP117907) Grants, the NSW Ministry of Health, and UNSW Sydney. The Kirby Institute receives funding from the Australian Government Department of Health and is affiliated with the Faculty of Medicine, UNSW Sydney. No pharmaceutical industry grants were received for this study.

Acknowledgments

The investigators on the NSW HIV Prevention Partnership Project are Andrew Grulich, Christine Selvey, Rebecca Guy, Garrett Prestage, Iryna Zablotska, Jo Holden, Tim Duck, Craig Cooper, Karen Price, Martin Holt, John de Wit, John Kaldor, Anthony Kelleher, and David Wilson. The project steering committee includes these investigators and Bill Whittaker, Phillip Keen, Prital Patel, Barbara Telfer, Steven Nigro, Cherie Power, and Nick Rose. The late Alan Brotherton, the late Adjunct Levinia Crooks, and the late David A. Cooper were also investigators on the project grant. The authors also want to acknowledge the numerous investigators of contributing studies of this cohort. We also thank Roger D. Kouyos, Huldrych F. Günthard, and Karin J. Metzner for their insightful comments.

Conflicts of Interest

F.J.L. has received educational grants and sponsorship from ViiV Healthcare, Gilead Sciences and MSD, speaker fees from AstraZeneca, MSD and the Australasian Society for HIV, Viral Hepatitis and Sexual Health Medicine (ASHM), and provided consultancy and served on advisory boards for Viiv Healthcare, MSD and Calimmune Australia. All other authors report no conflict of interest.

References

  1. UNAIDS. 90-90-90 An Ambitious Treatment Target to Help End the AIDS Epidemic; UNAIDS: Geneva, Switzerland, 2014. [Google Scholar]
  2. Insight Start Study Group; Lundgren, J.D.; Babiker, A.G.; Gordin, F.; Emery, S.; Grund, B.; Sharma, S.; Avihingsanon, A.; Cooper, D.A.; Fatkenheuer, G.; et al. Initiation of antiretroviral therapy in early asymptomatic HIV infection. N. Engl. J. Med. 2015, 373, 795–807. [Google Scholar] [PubMed]
  3. Vernazza, P.; Hirschel, B.; Bernasconi, E.; Flepp, M. Les personnes séropositives ne souffrant d’aucune autre MST et suivant un traitement antirétroviral efficace ne transmettent pas le VIH par voie sexuelle. Bulletin des Médecins Suisses 2008, 89, 165–169. [Google Scholar]
  4. Kirby Institute. Bloodborne Viral and Sexually Transmissible Infections in Aboriginal and Torres Strait Islander People: Annual Surveillance Report 2018; Kirby Institute, UNSW Sydney: Sydney, Australia, 2018. [Google Scholar]
  5. Keen, P.; Gray, R.T.; Telfer, B.; Guy, R.; Schmidt, H.M.; Whittaker, B.; Holden, J.; Holt, M.; Kelleher, A.; Wilson, D.; et al. The 2016 HIV diagnosis and care cascade in New South Wales, Australia: Meeting the UNAIDS 90-90-90 targets. J. Int. AIDS Soc. 2018, 21, e25109. [Google Scholar] [CrossRef]
  6. NSW Ministry of Health. NSW HIV Strategy 2016–2020 Quarter 4 and Annual 2018 Data Report; Centre for Population Health: North Parramatta, NSW, Australia, 2018.
  7. Grulich, A.E.; Guy, R.; Amin, J.; Jin, F.; Selvey, C.; Holden, J.; Schmidt, H.A.; Zablotska, I.; Price, K.; Whittaker, B.; et al. Population-level effectiveness of rapid, targeted, high-coverage roll-out of HIV pre-exposure prophylaxis in men who have sex with men: The EPIC-NSW prospective cohort study. Lancet HIV 2018, 5, e629–e637. [Google Scholar] [CrossRef]
  8. Hemelaar, J. The origin and diversity of the HIV-1 pandemic. Trends Mol. Med. 2012, 18, 182–192. [Google Scholar] [CrossRef] [PubMed]
  9. Kirby Institute. HIV, Viral Hepatitis and Sexually Transmissible Infections in Australia: Annual Surveillance Report 2018; Kirby Institute, UNSW Sydney: Sydney, Australia, 2018. [Google Scholar]
  10. Herring, B.L.; Ge, Y.C.; Wang, B.; Ratnamohan, M.; Zheng, F.; Cunningham, A.L.; Saksena, N.K.; Dwyer, D.E. Segregation of human immunodeficiency virus type 1 subtypes by risk factor in Australia. J. Clin. Microbiol. 2003, 41, 4600–4604. [Google Scholar] [CrossRef]
  11. REGA HIV-1 Subtyping Tool v3.0. Available online: http://dbpartners.stanford.edu:8080/RegaSubtyping/stanford-hiv/typingtool/ (accessed on 1 February 2019).
  12. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree: Computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 2009, 26, 1641–1650. [Google Scholar] [CrossRef]
  13. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In Proceedings of the Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; pp. 1–8. [Google Scholar]
  14. Turk, T.; Bachmann, N.; Kadelka, C.; Boni, J.; Yerly, S.; Aubert, V.; Klimkait, T.; Battegay, M.; Bernasconi, E.; Calmy, A.; et al. Assessing the danger of self-sustained HIV epidemics in heterosexuals by population based phylogenetic cluster analysis. eLife 2017, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kuraku, S.; Zmasek, C.M.; Nishimura, O.; Katoh, K. aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Res. 2013, 41, W22–W28. [Google Scholar] [CrossRef] [Green Version]
  16. Wensing, A.M.; Calvez, V.; Gunthard, H.F.; Johnson, V.A.; Paredes, R.; Pillay, D.; Shafer, R.W.; Richman, D.D. 2017 update of the drug resistance mutations in HIV-1. Top. Antivir. Med. 2017, 24, 132–133. [Google Scholar]
  17. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef]
  18. Anisimova, M.; Gascuel, O. Approximate likelihood ratio test for branchs: A fast, accurate and powerful alternative. Syst. Biol. 2006, 55, 539–552. [Google Scholar] [CrossRef]
  19. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017. [Google Scholar]
  20. Kusejko, K.; Kadelka, C.; Marzel, A.; Battegay, M.; Bernasconi, E.; Calmy, A.; Cavassini, M.; Hoffmann, M.; Boni, J.; Yerly, S.; et al. Inferring the age difference in HIV transmission pairs by applying phylogenetic methods on the HIV transmission network of the Swiss HIV Cohort Study. Virus Evol. 2018, 4. [Google Scholar] [CrossRef] [PubMed]
  21. Kouyos, R.D.; von Wyl, V.; Yerly, S.; Boni, J.; Taffe, P.; Shah, C.; Burgisser, P.; Klimkait, T.; Weber, R.; Hirschel, B.; et al. Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. J. Infect. Dis. 2010, 201, 1488–1497. [Google Scholar] [CrossRef]
  22. Ragonnet-Cronin, M.; Hodcroft, E.; Hue, S.; Fearnhill, E.; Delpech, V.; Brown, A.J.L.; Lycett, S.; UK HIV Drug Resistance Database. Automated analysis of phylogenetic clusters. BMC Bioinform. 2013, 14, 317. [Google Scholar] [CrossRef]
  23. Junqueira, D.M.; Sibisi, Z.; Wilkinson, E.; de Oliveira, T. Factors influencing HIV-1 phylogenetic clustering. Curr. Opin. HIV AIDS 2019. [Google Scholar] [CrossRef]
  24. Oberle, C.S.; Joos, B.; Rusert, P.; Campbell, N.K.; Beauparlant, D.; Kuster, H.; Weber, J.; Schenkel, C.D.; Scherrer, A.U.; Magnus, C.; et al. Tracing HIV-1 transmission: Envelope traits of HIV-1 transmitter and recipient pairs. Retrovirology 2016, 13, 62. [Google Scholar] [CrossRef]
  25. Vasylyeva, T.I.; du Plessis, L.; Pineda-Pena, A.C.; Kuhnert, D.; Lemey, P.; Vandamme, A.M.; Gomes, P.; Camacho, R.J.; Pybus, O.G.; Abecasis, A.B.; et al. Tracing the impact of public health interventions on HIV-1 transmission in Portugal using molecular epidemiology. J. Infect. Dis. 2019. [Google Scholar] [CrossRef] [PubMed]
  26. Chibo, D.; Birch, C. Increasing diversity of human immunodeficiency virus type 1 subtypes circulating in Australia. AIDS Res. Hum. Retrovir. 2012, 28, 578–583. [Google Scholar] [CrossRef]
  27. Castley, A.; Sawleshwarkar, S.; Varma, R.; Herring, B.; Thapa, K.; Dwyer, D.; Chibo, D.; Nguyen, N.; Hawke, K.; Ratcliff, R.; et al. A national study of the molecular epidemiology of HIV-1 in Australia 2005–2012. PLoS ONE 2017, 12, e0170601. [Google Scholar] [CrossRef] [PubMed]
  28. Chibo, D.; Kaye, M.; Birch, C. HIV transmissions during seroconversion contribute significantly to new infections in men who have sex with men in Australia. AIDS Res. Hum. Retrovir. 2012, 28, 460–464. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Changes in HIV-1 subtypes and circulating recombinant form 01_AE (CRF01_AE)_demographics over time. The proportion of subtypes is mapped across different years. Subtype B = yellow, subtype CRF01_AE = olive, subtype C = dark cyan, other non-recombinant subtypes = steel, and recombinants = blue.
Figure 1. Changes in HIV-1 subtypes and circulating recombinant form 01_AE (CRF01_AE)_demographics over time. The proportion of subtypes is mapped across different years. Subtype B = yellow, subtype CRF01_AE = olive, subtype C = dark cyan, other non-recombinant subtypes = steel, and recombinants = blue.
Viruses 11 00482 g001
Figure 2. Sampling size changes over time. Number of sequences mapped across different years (bottom). New South Wales (NSW)-specific clade duration is shown by plotting the duration from the oldest to the most recent sequence sampled within a clade (top). Clades are colored according to clade size that represents the number of sequences within each clade.
Figure 2. Sampling size changes over time. Number of sequences mapped across different years (bottom). New South Wales (NSW)-specific clade duration is shown by plotting the duration from the oldest to the most recent sequence sampled within a clade (top). Clades are colored according to clade size that represents the number of sequences within each clade.
Viruses 11 00482 g002
Figure 3. Discontinued and active ongoing NSW-specific clades. (A) Clade size for different clades is shown across the five sub-data sets. Discontinued clades are shown in orange with time period for the most recent sequence sampled shown with a red circle, sustained ongoing clades in blue, all other clades in black. Clade size for the analyzed time frame is indicated in the legend. (B) The number of singleton sequences is shown for the five different sub-data sets (dark blue). For each sub-data set new sequences belonging to an existing clade (blue) from the previous sub-data set, and new sequences forming a new clade (grey), are also shown. The x-axis defines the sequence sampling date grouped into the two-year time periods investigated.
Figure 3. Discontinued and active ongoing NSW-specific clades. (A) Clade size for different clades is shown across the five sub-data sets. Discontinued clades are shown in orange with time period for the most recent sequence sampled shown with a red circle, sustained ongoing clades in blue, all other clades in black. Clade size for the analyzed time frame is indicated in the legend. (B) The number of singleton sequences is shown for the five different sub-data sets (dark blue). For each sub-data set new sequences belonging to an existing clade (blue) from the previous sub-data set, and new sequences forming a new clade (grey), are also shown. The x-axis defines the sequence sampling date grouped into the two-year time periods investigated.
Viruses 11 00482 g003

Share and Cite

MDPI and ACS Style

Di Giallonardo, F.; Pinto, A.N.; Keen, P.; Shaik, A.; Carrera, A.; Salem, H.; Telfer, B.; Cooper, C.; Price, K.; Selvey, C.; et al. Limited Sustained Local Transmission of HIV-1 CRF01_AE in New South Wales, Australia. Viruses 2019, 11, 482. https://doi.org/10.3390/v11050482

AMA Style

Di Giallonardo F, Pinto AN, Keen P, Shaik A, Carrera A, Salem H, Telfer B, Cooper C, Price K, Selvey C, et al. Limited Sustained Local Transmission of HIV-1 CRF01_AE in New South Wales, Australia. Viruses. 2019; 11(5):482. https://doi.org/10.3390/v11050482

Chicago/Turabian Style

Di Giallonardo, Francesca, Angie N. Pinto, Phillip Keen, Ansari Shaik, Alex Carrera, Hanan Salem, Barbara Telfer, Craig Cooper, Karen Price, Christine Selvey, and et al. 2019. "Limited Sustained Local Transmission of HIV-1 CRF01_AE in New South Wales, Australia" Viruses 11, no. 5: 482. https://doi.org/10.3390/v11050482

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop