Skip to main content

Advertisement

Trends in prevalence of selected opportunistic infections associated with HIV/AIDS in Uganda

BMC Infectious Diseasesvolume 15, Article number: 187 (2015) | Download Citation

Abstract

Background

After more than a decade of establishing and expanding access to highly active antiretroviral therapy (HAART), empirical evidence on its impact on trends of opportunistic infections (OIs) associated with the deadly human immunodeficiency virus (HIV) in resource poor settings is scarce. The primary objective of this study was to assess the effect of HAART coverage on trends of five most common OIs in Uganda.

Methods

Observational data from January 2002 to December 2013 for 5972 HIV positive individuals attending the AIDS Support Organisation (TASO) HIV/AIDS care programme in Uganda were extracted and analysed. Trends were analysed using autoregressive moving average time series and mixed effects linear regression models adjusting for all available potential confounders.

Results

A total of 204,871 monthly medical reports were retrieved and analysed. Majority of the participants were female (73%) with a median age of 32 years (inter-quartile range 26–39). Overall, significant decreasing mean annual prevalence trends were observed for Mycobacterium tuberculosis, herpes zoster, genital ulcer and oral candidiasis (p < 0.05, X2trend). Non-significant declining mean annual prevalence trend was observed for cryptococcal meningitis (p = 0.181, X2trend). The largest impact of HAART was observed in Oral candidiasis and TB whose average annual prevalence reduced by 61% and 43% respectively following the introduction of HAART. Monthly series for TB, Herpes zoster and genital ulcers differed significantly by age and clinic but only genital ulcer series differed significantly by sex (p < 0.05, kruskal wallis). After controlling for the effects of age, sex and clinic (fixed) and monthly clustering (random effect) in a mixed effects linear regression model, all the five OIs showed a significant monthly change in prevalence (p < 0.001).

Conclusion

Overall, prevalence of most OIs declined especially after the introduction of HAART. However significant variations exist in the trends of different OIs in different geographical areas in Uganda. It is therefore important that site specific factors are properly identified to enable the development of targeted interventions.

Open Peer Review reports

Background

Since the outbreak of HIV in 1981, an estimated 39 million people worldwide have died and about 35 million are living with the deadly virus with Sub-Saharan Africa suffering the greatest brunt of the epidemic [1]. Opportunistic infections (OIs) remain the single main cause of ill-health and death among HIV-infected patients [2-4]. Research shows that about 90% of HIV-related morbidity and mortality are caused by opportunistic infections compared to 7% due to opportunistic cancers and 3% due to other causes [5]. However, this may have changed since the introduction of HAART in mid-1990s in developed countries [6-10]. HAART is known for effective suppression of systemic HIV viral load and immune restoration thereby reducing the frequency of opportunistic infections, deferring morbidity and mortality hence improving survival among HIV infected individuals [7,11,12].

Several developing countries are slowly scaling up access to HAART, amidst scarcity of resources and uncertainty for a sustained lifelong provision of treatment to an ever increasing number of eligible HIV patients [1]. By end of 2013, about 13million HIV patients had access to HAART globally with 9.2million from middle and low income countries [13]. In resource poor settings, HIV positive individuals usually access care and treatment with marked immune suppression associated with a higher risk of OIs whose spectrum and frequencies may vary over time and in different countries or even within the same country [3,14]. OIs lower the quality of life of persons living with HIV/AIDS (PLHA), increases stigma and limits one’s ability to work and are usually associated with high medical care costs. Opportunistic infections therefore have greatly contributed to poverty among those infected and affected by HIV hence may be an impediment to the attainment of the millennium development goals (MDGs) on health and poverty eradication in resource poor settings.

Previous studies in developed countries show varied results on the effect of HAART on opportunistic infections over time and in different geographical areas [11,15-17]. For example a study in the USA that evaluated annual trends for 13 most common AIDS-defining opportunistic infections by examining medical records in more than 90 hospitals and clinics in 9 US cities before HAART (1991–96) showed decreasing trends in 5 OIs (PCP, esophageal candidiasis, tuberculosis, herpes simplex and cryptosporidiosis) and an increasing trend in recurrent pneumonia [11]. The trends in the time of onset, spectrum and frequency of infections was found to be unique for different OIs and varied by level of immune-suppression [11].

In another study in the USA, opportunistic infection rates varied considerably among US-born, Mexican-born and central American-born Latinos in the era of HAART. U.S.-born Latino women were more likely than Central American born Latino women to develop an OI (hazard ratio = 2.9, CI: 1.3, 6.5). In a Poisson regression analysis, U.S.-born Latino men and women combined were at greater risk of Kaposi’s sarcoma (RR 2.9, 95% CIs: 1.1, 7.6, p = 0.03) and yet for esophageal candidiasis, there was no evidence of a change in rate between the three communities [18].

Another study in the USA that reviewed the trends in the epidemiology of opportunistic fungal infections associated with HIV/AIDS reported an increasing trend in the incidence of fungal infections (Aspergillus sp) attributed to increasing resistance to anti-fungal treatment and recommended documentation of epidemiological trends to gain more insights into the effectiveness of treatment strategies [19]. A study in Spain that examined temporal trends in the incidence of opportunistic infections (OIs) associated with AIDS in the period 1989–1997 showed a significant decreasing trend in esophageal candidiasis, pulmonary and extra-pulmonary tuberculosis, and cerebral toxoplasmosis. However an increasing but non-significant trend of MAC incidence over time [20]. In Brazil, trends of OIs among HIV-infected adults (>12 years old) at a national level in the period between 1980 and 1999 declined significantly but not homogeneously by regions, risk groups, education and sex [21]. In Northeast and Central-West regions, they showed an increasing trend for TB and Toxoplasmosis. TB had higher incidence among those with lower education (>8 years) while PCP and KS had higher incidence among those with 8+ years of education, despite having similar trends of decline.

Few studies have examined the effect of HAART on trends of OIs in the African settings. One of these few studies examined the trends in incidence rates of TB in HIV-negative adults in South Africa before HAART and reported an increasing trend of TB incidence (test for trend P = 0.17) [22]. Given the scale and speed at which HAART roll-out is taking place in sub-Saharan Africa and in view of the policy changes in HAART access over time [23], it is important to evaluate the impact of these costly HIV interventions in resource poor settings. Uganda is among African countries with the highest burden of HIV/AIDS [24]. Currently there are an estimated 1.5 million people living with HIV/AIDS with about 60% in need of HAART [13]. HAART rollout began in 2004 as part of the Global HAART roll out strategy and the national target is to provide HAART to 80% of the population in need by 2015 [25]. However, the impact of the increasing coverage of HAART on trends of OIs in the country has previously not been examined due to lack of sufficient data. The primary objective of this study was to evaluate the impact of HAART coverage on trends of five most common OIs over time and in different geographical areas in Uganda.

Methods

Study design

Observational data for 5972 HIV positive individuals were obtained from the AIDS support organisation (TASO) in Uganda and analysed. The period of study was from January 2002 to December 2013 including 2 years before and 10 years after HAART introduction in Uganda. The OIs of interest were cryptococcosis, tuberculosis, oral candida, herpes zoster and genital ulcer. These were chosen on the basis of being the commonest OIs among HIV+ individuals in Uganda and are easily diagnosed. The review period was categorised into three mutually exclusive time periods that signify important milestones in the prevention and treatment of opportunistic infections associated with HIV/AIDS in Uganda. First period (2 years) was from January 2002 to December 2003 in which HIV+ individuals had no access to HAART. Second period (5 years) was from January 2004 to December 2008 denoted as “Early HAART” in which HAART access was limited to severely ill HIV+ individuals (CD4 count <200 cells/μl of blood) as per World Health organisation 2006 HAART access guidelines [26]. Third period (5 years) was from January 2009 to December 2013 denoted as “Late HAART” in which a lot of experience in HAART administration had been gained and HAART access policy changed to include all HIV+ individuals with CD4 cell count ≤350 cells/ μl of blood as per World Health organisation 2010 HAART access guidelines [27].

Settings

The study was conducted in TASO, one of the oldest and largest HIV/AIDS care and treatment non-governmental organisation (NGO) in Uganda and sub-Saharan Africa [28]. TASO has 11 HIV clinics nationally recognised as centres of excellence with excellent laboratory practice and management supported by Centre for Disease Prevention and Control (CDC) [29]. With close to 3 decades of experience in HIV/AIDS care and treatment, TASO provided a good opportunity for assessing OI trends in real life programmatic settings in which HAART roll out is taking place in sub-Saharan Africa [30]. TASO ART programme started as part of the National ART roll-out programme in public health facilities in 2004. Being one of the largest ART providers in the country, TASO got support from different funders supporting ART programmes in sub-saharan Africa including the President’s Emergency Plan for AIDS Relief (PEPFAR) and the Global Fund to Fight AIDS, Tuberculosis and Malaria. Initially, access to HAART was based on the Ugandan ministry of Health and the World Health Organisation (WHO) 2006 guide lines i.e. WHO stage 3 or 4 illness or a CD4 cell count < 200 cells/μl for adults and adolescents and WHO stage III, advanced stage II or stage I with CD4 cell percentage less than 20% for those more than 18 months of age [26]. However, after 2010 access to HAART was expanded to include HIV patients with CD4 cell count ≤ 350 or WHO clinical stage3 or 4 irrespective of CD4 cell count as was recommended by Ugandan Ministry of Health and WHO [27,31]. Clients who are not eligible for HAART are usually offered cotrimoxazole or dapson prophylaxis. By end of 2013, a total of 91,218 were actively in care of which 58,051 clients (64%) were on HAART and 87,903 (96.4%) were on cotrimoxazole/dapson prophylaxis throughout the 11 HIV clinics across Uganda [29]. All services are free of charge including anti-retroviral drugs (ARVs) for those who are eligible [30]. Though TASO data system was not set up to answer any particular research question, but with assistance from Centre for Prevention and Disease Control (CDC) and other partners, TASO has been able to put up a large and up to date HIV/AIDS data base that can be used to inform and guide national policy on HIV/AIDS in Uganda [28].

Sampling and sample size

Three TASO HIV clinics were purposively selected basing on volume and quality of data necessary for analyzing trends of most common OIs as well as less common OIs and geographical representation. The HIV clinics selected were TASO Mulago HIV clinic in central Uganda, TASO Mbarara HIV clinic in south-western Uganda and Tororo HIV clinic in Eastern Uganda. All clients irrespective of age who were in care at the three selected HIV clinics by 1st January 2002 were included in the study.

Data collection

Data on OIs were obtained from the TASO electronic data bases at the three HIV clinics and TASO headquarters in Kampala. The data were collected by TASO medical staff following an established protocol for all TASO HIV clinics. In brief, each client was expected to attend the clinic at least once a month. At each clinic visit, data per client was collected on a standardized medical form detailing the client’s demographic information, clinical condition, medical history, OI diagnosis, ART use and level of adherence, prophylaxis use, any other treatment given and side effects/toxicities if any. Diagnosis of OIs was based on the Ugandan ministry of health guidelines: Sputum tuberculin positive test or a chest radiography and/or GeneXpert positive test for TB; A CRAG positive test for Cryptococal meningitis; persistent creamy white curd-like plagues or red patches on the tongue, palate or lining of the mouth for oral candida; progressive and painful genital ulceration for genital ulcers and painful vesicular skin blisters with a dermatomal distribution for herpes zoster. All medical data were compiled and entered into the TASO electronic data base by TASO data administrator using EPIINFO vs3 in Access format. Monthly medical data for each participant covering the period January 2002 to December 2003 were extracted by the data administrator, delinked from overt identifiers and then handed over for analysis.

Data management and analysis

Data management and analysis was done using STATA 12 (Stata Corp, Collage station, Texas, USA). Data were exported from epiinfo/access format to STATA 12 and cleaned. Data cleaning involved data reduction, deleting duplicates and missing data checks. To analyse trends, data on each OI was re-assembled and summarised by year and month. Then monthly and average annual prevalence for each OI were plotted. Monthly prevalence was calculated from the total number of clients recorded with an event in a month divided by the total number that attended the clinic in that month. Annual prevalence per OI was calculated from the mean of the monthly prevalences in a year.

To establish whether there is a trend, prevalence of each OI by month and calendar year were first tabulated and then a monthly time plot for prevalence of each OI was drawn. To filter out short-term fluctuations and random variation within monthly trends, we used the Box Jenkins moving average smoothing technique [32]. The moving average smoothing technique achieves this by replacing each element of the time series by n neighbouring elements, where n is the width of the smoothing window. We used a centred moving average including 3 observations before and 2 observations after the current observation inclusive, and then generated smoothed 144 monthly prevalence series per OI.

Annualized prevalence estimates were computed for each OI and X2 test for linear trend used to test for the significance of the trends. To measure the rate of change in monthly prevalence (b-coefficient) for each OI, we used mixed effects linear regression and modelled between monthly variability as a random effect. To test for the significance of the difference between trends, we used kruskal wallis test [33]. A scatter plot and histogram of the residuals were used to assess for assumption underlying linear regression (i.e. linearity, normal distribution and homogeneity of the variance). Monthly prevalence scatter plots that were non-linear (curvilinear) were first log-transformed before modelling the monthly rate of change in mixed effects linear regression.

Ethical considerations

The study obtained ethical approval from Makerere University School of Public Health Higher Degrees Research and Ethics committee and the Uganda National Council for Science and Technology. Informed consent from study participants was not required as this was routinely collected operational data and the above ethical committees waived the need for consent. However, written consent was obtained from TASO for conducting study and publication of findings with any accompanying images.

Results

Baseline characteristics

A total of 204,871 monthly medical reports from a cohort of 5972 HIV positive patients were retrieved and analysed. 73% (4301) were female with a median age of 32 years (inter-quartile range 26–39). More than half of the study participants were rural poor peasants (56%) with majority having only primary level education (50%) and 30% were widowed (Table 1) . By 2004, only 9% of these had access to HAART but by 2013, all (100%) who were still in care were on HAART. The overall median age at ART initiation was 44 years (inter-quartile range 37–50) but was significantly lower for women at 43 (inter-quartile range 37–49) than that of men at 45 (inter-quartile range 40–52) (p = 0.001). The median CD4 cell count at the start of HAART was 128 cells (inter-quartile range 55–190) and the most commonly used HAART regimes were zidovudine plus lamivudine and nevirapive (ZDV + 3TC+ NVP) (44%), followed by stavudine, lamivudine and nevirapine (d4T + 3TC + NVP) (40%), and other regimen combinations that includes efavirenz (EFV), Tenofovir (TDF) and Lopinavir/ritonavir, etc. were rarely used (20%) (Table 2).

Table 1 Baseline characteristics of the cohort at the start of the study, total and clinic-specific
Full size table
Table 2 Mean number of study participants who were in care and those on ART each year
Full size table

Pre-ART (2002–2003) retention was 79% (4725/5972) after 24 months of follow up. Post HAART (2004–2013) mean annual retention was 90% (2237/2473) with an average of 10% (236/2473) loss to follow up per year. In sensitivity analysis, the female gender of those who remained in care did not differ significantly from those who attended but did not return (X2 = 0.475, P = 0.491). The median age only differed by 1 year (X2 = 19, P < 0.001), with those who remained in care having median age of 32 yrs (inter-quartile rage 27–39) while those who were lost to follow up had a median age of 31 (inter-quartile rage 26–38). So we assumed cases that remained in care were a representative sample of all the cases and that missing clinic visits happened randomly.

Prevalence and trends of OIs

Overall monthly trends for each OI are shown in Figure 1. Best-fit regression lines, R2 and p-values for each OI are shown in Figure 2. Baseline characteristics at the commencement of HAART are shown in Table 2. Numbers of study participants accessing HAART over time are shown in Table 3. Average annual prevalence by each OI and statistical tests for trend are shown in Table 4. Average monthly changes in OI prevalence over time after controlling for age, sex and clinic as fixed effects and monthly clustering as a random effect are shown in Table 5. Figures 2, 3 and 4 show monthly prevalence trend comparisons by sex, age and HIV clinic. Genital ulcer and Herpes zoster monthly prevalence decreased linearly with time, whereas TB and oral candida monthly prevalence trends exhibited a curvilinear decrease with time, decreasing rapidly between 0-48months and thereafter levelled off in the following months.

Figure 1
figure1

Monthly prevalence trend for oral candida (orange line), Mycobacterium tuberculosis (redline), genital ulcer (black line), herpes zoster (green line) and cryptococcal meningitis (blue line) expressed as a proportion of HIV-positive persons diagnosed with a particular OI out of the total number who turned up for care per month for the period January 2002 to December 2013.

Full size image
Figure 2
figure2

Scatter plot and fitted regression line for monthly prevalence over time (in months) for each OI expressed as a proportion of HIV-positive persons diagnosed with a particular OI out of the total number who turned up for care per month for the period January 2002 to December 2013.

Full size image
Table 3 Background characteristics of the study participants at commencement of HAART, total and sex segregated
Full size table
Table 4 Mean annual OI prevalence (per 1000) in a cohort of HIV positive individuals before and after HAART in Uganda
Full size table
Table 5 Mixed effects linear regression analysis of monthly rate of change in the prevalence of each OI (per 1000) adjusted for fixed effects (age, sex and clinic) and random effects (monthly clustering)
Full size table
Figure 3
figure3

Monthly prevalence trend for each OI by sex: female (blue line), male (maroon line); expressed as a proportion of HIV-positive persons diagnosed with a particular OI out of the total number who turned up for care per month for the period January 2002 to December 2013.

Full size image
Figure 4
figure4

Monthly prevalence trend for each OI by age: <30 yrs (blackline), 30-39 yrs (greenline) & 40 yrs and above (redline); expressed as a proportion of HIV-positive persons diagnosed with a particular OI out of the total number who turned up for care per month for the period January 2002 to December 2013.

Full size image

Mycobacterium tuberculosis

TB monthly prevalence was curvilinear with a steep declining trend between 2003 and 2005 and thereafter dismal change in prevalence was observed (Figure 1). The mean annual prevalence dropped from 76.9/1000 persons at risk in 2002 to 6.6/1000 persons at risk in 2012 but slighly increased to 11.0/1000 persons at risk in 2013 (p < 0.0001, X2trend) (Table 4). The largest reduction occurred in 2005 (52%), but thereafter there was dismal change in prevalence (Table 4). TB monthly prevalence trends were not different by sex (X2 = 0.19, p = 0.66, kruskal wallis) (Figure 3) but differed significantly by age (X2 = 18.5, p = 0.0001, kruskal wallis) (Figure 4) and HIV clinic (X2 = 94.7,p = 0.0001, kruskal wallis) (Figure 5). After adjusting for age, sex and clinic as fixed effects and monthly clustering as a random effect in a mixed effects linear regression model, average monthly prevalence declined at a rate of 2% per month (p < 0.05). However the rate of decline differed significantly by age, sex and clinic (p < 0.05) (Table 5).

Figure 5
figure5

Monthly prevalence trend for each OI by HIV clinic: Mulago HIV clinic (orange line), Tororo HIV clinic (blackline), and Mbarara HIV clinic (blue line); expressed as a proportion of HIV-positive persons diagnosed with a particular OI out of the total number who turned up for care per month for the period January 2002 to December 2013.

Full size image

Oral candida

Oral candida monthly prevalence also showed a curvilinear trend with a steep declining trend between 2003 and 2005 and thereafter dismal change in prevalence was observed (Figure 1). The mean annual prevalence dropped from 173/1000 persons at risk in 2002 to 6.7/1000 persons at risk in 2013 (p < 0.0001, X2trend) (Table 4). The largest reduction in mean annual prevalence occurred in 2005 (61%) (Table 4). The monthly Oral candida prevalence trends were not significantly different by sex (X2 = 0.6, p = 0.44, kruskal wallis) (Figure 2) and age (X2 = 1.7,p = 0.42, kruskal wallis) (Figure 3) but differed significantly by clinic (Figure 4) with Mbarara HIV clinic (south western region) showing the highest prevalence trend, followed by Mulago HIV clinic (central region) and then Tororo HIV clinic (eastern region) (X2 = 63.2, p = 0.0001, kruskal wallis). After adjusting for age, sex and clinic as fixed effects and monthly clustering as a random effect in a mixed effects linear regression model, average monthly prevalence declined at a rate of 3% per month (p < 0.001) (Table 5). However the rate of decline did not differ significantly by age or sex (p > 0.05) but differed significantly by clinic (p < 0.05) (Table 5)

Herpes zoster

Herpes zoster monthly prevalence showed a consistent declining trend (Figure 1). The mean annual prevalence reduced from 13.4/1000 persons at risk in 2002 to 3.3/1000 persons at risk in 2013 (p = 0.0002, X2trend) (Table 4). Largest reduction in mean annual prevalence was observed in 2008 (53.4%) (Table 4). However, Herpes zoster monthly trends were not significantly different by sex (X2 = 0.08, p = 0.77, kruskal wallis) and clinic (X2 = 2.4, p = 0.29, kruskal wallis) but differed significantly by age (X2 = 7.2,p = 0.027, kruskal wallis). After adjusting for age, sex and clinic as fixed effects and monthly clustering as a random effect in a mixed effects linear regression model, monthly prevalence declined significantly at an average rate of 0.087 per month (p < 0.001). However this rate of decline did not differ significantly by age, sex and clinic (p > 0.05) (Table 5).

Genital ulcer

The prevalence of genital ulcer was relatively lower when comapred to Oral candida and TB. Monthly prevalence slightly but significantly reduced over the observation period (Figure 1). The mean annual prevalence dropped from 15.1/1000 persons at risk in 2002 to 8.9/1000 persons at risk in 2013 (p = 0.0029, X2trend) (Table 4). Its monthly prevalence series were significantly different by sex with women having a higher prevalence trend than men (X2 = 14.8, p = 0.0001, kruskal wallis) (Figure 3). The monthly series were also significantly different by age with HIV positive clients aged between 30-39 yrs having a relatively higher prevalence compared to those aged below 30 years and those aged 40 years and above (X2 = 10.6, p = 0.005, kruskal wallis) (Figure 4). The monthly series also differed significantly by HIV clinic with Mulago HIV clinic (central region) showing the highest prevalence trend, followed by Mbarara HIV clinic (south western region) and then lastly Tororo HIV clinic (eastern region) (X2 = 45.3, p = 0.0001) (Figure 5). After adjusting for age, sex and clinic as fixed effects and monthly clustering as a random effect, monthly prevalence declined significantly at an average rate of 0.078 per month. However the rate of decline did not differ significantly by age or sex or clinic (p > 0.05) (Table 5).

Crypotococcal meningitis

Crypotococcal meningitis prevalence was relatively the lowest in this cohort and showed a non-significant trend (Figure 1). The prevalence decreased between 2003 and 2004 but showed intermittent decreasing and increasing trends between 2005 and 2013 with a peak increase observed in 2010. The mean annual prevalence reduced from 5.5/1000 persons at risk in 2002 to 1.2/1000 persons at risk in 2009, then increased to 4.0/1000 persons at risk in 2010 and thereafter declined to 0.23/1000 persons at risk in 2013 (p = 0.181, X2trend). The monthly series were not significantly different by sex (X2 = 3.48, p = 0.062, kruskal wallis) and age (X2 = 3.7, p = 0.15, kruskal wallis) but were significantly different by HIV clinic (X2 = 15.2, p = 0.0005, Kruskal Wallis). After adjusting for age, sex and clinic as fixed effects and monthly clustering as a random effect, monthly prevalence declined significantly at an average rate of 0.3% per month. However the rate of decline did not differ significantly by age or sex or clinic (p > 0.05) (Table 5)

Discussion

On average globally, TB prevalence among HIV+ individuals slightly reduced from around 25% in 2007 to 21% in 2010 [34]. In the current study, we have observed a significant declining trend in mean annual prevalence of Mycobacterium tuberculosis. This is consistent with other previous studies, for example, a study in the USA that evaluated annual trends for 13 most common AIDS-defining opportunistic infections by examining medical records in more than 90 hospitals and clinics in 9 US cities before HAART (1991–96) and showed decreasing trends in tuberculosis [11]. Another study in Italy that examined temporal trends in the incidence of opportunistic infections (OIs) associated with AIDS in the period 1989–1997 showed significant decreasing trends in both pulmonary and extra-pulmonary tuberculosis [20]. However, differs from a cohort study in west London who investigated the incidence of 12 most frequent AIDS-defining illnesses in the pre-HAART and post-HAART time periods and found no significant reduction in the incidence of TB and other OIs but significant decrease in the incidence of PCP, Kaposis sarcoma and cryptosporidiosis [35].

Few studies though have examined the trends in prevalence of TB among HIV+ individuals in low income settings. One study in South Africa examined the trends in incidence rates of TB before HAART (1991–2000) reported a non-significant increasing trend of TB incidence (P = 0.17, X2trend) [22]. Another study in Malawi that examined trends in TB associated with HIV-infection between 1988 to 2001 showed incidence of active TB first increased up to 1990 and thereafter started to decline [36]. Uganda is among the 22 high-burden TB countries in the world [37] and previous studies in Uganda show that over 80% of the HIV-related morbidity and 30% of the HIV related death were due to TB [4,38]. TB co-infection with HIV has been associated with poor HAART prognostic outcomes [39]. A prospective cohort study that assessed the effect of HAART on TB incidence in Eastern Uganda showed TB reduced from 7.2% at baseline to 5.5% after 1.4 yrs of follow up [38]. Our findings are consistent with these previous findings and provide additional evidence of a significant HAART effect on TB prevalence in resource poor settings. However, the slight increase observed in 2013 implies that TB may never be completely eliminated since it is endemic in these settings and was even there before the advent of the HIV epidemic in Uganda.

Herpes zoster showed a significant declining trend though the effect of HAART seems to be relatively less compared to other OIs. Our findings are consistent with findings from a study in the USA that evaluated the effect of HAART on incidence trends of Herpes zoster from 1987 to 2011 and reported a significant reducing trend in annual incidence of Herpes zoster [40]. However our findings deviates from another study that followed HIV- infected patients from 1985 to 2003 in France found no significant difference in the incidence rate for Herpes zoster over time [41]. Another study in the USA reported an increasing trend in the prevalence of clinical Herpes zoster in the period from 1945 to 2008 [42]. Our findings provide additional evidence that increasing coverage of HAART is having a significant effect on Herpes zoster even in resource poor settings. However this effect on prevalence could also be attributed to increasing availability of other potent antiviral drugs like acyclovir.

Oral candida showed very sharp decline of 61% in 2005 probably because of the introduction of HAART in 2005 and increased availability of more potent systemic antifungal drugs like fluconazole. This is consistent with recent study findings in Spain which examined trends of candidiasis among HIV+ children (<15 years) for the period 1997–2008 and found significant decline in candidaisis diagnoses in this time period attributed to the effect of HAART . However, another study in the USA that reviewed the trends in the epidemiology of opportunistic fungal infections associated with HIV/AIDS showed a significant increasing trend in the incidence of fungal infections (Aspergillus sp) attributed to increasing resistance to anti-fungal treatment and recommended documentation of epidemiological trends to gain more insights into the effectiveness of treatment strategies [19]. Though a significant declining trend was observed, the OI has not been completely eliminated probably because it is highly endemic in Uganda or there could also be drug resistance which requires further investigations.

Genital ulcers probably caused by Herpes simplex virus type2 have previously been shown to be common among HIV positive patients [43-45]. A study in France found about 18% of the HIV infected adults were co-infected with herpes simplex virus type2 [46]. Herpes simplex virus-2 infection can be latent in the normal human body but become reactivated when the immune system is severely compromised [47]. In Uganda, a rural population study found the rate ratio for HSV-2 incidence was 3.69 in HIV-positive cases with genital ulcers compared to HIV-negative persons after adjusting for age and sex [43]. The study also found the prevalence of HSV-2 was much higher in women (71.5%) than in men (36.6%). Our study shows prevalence of genital ulcers declined slightly over time and the trends differed significantly by sex. A study in India reported an increasing prevalence trend of genital herpes between 2000 to 2004 but prevalence was reported higher in women compared to men [48]. Similar findings were also reported in another study in India [49]. Our findings provides additional evidence that increased availability of HAART and effective treatment particularly acyclovir could be responsible for the observed decline in the prevalence of genital ulcers in Uganda.

In the current study Cryptococcal meningitis showed a non-significant declining trend with a 50% reduction in 2004 probably due to the introduction of fluconazole around this time and highly active antiretroviral therapy in 2005 [50]. The overall declining trend is consistent with other previous studies in developed countries which also reported declining trends in the era of HAART. A study in Brazil [21] examined temporal trends in incidence of opportunistic infections (OI) among HIV-infected adults (>12 years old) at a national level for the period in 1980–1999. The study showed a declining trend for all OI including Cryptococcal meningitis among adult AIDS cases particularly after the introduction of HAART in 1996. In the USA, a population based surveillance for cryptococcosis conducted between 1992 and 2000 reported a declining trend in the incidence of cryptococcosis among HIV+ individuals attributed to widespread use of effective antiretroviral drugs [51]. Another study [19] that reviewed trends in epidemiology of opportunistic fungal infections showed that increasing use of effective antifungal prescriptions has reduced the frequency of invasive mycoses including cryptococcosis among HIV+ individuals in the US. Our study adduces more evidence that cryptococcosis is reducing among HIV+ individuals over time with increasing access to HAART and fluconazole in Uganda. However the slight increase in 2010 could be attributed to improvement in diagnosis following the introduction of the lateral flow cryptococcal Antigen (CrAg) rapid test strips [50,52]

The study also shows that complete elimination of opportunistic infections may not be possible partly because they are endemic and also as a consequence of gaps in timely access to treatment, poor adherence to HAART, inadequate staff and laboratory equipment, poverty and high patient attrition [53]. We also observed a mean annual attrition rate of 10% in the era of HAART whose cause we can only speculate. Previous studies have shown that despite the best efforts by care providers and donors to provide free lifelong ARVs, HIV patients would still drop out of programme care or withdraw from treatment for various reasons [54-57]. Sydney Rosen and colleagues in their systematic review on prevalence of loss to follow up among patients on ART, reported an average 40% prevalence of loss to follow up [56]. Some studies have also shown retention at 24 month in most ART programmes in Sub-Saharan Africa to be between 60-80% [56,58]. In Uganda, ART attrition rates after 12 months have been reported to reach up to 30% [53,59] while generally, among ART patients who are subject to attrition, i.e. who are lost to follow-up (LTFU), it is estimated that 30-60% are dead [53,59,60].

Lastly, missed clinic visits are inevitable over a lifelong HIV/AIDS care programme. However, our understanding of why people drop out of care and how to address this problem remains limited [61]. Norma Ware and colleagues conducted a large qualitative study among patients in HIV treatment programs in sub-Saharan Africa to investigate reasons for missed clinic visits and reported both intentional and unintentional reasons including cost of transport to the clinic, stigma and dissatisfaction with care, competing demands (both economic and social), finding an alternative source of care and discouragement due to perceived harsh treatment by care providers [57]. However, these factors are likely to vary from one geographical setting to another. Hence, it is crucial to correctly identify these factors so as to improve ART program retention and success in resource poor settings.

Limitations

This study had a number of limitations. First, data analysed was limited to what was available in the TASO electronic data base and so certain variables whose data was not captured could not be analysed. Being a retrospective cohort, there was no information available on the clients who never returned for care. So we did not have information on the connection to care or survival status of clients who were lost to follow up. However we assumed that these simply relocated to another clinic for care or died as had been reported in other previous studies [54,55,59,62,63]. The TASO clients whose medical records were analysed may have not been representative of all HIV+ individuals in Uganda which means generalisability could be limited but still this does not compromise the evidence we have adduced since most if not all beneficiaries of public health care programmes in Uganda are generally the same. There may also have been a possibility that some OIs went away undiagnosed or were not recorded during clinic visits as was reported in a previous study conducted in Uganda [64]. Though we believe this was minimized by selected OIs which were easy to diagnose and of great interest in the country given that the study by Kiragga et al. [64] reported that the level of under reporting was far less for TB and Cryptococal meningitis compared to other less important OIs. From sensitivity analysis, we established that cases who remained in care were not significantly different from those that were lost to follow up and so we conservatively assumed that cases who remained in care were a representative sample of all the cases in the original cohort and that missing clinic visits happened randomly. Even with these limitations, there is sufficient evidence that increasing HAART availability in resource poor settings is having a meaningful impact on prevalence of opportunistic infections.

Conclusions and recommendations

Overall, the data provides emperical evidence that the disease burden due to OIs has significantly reduced in the 12 years of follow up. The decrease in trends of individual OIs could be attribute to the increased access to HAART and other preventive prescriptions such as cotrimoxazole prophylaxis and effective therapeutic drugs over time. The study further shows that some opportunistic infections were still common though at reduced prevalence in spite of availability of effective treatment including HAART. This is probably because they are endemic or this could be a sign of resistance/treatment failure that needs to be further investigated. The differences in trends between OIs and geographical areas were likely to be due to differences in the level of exposure to infectious agents, social-economic status, drug resistance, immunity, nutrition, etc. Future research should explore the association of these factors with trends of individual OIs and knowledge of these factors could help in the design of taregeted interventions. The effect attributable to HAART alone or Cotrimoxazole prophylaxis and other therapeutic drugs was not established as this would require a comparison cohort which is exclusively on HAART alone which was not available. We recommend that future studies should examine this aspect. We also recommend that information on those lost to programme care should be routinely collected in order to have insight on the causes for drop out which will ultemately help to improve ART programmes and inform policy decisions on improvement of treatment outcomes in resource poor settings.

References

  1. 1.

    UNAIDS. Global report. UNAIDS report on the global AIDS epidemic 2013. Geneva: WHO; 2013.

  2. 2.

    Kaplan JE, Hu DJ, Holmes KK, Jaffe HW, Masur H, De Cock KM. Preventing opportunistic infections in human immunodeficiency virus-infected persons: implications for the developing world. Am J Trop Med Hyg. 1996;55(1):1–11.

  3. 3.

    Holmes BC, Losina E, Walensky RP, Yazdanpanah Y, Freedberg KA. Review of Human Immunodeficiency Virus Type 1–related opportunistic infections in Sub-Saharan Africa. Clin Infect Dis. 2003;36:652–62.

  4. 4.

    Mermin J, Were W, Ekwaru JP, Moore D, Dawning R, Behumbiize P, et al. Mortality in HIV-infected Ugandan adults receiving antiretroviral treatment and survival of their HIV-uninfected children: a prospective cohort study. Lancet. 2008;371:752–9.

  5. 5.

    Staine JG. AIDS up date 2007: An overview of acquired immune deficiency syndrome. Inc: New York McGraw-Hill co.; 2008.

  6. 6.

    Microsoft A, Ledergerber B, Katlama C, Kirk O, Reiss P, d’Arminio Monforte A, et al. Decline in the AIDS and death rates in the EuroSIDA study: an observational study. Lancet. 2003;362:22–9.

  7. 7.

    McNaghten AD, Hanson DL, Jones JL, Dworkin MS, Ward JW. Effects of antiretroviral therapy and opportunistic illness primary prophylaxis on survival after AIDS diagnosis. AIDS. 1999;13:1687–95.

  8. 8.

    Sansone GR, Frengley JD. Impact of HAART on causes of death of persons with late stage AIDS. J Urban Health. 2000;77(2):166–75.

  9. 9.

    Jones JL, Hanson DL, Dworkin MS, Alderton DL, Fleming PL, Kaplan JE, et al. Surveillance for AIDS-defining opportunistic illnesses (1992–97). MMWR CDC Surveill Summ. 1999;48(2):1–22.

  10. 10.

    Palella FJ, Delaney KM, Moorman AC, Loveless MO, Fuhrer J, Satten GA, et al. Declining morbidity and mortality among patients with advanced human immunodeficiency virus infection. HIV outpatient study investigators. N Engl J Med. 1998;338:853–60.

  11. 11.

    Jones JL, Hanson DL, Dworkin MS, Kaplan JE, Ward JW. Trends in AIDS-related opportunistic infections among Men Who Have Sex with Men and among injecting drug users, 1991–1996. J Infect Dis. 1998;178:114–20.

  12. 12.

    Kaplan JE, Hanson D, Dworkin MS, Frederick T, Bertolli J, Lindergren ML, et al. Epidemiology of human immunodeficiency virus-associated opportunistic infections in the united states in the era of highly active antiretroviral therapy. Clin Infect Dis. 2000;30 Suppl 1:S5–14.

  13. 13.

    UNAIDS. Access to Antiretroviral therapy in Africa: Satus report on progress towards 2015 targets. Geneva: UNAIDS; 2014.

  14. 14.

    Selik MR, Starcher ET, Curran JW. Opportunistic diseases reported in AIDS patients: Freqiencies, associations and trends. AIDS. 1987;1:175–82.

  15. 15.

    Pezzotti P, Serraino D, Rezza G, Maso LD, Vaccher E, Lepri AC, et al. The spectrum of AIDS-defining diseases: temporal trends in Italy prior to the use of highly active antiretroviral therapies, 1982–1996. Int J Epidemiol. 1999;28:975–81.

  16. 16.

    Nesheim SR, Kapogiannis BG, Soe MM, Sullivan KM, Abrams E, Farley J, et al. Trends in opportunistic infections in the pre- and post- HAART eras among HIV-infected children in the perinatal AIDS collaborative transmission study (1986–2004). Pediatrics. 2007;120:100–9.

  17. 17.

    Muñoz A, Schrager LK, Bacellar H, Speizer I, Vermund SH, Detels R, et al. Trends in the incidence of outcomes defining acquired immunodeficiency syndrome (AIDS) in the Multicenter AIDS Cohort Study: 1985–1991. Am J Epidemiol. 1993;137(4):423–38.

  18. 18.

    Wohl AR, Lu S, Turner J, Kovacs A, Witt M, Squires K, et al. Risk of opportunistic infections in the HAART era among HIV- infected Latinos born in the United States compared to Latinos born in Mexico and Central America. AIDS Patient Care STDs. 2003;17(6):267–75.

  19. 19.

    Singh N. Trends in the epidemiology of opportunistic fungal infections: predisposing factors and the impact of antimicrobial use practices. Clin Infect Dis. 2001;33:1692–6.

  20. 20.

    San-Andres FJ, Rubio R, Castilla J, Pulido F, Palao G, de Pedro I, et al. Incidence of acquired immunodeficiency syndrome–associated opportunistic diseases and the effect of treatment on a cohort of 1115 patients infected with human immunodeficiency virus, 1989–1997. Clin Infect Dis. 2003;36:1177–85.

  21. 21.

    Guimaraes MD, Castilho EA, Fonseca MG. Temporal trends in AIDS associated opportunistic infections in Brazil, 1980–1999. In: Int Conf AIDS, vol. 13. 2000. abstract MoPeB2286.

  22. 22.

    Corbett EL, Marston B, Churchyard GJ, De Cock KM. Tuberculosis in sub-Saharan Africa: Opportunities, challenges, and change in the era of anti-retroviral treatment. Lancet. 2006;367:926–37.

  23. 23.

    WHO. Consolidated guidelines on the use of antiretroviral drugs for treating and preventing HIV infection: recommendations for a public health approach. Geneva: World Health Organization; 2013.

  24. 24.

    UNAIDS. HIV and AIDS Uganda Country Progress Report 2013. Geneva: UNAIDS; 2014.

  25. 25.

    Ministry of Finance and Economic Development (MoFPED). Millennium Development Goals Report for Uganda 2010. Kampala; Uganda

  26. 26.

    WHO. Antiretroviral therapy for HIV infection in adults and adolescents; recommendations for a public health approach. Geneva: WHO; 2006.

  27. 27.

    WHO. Antiretroviral therapy for HIV infection in adults and adolescents: recommendations for a public health approach 2010 revision. Geneva: WHO; 2010.

  28. 28.

    Calvarese M, Bame S, Bakanume B. Historical analysis of AIDS patients in Uganda using innovative community clinical service: The AIDS Support Organisation (TASO). J Hum Soc Sci. 2007;1:1–13.

  29. 29.

    TASO. TASO Annual Report. Kampala, Uganda: TASO; 2013.

  30. 30.

    TASO Services and Programmes. [http://www.tasouganda.org/].

  31. 31.

    Ministry of Health Uganda. Antiretroviral treatment policy. Kampala: Uganda Ministry of Health; 2009.

  32. 32.

    Box, Jenkins, Reinsel, editor. Time Series Analysis, Forecasting and Control. 3rd ed. Englewood Clifs, NJ: Prentice Hall; 1994.

  33. 33.

    Stata corp. Nonparametric analysis: kruskal wallis rank test, vol. 2012. 2012th ed. College statation, Texas, USA: Stata corp LP; 2012.

  34. 34.

    WHO. Global Tuberculosis Control Report 2013. Geneva, Switzerland: WHO; 2013.

  35. 35.

    Ives NJ, Gazzard BG, Easternbrook PJ. The changing pattern of AIDS-defining illnesses with the introduction of highly active antiretroviral therapy (HAART) in a London clinic. J Infect. 2001;42:134–9.

  36. 36.

    Glynn JR, Crampin AC, Ngwira BMM, Mwaungulu FD, Mwafulirwa DT, Floyd S. Po¨ nnighaus JM, Warndorff D K, E.M. FP: Trends in tuberculosis and the influence of HIV infection in northern Malawi, 1988–2001. AIDS. 2004;18:1459–63.

  37. 37.

    WHO. Global tuberculosis control: surveillance, planning, financing: WHO report 2005. WHO/HTM/ TB/2005.349. Geneva: WHO; 2005.

  38. 38.

    Moore D, Liechty C, Ekwaru P, Were W, Mwima G, Solberg P, et al. Prevalence, incidence and mortality associated with tuberculosis in HIV-infected patients initiating antiretroviral therapy in rural Uganda. AIDS. 2007;21(6):713–9.

  39. 39.

    UNAIDS. Government of Uganda: UNGASS 2010 country progress report. Geneva: WHO; 2010.

  40. 40.

    Moanna A, Rimland D. Decreasing incidence of Herpes zoster in the HAART era. Clin Infect Dis. 2013;57(1):122–5.

  41. 41.

    Vanhems P, Voisin L, Gayet-Ageron A, Trepo C, Cotte L, Peyramond D, et al. The incidence of Herpes zoster is less likely than other opportunistic infectionss to be reduced by highly active antiretroviral therapy. J Acquir Immune Defic Syndr. 2005;38:111–3.

  42. 42.

    IDSA. Trends in Herpes zoster prevalence in a Populaion based cohort from 1945 to 2008. In: IDSA 48th Annual Meeting: 2010. Vancouver, Canada: IDSA; 2010.

  43. 43.

    Kamali A, Nunn AJ, Mulder DW, Dyck EV, Dobbins JG, Whitworth AG. Seroprevalence and incidence of genital ulcer infections in a rural Ugandan population. Sex Transm Infect. 1999;75:98–102.

  44. 44.

    Glynn RJ, Biraro S, Weiss HA. Herpes simplex virus type 2: a key role in HIV incidence. AIDS. 2009;23:1595–8.

  45. 45.

    Freeman EE, Weiss HA, Glynn JR, Cross PL, Whitworth JA, Hayes RJ. Herpes simplex virus 2 infection increases HIV acquisition in men and women: systematic review and meta-analysis of longitudinal studies. AIDS. 2006;20(1):73–83.

  46. 46.

    Jougla E, Pequignot F, Carbon C, Pavillon G, Eb M, Bourdais JP, et al. AIDS-related conditions: Study of a representative sample of 1203 patients deceased in 1992 in France. Int J Epidemiol. 1996;25:190–7.

  47. 47.

    Looker KJ, Garnett GP, Schmid GP. An estimate of the global prevalence and incidence of herpes simplex virus type 2 infection. In: Bulletin of the World Health Organization. Geneva: WHO; 2007.

  48. 48.

    Mendiratta V, Bhawna H. Genital herpes: changing patterns and current trends. Indian J Sex Transm Dis. 2006;27(1):6.

  49. 49.

    Kumar B, Sahoo B, Gupta S, Jain R. Raising incidence of genital herpes over two decades in a sexually transmitted disease clinic in North India. J Dermatol. 2002;29:74–8.

  50. 50.

    Jarvis JN, Lawn SD, Vogt M, Bangani N, Wood R, Harrison TS. Screening for cryptococcal antigenemia in patients accessing an antiretroviral treatment program in South Africa. Clin Infect Dis. 2009;48(7):856–62.

  51. 51.

    Warnock DW. Trends in the epidemiology of invasive Fungal infections. Jpn J Infect Dis. 2007;48:1–12.

  52. 52.

    WHO. Rapid Advice: Diagnosis, Prevention, and Management of Cryptococcal Disease in HIV-Infected Adults, Adolescents, and Children. Geneva: WHO; 2011.

  53. 53.

    WHO. Retention in HIV programmes: defining the challenges and identifying solutions: meeting report. Geneva: WHO; 2011.

  54. 54.

    Brinkhof MW, Dabis F, Myer L, Bangsberg DR, Boulle A, Nash D, et al. Early loss of HIV-infected patients on potent antiretroviral therapy programmes in lower-income countries. Bull World Health Organ. 2008;86:559–67.

  55. 55.

    Egger M, Spycher BD, Sidle J, Weigel R, Geng EH, Fox MP, et al. Correcting mortality for loss to follow-up: a nomogram applied to antiretroviral treatment programmes in sub-Saharan Africa. PLoS Med. 2011;8(1):e1000390.

  56. 56.

    Rosen S, Fox MP, Gill C. Patient retention in antiretroviral therapy programs in sub-Saharan Africa: a systematic review. PLoS Med. 2007;4:e298.

  57. 57.

    Ware CN, Wyatt MA, Geng HE, Kaaya SF, Agbaji OO, Muyindike WR, et al. Toward an understanding of disengagement from HIV treatment and care in Sub-Saharan Africa: a qualitative study. PLoS Med. 2013;10(1):e1001369.

  58. 58.

    Fox MP, Rosen S. Patient retention in antiretroviral therapy programs up to three years on treatment in sub-Saharan Africa, 2007–2009: systematic review. Trop Med Int Health. 2010;15(Sup1):1–15.

  59. 59.

    Geng EH, Glidden DV, Bwana MB, Musinguzi N, Emenyonu N, Muyindike W, et al. Retention in care and connection to care among HIV-infected patients on antiretroviral therapy in Africa: estimation via a sampling-based approach. PLoS One. 2011;6(7):e21797.

  60. 60.

    Nakiwogga-Muwanga A, Alamo-Talisuna S, Musaazi J, Kambugu A, Ssekawungu P, Katabira E, et al. Inadequate monitoring in advanced stages of disease with lack of supportive counseling increasesAttrition among patients on antiretroviral treatment at a large urban clinic in Uganda. J Int Assoc Prov AIDS Care (JIAPAC). 2014;13(6):547–54.

  61. 61.

    Mills E, Nabiryo C. Preventing antiretroviral treatment interruptions among HIV/AIDS patients in Africa. PLoS Med. 2013;10(1):e1001370.

  62. 62.

    Braitstein P, Brinkhof MW, Dabis F, Schechter M, Boulle A, Miotti P, et al. Mortality of HIV-1-infected patients in the first year of antiretroviral therapy: comparison between low-income and high-income countries. Lancet. 2006;367(9526):817–24.

  63. 63.

    Kiragga NA, Lok JJ, Musick SB, Bosch JR, Mwangi A, Wools-Kaloustian KK, et al. CD4 trajectory adjusting for dropout among HIV-positive patients receiving combination antiretroviral therapy in an East African HIV care centre. J Int AIDS Soc. 2014;17:18957.

  64. 64.

    Kiragga NA, Castelnuovo B, Schaefer P, Muwonge T, Easterbrook JP. Quality of data collection in a large HIVobservational clinic database in sub-Saharan Africa: implications for clinical research and audit of care. J Int AIDS Soc. 2011;14:3.

Download references

Acknowledgements

We are grateful to the AIDS Support Organisation (TASO) in Uganda for providing the data for the study. We would also like to extend our sincere appreciation to Mountains of the Moon University for the financial support.

Author information

Affiliations

  1. Department of Epidemiology and Biostatistics, School of Public Health, College of Health Sciences, Makerere University, Kampala, Uganda

    • John Rubaihayo
    • , Nazarius M Tumwesigye
    •  & Joseph Konde-Lule

  2. Department of Public Health, School of Health Sciences, Mountains of the Moon University, Fort Portal, Uganda

    • John Rubaihayo

Authors

  1. Search for John Rubaihayo in:

  2. Search for Nazarius M Tumwesigye in:

  3. Search for Joseph Konde-Lule in:

Corresponding author

Correspondence to John Rubaihayo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JR: Conceived the study, participated in the study design, data collection, analysis and manuscript writing; NMT: Participated in study design, data analysis and manuscript writing; JKL: Participated in study design, data analysis and manuscript writing. All authors read and approved the final manuscript

Rights and permissions

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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Human Immunodeficiency Virus
  • Herpes Zoster
  • Cryptococcal Meningitis
  • Human Immunodeficiency Virus Patient
  • Resource Poor Setting

BMC Infectious Diseases

ISSN: 1471-2334

Contact us