# An Atypical Age-Specific Pattern of Hepatocellular Carcinoma in Peru: A Threat for Andean Populations

Stephane Bertani, Pascal Pineau, Sebastian Loli, Julien Moura, Mirko Zimic, Eric Deharo, Eloy Ruiz

► To cite this version:

Stephane Bertani, Pascal Pineau, Sebastian Loli, Julien Moura, Mirko Zimic, et al. An Atypical Age-Specific Pattern of Hepatocellular Carcinoma in Peru: A Threat for Andean Populations. PLoS ONE, 2013, 8 (6), pp.e67756. 10.1371/journal.pone.0067756. hal-01356600

HAL Id: hal-01356600

https://hal.science/hal-01356600v1

Submitted on 26 Aug 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# An Atypical Age-Specific Pattern of Hepatocellular Carcinoma in Peru: A Threat for Andean Populations

Stéphane Bertani¹,²*§, Pascal Pineau³¶, Sebastian Loli⁴, Julien Moura⁵, Mirko Zimic⁴, Eric Deharo¹,², Eloy Ruiz⁶

¹Université de Toulouse 3, Unité de Pharmacochimie et Pharmacologie pour le Développement, UMR152, Toulouse, France, ²Institut de Recherche pour le Développement, Unité de Pharmacochimie et Pharmacologie pour le Développement, UMR152, Lima, Peru, ³Institut Pasteur, Institut National de la Santé et de la Recherche Médicale, Nuclear Organization and Oncogenesis Unit, U993, Paris, France, ⁴Universidad Peruana Cayetano Heredia, Facultad de Ciencias, Laboratorios de Investigación y Desarrollo, Unidad de Bioinformática y Biología Molecular, Lima, Peru, ⁵Institut de Recherche pour le Développement, Programme Andin de Recherche et de Formation sur la Vulnérabilité et les Risques en Milieu Urbain, Lima, Peru, ⁶Instituto Nacional de Enfermedades Neoplásicas, Departamento de Cirugía en Abdomen, Lima, Peru

## Abstract

**Background:** In South America, the highest incidence of primary liver cancer is observed in Peru. However, national estimations on hepatocellular carcinoma incidence and mortality are approximated using aggregated data from surrounding countries. Thus, there is a lack of tangible information from Peru that impairs an accurate description of the local incidence, presentation, and outcomes of hepatocellular carcinoma. The present study attempts to fill this gap and assesses the clinical epidemiology of hepatocellular carcinoma in this country.

**Methods:** A retrospective cohort study was conducted by analysing the medical charts of 1,541 patients with hepatocellular carcinoma admitted between 1997 and 2010 at the Peruvian national institute for cancer. The medical records including liver function, serologic status, and tumor pathology and stage were monitored. Statistical analyses were performed in order to characterize tumor presentation according to demographic features, risk factors, and regional origin.

**Results:** Surprisingly, the age distribution of the patient population displayed bimodality corresponding to two distinct age-based subpopulations. While an older group was in keeping with the age range observed for hepatocellular carcinoma around the world, a younger population displayed an abnormally juvenile mean age of 25.5 years old. In addition, each subpopulation displayed age-specific pathophysiological and clinical characteristics.

**Conclusions:** The analysis suggests two different age-specific natural histories of hepatocellular carcinoma in the Peruvian patient population. This otherwise unusual tumor process that is ongoing in younger patients leads to the hypothesis that there may be a Peru-endemic risk factor driving hepatocarcinogenesis in the local population.

**Citation:** Bertani S, Pineau P, Loli S, Moura J, Zimic M, et al. (2013) An Atypical Age-Specific Pattern of Hepatocellular Carcinoma in Peru: A Threat for Andean Populations. PLoS ONE 8(6): e67756. doi:10.1371/journal.pone.0067756

**Editor:** Jung Weon Lee, Seoul National University, Republic of Korea

**Received** January 29, 2013; **Accepted** May 22, 2013; **Published** June 28, 2013

**Copyright:** © 2013 Bertani et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** The authors have no support or funding to report.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: stephane.bertani@ird.fr

¶ These authors contributed equally to this work.

## Introduction

Hepatocellular carcinoma (HCC), the main form of primary liver cancer, is the sixth most common malignancy and the third leading cause of tumor-related death in the world [1,2]. Global clinical epidemiology of HCC defines a dominant patient profile corresponding grossly to males over 40 years old [3,4]. The incidence rate of HCC has doubled worldwide during the last two decades, with nearly 85% of the recorded cases occurring in developing countries [1,2,5]. The greatest burden of HCC is borne in sub-Saharan Africa and eastern Asia, where chronic infection with hepatitis B virus (HBV) is highly endemic [1–6].

Reference reviews on the worldwide burden of HCC have constantly neglected to include the epidemiology of the disease in South America creating a gap between the existing epidemiolog-

ical situation on the field and the global epidemics of HCC described in literature [2–5]. While the incidence rate of primary liver cancer in South America is considered low, the epidemiology of HCC on this continent displays intriguing characteristics. For example, the sex ratio of HCC is more balanced in South America than in any other regions of the world [1]. Nonetheless, with the exception of Brazil, a lack of information on primary liver cancers impairs an accurate description of HCC outcomes at both national and continental levels [1]. Incidence and mortality data for HCC in South America are often gathered from sparse cancer registries, and estimations are thus approximated using aggregated data from surrounding countries [1].

Peru is considered to have the highest incidence of primary liver cancer in South America [1]. Very few studies have been conducted concerning the clinical epidemiology of HCC among
the Peruvian population; therefore our knowledge of the disease in Peru is lacking. To address this issue, we performed a retrospective study on Peruvian patients diagnosed with HCC between 1997 and 2010. Specifically, this study was designed to assess the clinical epidemiology of Peruvian patient population with HCC in a large well-characterized cohort.

## Methods

### 2.1. Study Design

The current study was conducted retrospectively within a cohort assembled by analysing the medical records of 1,541 patients admitted for HCC at the Peruvian national hospital for cancer (INEN) between January 1997 and December 2010. Written consent was given by the patients for their information to be stored in the Department of Cancer Statistics and Epidemiology of INEN, and used for research. Human Subjects committees at both the INEN (Peru) and the Institut Pasteur (France) approved this study.

Under the Peruvian Ministry of Health, INEN is the health care institution in charge of the management of neoplastic diseases at the national level [7]. As a public hospital, INEN treats patients regardless of age, sex, ethnicity, place of residence, economic status, and health care coverage. The centralization of the Peruvian health care system means that INEN serves as a national hub for neoplasm management and handles a large ratio of cancer cases from across the nation, providing a valuable environment for assessing the clinical-epidemiological context of HCC in Peru [7].

The vital records (i.e. sex, age, birthplace, place of residence, social situation, and family health history), personal medical history, liver function, and physiological, biochemical, and immunological status of patients with HCC were monitored for the duration of their hospital stay. The overall patient population was 1,541, corresponding to the whole patient population with HCC admitted at INEN between 1997 and 2010. However, due to the context of limited resources in Peru, INEN could not routinely examine all clinical specimens for every biomedical analysis, notably for hepatitis infection and tumor grading [8]. The opportunity of each biomedical analysis was then at the own discretion of the practitioner taking into account the patient's best interests. The available cohorts for the different pathophysiological parameters analyzed from the initial 1,541 patient's number are presented (Table 1).

### 2.2. Hospital Procedure and Monitored Parameters

Individuals suspected to have malignant liver neoplasm were managed through the Department of Abdominal Surgery of INEN. The vital records and medical history of these patients were documented during their first consultation. Detection and preliminary clinical characterization of the tumor were carried out by both physical examination and non-invasive diagnostic imaging, i.e. computed axial tomography and abdominal ultrasonics. Cancer staging was classified according to the TNM and the Barcelona-Clinic Liver Cancer (BCLC) staging systems [9]. Pathologists assessed the type of cancer cells on hematoxylin-eosin-stained liver biopsy sections and scored the grade of tumor (G1-3) according to the recommendations of the American Joint Committee on Cancer [10]. Hepatoblastoma and intrahepatic cholangiocellular carcinoma cases were not included in the study.

Liver function tests were conducted by serum quantitation of gamma-glutamyl transpeptidase, alkaline phosphatase, total and direct bilirubin, albumin, and alanine and aspartate transaminases. Tumor marker alpha-fetoprotein (AFP) was monitored for its serum concentration (seroAFP) by radioimmunoassay (Roche)

during the patient's hospital stay. In parallel, infections with HBV and hepatitis C virus (HCV) were assessed from serum of HCC patients by electrogenerated chemiluminescence assay using antibodies against HBV surface antigen (HBsAg), anti-HBs antibody, anti-HBc IgG and IgM, HBe antigen, and anti-HCV antibody (all from Roche).

Patients with a Child-Pugh score of A were treated by anatomic liver resection, i.e. systematic removal of the tumoral liver segments confined by portal branches, to ensure tumor-free margins [11-13]. The surgical report recorded the size, location, and appearance of the tumor, as well as additional pathological observations (e.g. cirrhotic, multinodular, and/or metastatic status). Patient survival was followed-up prospectively by periodic medical consultations or phone interviews.

### 2.3. Statistical Analysis

Age distribution was analyzed using a Gaussian mixture model for detecting bimodality and calculating the bimodality index (BI) and the moment of mixture (MM) [14-16]. Description of variables was reported by median (inter-quartile range) or mean (95% confidence interval) depending on the normality of the distribution assessed by the Shapiro-Wilks test [17]. Bivariable comparisons were conducted by parametric or non-parametric tests depending on the normality of the distribution. A multivariable analysis was performed in order to identify the age-adjusted effect. Linear regressions were used to model the dispersion of normally distributed continuous variables. Simple exact and multinomial logistic regressions were used to model the outcomes with dichotomized age (i.e. ≤44 vs. >44) as the predictor after patients' geographical precedence and seasonal adjustments. Multidimensional outliers were assessed with the Hadi test and nested models were compared with the likelihood ratio test [18]. The tumor size was transformed with the square root function in order to make it normally distributed. Patient survival rates were analyzed using the Kaplan-Meier method and the log-rank test was used to compare survival distributions [19-21]. Statistical analyses were performed with a 5% significance level using Stata Statistical Software: Release 11 (College Station, TX: StataCorp LP) and SPSS 13.0 predictive analytics software (SPSS Inc.). Patient places of origin (birthplace) were spatially aggregated according to a standard deviation classification of mean age, and choropleth mapping was performed using Philcarto mapping software [22].

## Results

The INEN admitted an average of 112 HCC cases per year (Figure 1A). The mean age of the patients was 44.8±21.9 years old (Table 1). The age distribution of the HCC patients displayed bimodality with two frequency modes corresponding to two distinct age-based subpopulations (BI = 1.95) (Figure 1B,C and Table 1). The MM of the two modes corresponding to the in-between age was calculated at 44.8 years old (Figure 1B). The patient distribution for the two modes was 778 for the group younger or equal to age 44 (≤44 group), and 763 the group older than age 44 (>44 group) (Figure 1B,C and Table 1). The respective mean and median ages were 25.5±9.5 and 25 in the ≤44 group, and 64.5±10.4 and 65 in the >44 group (Figure 1B and Table 1). The overall sex ratio was 1.6, with 61.1% and 38.9% of men and women, respectively (Table 1). The age-specific sex ratios were significantly different with 2 and 1.2 in the ≤44 group and the >44 group, respectively (*P*<0.0001) (Figure 1C and Table 1).
<table><caption>**Table 1.** Patient characteristics at baseline.</caption><thead><tr><th></th><th></th><th>Number of patients<br>(≤44 and >44 groups)</th><th>Number of patients<br>(≤44 group)</th><th>Number of patients<br>(>44 group)</th><th>P-value</th></tr></thead><tbody><tr><td>Overall cohort</td><td></td><td>1,541 (100%)</td><td>778 (50.5%)</td><td>763 (49.5%)</td><td>-</td></tr><tr><td rowspan="3">Age</td><td>Mean±SEM</td><td>44.8±21.9</td><td>25.5±9.5</td><td>64.5±10.4</td><td>-</td></tr><tr><td>Median</td><td>44</td><td>25</td><td>65</td><td>-</td></tr><tr><td>Range</td><td>[0-89]</td><td>[0-44]</td><td>[45-89]</td><td>-</td></tr><tr><td rowspan="3">Sex</td><td>Male</td><td>942 (61.1%)</td><td>520 (66.8%)</td><td>422 (55.3%)</td><td>-</td></tr><tr><td>Female</td><td>599 (38.9%)</td><td>258 (33.2%)</td><td>341 (44.7%)</td><td>-</td></tr><tr><td>Total</td><td>1,541 (100%)</td><td>778 (100%)</td><td>763 (100%)</td><td>&lt;0.0001*</td></tr><tr><td rowspan="4">Tumor grade<sup>†</sup></td><td>G1</td><td>47 (13.7%)</td><td>22 (14.3%)</td><td>25 (13.2%)</td><td>-</td></tr><tr><td>G2</td><td>226 (65.7%)</td><td>86 (55.8%)</td><td>140 (73.7%)</td><td>-</td></tr><tr><td>G3</td><td>71 (20.6%)</td><td>46 (29.9%)</td><td>25 (13.2%)</td><td>-</td></tr><tr><td>Total</td><td>344 (100%)</td><td>154 (100%)</td><td>190 (100%)</td><td>0.037**</td></tr><tr><td rowspan="3">Cirrhosis</td><td>Yes</td><td>163 (11%)</td><td>38 (5.1%)</td><td>125 (17.1%)</td><td>-</td></tr><tr><td>No</td><td>1,316 (89%)</td><td>709 (94.9%)</td><td>607 (82.9%)</td><td>-</td></tr><tr><td>Total</td><td>1,479 (100%)</td><td>747 (100%)</td><td>732 (100%)</td><td>&lt;0.0001*</td></tr><tr><td rowspan="3">Hepatitis B</td><td>Yes</td><td>453 (50.1%)</td><td>365 (71%)</td><td>88 (22.5%)</td><td>-</td></tr><tr><td>No</td><td>452 (49.9%)</td><td>149 (29%)</td><td>303 (77.5%)</td><td>-</td></tr><tr><td>Total</td><td>905 (100%)</td><td>514 (100%)</td><td>391 (100%)</td><td>&lt;0.0001*</td></tr><tr><td rowspan="3">Hepatitis C</td><td>Yes</td><td>29 (4.7%)</td><td>9 (2.7%)</td><td>20 (6.8%)</td><td>-</td></tr><tr><td>No</td><td>593 (95.3%)</td><td>319 (97.3%)</td><td>273 (93.2%)</td><td>-</td></tr><tr><td>Total</td><td>622 (100%)</td><td>328 (100%)</td><td>293 (100%)</td><td>0.021*</td></tr><tr><td rowspan="3">Hepatitis B+C</td><td>Yes</td><td>8 (1.4%)</td><td>5 (1.6%)</td><td>3 (1.1%)</td><td>-</td></tr><tr><td>No</td><td>583 (98.6%)</td><td>301 (98.4%)</td><td>282 (98.9%)</td><td>-</td></tr><tr><td>Total</td><td>591 (100%)</td><td>306 (100%)</td><td>285 (100%)</td><td>NS*</td></tr><tr><td rowspan="3">Nodule(s)</td><td>=1</td><td>450 (41.4%)</td><td>246 (38.3%)</td><td>204 (45.8%)</td><td>-</td></tr><tr><td>&gt;1</td><td>638 (58.6%)</td><td>397 (61.7%)</td><td>241 (54.2%)</td><td>-</td></tr><tr><td>Total</td><td>1,088 (100%)</td><td>643 (100%)</td><td>445 (100%)</td><td>0.015*</td></tr><tr><td rowspan="3">Metastasis</td><td>Yes</td><td>430 (30.1%)</td><td>244 (33.3%)</td><td>186 (26.7%)</td><td>-</td></tr><tr><td>No</td><td>1,000 (69.9%)</td><td>489 (66.7%)</td><td>511 (73.3%)</td><td>-</td></tr><tr><td>Total</td><td>1,430 (100%)</td><td>733 (100%)</td><td>697 (100%)</td><td>0.007*</td></tr><tr><td rowspan="3">Operation</td><td>Yes</td><td>298 (19.3%)</td><td>159 (20.4%)</td><td>139 (18.2%)</td><td>-</td></tr><tr><td>No</td><td>1,243 (80.7%)</td><td>619 (79.6%)</td><td>624 (81.8%)</td><td>-</td></tr><tr><td>Total</td><td>1,541 (100%)</td><td>778 (100%)</td><td>763 (100%)</td><td>NS*</td></tr><tr><td rowspan="3">Recurrence<sup>‡</sup></td><td>Yes</td><td>132 (44.7%)</td><td>76 (48.7%)</td><td>56 (40.3%)</td><td>-</td></tr><tr><td>No</td><td>163 (55.3%)</td><td>80 (51.3%)</td><td>83 (59.7%)</td><td>-</td></tr><tr><td>Total</td><td>295 (100%)</td><td>156 (100%)</td><td>139 (100%)</td><td>NS*</td></tr></tbody></table>

≤44 and >44 groups are defined as patients below or equal age 44 and those above age 44, respectively. Mean values are presented with ± standard errors of the means (SEM). Percentages are expressed as ratio of the total patient population for the considered parameter. P-values were obtained by exact simple (*) or multinomial (**) logistic regression analysis, modelling the specific outcomes with age (≤44 groups vs. >44 group) after geographical and seasonal adjustments. NS: not significant.

†Tumor grades were defined as follow: well differentiated (G1), moderately differentiated (G2), and poorly differentiated (G3).

‡The data on recurrence are presented for patients who developed a new HCC within the following 12 months after surgical intervention. doi:10.1371/journal.pone.0067756.t001

Patients were hospitalized with an average tumor size of 12.3±5.6 cm-diameter. The age-specific average tumor size was 13.1±5.7 and 11.4±5.4 cm-diameter for the ≤44 and >44 groups, respectively. The tumor size was transformed with the square root function in order to make it normally distributed for linear regression analysis. ≤44 and >44 groups displayed significant difference in tumor dimension (*P*<0.0001). A higher prevalence of distant metastases was recorded in the ≤44 group (*P*=0.007), while multinodular HCCs were preponderant in the >44 group (*P*=0.015) (Table 1). A significant ratio of undiffer-

entiated or poorly differentiated tumors with high grades (G3) among the ≤44 group compared to the >44 group (*P*=0.037) (Table 1). The tumor size standardization demonstrated a drastic 6.7 fold elevation of AFP in the ≤44 group with 29,124±3,382 ng/ml/tumor-cm, compared to the 4,331±311 ng/ml/tumor-cm of the >44 group (*P*<0.0001) (Figure 2A). Although it remained remarkably low in both subpopulations, the prevalence of associated cirrhosis was significantly higher in the >44 group, showing 17.1% vs. 5.1% in the ≤44 group (*P*<0.0001) (Figure 2B and Table 1). Linear regression
analyses on the quantitation of gamma-glutamyl transpeptidase, alkaline phosphatase, bilirubin, albumin, and transaminases did not demonstrate neither deviant values nor significant difference between the ≤44 and >44 groups (*P*>0.05).

While the rate of HBV infection was substantially higher in the younger population (*P*<0.0001), the HCV infection rate was significantly higher in the older one (*P*=0.021) (Figure 2B and Table 1). Logistic regression analyses showed that cancer history in both nuclear and extended families was not a discriminating trait between the two age-based subpopulations (*P*=0.123). Overall, 46.2% of patients did not present any associated risk factor, such as cirrhosis or hepatitis (Figure 2B).

Among the 298 patients selected for surgery, 44.7% developed a recurrence within 12 months following anatomic liver resection (Table 1). The survival rate at 5 years following the intervention was not significantly different between the two subpopulations, with 25.9% and 19.9% for the ≤44 and >44 groups, respectively (*P*=0.104) (Figure 2C-E and Table 1).

Mapping patients' place of origin revealed geographic disparities in age-based HCC development (Figure 3). While the youngest mean ages were recorded in the inland Andean regions with the lowest value in Apurimac region [31.6±17.9 years old (Figure 3, region number 3)], the oldest mean ages were monitored in the coastal and southern regions, culminating in La Libertad region [61.6±19.3 years old (Figure 3, region number 13)].

## Discussion

This study sought to characterize the clinical-epidemiological situation of HCC in Peru, the country in South America with highest incidence of primary liver cancers [1]. The study was conducted by retrospectively analysing the clinical records of patients with HCC hospitalized at INEN between 1997 and 2010. During this period, INEN admitted over 100 cases of HCC per year, about 25% of the national incident cases of primary liver cancer according to the GLOBOCAN 2008 estimates (Figure 1A and Table 1) [1]. Patients hospitalized with HCC were mostly low- and medium-income patients with poor health care coverage. Positioned at a central hub for neoplasm treatment in Peru, INEN received patients from all 25 administrative regions of the country providing an *ad rem* nationwide sampling perspective (Figure 3).

It is generally admitted that HCC occurs mostly in men over 40 years old [1,3,5]. In high incidence areas of West Africa and East Asia, age-specific incidence distributions are unimodal and increase after 40 years of age to culminate around 55 [5]. In our study, the age-specific distribution of the patients exhibited an astonishing bimodal pattern corresponding to two distinct age-based subpopulations (Figure 1B,C). This unusual feature prompted us to analyze subsequently the cohort according to age by dividing the series in two distinct subsets below and above age 44, which corresponds to the MM of the bimodal Gaussian distribution. Between 1997 and 2010, the proportions of these two age-based subpopulations were similar with the ≤44 and >44 groups representing respectively 50.5% and 49.5% of the overall patient population (Table 1). A careful examination of the ≤44 and >44 group distributions on an annual basis did not show any noticeable evolution of age repartition over the period; indicating that the presence of young HCC patients is a constitutive feature of the disease in Peru. While the older group was in keeping with the age range observed over the world, the ≤44 group with a mean age of 25 corresponded to an epidemiologically deviant population of patients with regard to the world standard of HCC development [1-6]. A study conducted in Mozambique in 1985 showed a significant HCC incidence among young patients,

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Distribution of the patients admitted at INEN with HCC from January 1997 to December 2010. Histograms show the distribution of patients over the period in function of their diagnosis date (A), age at the time of diagnosis (B,C), and gender (C). X-axes display date (annual interval) (A) and age (3-year interval) (B,C); Y-axes display fraction (A,B) and frequency (Freq) (C) of the total number of patients (N = 1,541). Straight line represents the histogram curve fitting with a Gaussian mixture function (BI = 1.95; MM: 44.8) (B). Female and male patients are represented in light and dark greys, respectively (C). doi:10.1371/journal.pone.0067756.g001</figcaption></figure>
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

E

<table><thead><tr><th></th><th>OVERALL</th><th>&le;44 Group</th><th>&gt;44 Group</th></tr></thead><tbody><tr><td>Median survival (months)</td><td>21.2 (16.4-36)</td><td>22.3 (12.7-31.8)</td><td>17.9 (7.7-28)</td></tr><tr><td>1-year survival rate (%)</td><td>61.2</td><td>66.6</td><td>53.9</td></tr><tr><td>3-year survival rate (%)</td><td>35.5</td><td>38.5</td><td>31.5</td></tr><tr><td>5-year survival rate (%)</td><td>23.4</td><td>25.9</td><td>19.9</td></tr></tbody></table>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Pathophysiological parameters associated with HCC in the two age-based patient populations. (A) Box plots show quartile distribution of seroAFP in ng/ml in ≤44 group (grey box; N=631) and >44 group (black box; N=558). \*P<0.0001 vs. >44 group (i.e. adjusted multiple linear regression analysis). (B) Estimated fractions (%) of HCC attributable to HBV and HCV etiologies (VirEtiol) or non-viral cirrhosis (CIR), or non-attributable to bona fide risk factor (ND) in ≤44 females (hatched grey) and males (solid grey) and >44 females (hatched black) and males (solid black). (C) Overall Kaplan-Meier survival curve for patients selected for anatomic liver resection (N = 298). (D) Age-based Kaplan-Meier survival curves for the ≤44 group (grey curve; N = 156) and the >44 group (black curve; N = 139). (E) Kaplan-Meier survival rates (%) monitored for the 2 age-based populations at 1, 3, and 5 years following hepatectomy. Median survival is presented as month (upper - lower interquartile range). doi:10.1371/journal.pone.0067756.g002</figcaption></figure>

however they were included within a unimodal age distribution and a clear distinction between two age-based subpopulations was not established [23].

In Peru, a large majority of patients were diagnosed with a tumor exceeding 10 cm in diameter. The consistent detection of advanced HCC may reflect a restricted access to health care in the country [8]. However, this situation is most likely due to several features intrinsic to the type of liver disease that affects Peruvian patients. Amazingly, proportions of HCC patients with and without liver cirrhosis in Peru (11% and 89%, respectively) are inverted when compared to the European and North American standards, in which rates of up to 80-90% of cirrhosis-associated HCC are observed (Table 1) [3-6]. In addition, almost half the

patients lack a bona fide risk factor (Figure 2B). Both characteristics imply that in most Peruvian cases, HCC is diagnosed without any prodromes or long-lasting liver disease enabling the prediction of an ongoing tumorigenic process. As Peruvian HCC develops essentially in healthy liver, the clinical symptoms commonly associated with an advanced malignant tumor represent the inaugural signs of hepatic neoplasia in Peru. Such a situation is exceptional in Asia, Europe, and North America, whereas it is the presentation of a minority of HCC cases in Africa [2-6,9].

Current staging systems, e.g. the TNM and BCLC classifications, seem to be inadequate for the Peruvian context of HCC, as they have been validated by surgically oriented groups from the northern hemisphere treating patients in an utterly different
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Regional disparities in mean age of Peruvian patients with HCC. The map was generated using the incidence rates of patients admitted at INEN with HCC between 1997 and 2010. The regions are classified as follows: (1) Amazonas (N=20; mean age: 36.2±20.3); (2) Ancash (N=117; 49.7±22.3); (3) Apurimac (N=121; 31.6±17.9); (4) Arequipa (N=39; 56.6±22.3); (5) Ayacucho (N=161; 39.5±20.6); (6) Cajamarca (N=36; 50.2±19.4); (7) Callao (N=21; 52.3±20.7); (8) Cuzco (N=96; 33.8±16.6); (9) Huancavelica (N=32; 45.3±23.9); (10) Huanuco (N=68; 37.2±20.8); (11) Ica (N=43; 56.1±19.3); (12) Junin (N=180; 37.2±19.9); (13) La Libertad (N=68; 61.6±19.3); (14) Lambayeque (N=53; 56.8±19.8); (15) Lima (N=289; 48±21.5); (16) Loreto (N=21; 42.7±20.8); (17) Madre de Dios (N=2; ND); (18) Moquegua (N=5; ND); (19) Pasco (N=29; 47.3±19.6); (20) Piura (N=64; 58.2±18.5); (21) Puno (N=18; 55.4±13.7); (22) San Martin (N=26; 43.3±19.6); (23) Tacna (N=4; ND); (24) Tumbes (N=3; ND); and (25) Ucayali (N=25; 38.2±18.8). Mean age intervals (years) are established as follows: [31–38] (dark blue); [39–44] (light blue); [45–48] (yellow); [49–52] (pink); and [53–62] (red). ND: Mean age was not determined because the cohort of patients admitted at INEN during the period and originating from these regions was too small. Administrative regional limits from the Peruvian National Institute of Statistics and Informatics (INEI).
doi:10.1371/journal.pone.0067756.g003</figcaption></figure>

pathophysiological context [9,24,25]. The generally late detection of HCC in Peru resulted in merely 20% of the ≤44 patients being selected for tumor resection over the period (Table 1). The aggressiveness of HCC in younger patients was illustrated by the higher percentage of early recurrence and lower 5-year survival rate among ≤44 group of patients which underwent resection of large HCC, as compared to previous studies (Figure 2D,E and

Table 1) [26–30]. Such astounding figures highlight the urgent need for a population-wide screening program of HCC in Peru.

Interestingly, the ≤44 and >44 groups displayed age-specific pathophysiological characteristics (Figure 2A,B and Table 1). The ratios of affected men to affected women were significantly higher in the younger population than in the older one with 2 and 1.2, respectively (Figure 1C and Table 1). Withal, sex ratio in both subpopulations appeared to be more balanced than those observed elsewhere in the world (up to 7 in selected populations) [1,26]. Besides, metastasis-associated HCC occurred significantly more in the younger population and development of intrahepatic multinudolar tumors took place preferentially in the older group (Table 1). Albeit surprisingly low, the rate of cirrhosis was higher in the >44 group than in the ≤44 group with 17% and 5%, respectively (Figure 2B and Table 1) [3–6]. The biological base of such situations remains unexplained and it is difficult to determine whether the outcome is due to an explosive tumor process that does not require chronic insult of liver tissue to develop. Otherwise, it may be the consequence of a particular resistance of the Peruvian populations to cirrhotogenic conditions.

The seroAFP was exceedingly high in the ≤44 group compared with the >44 group (Figure 2A). The mean of AFP concentration for the ≤44 group reached 332,000 ng per millilitre, while AFP biosynthesis is usually augmented during HCC up to hundreds or thousands ng per millilitre [31–33]. AFP is transiently produced during embryonic development in the establishment of the endodermal cell lineages (e.g. hepatocytes) and its expression is reactivated during HCC development [31–37]. Carcinogenesis shares underlying molecular mechanisms with cellular differentiation [38–41]. Homeotic genes that function in cell specialization by regulating differentiation marker expression such as AFP are misregulated in cancer [42–45]. The tremendously high seroAFP monitored in the ≤44 group therefore supports the hypothesis that aberrant activity of the genes involved in the cell identity and homeostasis drives hepatocarcinogenesis in this patient subpopulation. This hypothesis is reinforced by the predominance of undifferentiated and poorly differentiated hepatic cancer cells encountered in the ≤44 group (Table 1).

Multiple factors can be responsible of HCC development, including genetics, infectious agents, comorbid conditions, and environmental exposures [6,46–51]. Each of these factors has a specific impact on the development of liver cancer and influences the natural course of the neoplasm, both at cellular and pathophysiological levels [33,51]. As in western Africa and eastern Asia, the major of risk factors identified in Peru is the infection with HBV (Figure 2B) [3–6]. A large majority of HBV strains found in Peru belongs to subgenotype F1 and present peculiar mutations and deletions in the viral polymerase and PreS regions, respectively [52]. However, the tumorigenic potential of these strains remains unknown [52]. Classically, the age at which HCC develops is closely related to the age of acquisition of infection with hepatitis and to the rate of viral replication [5]. In a case of an early infection with HBV and an active viral replication (e.g. East Asia), the age-specific distribution of HCC starts to increase at the age of 50 without reaching a plateau [5]. In the case of a later infection with HBV and a rapid decline of viral replication after adolescence (e.g. West Africa), the age-specific distribution of HCC represents a single peak starting between the ages of 45 and plateauing thereafter [5,23]. In cases where HBV infection occurs in adulthood, the HCC is rarely developed before the age of 50 years and the highest age-specific incidence rates are observed in people over the age of 75 years [5]. None of these age distribution patterns were observed for the Peruvian patients admitted at INEN (Figure 1B). Moreover, the HCC etiology appears to be
different between the ≤44 and >44 groups, as the percentage of people over age 44 infected with HBV was significantly lower than in the ≤44 group (Figure 2B and Table 1). In addition, 50% of the patients tested for hepatitis B infection were serologically negative for HBsAg (Figure 2B and Table 1). Thus, while the role of HBV infection as a major etiology of HCC development in Peru is certain, it is not always present and appears, therefore, to be insufficient in and of itself to explain the peculiarities of the explosive tumor process as it is currently observed in Peru.

The particularity of the HCC epidemic is reinforced by the geographical mapping of the patients according to their place of origin. This process resulted in a striking geographical pattern in which the younger patients were concentrated mostly in the inland regions of the Andes, while the coastal regions provided HCC cases with a more elevated age (Figure 3). As a consequence, geographically neighbouring regions could display a 25 years difference with regard to the mean age of HCC patients, as observed for example in both Apurimac/Arequipa and Huanuco/La Libertad tandems (Figure 3; region number 3/4 and 10/13, respectively). At the moment, environmental or ethnogeographical differences that may explain such dramatic variations in tumor presentation remain unknown. In order to confirm this data, the eventuality of a selection effect has to be discarded with a prospective field study. Indeed, the lack of cancer surveillance at the national level and the deficit of information about the patients' clinical pathway prior to arrival at INEN prevent the accurate evaluation of the real proportions of the ≤44 and >44 HCC populations in Peru. Uncharacterized biasing parameter(s) may result in enhancing the number of ≤44 patients hospitalized at INEN. However, while such selection effect could have impacted quantitatively the perception of the HCC context in our study, the

age-specific HCC pathophysiology in the ≤44 and >44 groups clearly marked the distinction between these two subpopulations, suggesting two different age-specific natural histories of HCC in the Peruvian patient population.

The bimodal age-based distribution, the young age of a distinct subgroup, the low rate of associated cirrhosis, the high seroAFP, and the geographical clustering of the cases by age delineate an unusual tumoral process ongoing in the Peruvian population. These findings support the hypothesis that a singular risk factor, environmental, infectious, or genetic, and possibly endemic to the region is driving liver tumorigenesis in a significant proportion of Peruvian HCC patients. Research programs embracing the various aspects of this intriguing disease are required to decipher the biological process ongoing in young Peruvian HCC patients with the aim to generate novel prediction and prevention tools and assays for the local public health authorities.

## Acknowledgments

The authors are grateful to Macarena Farías Loza for her leadership in aggregating the medical records; Agnès Marchio for critical discussions; Anastasia Brion, Brian Gadd, and David D. Parker Jr. for their suggestions in preparing the manuscript. SB and ED thank the Ando-Amazonian International Mix Laboratory of Life Chemistry (LMI-LAVI), IRD-UPCH, for its logistic support. ER thanks the Medical Laboratory of INEN for its logistic support.

## Author Contributions

Conceived and designed the experiments: SB PP ED ER. Performed the experiments: SB PP ER. Analyzed the data: SB PP SL JM MZ ER. Contributed reagents/materials/analysis tools: JM. Wrote the paper: SB PP ED.

## References

1. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, et al. (2008) Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer 127: 2893-2917.

2. Soerjomataram I, Lortet-Tieulent J, Parkin DM, Ferlay J, Mathers C, et al. (2012) Global burden of cancer in 2008: a systematic analysis of disability-adjusted life-years in 12 world regions. Lancet 380: 1840-1850.

3. El-Serag HB (2011) Hepatocellular carcinoma. N Engl J Med 365: 1118-1127.

4. Forner A, Llovet JM, Bruix J (2012) Hepatocellular carcinoma. Lancet 379: 1245-1255.

5. Yang JD, Roberts LR (2010) Hepatocellular carcinoma: A global view. Nat Rev Gastroenterol Hepatol 7: 448-458.

6. de Martel C, Ferlay J, Franceschi S, Vignat J, Bray F, et al. (2012) Global burden of cancers attributable to infections in 2008: a review and synthetic analysis. Lancet Oncol 13: 607-615.

7. Instituto Nacional de Enfermedades Neoplásicas (2011) Plan estratégico institucional 2011-2015. Resolución No. 260-2011-/INEN. 62 p. Spanish.

8. Robinson HM, Hort K (2012) Non-communicable diseases and health systems reform in low-and-middle-income countries. Pac Health Dialog 18: 179-190.

9. Forner A, Reig ME, de Lope CR, Bruix J (2010) Current strategy for staging and treatment: the BCLC uptake and future prospects. Semin Liver Dis 30: 61-74.

10. American Joint Committee on Cancer (2002) AJCC Cancer Staging Manual. 6th ed. New York, NY: Springer. 421 p.

11. Schneider PD (2004) Preoperative assessment of liver function. Surg Clin North Am 84: 355-373.

12. Hasegawa K, Kokudo N, Makuuchi M (2007) Surgical management of hepatocellular carcinoma. Liver resection and liver transplantation. Saudi Med 28: 1171-1179.

13. Tang YH, Wen TF, Chen X (2012) Resection margin in hepatectomy for hepatocellular carcinoma: a systematic review. Hepato-gastroenterology 59: 1393-1397.

14. Teschendorff AE, Naderi A, Barbosa-Morais NL, Caldas C (2006) PACK: profile analysis using clustering and kurtosis to find molecular classifiers in cancer. Bioinformatics 22: 2269-2275.

15. Ahmed SE, Goria MN, Hussein A (2008) Gamma mixture: bimodality, inflexions and L-moments. Commun Stat-Theor M 37: 1147-1161.

16. Wang J, Wen S, Symmans WF, Pusztai L, Coombes KR (2009) The bimodality index: a criterion for discovering and ranking bimodal signatures from cancer gene expression profiling data. Cancer Inform 7: 199-216.

17. Shapiro SS, Wilk MB (1965) An analysis of variance test for normality (complete samples). Biometrika 52: 591-611.

18. Hadi AS, Simonoff,JS (1993) Procedures for the identification of multiple outliers in linear models. J Am Stat Assoc 88: 1264-1272.

19. Kaplan EL, Meier P (1958) Nonparametric estimation from incomplete observations. J Am Stat Assoc 53: 457-481.

20. Mantel N (1966) Evaluation of survival data and two new order statistics arising in its consideration. Cancer Chemother Rep 50: 163-170.

21. Peto R, Peto J (1972) Asymptotically efficient rank invariant test procedures. J R Statist Soc A 135: 185-207.

22. Waniez P (2011) Université Bordeaux Segalen, UFR Sciences et Modélisation. Available: http://philcarto.free.fr. Accessed 2013 Jan 28. French.

23. Van Rensburg SJ, Cook-Mozaffari P, Van Schalkwyk DJ, Van der Watt JJ, Vincent TJ, et al. (1985) Hepatocellular carcinoma and dietary aflatoxin in Mozambique and Transkei. Brit J Cancer 51: 713-726.

24. Marrero JA (2006) Staging systems for hepatocellular carcinoma: should we all use the BCLC system? J Hepatol 44: 630-632.

25. Villanueva A, Hernandez-Gea V, Llovet JM (2012) Medical therapies for hepatocellular carcinoma: a critical view of the evidence. Nat Rev Gastroenterol Hepatol 10: 34-42.

26. Parkin DM, Bray F, Ferlay J, Pisani P (2005) Global cancer statistics, 2002. CA Cancer J Clin 55: 74-108.

27. Sherman M (2008) Recurrence of hepatocellular carcinoma. N Engl J Med 359: 2045-2047.

28. Choi GH, Han DH, Kim DH, Choi SB, Kang CM, et al. (2009) Outcome after curative resection for a huge (>or = 10cm) hepatocellular carcinoma and prognostic significance of gross tumor classification. Am. J Surg 198: 693-701.

29. Yamashita Y, Taketomi A, Shirabe K, Aishima S, Tsuijita E, et al. (2011) Outcomes of hepatic resection for huge hepatocellular carcinoma (≥10 cm in diameter). J Surg Oncol 104: 292-298.

30. Allemann P, Demartines N, Bouzourene H, Tempia A, Halkic N (2012) Long-term outcome after liver resection for hepatocellular carcinoma larger than 10 cm. World J Surg 36: 2684-2691.

31. Yang SH, Suh KS, Lee HW, Cho EH, Cho JY, et al. (2007) A revised scoring system utilizing serum alphafetoprotein levels to expand candidates for living donor transplantation in hepatocellular carcinoma. Surgery 141: 598-609.

32. Sherman M (2010) The resurrection of alphafetoprotein. J Hepatol 52: 939-940.

33. Hiotis SP, Rahbari NN, Villanueva GA, Klegar E, Luan W, et al. (2012) Hepatitis B vs. hepatitis C infection on viral hepatitis-associated hepatocellular carcinoma. BMC Gastroenterol 12: 64.

34. Dzidek MA, Andrews GK (1983) Tissue specificity of alpha-fetoprotein messenger RNA expression during mouse embryogenesis. EMBO J 2: 549-554.
35. Tyner AL, Godbout R, Compton RS, Tilghman SM (1990) The ontogeny of alpha-fetoprotein gene expression in the mouse gastrointestinal tract. J Cell Biol 110: 915-927.

36. Shiojiri N, Lemire JM, Fausto N (1991) Cell lineages and oval cell progenitors in rat liver development. Cancer Res 51: 2611-2620.

37. Gualdi R, Bossard P, Zheng M, Hamada Y, Coleman JR, et al. (1996) Hepatic specification of the gut endoderm *in vitro*: cell signalling and transcriptional control. Gene Dev 10: 1670-1682.

38. Reya T, Morrison SJ, Clarke MF, Weissman IL (2001) Stem cells, cancer, and cancer stem cells. Nature 414: 105-111.

39. Reya T, Clevers H (2005) Wnt signalling in stem cells and cancer. Nature 434: 843-850.

40. Müller FJ, Laurent LC, Kostka D, Ulitsky I, Williams R, et al. (2008) Regulatory networks define phenotypic classes of human stem cell lines. Nature 455: 401-405.

41. Wong DJ, Liu H, Ridky TW, Cassarino D, Segal E, et al. (2008) Module map of stem cell genes guides creation of epithelial cancer stem cells. Cell Stem Cell 2: 333-344.

42. Krumlauf R (1994) Hox genes in vertebrate development. Cell 78: 191-201.

43. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, et al. (2010) Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature 464: 1071-1076.

44. Bertani S, Sauer S, Bolotin E, Sauer F (2011) The noncoding RNA *Mistral* activates *Hoxa6* and *Hoxa7* expression and stem cell differentiation by recruiting MLL1 to chromatin. Mol Cell 43: 1040-1046.

45. Tabuse M, Ohta S, Ohasi Y, Fukaya R, Misawa A, et al. (2011) Functional analysis of HOXD9 in human glioma cancer stem cells. Mol Cancer 10: 60.

46. Liaw YF, Chu CM (2009) Hepatitis B virus infection. Lancet 373: 582-592.

47. Averhoff FM, Glass N, Holtzman D (2012) Global burden of hepatitis C: considerations for heathcare providers in United States. Clin Infect Dis (Suppl 1): 10-15.

48. Calle EE, Rodriguez C, Walker-Thurmond K, Thun MJ (2003) Overweight, obesity, and mortality from cancer in a prospectively studied cohort of US adults. N Engl J Med 348: 1625-1638.

49. Schuppan D, Afdhal NH (2008) Liver cirrhosis. Lancet 371: 838-851.

50. Mack TM (2001) The pan-asian paan problem. Lancet 357: 1638-1639.

51. Van Thiel DH, Ramadori G (2011) Non-viral causes of hepatocellular carcinoma. J Gastrointest Cancer 42: 191-194.

52. von Meltzer M, Vasquez S, Sun J, Wendt UC, May A, et al. (2008) A new clade of hepatitis B virus subgenotype F1 from Peru with unusual properties. Virus Genes 37: 225-230.