Characterizing trends in HIV infection among men who have sex with men in Australia by birth cohorts: results from a modified back-projection method
© Wand et al. 2009
Received: 6 April 2009
Accepted: 18 September 2009
Published: 18 September 2009
We set out to estimate historical trends in HIV incidence in Australian men who have sex with men with respect to age at infection and birth cohort.
A modified back-projection technique is applied to data from the HIV/AIDS Surveillance System in Australia, including " newly diagnosed HIV infections ", " newly acquired HIV infections " and " AIDS diagnoses ", to estimate trends in HIV incidence over both calendar time and age at infection.
Our results demonstrate that since 2000, there has been an increase in new HIV infections in Australian men who have sex with men across all age groups. The estimated mean age at infection increased from ~35 years in 2000 to ~37 years in 2007. When the epidemic peaked in the mid 1980s, the majority of the infections (56%) occurred among men aged 30 years and younger; 30% occurred in ages 31 to 40 years; and only ~14% of them were attributed to the group who were older than 40 years of age. In 2007, the proportion of infections occurring in persons 40 years or older doubled to 31% compared to the mid 1980s, while the proportion of infections attributed to the group younger than 30 years of age decreased to 36%.
The distribution of HIV incidence for birth cohorts by infection year suggests that the HIV epidemic continues to affect older homosexual men as much as, if not more than, younger men. The results are useful for evaluating the impact of the epidemic across successive birth cohorts and study trends among the age groups most at risk.
After a steady decline since the mid 1980s, there is now growing evidence that HIV infection has been increasing in parts of the developed world . Despite apparent successes in the past, the most current data  indicate an ~40% increase in new HIV diagnoses in 2007 from the year 2000 among Australian men who have sex with men (MSM).
Increases over time in the epidemic could potentially reflect the risk profile of new generations of MSM as they become sexually active. Therefore, describing the trends in HIV incidence by year of infection for successive birth cohorts can potentially provide a comprehensive picture of the epidemic. Comparing incidence estimates obtained from these sub-populations (e.g., "youngest" versus "oldest") can be fundamental in providing indications of the epidemic trends in terms of demography for informing public health responses.
Australia established an HIV surveillance system in the early 1980s, whereby all new HIV diagnoses are reported. Since 1991, further surveillance has been supplemented by national notification of HIV diagnoses with evidence of newly acquired HIV infection, defined as new HIV diagnoses with either a previous negative HIV test within 12 months, or with evidence of a recent seroconversion illness. Although these data are indicative of trends in the HIV epidemic, they cannot be used directly to estimate HIV incidence (number of new infections per year). Accurate estimates of HIV incidence by population subgroup are required to determine trends in the epidemic and to evaluate the groups most at risk for acquiring HIV.
Methods based on back projection  have historically been used to estimate HIV incidence from AIDS surveillance data based on an assumed probability distribution for the incubation period, that is, the time from HIV infection to AIDS. The availability of effective antiretroviral therapies since 1997 has altered the distribution of the incubation period in ways that are difficult to quantify. The current study used a modified back-projection methodology to estimate the annual number of HIV infections by using the HIV/AIDS surveillance data to reconstruct the HIV epidemic among successive birth cohorts since the onset of the epidemic.
This method allows us to obtain estimates of the past and current incidence trends in a population over each calendar year and with respect to age at infection. Similar approaches have been used to estimate age-specific historical trends in the past [4–7], but the methodology used in this study extends previous methods and makes maximal use of available HIV/AIDS data sources in Australia's surveillance databases, including " newly diagnosed HIV infections ", " newly acquired HIV infections " and " AIDS diagnoses ", to estimate trends in HIV incidence.
Since there is no established statistical model to link HIV incidence to HIV diagnosis with respect to HIV testing patterns, the current methodology assumed that if an individual was infected before, or in, a certain year, it was more likely that this individual sought an HIV diagnostic test at the onset of clinical symptoms. However, as HIV testing became more available and promoted, individuals infected in later years tended to be more likely to seek testing independent of the onset of clinical symptoms. Unlike the similar models used in the literature  where both "time at HIV diagnoses" and "time at AIDS diagnoses" data have almost equal roles, in our approach, "time of HIV infection" is the predominant variable, whereas "time at AIDS diagnosis" data is less influential.
This methodology is applied to Australia's HIV/AIDS National Surveillance data to estimate the number of HIV infections among MSM and to evaluate the impact of the year of birth and age at infection. It is also used to assess whether incremental increases in age at diagnoses  in recent years truly reflect increases in age at infection.
In Australia, HIV transmission is monitored through the notification of cases of newly diagnosed HIV infection, including cases with evidence of newly acquired HIV infection (which is defined as HIV infection with evidence of a prior negative test, a diagnosis of primary HIV infection, or an indeterminate western blot within 12 months of HIV diagnosis). There are potentially three data sources of HIV surveillance data available in each calendar year: new HIV-positive diagnoses by year of diagnosis; newly acquired HIV infection (recent infections among new HIV diagnoses); and new AIDS cases (based on physicians' reporting on diagnoses of clinical events subject to AIDS Case Definition) .
Modified back-projection method by birth cohorts
The back-projection method was originally proposed by Brookmeyer and Gail and used in western countries in the late 1980s and early 1990s to estimate trends in HIV infections based on reported AIDS diagnoses . Later studies modified the back-projection methodology to account for the increased incubation period due to treatment [8, 9].
The method used here differs from similar approaches in the literature in that it does not require data linkage between the HIV and AIDS diagnostic registries. It is based on a parametric formulation of the duration of time between the time of acquisition of HIV infection and the time of earliest diagnosis of HIV infection from enhanced HIV surveillance systems or from laboratory confirmed testing. Many factors may influence this distribution (e.g., awareness of recent exposure due to onset of symptoms related to the HIV progression process, random detection due to accidental circumferences, or frequent testing), but in this study we consider two testing "forces" as described below.
Sub-model 1: HIV testing during asymptomatic infection
for the probability of testing x = 0, 1, 2,... years after infection.
Sub-model 2: HIV testing driven by clinical symptoms at late stage of HIV progression
The Weibull distribution has the property that the hazard increases with time from infection, which intuitively would mirror the risk of progression to HIV-related symptoms in untreated HIV infection.
Overall rate of progression to HIV diagnosis
where , t ≥ t 0 is a mixing function with φ = (π, δ, γ); π represents the proportion of infected individuals who were not tested because of clinical symptoms; δ determines the overall shape of the curves; and γ denotes the rate of increase in infection at time t.
When HIV testing just became available at t = t 0, , as m t (φ) increases with time t at rate γ to a saturation level π = limt → ∞ m t (φ) We assume that there will be a proportion 1 - π of infected individuals who are driven by clinical symptoms to be tested (as specified by Sub-model 2).
The HIV incidence curve was then reconstructed by combining two back-projection estimated HIV incidence curves from AIDS diagnostic data (up to 1994, prior to which effective antiretroviral treatment was not available) and HIV diagnostic data using the combined progression rate distribution.
This model uses three key parameters: the proportion of infected individuals who are "late testers" (from Sub-model 2); the rate of increase in incidence over time; and a shape parameter. The final estimates for HIV incidence were produced using an optimization method based on Nelder-Mead, quasi-Newton and conjugate-gradient algorithms written in R-language.
Sensitivity analyses were performed to assess the robustness of the hypothetical HIV incidence curve for a variety of shape parameters. Figure 1a presents the hazard function h(x | t; φ) as defined by m t (φ) at the values of 0.2, 0.4, 0.6 and 0.8 for the shape parameter.
The properties of the parameters used in this approach were also studied with a series of simulation studies. We predicted the mean numbers of HIV diagnoses over time using a variety of parameter estimates. We used the log-logistic distribution with mean = 10 years and a shape parameter 3.08 for the HIV to AIDS incubation period to predict the mean numbers of new AIDS diagnoses up to the time of the availability of effective antiretroviral therapies in 1997. We set the testing behaviour parameter to 0.5 per year before we predicted the proportions of recent infections.
For each year, we generated 200 random numbers for each time point from a mixed-Poisson distribution around the predicted mean values with an over-dispersion factor to account for the extra-Poisson variation. The resulting 2003 sets of simulated data from all parameter combinations contained: (i) the primary data for HIV diagnoses; (ii) the AIDS diagnostics data to adjust for the ramp-up period; and (iii) the auxiliary data for proportions of recent infections among newly diagnosed cases. Results are presented in Figure 1b.
Application to the Australian birth cohorts
To explore the possible effects of birth cohorts on annual HIV incidence over time, birth years were grouped in 10 five-year birth cohorts; in addition, men born after 1980 or men born before 1940 were grouped together. Changes in the age distribution of HIV infections over time were examined by considering "age at diagnosis" and "age at infection". Age at diagnosis was obtained directly from available data on the date of person's new HIV diagnosis and year of birth.
Following a long-term decline, the annual number of new HIV diagnoses, according to data from Australia's HIV/AIDS surveillance system, has recently been increasing: from 557 cases in 2000 to 774 in 2007. Among newly diagnosed HIV cases, an increasing number were in people who had acquired HIV infection within the previous year. A total of ~21,000 cases of HIV infection were diagnosed by December 2007 among homosexually active men in Australia, which accounted for more than 75% of all reported cases. Sensitivity analyses revealed that the magnitudes and shape of the estimated HIV curves were slightly sensitive to assumptions about the parameters in the models (Figures 1a and 1b). However, the primary qualitative conclusion of the analysis was very robust.
Patterns of "age at diagnosis" and "age at infection"
Consistent with overall patterns, it is estimated that all age groups experienced a sharp decline in HIV incidence from 1985 until the mid 1990s, and then HIV incidence levelled off between 1995 and 2000 (Figure 2). But since then, there has been a steady increase in incidence across all age groups, with the older age groups experiencing the greatest magnitude and earliest nadir in incidence.
Summary of age at infection and age at diagnosis by calendar year
age ± se 1
age ± se 1
31.21 ± 0.17
32.31 ± 0.65
29.80 ± 0.18
32.44 ± 0.68
29.64 ± 0.19
33.86 ± 0.20
29.29 ± 0.22
34.61 ± 0.24
28.34 ± 0.24
34.43 ± 0.25
27.42 ± 0.26
34.03 ± 0.28
27.11 ± 0.28
34.41 ± 0.28
27.23 ± 0.29
35.58 ± 0.32
27.67 ± 0.30
34.65 ± 0.34
28.59 ± 0.31
35.49 ± 0.36
29.92 ± 0.31
36.00 ± 0.37
31.51 ± 0.31
36.19 ± 0.38
33.10 ± 0.32
37.00 ± 0.41
34.12 ± 0.31
37.60 ± 0.46
34.68 ± 0.30
36.90 ± 0.45
35.10 ± 0.29
37.37 ± 0.43
35.72 ± 0.27
36.77 ± 0.42
36.57 ± 0.26
37.61 ± 0.41
37.19 ± 0.25
38.54 ± 0.42
37.39 ± 0.25
38.15 ± 0.40
37.50 ± 0.23
38.52 ± 0.40
37.72 ± 0.19
38.77 ± 0.39
37.81 ± 0.27
38.53 ± 0.38
Overall, it was found that the average age at infection has tended to increase over time since the early 1990s, from ~27 (standard error, se = 0.26) years in 1990, to ~35 (se = 0.29) years in 2000, and ~38 (se = 0.27) years in 2007. Such a trend is suggestive of an increase over time in HIV incidence among older individuals relative to younger individuals. This is reflected by a continual increase in the average age of HIV diagnosis (Table 1).
Although it is expected that the average age at infection would be less than the average age at diagnosis, this difference was pronounced between 1990 and 1999 (approximately six to seven years of difference on average). The difference between the average infection and diagnosis ages is a surrogate indicator for the average time between infection and diagnosis. The difference in average ages of infection and diagnosis has decreased significantly over time, consistent with increases in HIV testing rates . Currently, it is estimated that there is ~1 year difference between average infection and diagnosis ages.
Patterns of new HIV infections by birth cohorts
HIV incidence was also estimated for each birth cohort and plotted as a smoothed function over calendar years in Figure 2 (birth cohorts before 1954 were not shown). The annual incidence among MSM who were born before 1965 peaked around 1985. Men in this age group were highly sexually active at the beginning of the HIV epidemic. During these peak years, people who were born after 1965 were younger than 19 years old and were at low risk.
However, this younger group subsequently contributed to a rising HIV infection trend after 1985. It is clear that the younger age groups experienced delays in the onset of HIV infections (until they were sufficiently sexually active as expected). But younger age groups also experience longer incidence levels; that is, as age decreases, progressively the duration of the peak is longer and, in fact, the first peak in incidence has now been reached for MSM born after 1970. Figures also indicate that the annual HIV incidence of new infections appears to have reached a plateau between 1995 and 2000 in all birth cohorts, but since 2000, HIV incidence has been steadily increasing in all birth cohorts.
Estimates of past and current incidence of HIV are important for assessing the impact of public health prevention strategies and need to be updated regularly to describe the changing face of the epidemic. Particularly, monitoring new HIV infections among various age groups may provide one of the sources of information on the spread of HIV that can contribute to a better understanding of the HIV epidemic at the national level. Results based on age at infection are useful to assess if the average age of infection is increasing over time, whereas results based on birth cohort are important for understanding demographic trends in HIV infections. Together, they provide considerable insight into the epidemic to inform HIV prevention programmes.
To our knowledge, this is the first study to investigate the impact of age and birth cohorts on an HIV epidemic by using advanced back-projection methodology.
As the epidemic has matured, the primary route of HIV transmission in Australia has remained unprotected homosexual contact. Annual HIV diagnoses and the incidence of HIV infection in Australia during the 1990s fell to levels below the peak of the mid 1980s. Although this decline was seen in all age groups, it was more pronounced in homosexual men older than 30 years, with more modest reductions in age groups younger than 30 years .
Using a novel adaptation of a back-projection method, we estimated that the annual incidence of HIV infections among men who have sex with men in Australia has been gradually increasing since the late 1990s, from ~500 new infections per year to 874 new infections in 2007. An increase occurred among all age groups.
These analyses may help in providing greater understanding of the dynamics of the HIV epidemic, based on high-quality surveillance data, and provide reasonably reliable estimates of HIV infection. Our improved methodology has allowed quantitative assessments of the HIV epidemic by birth cohorts, thus providing a sound basis for informing targeted public health policy. This method could also be easily applied to other settings by using the publicly available software written in R-language. Further technical and methodological documents are available upon request email@example.com
The National Centre in HIV Epidemiology and Clinical Research (NCHECR) is funded by the Australian Government's Department of Health and Ageing, and is affiliated with the Faculty of Medicine, University of New South Wales, Sydney, NSW. Its work is overseen by the Ministerial Advisory Committee on AIDS, Sexual Health and Hepatitis. The NCHECR Surveillance Program is a collaborating unit of the Australian Institute of Health and Welfare.
- Grulich AE, Kaldor JM, Grulich AE, Kaldor JM: Trends in HIV incidence in homosexual men in developed countries. Sexual Health 2008, 5: 113–118.View ArticlePubMedGoogle Scholar
- NCHECR: HIV/AIDS, viral hepatitis and sexually transmissible infections in Australia Annual Surveillance Report. Sydney, NSW: National Centre in HIV Epidemiology and Clinical Research (NCHECR), The University of New South Wales; 2008.Google Scholar
- Brookmeyer R, Gail MH: Minimum size of the acquired immunodeficiency syndrome (AIDS) epidemic in the United States. Lancet 1986, 2: 1320–1322.View ArticlePubMedGoogle Scholar
- Becker NG, Watson LF, Carlin JB: A method of non-parametric back-projection and its application to AIDS data. Statistics in Medicine 1991, 10: 1527–1542.View ArticlePubMedGoogle Scholar
- Becker NG, Watson LF, Marschner IC, Motika M, Newstead SV, Carlin JB: Assessing the extent of the Australian HIV epidemic from AIDS surveillance data. Australian Journal of Public Health 1993, 17: 226–231.View ArticlePubMedGoogle Scholar
- Rosenberg PS: Backcalculation models of age-specific HIV incidence rates. Statistics in Medicine 1994, 13: 1975–1990.View ArticlePubMedGoogle Scholar
- Law MG, Rosenberg PS, McDonald A, Kaldor JM: Age-specific HIV incidence among homosexually active men in Australia. Med J Aust 1996, 164: 715–718.PubMedGoogle Scholar
- Muñoz A, Hoover DR: Use of cohort studies for evaluating AIDS therapies. In AIDS Clinical Trials. New York: Wiley; 1995:423–446.Google Scholar
- Cui J, Becker NG: Estimating HIV incidence using dates of both HIV and AIDS diagnoses. Statistics in Medicine 2000, 19: 1165–1177.View ArticlePubMedGoogle Scholar
- Brookmeyer R, Gail MH: Generalized back-calculation: extension to account for nonstationary incubation distributions. In AIDS Epidemiology: A Quantitive Approach. New York: Oxford University Press; 1994:219.Google Scholar
- NSW, VIC and QLD Gay Periodic Surveys Sydney, Australia: National Centre in HIV Social Research, University of New South Wales; 1999.Google Scholar
- Law MG, Rosenberg PS, McDonald A, Kaldor JM: Age-specific HIV incidence among homosexually active men in Australia. Medical Journal of Australia 1996, 164: 715–718.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.