Volume 28, Number 2—February 2022
Research
Comparative Effectiveness of Coronavirus Vaccine in Preventing Breakthrough Infections among Vaccinated Persons Infected with Delta and Alpha Variants
, Eduardo Freire Rodrigues1, Vítor Borges, João P. Gomes, Carlos Sousa, José P. Almeida, André Peralta-Santos, Baltazar Nunes, and PT-COVID-19 Group2
Abstract
We developed a case–case study to compare mRNA vaccine effectiveness against Delta versus Alpha coronavirus variants. We used data on 2,097 case-patients with PCR-positive severe acute respiratory syndrome coronavirus 2 infections reported in Portugal during May–July 2021. We estimated the odds of vaccine breakthrough infection in Delta-infected versus Alpha-infected patients by using conditional logistic regression adjusted for age group and sex and matched by the week of diagnosis. We compared reverse-transcription PCR cycle threshold values by vaccination status and variant as an indirect measure of viral load. We found significantly higher odds of vaccine breakthrough infection in Delta-infected patients than in Alpha-infected patients (odds ratio 1.96 [95% CI 1.22–3.14]), suggesting lower effectiveness of the mRNA vaccines in preventing infection with the Delta variant. We estimated lower mean cycle threshold values for the Delta cases (mean difference −2.10 [95% CI −2.74 to −1.47]), suggesting higher infectiousness than the Alpha variant.
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.617.2 lineage, also known as Delta variant of concern (VOC), first sequenced in India in December 2020, was identified in Portugal in late April and quickly became dominant, reaching 90% of all sequenced cases in late June 2021 (epidemiologic week 26), just 2 months after it was first identified (1). Available evidence suggests that this VOC is associated with higher transmissibility, higher risk for hospitalization, and reduced antibody neutralization compared with other VOCs (2,3).
Vaccination is the primary pharmacologic measure to control the transmission of SARS-CoV-2 and mitigate its effect on hospitalizations and mortality rates. In Portugal, vaccination was initiated in late December 2020 for those at higher risk for severe disease or exposure and since February 2021 has been rolled out by descending age criteria. By week 26 (June 28–July 4), 36% of the population was fully vaccinated, and 56% had started or completed vaccination (Appendix Figure 1), most (75%) with mRNA vaccines (BNT162b2 [Pfizer-BioNTech, https://www.pfizer.com] or mRNA-1273 [Moderna, https://www.modernatx.com]) administrated with a 28-day dose interval (4).
Early reports of vaccine effectiveness indicate a high protection for mRNA vaccines against infection and disease (5,6) and a reduced viral load in the vaccinated case-patients (7,8). However, reports of vaccine effectiveness against Delta have shown decreased protection of the vaccines compared with the Alpha variant (2,9). Validating this potential reduction of vaccine effectiveness against the Delta VOC is critical to inform further public health measures, particularly as the variant becomes globally dominant.
We aimed to provide a measure of comparative effectiveness of mRNA vaccines (BNT162b2 and mRNA-1273) against B.1.617.2 (Delta) versus B.1.1.7 (Alpha) VOCs, using a case–case study design. As a secondary objective, we compared reverse-transcription PCR (RT-PCR) cycle threshold (Ct) values between vaccine status for Alpha and Delta variants as an indirect measure of viral load and, thus, transmissibility of the vaccine breakthrough cases for both variants.
Study Design
We developed an observational case–case study (10) comparing odds of vaccination (partial or complete) between RT-PCR–positive cases (symptomatic or asymptomatic) classified as infected with Delta versus Alpha VOCs. The study period was May 17–July 4, 2021 (epidemiologic weeks 20–26), to cover the period of VOC replacement in Portugal, from the Alpha (84.8%, week 19) to Delta dominance (96.1%, week 27) (1). Our analysis included persons with data on whole-genome sequencing (WGS) or spike (S) gene target failure (SGTF) who were >40 years of age and eligible for vaccination during the study period. Persons for whom data on national health registry number, age, sex, or diagnosis date were missing, and those vaccinated with Ad26.COV2-S (Johnson & Johnson/Janssen, https://www.janssen.com) or ChAdOx1 nCoV-19 (AstraZeneca, https://www.astrazeneca.com) vaccines were excluded from the study.
To indirectly infer the level of infectiousness of case-patients according to vaccination status and VOC type, we performed a secondary analysis by comparing the paired means of RT-PCR Ct values for nucleocapsid and open reading frame 1ab genes, by using data from a single large laboratory (Unilabs, https://unilabs.com). Lower Ct values reflect a reduced number of RT-PCR cycles required for amplification of SARS-CoV-2 RNA and, therefore, a higher number of virus copies within the sample. As such, studies have used Ct values to estimate viral load or viral shedding (8,11).
Data Sources
SARS-CoV-2 Cases
RT-PCR testing for SARS-CoV-2 in Portugal is done by hospitals as well as by public and private laboratories, and is available free of charge to anyone with symptoms consistent with coronavirus disease (COVID-19) (12). Laboratory-confirmed cases are reported to the mandatory National Epidemiologic Surveillance Information System (Sistema Nacional de Vigilância Epidemiológica, https://www.spms.min-saude.pt/2020/07/sinave-2). For this study, each notifying laboratory selected a random subset of RT-PCR–positive nasopharyngeal samples collected during the study period to be sent to the National SARS-CoV-2 Genomic Surveillance Network (1,13) and, thus, to be included in the study. We also included samples from a private molecular biology laboratory (Unilabs) with nationwide coverage that routinely performs analysis on SGTF. We collected RT-PCR Ct values as an indirect measure of viral load (11). We removed duplicate records on the basis of national health register numbers, maintaining only the first collected sample.
Variant Classification
We classified SARS-CoV-2 variants by viral WGS or inferred by SGTF data. For nonsequenced samples, we considered S-positive specimens (with amplification of structural gene) as Delta and SGTF samples as Alpha by using the TaqPath COVID-19 CE IVD RT-PCR Kit (Thermo Fisher Scientific, https://www.thermofisher.com) that targets 3 genes (structural, nucleocapsid, and open reading frame 1ab), performed according to the manufacturer’s specifications, as described elsewhere (13).
Vaccination Status, Demographics, and Data Linkage
We obtained COVID-19 vaccination status through the electronic national vaccination register (https://www.sns.gov.pt/monitorizacao-do-sns/vacinas-covid-19). We classified vaccination exposure as no register of vaccine administration before diagnosis (i.e., unvaccinated); SARS-CoV-2 infection diagnosis <14 days after first dose mRNA vaccination (1 dose [<14 days]); SARS-CoV-2 infection diagnosis >14 days after first dose or <14 days after second dose (1 dose [>14 days] or 2 doses [<14 days]) (i.e., partial vaccination); and >14 days after second dose of the mRNA vaccine (2 doses [>14 days]) (i.e., complete vaccination). Information about age, sex, and date at diagnosis was routinely collected by National Epidemiologic Surveillance Information System. We performed a deterministic record linkage to join all data sources, namely on vaccination status, outcomes (VOC classification), and other covariates (e.g., age group and sex) and to remove duplicate data from the dataset.
Statistical Analysis
We compared characteristics of delta and alpha SARS-CoV-2 case-patients by using the χ2 test. We considered Delta-infected case-patients as case-patients of interest and Alpha-infected case-patients as the reference group. We used conditional logistic regression matched by the week of diagnosis and adjusted for age group and sex to estimate confounder-adjusted odds of having been infected by SARS-CoV-2 and vaccinated among Delta case-patients compared with Alpha case-patients. These covariates might be associated with the probability of having been vaccinated and being exposed to the virus and type of variant.
In our analysis, odds ratio (OR) = 1 indicates no difference in odds of having been infected by SARS-CoV-2 and vaccinated and, thus, a proxy of no difference between mRNA vaccine effectiveness against the Delta versus Alpha VOC. OR >1 indicates a higher odds of having been infected by SARS-CoV-2 and vaccinated, thus lower vaccine effectiveness against the Delta versus Alpha VOC, whereas OR <1 indicates a lower odds of having been infected by SARS-CoV-2 and vaccinated among Delta case-patients and a higher vaccine effectiveness against the Delta versus Alpha VOC (Appendix).
We stratified mean and SD Ct values for Alpha-infected and Delta-infected case-patients on the basis of vaccination status. We evaluated differences between mean Ct values by vaccination status and VOC by fitting a linear multiple regression model with Ct values as outcome, adjusting for sex, age group, and week of case diagnosis. We included an interaction term between vaccination status and VOC type in the regression model to determine whether the effect of vaccination status on Ct values differed between Delta and Alpha.
Sensitivity Analysis
To assess the change of the sampling strategy for WGS from a monthly to weekly basis, which occurred on week 21, we restricted our analysis to weeks 22–26. In addition, to assess the bias of misclassification error associated with the SGTF method (particularly in the early weeks of the study period, when overall prevalence of the Delta variant was lower and SGTF sensitivity may also have been lower), we analyzed samples identified exclusively through WGS during weeks 22–26. Finally, to address whether having been infected and vaccinated was associated with lower infectiousness in any of the studied VOCs, we restricted analysis to samples with Ct values <25 (8).
Ethics Considerations
Genomic surveillance of SARS-CoV-2 in Portugal is regulated by Assistant Secretary of State and Health Executive Order no. 331/2021, issued on January 11, 2021. The research on genomic epidemiology of SARS-COV-2 received the clearance of the Ethics Committee of Instituto Nacional de Saúde Doutor Ricardo Jorge on March 30, 2021.
Main Analysis
A total of 22,784 SARS-CoV-2–positive cases were reported in Portugal during May 17–July 4, 2021, among persons >40 years of age. Of 2,097 cases included in the analysis, 966 (46.1%) were variant-classified with WGS and 1,131 (53.9%) with SGTF. During the study period, 94.7% (827/873) of the S-positive sequenced samples were confirmed as Delta and 96.9% (372/384) of SGTF samples were classified as Alpha through WGS, thus indicating that the SGTF-derived VOC classification was robust.
Among Delta case-patients, we observed a higher proportion of persons >70 years of age (p<0.001) (Table 1), and a higher proportion of vaccinated persons (p<0.001) than among the Alpha case-patients. We report a statistically significant higher odds of being partially vaccinated (OR 1.70 [95% CI 1.18–2.47]) or completely vaccinated (OR 1.96 [95% CI 1.22–3.14]) among the Delta case-patients than among the Alpha case-patients, suggesting lower mRNA vaccine effectiveness for the Delta variant (Table 2). After adjustment for age group and sex, similar estimated ORs were observed for the complete vaccination scheme (OR 1.96 [95% CI 1.43–2.69]) or for partial vaccination (OR 1.81 [95% CI 1.37–2.39]).
Secondary Analysis
Figure
Figure. Distribution of Ct paired mean (reverse-transcription PCR Ctvalues for nucleocapsid and open reading frame 1ab genes) by coronavirus variant and vaccination status, Portugal, epidemiologic weeks 20–26...
We observed statistically significant higher mean Ct values among those with complete vaccination compared with unvaccinated case-patients for Delta (17.7 vs. 16.5 ) as well as for Alpha (21.8 vs. 18.4) (Table 3; Figure), suggesting lower viral loads in vaccinated compared with unvaccinated case-patients for both VOCs. Although the Alpha variant cases had statistically significant confounder-adjusted Ct values mean difference (MD) of 4.49 (95% CI 2.07–6.91) after complete vaccination, representing an increase of Ct values (lower infectiousness), the Delta variant cases showed only about half of that increase, with a statistically significant confounder-adjusted Ct value mean difference point estimate of 2.24 (95% CI 0.85–3.64) between unvaccinated and fully vaccinated breakthrough case-patients. For partial vaccination, statistically significant differences in mean Ct values were observed for Alpha (MD 1.87 [95% CI 0.2–3.53]) but not for Delta cases (MD −0.15 [95% CI −0.99 to 0.96]), suggesting similar viral load between unvaccinated and partially vaccinated Delta case-patients.
The confounder-adjusted paired-mean difference in Ct values between Delta and Alpha was statistically significant for unvaccinated persons (MD –1.66 [95% CI –2.37 to –0.95]). For partially (MD –1.88 [95% CI –3.77 to 0.003]) and completely vaccinated persons (MD –2.24 [95% CI –4.8 to 0.32]) we did not observed statistically significant mean differences in Ct values between delta and alpha.
Sensitivity Analysis
In a sensitivity analysis (Table 4), restricted to weeks 22–26, confounder-adjusted OR estimates for having been infected and completely vaccinated (OR 1.78 [95% CI 1.01–3.13]) and partially vaccinated (OR 1.91 [95% CI 1.22–3.00]) remained similar and statistically significant. Restricting the analysis to the cases identified through WGS, we observed a drop in the adjusted OR point estimate of complete vaccination (OR 1.48 [95% CI 0.75–2.93]) with a loss of statistical significance, but for partial vaccination, estimates yielded increased significant results (OR 2.5 [95% CI 1.23–5.08]). Restricting to cases with Ct values <25 (n = 1,363), we observed an increase of the OR of infection and vaccination against Delta versus Alpha variants both for partial (OR 2.47 [95% CI 1.48–4.12]) and complete vaccination (OR 2.42 [95% CI 1.06–5.51).
We observed statistically significant higher odds of having been infected and vaccinated (vaccine infection breakthrough) among Delta-infected versus Alpha-infected case-patients, suggesting a lower mRNA vaccines effectiveness for SARS-CoV-2 infection with the Delta VOC. The findings were consistent for both complete and partial vaccination. Delta breakthrough case-patients have a higher viral load (lower Ct values) compared with Alpha breakthrough case-patients.
The OR estimates for complete vaccination (vaccine infection breakthrough) (1.96) are in line with findings of test-negative design studies on vaccine effectiveness from Scotland and England (2,9) based only on SGTF or mixed SGTF and WGS methodology for variant identification. Those studies reported a 5.9 percentage point reduction of BNT162b2 vaccine effectiveness against the Delta VOC compared with 13.0 percentage points for the alpha VOC for complete vaccination (2,9), with nonoverlapping CIs for vaccine effectiveness estimates.
For partial vaccination, our results indicated statistically significant lower mRNA vaccine effectiveness against the Delta VOC (OR = 1.8), supporting the need to promptly complete vaccination schedules to account for swiftly reduced effectiveness against this variant. This result is not in line with previous research conducted in England (9). Several factors may explain the differences between that study (9) and our work: the target population is different (persons >16 years of age vs. >40 years of age); differences exist in the Alpha and Delta relative frequencies during the study period; differences exist in vaccination calendar and time between doses administration, with England having a higher proportion of persons exposed to 2 doses because of an earlier campaign roll out and Portugal having shorter time of exposure to a single dose; and a sample selection methodology for WGS.
On our secondary analysis, we observed lower Ct values (indicative of higher viral loads) among Delta compared with Alpha case-patients (MD −2.10 [95% CI −2.74 to −1.47]). Furthermore, although complete vaccination increases Ct values (thus reducing estimated viral loads) in Alpha case-patients by 4.49 (95% CI 2.07–6.91), for the Delta variant, complete vaccination had a much lower increase of Ct values (2.24 [95% CI 0.85–3.64]), only half of the difference observed for Alpha case-patients. These findings are consistent with vaccine infection breakthrough cases; Delta variant case-patients have higher infectiousness than alpha variant case-patients. Our findings were similar to results of the studies performed in Israel, when Alpha variant was the dominant variant, that also found a 5.09 (95% CI 2.8–7.4) increase of nucleocapsid gene Ct values (8) between completely vaccinated and unvaccinated cases and 1.51 for those partially vaccinated with BNT162b2 (7). Those findings are consistent with our results for the Alpha variant.
Our results were robust to variations in the sampling strategy, including changing the weeks of diagnosis included. Restricting analysis to samples with higher viral loads (selected on the basis of low Ct values [Ct <25]), we observed an even higher OR of vaccine infection breakthrough between Delta and Alpha variant cases, which supports the relevance of the Delta relative vaccine effectiveness reduction effect in the transmission of SARS-CoV-2, both for partial (OR 2.47 [95% CI 1.48–4.12]) and complete vaccination (OR 2.42 [95% CI 1.06–5.51]). The restriction to WGS samples only, however, resulted in a loss in a sample size and lack of statistical power to detect differences for complete vaccination.
This study has several limitations; in particular, possible misclassification and selection and confounding bias should be acknowledged. Although we observed a high positive predictive value (94.7%) of non-SGTF data to identify Delta cases within the study sample, a misclassification error may have led to underestimating the reported effect. However, SGTF methodology has previously shown good classification accuracy to identify B.1.1.7 in Portugal (13) and has more recently been shown to distinguish between Alpha and Delta VOCs in Scotland and England (2,9) and may be highly useful when a large-scale testing strategy is in place and electronic vaccination registers are used. In addition, information bias caused by misclassification of vaccination status can arise from delays in registering vaccination status and diagnosis. However, data extraction was performed 3 weeks after the end of the study period to minimize bias (e.g., vaccination centers are instructed to register vaccinations up to 24 hours after administration). Our sampling strategy resulted in some selection bias. Overall, 22,784 cases were identified in Portugal during the study period among those >40 years of age, and their age distribution was different from the study sample (Appendix Figure 4), skewed toward younger ages in the study sample, possibly because our sample was collected mainly through ambulatory laboratories, whereas older persons (>80 years of age) are expected to be more frequently diagnosed by hospital laboratories. This result could bias our estimates if the reduction of vaccine effectiveness between Delta and Alpha variant is age-dependent. For example, if the vaccine effectiveness reduction is higher among older persons, our results could be underestimated. Moreover, Delta cases occurred more frequently among older participants, and older participants were the first to be vaccinated in Portugal (4) and, thus, had a longer exposure time after their second dose. With time, a waning of the vaccine effect can occur, possibly contributing to the observed differences among persons who were completely vaccinated.
After adjusting for confounding, we did not observe a substantial change in OR estimates. Although we cannot exclude residual confounding bias, given that other factors not accounted for might be associated with the probability of exposure to the virus and of being vaccinated (e.g., health and social services worker status, ethnicity, and education). However, we found that other studies (9) that adjusted for these potential confounding factors not accounted in our study did not observe a substantial difference between crude and confounder-adjusted estimates. Hence, if any residual confounding occurred, it could have a small effect on the effect estimates.
The results must always be interpreted in context because they do not provide evidence to question the benefits of the mRNA vaccines to individual health, such as reducing symptoms, disease severity, or the impact on health services capacity. We reported odds of vaccine breakthrough between Delta and Alpha VOCs, which can be interpreted as a measure of the relative vaccine effectiveness. Although a case–case study design does not provide a direct measure of effectiveness against a specific VOC, it may be useful for rapidly detect changes in vaccine effectiveness in the context of novel VOC emergence, providing substantial evidence for further public health measures to control the transmission of SARS-CoV-2.
Overall, we found significantly higher odds of vaccination in Delta case-patients than in Alpha case-patients, suggesting possible lower effectiveness of the mRNA vaccines in preventing infection with the Delta VOC. Case–case design has proven to be helpful to compare vaccine effectiveness for SARS-CoV-2 VOCs because of its quick implementation and valuable insights in the context of frequent and swift VOC emergence. These findings can help decision-makers as they consider applying or lifting of control measures and adjusting vaccine roll-out depending on the predominance of the Delta variant and levels of partial and complete mRNA vaccination coverage.
Ms. Kislaya is a biostatistician in the Department of Epidemiology at the Instituto Nacional de Saúde Doutor Ricardo Jorge, Lisbon, Portugal. Her primary research interests include vaccine effectiveness and population-level impact of public health interventions.
Members of the PT-COVID-19 Group: Ausenda Machado, Ana Paula Rodrigues, Carlos Matias Dias, Rita Ferreira, Joana Isidro, Miguel Pinto, Sílvia Duarte, Luís Vieira, Pedro Pinto Leite, Joana Moreno, António Maia Gonçalves, and Luís Menezes.
Acknowledgment
E.F. Rodrigues is a board member and owns stocks at UPHILL, a software company that provides digital training solutions and has customers in the pharmaceutical sector (e.g., Pfizer). No business was conducted between the 2 entities regarding mRNA vaccine products, similar products, or anything pertaining to this study.
References
- Instituto Nacional de Saúde Doutor Ricardo Jorge. Genetic diversity of the novel coronavirus SARS-CoV-2 (COVID-19) in Portugal. 2021 [cited 2021 Aug 13]. https://insaflu.insa.pt/covid19
- Sheikh A, McMenamin J, Taylor B, Robertson C; Public Health Scotland and the EAVE II Collaborators. SARS-CoV-2 Delta VOC in Scotland: demographics, risk of hospital admission, and vaccine effectiveness. Lancet. 2021;397:2461–2. DOIPubMedGoogle Scholar
- Wall EC, Wu M, Harvey R, Kelly G, Warchal S, Sawyer C, et al. Neutralising antibody activity against SARS-CoV-2 VOCs B.1.617.2 and B.1.351 by BNT162b2 vaccination. Lancet. 2021;397:2331–3. DOIPubMedGoogle Scholar
- Direção-Geral da Saúde. Norma circular no. 002/2021. Vaccination campaign against COVID-19 [in Portuguese]. 2021 [cited 2021 Aug 13]. https://covid19.min-saude.pt/wp-content/uploads/2021/10/Norma_002_2021_act_15_10_2021.pdf
- Haas EJ, Angulo FJ, McLaughlin JM, Anis E, Singer SR, Khan F, et al. Impact and effectiveness of mRNA BNT162b2 vaccine against SARS-CoV-2 infections and COVID-19 cases, hospitalisations, and deaths following a nationwide vaccination campaign in Israel: an observational study using national surveillance data. Lancet. 2021;397:1819–29. DOIPubMedGoogle Scholar
- Kissling E, Hooiveld M, Sandonis Martín V, Martínez-Baz I, William N, Vilcu A-M, et al. I-MOVE-COVID-19 primary care study team; I-MOVE-COVID-19 primary care study team. Vaccine effectiveness against symptomatic SARS-CoV-2 infection in adults aged 65 years and older in primary care: I-MOVE-COVID-19 Project, Europe, December 2020 to May 2021. Euro Surveill. 2021;26:
2100670 . DOIGoogle Scholar - Levine-Tiefenbrun M, Yelin I, Katz R, Herzel E, Golan Z, Schreiber L, et al. Initial report of decreased SARS-CoV-2 viral load after inoculation with the BNT162b2 vaccine. Nat Med. 2021;27:790–2. DOIPubMedGoogle Scholar
- Regev-Yochay G, Amit S, Bergwerk M, Lipsitch M, Leshem E, Kahn R, et al. Decreased infectivity following BNT162b2 vaccination: A prospective cohort study in Israel. Lancet Reg Health Eur. 2021;7:
100150 . DOIPubMedGoogle Scholar - Lopez Bernal J, Andrews N, Gower C, Gallagher E, Simmons R, Thelwall S, et al. Effectiveness of Covid-19 vaccines against the B.1.617.2 (delta) variant. N Engl J Med. 2021;385:585–94. DOIPubMedGoogle Scholar
- Pogreba-Brown K, Austhof E, Ellingson K. Methodology minute: An overview of the case-case study design and its applications in infection prevention. Am J Infect Control. 2020;48:342–4. DOIPubMedGoogle Scholar
- Rabaan AA, Tirupathi R, Sule AA, Aldali J, Mutair AA, Alhumaid S, et al. Viral dynamics and real-time RT-PCR Ct values correlation with disease severity in COVID-19. Diagnostics (Basel). 2021;11:1091. DOIPubMedGoogle Scholar
- Direção-Geral de Saúde. Norma circular no. 004/2020. Approach to suspected and confirmed COVID-19 cases [in Portuguese]. 2020 [cited 2021 Aug 13]. https://covid19.min-saude.pt/wp-content/uploads/2021/04/Norma_004_2020_act_19_04_2021.pdf
- Borges V, Sousa C, Menezes L, Gonçalves AM, Picão M, Almeida JP, et al. Tracking SARS-CoV-2 lineage B.1.1.7 dissemination: insights from nationwide spike gene target failure (SGTF) and spike gene late detection (SGTL) data, Portugal, week 49 2020 to week 3 2021. Euro Surveill. 2021;26:1–6. DOIPubMedGoogle Scholar
Figure
Tables
Cite This ArticleOriginal Publication Date: December 07, 2021
1These first authors contributed equally to this article.
2Members of the PT-Covid-19 Group are listed at the end of this article.
Table of Contents – Volume 28, Number 2—February 2022
| EID Search Options |
|---|
Advanced Article Search – Search articles by author and/or keyword.
|
Articles by Country Search – Search articles by the topic country.
|
Article Type Search – Search articles by article type and issue.
|
Please use the form below to submit correspondence to the authors or contact them at the following address:
Irina Kislaya, Instituto Nacional de Saúde Doutor Ricardo Jorge, Av Padre Cruz, 1649-016 Lisbon, Portugal
Top