Background:

The potential effect of alcohol or tea intake on the risk of nasopharyngeal carcinoma (NPC) remains controversial.

Methods:

In a population-based case–control study in southern China, we assessed alcohol or tea intake from 2,441 histopathologically confirmed NPC cases and 2,546 controls. We calculated mean daily ethanol (g/day) and tea intake (mL/day). Fully adjusted ORs with 95% confidence intervals (CI) were estimated using logistic regression; potential dose–response trends were evaluated using restricted cubic spline analysis.

Results:

Compared with nondrinkers, no significantly increased NPC risk in men was observed among current alcohol drinkers overall (OR, 1.08; 95% CI, 0.93–1.25), nor among current heavy drinkers (OR for ≥90 g/day ethanol vs. none, 1.32; 95% CI, 0.95–1.84) or former alcohol drinkers. Current tea drinking was associated with a decreased NPC risk (OR, 0.73; 95% CI, 0.64–0.84). Compared with never drinkers, those with the low first three quintiles of mean daily current intake of tea were at significantly lower NPC risk (OR, 0.53, 0.68, and 0.65, respectively), but not significant for the next two quintiles. Current daily tea intake had a significant nonlinear dose–response relation with NPC risk.

Conclusions:

Our study suggests no significant association between alcohol and NPC risk. Tea drinking may moderately reduce NPC risk, but the lack of a monotonic dose–response association complicates causal inference.

Impact:

Tea drinking might be a healthy habit for preventing NPC. More studies on biological mechanisms that may link tea with NPC risk are needed.

The incidence of nasopharyngeal carcinoma (NPC) varies substantially by geographic location. Recent data from the International Agency for Research on Cancer showed that 82.9% of incident NPC cases in 2018 occurred in East Asia, South-East Asia, and South-Central Asia (1). Studies especially among high-risk areas that can clarify the associations of modifiable risk factors with the risk of NPC are of potential importance for primary cancer prevention (2–6). However, so far there are still few well-established modifiable risk factors for NPC, except for smoking, poor oral health, and perhaps some dietary components (7–9).

Previously reported associations between alcohol or tea intake and the risk of NPC remain inconsistent. Most case–control studies reported no significant association with alcohol intake (10–15), whereas two systematic reviews and some other studies showed that alcohol drinking had a J-shaped dose–response trend: light or moderate alcohol drinkers exhibited a significant or nonsignificant decreased NPC risk (16), while heavy drinkers had an increased risk (17–22). Several studies investigated the potential association between tea drinking and the risk of NPC: studies in Guangdong and Taiwan or in Los Angeles (14, 17, 23, 24) reported that tea drinkers had a decreased NPC risk, whereas one study in Turkey found an increased risk among tea drinkers who consumed <10 glasses per day (10).

However, most of these studies lacked detailed evaluation of intake of alcohol and tea, or they did not collect precise information on different types of alcoholic beverages or tea (2, 10, 19). This information might be important for understanding the dose–response pattern of a potential causal association between consumption of ethanol or tea and NPC risk. To this end, we took advantage of our NPC Genes, Environment, and Epstein-Barr virus (NPCGEE) study, a large population-based case–control study conducted in southern China, to investigate associations between detailed alcohol and tea drinking over time and NPC risk.

Study design and population

NPCGEE is a population-based case–control study conducted in three areas of southern China. The rationale and design of the NPCGEE study have been described in detail elsewhere (8, 9, 25, 26). Briefly, between 2010 and 2014, we enrolled a total of 2,554 incident NPC cases and 2,648 population-based controls. All histopathologically diagnosed cases of NPC were ascertained by a rapid reporting system. Controls were randomly selected every 6 to 12 months from total population registries covering three geographic areas, with frequency matching to the 5-year age and sex distribution of the cases by geographic region. The overall participation rates for cases and controls were 83.8% and 82.7%, respectively. Inclusion criteria for both cases and controls were age 20 to 74 years, current residence in the study area during the recruitment period, and no previous malignant disease or congenital or acquired immunodeficiency. All eligible participants were interviewed in-person by trained interviewers using an electronic questionnaire to collect data on individual lifestyle factors.

Institutional review boards or ethics vetting boards from all study centers approved this study. Eligible subjects granted written or oral informed consent.

All eligible cases and controls with available alcohol drinking and tea drinking information were included in this analyses. We excluded subjects with data lost during the upload process (n = 18), missing information on alcohol drinking or tea drinking (n = 57), and those with interview of poor quality or missing on confounding factors listed in Table 1 (n = 140). The final analytic population included 2,441 NPC cases and 2,546 controls (Supplementary Fig. S1).

Table 1.

Characteristics of cases of nasopharyngeal carcinoma and control subjects, overall and stratified by sex.

OverallMenWomen
CharacteristicsCases n = 2,441 n (%)Controls n = 2,546 n (%)PaCases n = 1,785 n (%)Controls n = 1,869 n (%)PaCases n = 656 n (%)Controls n = 677 n (%)Pa
Residential area   0.418   0.418   0.935 
 Zhaoqing 1,227 (50.3) 1,284 (50.4)  889 (49.8) 933 (49.9)  338 (51.5) 351 (51.9)  
 Wuzhou 665 (27.2) 658 (25.8)  488 (27.3) 481 (25.7)  177 (27.0) 177 (26.1)  
 Guiping/Pingnan 549 (22.5) 604 (23.7)  408 (22.9) 455 (24.3)  141 (21.5) 149 (22.0)  
Age at interview, years   0.003   0.018   0.202 
 20–39 491 (20.1) 446 (17.5)  345 (19.3) 309 (16.5)  146 (22.3) 137 (20.2)  
 40–49 878 (36.0) 869 (34.1)  622 (34.9) 620 (33.2)  256 (39.0) 249 (36.8)  
 50–59 662 (27.1) 715 (28.1)  506 (28.4) 554 (29.6)  156 (23.8) 161 (23.8)  
 60–74 410 (16.8) 516 (20.3)  312 (17.5) 386 (20.7)  98 (14.9) 130 (19.2)  
Educational level, years   0.002   0.018   0.074 
 ≤6 1,880 (39.6) 967 (35.9)  574 (32.2) 539 (28.8)  393 (59.9) 374 (55.2)  
 7–9 1,996 (40.0) 977 (40.0)  788 (44.2) 819 (43.8)  189 (28.8) 200 (29.5)  
 ≥10 1,111 (20.4) 497 (24.1)  423 (23.7) 511 (27.3)  74 (11.3) 103 (15.2)  
Current housing typeb   <0.001   0.029   <0.001 
 Building 1,766 (72.4) 1,985 (78.0)  1,301 (72.9) 1,421 (76.0)  465 (70.9) 564 (83.3)  
 Cottage/boat 675 (27.7) 561 (22.0)  484 (27.1) 448 (24.0)  191 (29.1) 113 (16.7)  
Current occupation   <0.001   <0.001   0.035 
 Unemployed/farmer 903 (37.0) 1,063 (41.8)  546 (30.6) 668 (35.7)  357 (54.4) 395 (58.4)  
 Blue-collar 990 (40.6) 882 (34.6)  815 (45.7) 742 (39.7)  175 (26.7) 140 (20.7)  
 White-collar 336 (13.8) 410 (16.1)  266 (14.9) 318 (17.0)  70 (10.7) 92 (13.6)  
 Other/unknown 212 (8.7) 191 (7.5)  158 (8.9) 141 (7.5)  54 (8.2) 50 (7.4)  
Cigarette smoking   0.133   0.020   0.412 
 Never 1,093 (44.8) 1,194 (46.9)  450 (25.2) 535 (28.6)  643 (98.0) 659 (97.3)  
 Ever 1,348 (55.2) 1,352 (53.1)  1,335 (74.8) 1,334 (71.4)  13 (2.0) 18 (2.7)  
Salt-preserved fish consumption in 2000–2002  0.020   0.268   0.016 
 ≤Yearly 1,787 (73.2) 1,873 (73.6)  1,296 (72.6) 1,359 (72.7)  491 (74.9) 514 (75.9)  
 Monthly 469 (19.2) 527 (20.7)  352 (19.7) 390 (20.9)  117 (17.8) 137 (20.2)  
 ≥Weekly 185 (7.6) 146 (5.7)  137 (7.7) 120 (6.4)  48 (7.3) 26 (3.8)  
Nasopharyngeal carcinoma among first-degree relatives  <0.001   <0.001   <0.001 
 No 2,156 (88.3) 2,457 (96.5)  1,577 (88.4) 1,809 (96.8)  579 (88.3) 648 (95.7)  
 Yes 264 (10.8) 65 (2.6)  192 (10.8) 46 (2.5)  72 (11.0) 19 (2.8)  
 Unknown 21 (0.9) 24 (0.9)  16 (0.9) 14 (0.8)  5 (0.8) 10 (1.5)  
Body mass index at age 20 years, kg/m2   0.001   0.006   0.044 
 <18.5 394 (16.1) 515 (20.2)  275 (15.4) 357 (19.1)  119 (18.1) 158 (23.3)  
 18.5–22.9 1,763 (72.2) 1,739 (68.3)  1,326 (74.3) 1,306 (69.9)  437 (66.6) 433 (64.0)  
 ≥23.0 284 (11.6) 292 (11.5)  184 (10.3) 206 (11.0)  100 (15.2) 86 (12.7)  
Vegetable consumption in 2000–2002, g/day  <0.001   <0.001   <0.001 
 <252.0 682 (27.9) 715 (28.1)  505 (28.3) 510 (27.3)  177 (27.0) 205 (30.3)  
 252.0–335.9 617 (25.3) 829 (32.6)  450 (25.2) 588 (31.5)  167 (25.5) 241 (35.6)  
 ≥336.0 1,142 (46.8) 1,002 (39.4)  830 (46.5) 771 (41.3)  312 (47.6) 231 (34.1)  
Frequency of tooth brushing, times/day   <0.001   <0.001   <0.001 
 ≤1 1,626 (66.6) 1,408 (55.3)  1,233 (69.1) 1,072 (57.4)  393 (59.9) 336 (49.6)  
 2 733 (30.0) 1,036 (40.7)  497 (27.8) 726 (38.8)  236 (36.0) 310 (45.8)  
 ≥3 82 (3.4) 102 (4.0)  55 (3.1) 71 (3.8)  27 (4.1) 31 (4.6)  
Frequency of herbal soup consumption   <0.001   <0.001   0.290 
 ≤Yearly 671 (27.5) 585 (23.0)  498 (27.9) 430 (23.0)  173 (26.4) 155 (22.9)  
 Monthly 961 (39.4) 1,122 (44.1)  689 (38.6) 819 (43.8)  272 (41.5) 303 (44.8)  
 ≥Weekly 809 (33.1) 839 (33.0)  598 (33.5) 620 (33.2)  211 (32.2) 219 (32.4)  
OverallMenWomen
CharacteristicsCases n = 2,441 n (%)Controls n = 2,546 n (%)PaCases n = 1,785 n (%)Controls n = 1,869 n (%)PaCases n = 656 n (%)Controls n = 677 n (%)Pa
Residential area   0.418   0.418   0.935 
 Zhaoqing 1,227 (50.3) 1,284 (50.4)  889 (49.8) 933 (49.9)  338 (51.5) 351 (51.9)  
 Wuzhou 665 (27.2) 658 (25.8)  488 (27.3) 481 (25.7)  177 (27.0) 177 (26.1)  
 Guiping/Pingnan 549 (22.5) 604 (23.7)  408 (22.9) 455 (24.3)  141 (21.5) 149 (22.0)  
Age at interview, years   0.003   0.018   0.202 
 20–39 491 (20.1) 446 (17.5)  345 (19.3) 309 (16.5)  146 (22.3) 137 (20.2)  
 40–49 878 (36.0) 869 (34.1)  622 (34.9) 620 (33.2)  256 (39.0) 249 (36.8)  
 50–59 662 (27.1) 715 (28.1)  506 (28.4) 554 (29.6)  156 (23.8) 161 (23.8)  
 60–74 410 (16.8) 516 (20.3)  312 (17.5) 386 (20.7)  98 (14.9) 130 (19.2)  
Educational level, years   0.002   0.018   0.074 
 ≤6 1,880 (39.6) 967 (35.9)  574 (32.2) 539 (28.8)  393 (59.9) 374 (55.2)  
 7–9 1,996 (40.0) 977 (40.0)  788 (44.2) 819 (43.8)  189 (28.8) 200 (29.5)  
 ≥10 1,111 (20.4) 497 (24.1)  423 (23.7) 511 (27.3)  74 (11.3) 103 (15.2)  
Current housing typeb   <0.001   0.029   <0.001 
 Building 1,766 (72.4) 1,985 (78.0)  1,301 (72.9) 1,421 (76.0)  465 (70.9) 564 (83.3)  
 Cottage/boat 675 (27.7) 561 (22.0)  484 (27.1) 448 (24.0)  191 (29.1) 113 (16.7)  
Current occupation   <0.001   <0.001   0.035 
 Unemployed/farmer 903 (37.0) 1,063 (41.8)  546 (30.6) 668 (35.7)  357 (54.4) 395 (58.4)  
 Blue-collar 990 (40.6) 882 (34.6)  815 (45.7) 742 (39.7)  175 (26.7) 140 (20.7)  
 White-collar 336 (13.8) 410 (16.1)  266 (14.9) 318 (17.0)  70 (10.7) 92 (13.6)  
 Other/unknown 212 (8.7) 191 (7.5)  158 (8.9) 141 (7.5)  54 (8.2) 50 (7.4)  
Cigarette smoking   0.133   0.020   0.412 
 Never 1,093 (44.8) 1,194 (46.9)  450 (25.2) 535 (28.6)  643 (98.0) 659 (97.3)  
 Ever 1,348 (55.2) 1,352 (53.1)  1,335 (74.8) 1,334 (71.4)  13 (2.0) 18 (2.7)  
Salt-preserved fish consumption in 2000–2002  0.020   0.268   0.016 
 ≤Yearly 1,787 (73.2) 1,873 (73.6)  1,296 (72.6) 1,359 (72.7)  491 (74.9) 514 (75.9)  
 Monthly 469 (19.2) 527 (20.7)  352 (19.7) 390 (20.9)  117 (17.8) 137 (20.2)  
 ≥Weekly 185 (7.6) 146 (5.7)  137 (7.7) 120 (6.4)  48 (7.3) 26 (3.8)  
Nasopharyngeal carcinoma among first-degree relatives  <0.001   <0.001   <0.001 
 No 2,156 (88.3) 2,457 (96.5)  1,577 (88.4) 1,809 (96.8)  579 (88.3) 648 (95.7)  
 Yes 264 (10.8) 65 (2.6)  192 (10.8) 46 (2.5)  72 (11.0) 19 (2.8)  
 Unknown 21 (0.9) 24 (0.9)  16 (0.9) 14 (0.8)  5 (0.8) 10 (1.5)  
Body mass index at age 20 years, kg/m2   0.001   0.006   0.044 
 <18.5 394 (16.1) 515 (20.2)  275 (15.4) 357 (19.1)  119 (18.1) 158 (23.3)  
 18.5–22.9 1,763 (72.2) 1,739 (68.3)  1,326 (74.3) 1,306 (69.9)  437 (66.6) 433 (64.0)  
 ≥23.0 284 (11.6) 292 (11.5)  184 (10.3) 206 (11.0)  100 (15.2) 86 (12.7)  
Vegetable consumption in 2000–2002, g/day  <0.001   <0.001   <0.001 
 <252.0 682 (27.9) 715 (28.1)  505 (28.3) 510 (27.3)  177 (27.0) 205 (30.3)  
 252.0–335.9 617 (25.3) 829 (32.6)  450 (25.2) 588 (31.5)  167 (25.5) 241 (35.6)  
 ≥336.0 1,142 (46.8) 1,002 (39.4)  830 (46.5) 771 (41.3)  312 (47.6) 231 (34.1)  
Frequency of tooth brushing, times/day   <0.001   <0.001   <0.001 
 ≤1 1,626 (66.6) 1,408 (55.3)  1,233 (69.1) 1,072 (57.4)  393 (59.9) 336 (49.6)  
 2 733 (30.0) 1,036 (40.7)  497 (27.8) 726 (38.8)  236 (36.0) 310 (45.8)  
 ≥3 82 (3.4) 102 (4.0)  55 (3.1) 71 (3.8)  27 (4.1) 31 (4.6)  
Frequency of herbal soup consumption   <0.001   <0.001   0.290 
 ≤Yearly 671 (27.5) 585 (23.0)  498 (27.9) 430 (23.0)  173 (26.4) 155 (22.9)  
 Monthly 961 (39.4) 1,122 (44.1)  689 (38.6) 819 (43.8)  272 (41.5) 303 (44.8)  
 ≥Weekly 809 (33.1) 839 (33.0)  598 (33.5) 620 (33.2)  211 (32.2) 219 (32.4)  

aP values were calculated using χ2 test.

bCurrent housing type includes building (concrete structure), or cottage (clay brick structure) or boat.

Data collection

Alcohol consumption

All participants were first asked whether they had ever consumed alcohol at least once a week for 6 months. Those who responded “yes” further reported their drinking patterns. Information was collected for seven types of alcoholic beverages: beer, yellow rice wine, fruit wine, moderate distillate spirits, strong distillate spirits, and imported liquors. The corresponding ages at starting and stopping (if applicable), cups per day (for beer, yellow rice wine, fruit wine, and other wines) or liang per day (for distillate spirits and imported liquor), and weekly and monthly frequency of intake were ascertained. One cup was defined as 200 mL and one liang (Chinese unit) as 50 mL.

We assumed an ethanol content of 4.3 g/100 mL for beer, 10.35 g/100 mL for yellow rice wine, 10.2 g/100 mL for fruit wine or wines, 30.8 g/100 mL for moderate distillate spirits, 37.2 g/100 mL for strong distillate spirits, and 42.6 g/100 mL for imported liquors, based on China Food Composition Table (27). For each specific beverage and for combined alcoholic beverages, mean daily ethanol intake (g/day) was calculated by multiplying ethanol content (g/mL), and volume of alcohol intake (mL), weighted by lifetime days of intake.

Tea consumption

Similarly, habitual tea drinkers were defined as those who ever drank tea at least once a week for 6 months. Tea drinkers further reported their drinking habits for five types of tea: black tea, green tea, jasmine tea, oolong tea, and pu'er tea. For each type of tea, the collected information included serving size [small cup (10 mL), medium cup (50 mL), or large cup (300 mL)], ages at starting and stopping (if applicable), cups per day, temperature, and concentration. Temperature of tea was defined according to the waiting time before drinking, for example, “hot” referred to drinking boiled tea immediately, whereas “cold” referred to drinking tea after more than 30 minutes or as a chilled drink. The concentration for tea intake (strong or not) was recorded according to the subjective taste of individual subject.

Lifetime volume (mL) of intake of each specific type of tea was calculated by multiplying cups per day by serving size and by days of drinking each specific type of tea, and mean daily tea intake for all kinds of tea combined was calculated by summing the cumulative lifetime volume of intake of each specific type of tea, divided by days of lifetime exposure.

Quality control

Two of the most challenging aspects of conducting a valid population-based case–control study were to ensure complete case ascertainment and unbiased control selection. Therefore, we first chose a study setting with a relatively high incidence of NPC and a relatively low rate of residential mobility. Then we developed a rapid case ascertainment system involving a network of physicians who diagnosed and/or treated NPC. Frequency-matched controls were randomly selected from a computerized, continuously updated total population registry to ensure the representativeness of controls, and we used multiple ways to contact the selected controls to improve the participation rate. In addition, we used several approaches to reduce the risk of recall bias due to the use of self-reported, retrospectively collected data. First, we developed an electronic questionnaire with built-in logic checking and several rounds of pilot testing. Second, we used an interviewer manual that described standard survey techniques to be implemented for all participants, common errors to avoid, and pictures of serving sizes for food and beverage items, including alcohol and tea drinking. Third, interviews were audiotaped for quality control, and subjects were recontacted, if necessary, to clarify responses. The detailed quality-control approaches of the NPCGEE study have been described previously (26).

Statistical analysis

Characteristics among NPC cases and population controls were compared using the χ2 test for categorical variables. To accommodate different patterns of alcohol and tea drinking between men and women, we evaluated associations stratified by sex. We used adjusted unconditional logistic regression models to estimate ORs and their corresponding 95% confidence intervals (CI). Those who never habitually drank alcohol (i.e., habitual nondrinkers) were used as the reference category for analyses of alcohol, and those who never habitually drank tea were the reference category for analyses of tea.

For both analyses, former drinkers were defined as ever drinkers with >2 years since cessation. Intake of total ethanol was categorized as 0.1–9.9, 10.0–19.9, 20.0–49.9, 50.0–89.9, and ≥90.0 g/day, according to cut-off points commonly used for ethanol intake. Intake of total tea was categorized as ≤150.0, 150.1–500.0, 500.1–900.0, 900.1–1800.0, or >1800.0 mL/day according to quintiles among control subjects. Duration of alcohol drinking or tea drinking was categorized into quartiles among male controls and all controls, respectively. For the analysis of specific tea types, black and pu'er teas (fermented teas) were combined and green and jasmine teas (nonfermented teas) were combined.

Potential confounders were selected on the basis of known or suspected risk factors for NPC that are plausibly related to patterns of tea or alcohol intake. Minimally adjusted models included age at diagnosis/interview, sex, and residential area. Fully adjusted models were additionally adjusted for education level, current housing type, current occupation, cigarette smoking, body mass index (BMI) at age 20 years, salt-preserved fish consumption during 2000–2002, NPC among first-degree relatives, vegetable consumption in 2000–2002, frequency of tooth brushing, and frequency of herbal soup consumption.

Tests of linear trend with NPC risk were performed using the median value of each category in the logistic regression models. We also used restricted cubic splines with three knots as a nonparametric approach to examine the potential nonlinear relation between ethanol consumption or tea intake and risk of NPC (28). Overall and nonlinear associations in the nonparametric regression analysis were tested using the Wald χ2 test. We also evaluated the presence of effect modification by age, sex, socioeconomic status indicators, and current smoking status by evaluating stratified results.

Analyses were performed with SAS version 9.4 (SAS Institute). The two-sided significance level for all statistical tests was 0.05.

Baseline characteristics of study subjects

Table 1 presents the distribution of characteristics among 2,441 NPC cases and 2,546 controls, overall and stratified by sex. Both among overall participants and stratified by sex, in comparison with controls, cases tended to be younger, to have a lower education level, to live in a cottage (clay brick structure) or boat, to be blue-collar workers, to have higher BMI at age 20 years, to be more likely to have a family history of NPC, to brush teeth less frequently, and to have lower consumption of herbal soup. Male cases were more likely to be ever smokers.

Risk of NPC in association with alcohol drinking

Table 2 shows the ORs and 95% CIs for associations between alcohol drinking status and risk of NPC among subjects combined and stratified by sex. Similar estimates were found in minimally and fully adjusted models. Compared with nondrinkers of alcohol, former drinkers had a nonsignificantly increased risk of NPC (OR, 1.31; 95% CI, 0.99–1.74), and no excess risk was found among current drinkers (OR, 1.08; 95% CI, 0.93–1.25). We conducted further analyses for heavy ethanol intake and found no excess NPC risk for alcohol drinkers of 90.0–119.9 g/day (OR, 1.09; 95% CI, 0.68–1.73), but a significantly 1.57-fold (95% CI, 1.01–2.46) increased risk for alcohol drinkers of ≥120.0 g/day. No significant associations were found for former or current alcohol drinkers compared with nondrinkers after stratification by sex. Among women, only 5.5% of NPC cases and 5.1% of control subjects reported ever habitual consumption of alcohol (Table 2). Alcohol intake among women was not analyzed in further detail due to the low frequency of ever consumption.

Table 2.

ORs and corresponding 95% CIs for risk of NPC in association with habitual alcohol or tea drinking status, stratified by sex.

Cases n = 2,441 n (%)Controls n = 2,546 n (%)Minimally adjusted OR (95% CI)aFully adjusted OR (95% CI)b
Alcohol drinking among women and men 
 Never 1,686 (69.1) 1,791 (70.4) 1.00 (reference) 1.00 (reference) 
 Former drinkers 130 (5.3) 107 (4.2) 1.38 (1.05–1.80) 1.31 (0.99–1.74) 
 Current drinkers 625 (25.6) 648 (25.5) 1.04 (0.91–1.19) 1.08 (0.93–1.25) 
Among women 
 Never 620 (94.5) 643 (95.0) 1.00 (reference) 1.00 (reference) 
 Former drinkers 12 (1.8) 10 (1.5) 1.35 (0.57–3.16) 1.72 (0.70–4.26) 
 Current drinkers 24 (3.7) 24 (3.6) 1.04 (0.59–1.86) 0.94 (0.51–1.73) 
Among men 
 Never 1,066 (59.7) 1,148 (61.4) 1.00 (reference) 1.00 (reference) 
 Former drinkers 118 (6.6) 97 (5.2) 1.38 (1.04–1.84) 1.29 (0.95–1.74) 
 Current drinkers 601 (33.7) 624 (33.4) 1.04 (0.91–1.20) 1.08 (0.93–1.25) 
Tea drinking among women and men 
 Never 1,573 (64.4) 1,490 (58.5) 1.00 (reference) 1.00 (reference) 
 Former drinkers 68 (2.8) 66 (2.6) 0.97 (0.69–1.38) 0.99 (0.68–1.43) 
 Current drinkers 800 (32.8) 990 (38.9) 0.74 (0.65–0.84) 0.73 (0.64–0.84) 
Among women 
 Never 562 (85.7) 553 (81.7) 1.00 (reference) 1.00 (reference) 
 Former drinkers 12 (1.8) 10 (1.5) 1.19 (0.51–2.77) 1.46 (0.59–3.63) 
 Current drinkers 82 (12.5) 114 (16.8) 0.71 (0.52–0.97) 0.70 (0.50–0.97) 
Among men 
 Never 1,011 (56.6) 937 (50.1) 1.00 (reference) 1.00 (reference) 
 Former drinkers 56 (3.1) 56 (3.0) 0.93 (0.64–1.37) 0.93 (0.62–1.40) 
 Current drinkers 718 (40.2) 876 (46.9) 0.74 (0.64–0.85) 0.74 (0.63–0.86) 
Cases n = 2,441 n (%)Controls n = 2,546 n (%)Minimally adjusted OR (95% CI)aFully adjusted OR (95% CI)b
Alcohol drinking among women and men 
 Never 1,686 (69.1) 1,791 (70.4) 1.00 (reference) 1.00 (reference) 
 Former drinkers 130 (5.3) 107 (4.2) 1.38 (1.05–1.80) 1.31 (0.99–1.74) 
 Current drinkers 625 (25.6) 648 (25.5) 1.04 (0.91–1.19) 1.08 (0.93–1.25) 
Among women 
 Never 620 (94.5) 643 (95.0) 1.00 (reference) 1.00 (reference) 
 Former drinkers 12 (1.8) 10 (1.5) 1.35 (0.57–3.16) 1.72 (0.70–4.26) 
 Current drinkers 24 (3.7) 24 (3.6) 1.04 (0.59–1.86) 0.94 (0.51–1.73) 
Among men 
 Never 1,066 (59.7) 1,148 (61.4) 1.00 (reference) 1.00 (reference) 
 Former drinkers 118 (6.6) 97 (5.2) 1.38 (1.04–1.84) 1.29 (0.95–1.74) 
 Current drinkers 601 (33.7) 624 (33.4) 1.04 (0.91–1.20) 1.08 (0.93–1.25) 
Tea drinking among women and men 
 Never 1,573 (64.4) 1,490 (58.5) 1.00 (reference) 1.00 (reference) 
 Former drinkers 68 (2.8) 66 (2.6) 0.97 (0.69–1.38) 0.99 (0.68–1.43) 
 Current drinkers 800 (32.8) 990 (38.9) 0.74 (0.65–0.84) 0.73 (0.64–0.84) 
Among women 
 Never 562 (85.7) 553 (81.7) 1.00 (reference) 1.00 (reference) 
 Former drinkers 12 (1.8) 10 (1.5) 1.19 (0.51–2.77) 1.46 (0.59–3.63) 
 Current drinkers 82 (12.5) 114 (16.8) 0.71 (0.52–0.97) 0.70 (0.50–0.97) 
Among men 
 Never 1,011 (56.6) 937 (50.1) 1.00 (reference) 1.00 (reference) 
 Former drinkers 56 (3.1) 56 (3.0) 0.93 (0.64–1.37) 0.93 (0.62–1.40) 
 Current drinkers 718 (40.2) 876 (46.9) 0.74 (0.64–0.85) 0.74 (0.63–0.86) 

aORs and 95% CIs were calculated by logistic regression models adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), and sex (male or female).

bORs and 95% CIs were calculated by logistic regression models adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), sex (male or female), education level (≤6, 7–9, or ≥10 years), current housing type [building (concrete structure) or boat/cottage (clay brick structure)], current occupation (unemployed/farmer, blue-collar, white-collar, or other/unknown), current smoking status (ever or never), tea drinking status (never, former, or current, for alcohol analysis only), alcohol drinking status (never, former, or current, for tea analysis only), body mass index at age 20 years (<18.5, 18.5–22.9, or ≥23.0 kg/m2), salt-preserved fish consumption during 2000–2002 (≤yearly, monthly, or ≥weekly), nasopharyngeal carcinoma among first-degree relatives (yes, no, or unknown), vegetable consumption in 2000–2002 (<252.0, 252.0–335.9, or ≥336.0 g/day), frequency of tooth brushing (≤1, 2, or ≥3 times/day), and frequency of herbal soup consumption (≤yearly, monthly, or ≥weekly).

Table 3 shows the ORs and 95% CIs for associations between current alcohol drinking and risk of NPC among men. No clear exposure-response trends were seen for duration of alcohol drinking or intake of specific alcoholic beverages (Ptrend > 0.05). Nonparametric regression curves indicated a nonsignificant association between current daily ethanol intake and NPC risk in men (P for nonlinear association = 0.200; Fig. 1A).

Table 3.

ORs and corresponding 95% CIs for risk of NPC in association with current habitual alcohol consumption among men.

Cases n = 1,667 n (%)Controls n = 1,772 n (%)Minimally adjusted OR (95% CI)aFully adjusted OR (95% CI)b
Total ethanol consumption, g/day 
 Neverc 1,066 (64.0) 1,148 (64.8) 1.00 (reference) 1.00 (reference) 
 0.1–9.9 167 (10.0) 190 (10.7) 0.93 (0.74–1.17) 1.04 (0.82–1.31) 
 10.0–19.9 108 (6.5) 114 (6.4) 1.02 (0.78–1.35) 1.13 (0.85–1.51) 
 20.0–49.9 137 (8.2) 153 (8.6) 0.98 (0.76–1.25) 0.95 (0.73–1.24) 
 50.0–89.9 96 (5.8) 91 (5.1) 1.16 (0.86–1.57) 1.14 (0.83–1.57) 
 ≥90.0 93 (5.6) 76 (4.3) 1.35 (0.98–1.85) 1.32 (0.95–1.84) 
Ptrend   0.043 0.102 
Duration of alcohol drinking, years 
 <14.0 133 (8.0) 140 (7.9) 0.96 (0.75–1.24) 1.03 (0.79–1.35) 
 14.0–22.9 175 (10.5) 158 (8.9) 1.14 (0.91–1.45) 1.23 (0.96–1.58) 
 23.0–31.9 157 (9.4) 152 (8.6) 1.13 (0.89–1.44) 1.17 (0.91–1.51) 
 ≥32.0 136 (8.2) 174 (9.8) 0.94 (0.73–1.2) 0.91 (0.69–1.18) 
Ptrend   0.675 0.628 
Type of alcoholic beverage 
 Beer 218 (17.0) 195 (14.5) 1.14 (0.92–1.42) 1.23 (0.98–1.55) 
 Winesd 41 (3.7) 29 (2.5) 1.50 (0.92–2.43) 1.88 (1.13–3.13) 
 Moderate distillate spirits 464 (30.3) 489 (29.9) 1.05 (0.90–1.22) 1.08 (0.91–1.27) 
 Strong distillate spirits/liquors 140 (11.6) 110 (8.7) 1.35 (1.03–1.76) 1.56 (1.18–2.08) 
Beer intake, g/day 
 0.1–9.9 141 (11.0) 134 (10.0) 1.09 (0.84–1.40) 1.23 (0.94–1.61) 
 ≥10.0 77 (6.0) 61 (4.5) 1.27 (0.89–1.80) 1.25 (0.86–1.81) 
Ptrend   0.167 0.189 
Moderate distillate spirits intake, g/day 
 0.1–9.9 138 (9.0) 146 (8.9) 1.03 (0.80–1.32) 1.16 (0.90–1.51) 
 ≥10.0 326 (21.3) 343 (21.0) 1.06 (0.89–1.26) 1.04 (0.86–1.26) 
Ptrend   0.539 0.758 
Strong distillate spirits and liquors intake, g/day 
 0.1–9.9 58 (4.8) 43 (3.4) 1.41 (0.94–2.11) 1.88 (1.23–2.88) 
 ≥10.0 82 (6.8) 67 (5.3) 1.31 (0.94–1.84) 1.37 (0.96–1.97) 
Ptrend   0.106 0.077 
Cases n = 1,667 n (%)Controls n = 1,772 n (%)Minimally adjusted OR (95% CI)aFully adjusted OR (95% CI)b
Total ethanol consumption, g/day 
 Neverc 1,066 (64.0) 1,148 (64.8) 1.00 (reference) 1.00 (reference) 
 0.1–9.9 167 (10.0) 190 (10.7) 0.93 (0.74–1.17) 1.04 (0.82–1.31) 
 10.0–19.9 108 (6.5) 114 (6.4) 1.02 (0.78–1.35) 1.13 (0.85–1.51) 
 20.0–49.9 137 (8.2) 153 (8.6) 0.98 (0.76–1.25) 0.95 (0.73–1.24) 
 50.0–89.9 96 (5.8) 91 (5.1) 1.16 (0.86–1.57) 1.14 (0.83–1.57) 
 ≥90.0 93 (5.6) 76 (4.3) 1.35 (0.98–1.85) 1.32 (0.95–1.84) 
Ptrend   0.043 0.102 
Duration of alcohol drinking, years 
 <14.0 133 (8.0) 140 (7.9) 0.96 (0.75–1.24) 1.03 (0.79–1.35) 
 14.0–22.9 175 (10.5) 158 (8.9) 1.14 (0.91–1.45) 1.23 (0.96–1.58) 
 23.0–31.9 157 (9.4) 152 (8.6) 1.13 (0.89–1.44) 1.17 (0.91–1.51) 
 ≥32.0 136 (8.2) 174 (9.8) 0.94 (0.73–1.2) 0.91 (0.69–1.18) 
Ptrend   0.675 0.628 
Type of alcoholic beverage 
 Beer 218 (17.0) 195 (14.5) 1.14 (0.92–1.42) 1.23 (0.98–1.55) 
 Winesd 41 (3.7) 29 (2.5) 1.50 (0.92–2.43) 1.88 (1.13–3.13) 
 Moderate distillate spirits 464 (30.3) 489 (29.9) 1.05 (0.90–1.22) 1.08 (0.91–1.27) 
 Strong distillate spirits/liquors 140 (11.6) 110 (8.7) 1.35 (1.03–1.76) 1.56 (1.18–2.08) 
Beer intake, g/day 
 0.1–9.9 141 (11.0) 134 (10.0) 1.09 (0.84–1.40) 1.23 (0.94–1.61) 
 ≥10.0 77 (6.0) 61 (4.5) 1.27 (0.89–1.80) 1.25 (0.86–1.81) 
Ptrend   0.167 0.189 
Moderate distillate spirits intake, g/day 
 0.1–9.9 138 (9.0) 146 (8.9) 1.03 (0.80–1.32) 1.16 (0.90–1.51) 
 ≥10.0 326 (21.3) 343 (21.0) 1.06 (0.89–1.26) 1.04 (0.86–1.26) 
Ptrend   0.539 0.758 
Strong distillate spirits and liquors intake, g/day 
 0.1–9.9 58 (4.8) 43 (3.4) 1.41 (0.94–2.11) 1.88 (1.23–2.88) 
 ≥10.0 82 (6.8) 67 (5.3) 1.31 (0.94–1.84) 1.37 (0.96–1.97) 
Ptrend   0.106 0.077 

aORs and 95% CIs were calculated by logistic regression models adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), and sex (male or female). Tests of linear trend with risk of nasopharyngeal carcinoma were performed using the median value of each category in the logistic regression models.

bORs and 95% CIs were calculated by logistic regression models adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), sex (male or female), education level (≤6, 7–9, or ≥10 years), current housing type [building (concrete structure) or boat/cottage (clay brick structure)], current occupation (unemployed/farmer, blue-collar, white-collar, or other/unknown), current smoking status (ever or never), tea drinking status (never, former, or current), body mass index at age 20 years (<18.5, 18.5–22.9, or ≥23.0 kg/m2), salt-preserved fish consumption during 2000–2002 (≤yearly, monthly, or ≥weekly), nasopharyngeal carcinoma among first-degree relatives (yes, no, or unknown), vegetable consumption in 2000–2002 (<252.0, 252.0–335.9, or ≥336.0 g/day), frequency of tooth brushing (≤1, 2, or ≥3 times/day), and frequency of herbal soup consumption (≤yearly, monthly, or ≥weekly). Tests of linear trend with risk of nasopharyngeal carcinoma were performed using the median value of each category in the logistic regression models.

cNever habitual drinkers of alcohol were the reference group for all comparisons.

dWines included any type of the yellow rice wine, wines, or fruit wine.

Figure 1.

Natural logarithm of the OR and 95% CI for NPC in relation to current daily mean ethanol intake and tea consumption among habitual drinkers. The restricted cubic spline analysis with three knots for daily current ethanol intake was conducted among men (A; 1,667 cases and 1,772 controls of male nondrinkers and current alcohol drinkers), and for daily current tea intake among men and women (B; 2,373 cases and 2,480 controls of male and female nondrinkers and current tea drinkers). The restricted cubic spline analysis was adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), sex (male or female), education level (≤6, 7–9, or ≥10 years), current housing type [building (concrete structure) or boat/cottage (clay brick structure)], current occupation (unemployed/farmer, blue-collar, white-collar, or other/unknown), current smoking status (ever or never), tea drinking status (never, former, or current, for alcohol analysis only), alcohol drinking status (never, former, or current, for tea analysis only), BMI at age 20 years (<18.5, 18.5–22.9, or ≥23.0 kg/m2), salt-preserved fish consumption during 2000–2002 (≤yearly, monthly, or ≥weekly), NPC among first-degree relatives (yes, no, or unknown), vegetable consumption in 2000–2002 (<252.0, 252.0–335.9, or ≥336.0 g/day), frequency of tooth brushing (≤1, 2, or ≥3 times/day), and frequency of herbal soup consumption (≤yearly, monthly, or ≥weekly).

Figure 1.

Natural logarithm of the OR and 95% CI for NPC in relation to current daily mean ethanol intake and tea consumption among habitual drinkers. The restricted cubic spline analysis with three knots for daily current ethanol intake was conducted among men (A; 1,667 cases and 1,772 controls of male nondrinkers and current alcohol drinkers), and for daily current tea intake among men and women (B; 2,373 cases and 2,480 controls of male and female nondrinkers and current tea drinkers). The restricted cubic spline analysis was adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), sex (male or female), education level (≤6, 7–9, or ≥10 years), current housing type [building (concrete structure) or boat/cottage (clay brick structure)], current occupation (unemployed/farmer, blue-collar, white-collar, or other/unknown), current smoking status (ever or never), tea drinking status (never, former, or current, for alcohol analysis only), alcohol drinking status (never, former, or current, for tea analysis only), BMI at age 20 years (<18.5, 18.5–22.9, or ≥23.0 kg/m2), salt-preserved fish consumption during 2000–2002 (≤yearly, monthly, or ≥weekly), NPC among first-degree relatives (yes, no, or unknown), vegetable consumption in 2000–2002 (<252.0, 252.0–335.9, or ≥336.0 g/day), frequency of tooth brushing (≤1, 2, or ≥3 times/day), and frequency of herbal soup consumption (≤yearly, monthly, or ≥weekly).

Close modal

In contrast, among former drinkers, there was no trend either with ethanol intake or with duration of drinking in fully adjusted models (Supplementary Table S1).

Risk of NPC in association with tea intake

Table 2 also shows the ORs and 95% CIs for associations between tea drinking status and risk of NPC among overall subjects and stratified by sex. Compared with the nondrinkers, current tea drinking was significantly associated with a 27% decrease in risk of NPC (95% CI, 16%–36%) in the fully adjusted model. Similar results were found after stratification by sex, although men were more likely than women to have ever had a tea drinking habit (Table 2).

Table 4 lists the ORs and 95% CIs for NPC risk in relation to current tea drinking among overall participants. Mean daily current intake of tea (mL/day) showed no clear pattern of association with NPC risk; the ORs were significantly lower than 1 for first three categories, compared with never drinkers, but not significant for the next two categories. Nonparametric regression curves by restricted cubic spline analysis with three knots suggested a significant nonlinear dose–response relation between volume of current daily tea intake and NPC risk (P for overall association = 0.038; P for nonlinear association = 0.033; Fig. 1B).

Table 4.

ORs and corresponding 95% CIs for risk of NPC in association with current habitual tea consumption among men and women.

Cases n = 2,373 n (%)Controls n = 2,480 n (%)Minimally adjusted OR (95% CI)aFully adjusted OR (95% CI)b
Tea consumption, mL/day 
 Neverc 1,573 (66.3) 1,490 (60.1) 1.00 (reference) 1.00 (reference) 
 ≤150.0 117 (4.9) 205 (8.3) 0.52 (0.40–0.66) 0.53 (0.41–0.68) 
 150.1–500.0 163 (6.9) 220 (8.9) 0.68 (0.55–0.85) 0.68 (0.54–0.85) 
 500.1–900.0 132 (5.6) 191 (7.7) 0.65 (0.51–0.82) 0.65 (0.51–0.84) 
 900.1–1,800.0 198 (8.3) 180 (7.3) 1.03 (0.82–1.28) 1.03 (0.81–1.30) 
 >1,800.0 190 (8.0) 194 (7.8) 0.93 (0.74–1.16) 0.89 (0.70–1.13) 
Ptrend   0.918 0.999 
Years of tea drinking 
 <12.0 159 (6.7) 242 (9.8) 0.58 (0.47–0.72) 0.58 (0.46–0.72) 
 12.0–22.9 216 (9.1) 226 (9.1) 0.84 (0.68–1.03) 0.88 (0.71–1.10) 
 23.0–32.9 208 (8.8) 269 (10.9) 0.70 (0.57–0.86) 0.68 (0.55–0.85) 
 ≥33.0 217 (9.1) 253 (10.2) 0.86 (0.70–1.07) 0.84 (0.67–1.05) 
Ptrend   0.010 0.011 
Tea temperature 
 Hot 51 (2.2) 77 (3.1) 0.60 (0.42–0.86) 0.61 (0.41–0.89) 
 Warm 168 (7.1) 251 (10.1) 0.61 (0.49–0.76) 0.60 (0.48–0.75) 
 Mild 527 (22.2) 622 (25.1) 0.77 (0.67–0.90) 0.78 (0.67–0.91) 
 Cold 54 (2.3) 40 (1.6) 1.27 (0.83–1.92) 1.11 (0.72–1.72) 
Strong tea drinking 
 Yes 336 (14.2) 329 (13.3) 0.94 (0.79–1.12) 0.89 (0.74–1.08) 
 No 464 (19.6) 661 (26.7) 0.64 (0.56–0.75) 0.66 (0.57–0.77) 
Tea type 
 Black/Pu'er tea 235 (13.0) 288 (16.2) 0.74 (0.61–0.90) 0.79 (0.64–0.97) 
 Green/Jasmine tea 681 (30.2) 833 (35.9) 0.75 (0.66–0.86) 0.75 (0.65–0.87) 
 Oolong tea 63 (3.9) 44 (2.9) 1.29 (0.87–1.91) 1.40 (0.92–2.12) 
Cases n = 2,373 n (%)Controls n = 2,480 n (%)Minimally adjusted OR (95% CI)aFully adjusted OR (95% CI)b
Tea consumption, mL/day 
 Neverc 1,573 (66.3) 1,490 (60.1) 1.00 (reference) 1.00 (reference) 
 ≤150.0 117 (4.9) 205 (8.3) 0.52 (0.40–0.66) 0.53 (0.41–0.68) 
 150.1–500.0 163 (6.9) 220 (8.9) 0.68 (0.55–0.85) 0.68 (0.54–0.85) 
 500.1–900.0 132 (5.6) 191 (7.7) 0.65 (0.51–0.82) 0.65 (0.51–0.84) 
 900.1–1,800.0 198 (8.3) 180 (7.3) 1.03 (0.82–1.28) 1.03 (0.81–1.30) 
 >1,800.0 190 (8.0) 194 (7.8) 0.93 (0.74–1.16) 0.89 (0.70–1.13) 
Ptrend   0.918 0.999 
Years of tea drinking 
 <12.0 159 (6.7) 242 (9.8) 0.58 (0.47–0.72) 0.58 (0.46–0.72) 
 12.0–22.9 216 (9.1) 226 (9.1) 0.84 (0.68–1.03) 0.88 (0.71–1.10) 
 23.0–32.9 208 (8.8) 269 (10.9) 0.70 (0.57–0.86) 0.68 (0.55–0.85) 
 ≥33.0 217 (9.1) 253 (10.2) 0.86 (0.70–1.07) 0.84 (0.67–1.05) 
Ptrend   0.010 0.011 
Tea temperature 
 Hot 51 (2.2) 77 (3.1) 0.60 (0.42–0.86) 0.61 (0.41–0.89) 
 Warm 168 (7.1) 251 (10.1) 0.61 (0.49–0.76) 0.60 (0.48–0.75) 
 Mild 527 (22.2) 622 (25.1) 0.77 (0.67–0.90) 0.78 (0.67–0.91) 
 Cold 54 (2.3) 40 (1.6) 1.27 (0.83–1.92) 1.11 (0.72–1.72) 
Strong tea drinking 
 Yes 336 (14.2) 329 (13.3) 0.94 (0.79–1.12) 0.89 (0.74–1.08) 
 No 464 (19.6) 661 (26.7) 0.64 (0.56–0.75) 0.66 (0.57–0.77) 
Tea type 
 Black/Pu'er tea 235 (13.0) 288 (16.2) 0.74 (0.61–0.90) 0.79 (0.64–0.97) 
 Green/Jasmine tea 681 (30.2) 833 (35.9) 0.75 (0.66–0.86) 0.75 (0.65–0.87) 
 Oolong tea 63 (3.9) 44 (2.9) 1.29 (0.87–1.91) 1.40 (0.92–2.12) 

aORs and 95% CIs were calculated by logistic regression models adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), and sex (male or female). Tests of linear trend with risk of nasopharyngeal carcinoma were performed using the median value of each category in the logistic regression models.

bORs and 95% CIs were calculated by logistic regression models adjusted for age (20–39, 40–49, 50–59, or 60–74 years), area of residence (Zhaoqing, Wuzhou, or Guiping/Pingnan), sex (male or female), education level (≤6, 7–9, or ≥10 years), current housing type [building (concrete structure) or boat/cottage (clay brick structure)], current occupation (unemployed/farmer, blue-collar, white-collar, or other/unknown), current smoking status (ever or never), alcohol drinking status (never, former, or current), body mass index at age 20 years (<18.5, 18.5–22.9, or ≥23.0 kg/m2), salt-preserved fish consumption during 2000–2002 (≤yearly, monthly, or ≥weekly), nasopharyngeal carcinoma among first-degree relatives (yes, no, or unknown), vegetable consumption in 2000–2002 (<252.0, 252.0–335.9, or ≥336.0 g/day), frequency of tooth brushing (≤1, 2, or ≥3 times/day) and frequency of herbal soup consumption (≤yearly, monthly, or ≥weekly). Tests of linear trend with risk of nasopharyngeal carcinoma were performed using the median value of each category in the logistic regression models.

cNever habitual drinkers of tea were the reference group for all comparisons.

Analyses by tea temperature, concentration of tea, and types of tea revealed that those who consumed warm or hot tea, those who did not drink strong tea, and those who drank black and green tea tended to have a decreased NPC risk (Table 4). We found no significant interaction between tea drinking and age, sex, educational level, current housing type, and current occupation (Supplementary Table S2).

To evaluate whether tea cup size played an effect-modification role, we conducted further analyses stratified by tea cup size. We found a linear decreasing trend between daily tea intake and risk of NPC among tea drinkers using small tea cups (Supplementary Fig. 2A), whereas tea drinkers using medium or large tea cups had attenuated or no decreased NPC risk (Supplementary Fig. S2B and S2C).

To our knowledge, this is among the first endeavors to assess the association of detailed alcohol and tea consumption with the risk of NPC in a large population-based study. We found no increased risk of NPC among current alcohol drinkers, and analyses by daily consumption, duration, and type of alcoholic beverage provided no compelling evidence for a causal association. In contrast, tea intake was associated with a decreased risk, although the lack of a monotonic dose–response trend in relation to amount or duration of tea drinking complicates causal inference and suggests that residual confounding or bias may play a role.

Consistent with our findings, most studies reported no association between overall alcohol drinking and NPC risk (10–16, 18, 21, 29). A J-shaped dose–response relationship was found in a meta-analysis based on six case–control studies (19), and results from a subsequent case–control study indicated that light or moderate alcohol intake was associated with a reduced risk, while heavy consumption increased risk (17). Although we found a moderately increased risk for heavy ethanol intake of ≥120.0 g/day, caution is needed when interpreting this result, as the number of heavy alcohol drinkers was small. Another recent meta-analysis of two cohort studies and population-based case–control studies showed no significant association between alcohol consumption and NPC (22). Studies conducted in other high-risk areas for NPC, such as Hong Kong, Taiwan, Malaysia, Thailand, and Singapore, as well as some low-risk areas, showed similar null results (2–4, 6, 24, 30, 31). Most of these studies, however, had a small sample size and could not calculate specific ethanol intake.

Several previous studies examined four types of alcoholic beverages in association with NPC risk, but reported inconsistent associations for beer or spirits (13, 17, 24, 32), and consistently no increased risk for wine or Chinese rice wine (17, 24, 33). Most of these studies were unable to evaluate dose–response trends for each type of alcoholic beverage. We found an increased risk of NPC among current alcohol drinkers with wines or low intake of strong distillate spirits/liquors (<10 g/day), but these associations were attenuated and not significant after additional adjustment for dietary pattern (Supplementary Table S3), suggesting that the positive association could be partly explained by dietary pattern or other unknown residual confounding.

We found a significant inverse association between current intake of tea, especially black/pu'er tea, and risk of NPC; this decreased risk was confined to drinkers of hot/warm/mild tea or weak tea. Consistent with our results, four case–control studies also reported a significantly decreased risk of NPC associated with overall tea drinking (10, 14, 17, 23, 24); this association differed by specific tea type. One study reported a reduced risk among green tea drinkers (29); another found an inverse association among drinkers of any of six tea types (17). The inconsistencies may arise due to different patterns of tea drinking or tea manufacturing across areas, or they may be attributable to chance or residual confounding. To our knowledge, no prospective cohort studies have explored the relationship between tea drinking and NPC risk.

Our results showed that there seemed a nonlinear relationship between daily tea intake and NPC risk after adjusting for most confounders. Misclassification or differential recall bias may partly contribute to this nonlinear relationship, although exposure misclassification often leads to linearization of apparent exposure-response trends, rather than nonlinearization (34). However, we cannot eliminate the influence of residual confounding, even after adjusting for numerous confounders. Another possible reason for our findings may be that different tea drinking patterns in the study area are correlated with the use of characteristic tea cups of distinct sizes. Those with low daily tea intake were more likely to use small tea cups and drink Kungfu tea, which was typically consumed at a higher concentration with a short brewing time. Those with high daily tea intake more often used medium or large tea cups, in which tea was brewed for a relatively long time in a large teapot and, potentially allowing for oxidation of protective substances such as tea polyphenols. Although we did not find a significant association between highly concentrated tea drinking and NPC risk, this characteristic was subjectively reported in our study. We found a significant effect modification by tea cup size, which support this interpretation; however, these results were based on small numbers and require independent confirmation.

A protective effect of tea drinking is, nevertheless, biologically plausible. In animal studies, the epigallocatechin gallate, which can be extracted from tea, may reactivate methylation-silenced genes in cancer cell lines and inhibit cancer growth (35, 36). In addition, the tea polyphenol (-)-epigallocatechin-3-gallate has potent antimicrobial activity against viruses and block EBV infection–induced cytokine expression and EBV-induced B-lymphocyte transformation (37). These catechins and others are found in multiple tea types at various concentrations; for example, epigallocatechin-3-gallate levels in oolong tea were found to be lower than in other tea types, although findings may be method dependent (38).

Strengths of our study include its large size and strictly population-based design in an NPC-endemic area, with high and equal participation rate among eligible cases and controls. Potential but nonquantifiable limitations include recall bias and reverse causality. We undertook, however, extensive efforts to reduce or eliminate such bias through numerous approaches described in detail in earlier publications (7, 8, 25, 26). Furthermore, alcohol and tea are not well-established risk or protective factor for NPC. Thus, the influence of recall bias might not be a big concern. We did not collect information about the specific reasons for quitting alcohol or tea use among former drinkers, which hampers assessment of reverse causality. However, we defined former ethanol drinkers as those who had stopped for >2 years, thus symptoms of NPC might not be the major reasons for them to stop drinking ethanol. We were unable to collect detailed information on the ethanol content of individual alcoholic beverages; instead, we used the mean or most commonly used concentrations (39, 40) to calculate cumulative ethanol consumption. Finally, we were unable to distinguish between patterns of regular drinking and episodic drinking of either alcohol or tea.

In conclusion, a statistical association between alcohol drinking and the risk of NPC is more likely nonexisting. Tea drinking may be associated with moderately reduced NPC risk. However, these conclusions warrant further investigation and should be interpreted with caution for some residual confounding.

E.T. Chang reports grants from NCI during the conduct of the study and other from Exponent outside the submitted work. I. Ernberg reports grants from Swedish Cancer Society during the conduct of the study. No disclosures were reported by the other authors.

The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies.

R. Feng: Conceptualization, data curation, formal analysis, visualization, writing–original draft. E.T. Chang: data curation, writing–review and editing. Q. Liu: Conceptualization, data curation. Y. Cai: Resources, data curation. Z. Zhang: Conceptualization, data curation, funding acquisition. G. Chen: Conceptualization, data curation. Q.-H. Huang: Conceptualization, resources, data curation. S.-H. Xie: Resources, data curation, project administration. S.-M. Cao: Conceptualization, resources, data curation. Y. Zhang: Conceptualization, resources, data curation. J.-P. Yun: Resources, data curation. W.-H. Jia: Conceptualization, resources, data curation, methodology. Y. Zheng: Conceptualization, resources, data curation. J. Liao: Conceptualization, resources, data curation, investigation. Y. Chen: Resources, data curation, investigation. T. Huang: Resources, data curation, formal analysis, investigation. L. Lin: Conceptualization, resources, data curation. I. Ernberg: Conceptualization, resources, data curation, methodology. G. Huang: Conceptualization, resources, data curation, supervision, methodology. Y.-X. Zeng: Conceptualization, resources, data curation, supervision, methodology. H.-O. Adami: Conceptualization, resources, data curation, supervision, funding acquisition, methodology, writing–review and editing. W. Ye: Conceptualization, resources, data curation, supervision, methodology, project administration, writing–review and editing.

The NCI at the U.S. NIH [grant number R01 CA115873; principal investigator (PI): H.-O. Adami; co-PIs: Y. Zeng, Y.-X. Zeng]; the Swedish Research Council (2015-02625, 2015-06268, 2017-05814, to W. Ye); and the Karolinska Institutet Distinguished Professor Award (Dnr: 2368/10-221, to H.-O. Adami). The work in the Guiping/Pingnan area was supported by grants from the New Century Excellent Talents in University (no. NCET-12-0654, to Z. Zhang), National Basic Research Program of China (no. 2011CB504300, to G. Huang), and Guangxi Natural Science Foundation (2013GXNSFGA 019002, to Z. Zhang). R. Feng was partly supported by the Postdoctoral Research Foundation of China (2018M640874).

The authors sincerely thank their deceased colleague, Yi Zeng (Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention), who played a pivotal role during the planning and conduct of this study. They also thank the members of the External Advisory Board of the NPC Genes, Environment, and EBV Study, including Curtis Harris (U.S. NCI), Xihong Lin (Harvard T.H. Chan School of Public Health), Allan Hildesheim (U.S. NCI), Wei-Cheng You (Peking University), and You-Lin Qiao (Chinese Academy of Medical Sciences) for their guidance of the overall case–control study.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Bray
F
,
Ferlay
J
,
Soerjomataram
I
,
Siegel
RL
,
Torre
LA
,
Jemal
A
. 
Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries
.
CA Cancer J Clin
2018
;
68
:
394
424
.
2.
Armstrong
RW
,
Imrey
PB
,
Lye
MS
,
Armstrong
MJ
,
Yu
MC
,
Sani
S
. 
Nasopharyngeal carcinoma in Malaysian Chinese: salted fish and other dietary exposures
.
Int J Cancer
1998
;
77
:
228
35
.
3.
Ekburanawat
W
,
Ekpanyaskul
C
,
Brennan
P
,
Kanka
C
,
Tepsuwan
K
,
Temiyastith
S
, et al
Evaluation of non-viral risk factors for nasopharyngeal carcinoma in thailand: results from a case-control study
.
Asian Pac J Cancer Prev
2010
;
11
:
929
32
.
4.
Mai
ZM
,
Lin
JH
,
Ngan
RK
,
Kwong
DL
,
Ng
WT
,
Ng
AW
, et al
Milk consumption across life periods in relation to lower risk of nasopharyngeal carcinoma: a multicentre case-control study
.
Front Oncol
2019
;
9
:
253
.
5.
Mugheri
MH
,
Amur
SA
,
Channa
NA
,
Soomro
NA
,
Khuhro
Q
,
Paras
M
. 
Incidence of cancer at Liaquat University of Medical and Health Sciences Hospital, Jamshoro from 2010–2016: a retrospective study
.
Rawal Med J
2019
;
44
:
482
5
.
6.
Yong
SK
,
Ha
TC
,
Yeo
MC
,
Gaborieau
V
,
McKay
JD
,
Wee
J
. 
Associations of lifestyle and diet with the risk of nasopharyngeal carcinoma in Singapore: a case-control study
.
Chin J Cancer
2017
;
36
:
3
.
7.
Liu
Z
,
Chang
ET
,
Liu
Q
,
Cai
Y
,
Zhang
Z
,
Chen
G
, et al
Oral hygiene and risk of nasopharyngeal carcinoma-a population-based case-control study in China
.
Cancer Epidemiol Biomarkers Prev
2016
;
25
:
1201
7
.
8.
Chang
ET
,
Liu
Z
,
Hildesheim
A
,
Liu
Q
,
Cai
Y
,
Zhang
Z
, et al
Active and passive smoking and risk of nasopharyngeal carcinoma: a population-based case-control study in southern China
.
Am J Epidemiol
2017
;
185
:
1272
80
.
9.
Barrett
D
,
Ploner
A
,
Chang
ET
,
Liu
Z
,
Zhang
CX
,
Liu
Q
, et al
Past and recent salted fish and preserved food intakes are weakly associated with nasopharyngeal carcinoma risk in adults in southern China
.
J Nutr
2019
;
149
:
1596
605
.
10.
Turkoz
FP
,
Celenkoglu
G
,
Dogu
GG
,
Kalender
ME
,
Coskun
U
,
Alkis
N
, et al
Risk factors of nasopharyngeal carcinoma in Turkey-an epidemiological survey of the anatolian society of medical oncology
.
Asian Pac J Cancer Prev
2011
;
12
:
3017
21
.
11.
Sriamporn
S
,
Vatanasapt
V
,
Pisani
P
,
Yongchaiyudha
S
,
Rungpitarangsri
V
. 
Environmental risk factors for nasopharyngeal carcinoma: a case-control study in northeastern Thailand
.
Cancer Epidemiol Biomarkers Prev
1992
;
1
:
345
8
.
12.
Chen
CJ
,
Liang
KY
,
Chang
YS
,
Wang
YF
,
Hsieh
T
,
Hsu
MM
, et al
Multiple risk factors of nasopharyngeal carcinoma: Epstein-Barr virus, malarial infection, cigarette smoking and familial tendency
.
Anticancer Res
1990
;
10
:
547
53
.
13.
Armstrong
RW
,
Armstrong
MJ
,
Yu
MC
,
Henderson
BE
. 
Salted fish and inhalants as risk factors for nasopharyngeal carcinoma in Malaysian Chinese
.
Cancer Res
1983
;
43
:
2967
70
.
14.
Xu
FH
,
Xiong
D
,
Xu
YF
,
Cao
SM
,
Xue
WQ
,
Qin
HD
, et al
An epidemiological and molecular study of the relationship between smoking, risk of nasopharyngeal carcinoma, and Epstein-Barr virus activation
.
J Natl Cancer Inst
2012
;
104
:
1396
410
.
15.
Fachiroh
J
,
Sangrajrang
S
,
Johansson
M
,
Renard
H
,
Gaborieau
V
,
Chabrier
A
, et al
Tobacco consumption and genetic susceptibility to nasopharyngeal carcinoma (NPC) in Thailand
.
Cancer Causes Control
2012
;
23
:
1995
2002
.
16.
Yuan
JM
,
Wang
XL
,
Xiang
YB
,
Gao
YT
,
Ross
RK
,
Yu
MC
. 
Non-dietary risk factors for nasopharyngeal carcinoma in Shanghai, China
.
Int J Cancer
2000
;
85
:
364
9
.
17.
Ruan
HL
,
Xu
FH
,
Liu
WS
,
Feng
QS
,
Chen
LZ
,
Zeng
YX
, et al
Alcohol and tea consumption in relation to the risk of nasopharyngeal carcinoma in Guangdong, China
.
Front Med China
2010
;
4
:
448
56
.
18.
Ji
X
,
Zhang
W
,
Xie
C
,
Wang
B
,
Zhang
G
,
Zhou
F
. 
Nasopharyngeal carcinoma risk by histologic type in central China: impact of smoking, alcohol and family history
.
Int J Cancer
2011
;
129
:
724
32
.
19.
Chen
L
,
Gallicchio
L
,
Boyd-Lindsley
K
,
Tao
XG
,
Robinson
KA
,
Lam
TK
, et al
Alcohol consumption and the risk of nasopharyngeal carcinoma: a systematic review
.
Nutr Cancer
2009
;
61
:
1
15
.
20.
Nam
JM
,
McLaughlin
JK
,
Blot
WJ
. 
Cigarette smoking, alcohol, and nasopharyngeal carcinoma: a case-control study among U.S. whites
.
J Natl Cancer Inst
1992
;
84
:
619
22
.
21.
Friborg
JT
,
Yuan
JM
,
Wang
R
,
Koh
WP
,
Lee
HP
,
Yu
MC
. 
A prospective study of tobacco and alcohol use as risk factors for pharyngeal carcinomas in Singapore Chinese
.
Cancer
2007
;
109
:
1183
91
.
22.
Du
T
,
Chen
K
,
Zheng
S
,
Bao
M
,
Huang
Y
,
Wu
K
. 
Association between alcohol consumption and risk of nasopharyngeal carcinoma: a comprehensive meta-analysis of epidemiological studies
.
Alcohol Clin Exp Res
2019
;
43
:
2262
73
.
23.
Hsu
WL
,
Pan
WH
,
Chien
YC
,
Yu
KJ
,
Cheng
YJ
,
Chen
JY
, et al
Lowered risk of nasopharyngeal carcinoma and intake of plant vitamin, fresh fish, green tea and coffee: a case-control study in Taiwan
.
PLoS One
2012
;
7
:
e41779
.
24.
Henderson
BE
,
Louie
E
,
SooHoo Jing
J
,
Buell
P
,
Gardner
MB
. 
Risk factors associated with nasopharyngeal carcinoma
.
N Engl J Med
1976
;
295
:
1101
6
.
25.
Liu
Z
,
Chang
ET
,
Liu
Q
,
Cai
Y
,
Zhang
Z
,
Chen
G
, et al
Quantification of familial risk of nasopharyngeal carcinoma in a high-incidence area
.
Cancer
2017
;
123
:
2716
25
.
26.
Ye
W
,
Chang
ET
,
Liu
Z
,
Liu
Q
,
Cai
Y
,
Zhe
Z
, et al
Development of a population-based cancer case-control study in southern china
.
Oncotarget
2017
;
8
:
87073
85
.
27.
Yang
YX
,
He
M
,
Pan
XC
.
Institute of Nutrition and Food Safety China
.
China food composition table 2004 book 1
, 2nd ed.
Beijing, China
:
Peking University Medical Press
; 
2004
.
28.
Desquilbet
L
,
Mariotti
F
. 
Dose-response analyses using restricted cubic spline functions in public health research
.
Stat Med
2010
;
29
:
1037
57
.
29.
Cheng
YJ
,
Hildesheim
A
,
Hsu
MM
,
Chen
IH
,
Brinton
LA
,
Levine
PH
, et al
Cigarette smoking, alcohol consumption and risk of nasopharyngeal carcinoma in Taiwan
.
Cancer Causes Control
1999
;
10
:
201
7
.
30.
Huang
HC
,
Tantoh
DM
,
Hsu
SY
,
Nfor
ON
,
Frank
CL
,
Lung
CC
, et al
Association between coarse particulate matter (PM10–2.5) and nasopharyngeal carcinoma among Taiwanese men
.
J Investig Med
2020
;
68
:
419
24
.
31.
Polesel
J
,
Serraino
D
,
Negri
E
,
Barzan
L
,
Vaccher
E
,
Montella
M
, et al
Consumption of fruit, vegetables, and other food groups and the risk of nasopharyngeal carcinoma
.
Cancer Causes Control
2013
;
24
:
1157
65
.
32.
Mabuchi
K
,
Bross
DS
,
Kessler
II
. 
Cigarette smoking and nasopharyngeal carcinoma
.
Cancer
1985
;
55
:
2874
6
.
33.
Geser
A
,
Charnay
N
,
Day
NE
,
Dethe
G
,
Ho
HC
. 
Environmental factors in the etiology of nasopharyngeal carcinoma: report on a case-control study in Hong Kong
.
IARC Sci Publ
1978
;
20
:
213
29
.
34.
Rhomberg
LR
,
Chandalia
JK
,
Long
CM
,
Goodman
JE
. 
Measurement error in environmental epidemiology and the shape of exposure-response curves
.
Crit Rev Toxicol
2011
;
41
:
651
71
.
35.
Fang
CY
,
Wu
CC
,
Hsu
HY
,
Chuang
HY
,
Huang
SY
,
Tsai
CH
, et al
EGCG inhibits proliferation, invasiveness and tumor growth by up-regulation of adhesion molecules, suppression of gelatinases activity, and induction of apoptosis in nasopharyngeal carcinoma cells
.
Int J Mol Sci
2015
;
16
:
2530
58
.
36.
Yang
CS
,
Wang
X
,
Lu
G
,
Picinich
SC
. 
Cancer prevention by tea: animal studies, molecular mechanisms and human relevance
.
Nat Rev Cancer
2009
;
9
:
429
39
.
37.
Liu
S
,
Li
H
,
Chen
L
,
Yang
L
,
Li
L
,
Tao
Y
, et al
(-)-Epigallocatechin-3-gallate inhibition of Epstein-Barr virus spontaneous lytic infection involves ERK1/2 and PI3-K/Akt signaling in EBV-positive cells
.
Carcinogenesis
2013
;
34
:
627
37
.
38.
Lin
YS
,
Tsai
YJ
,
Tsay
JS
,
Lin
JK
. 
Factors affecting the levels of tea polyphenols and caffeine in tea leaves
.
J Agric Food Chem
2003
;
26
;
51
:
1864
73
.
39.
Gong
Z
,
Kristal
AR
,
Schenk
JM
,
Tangen
CM
,
Goodman
PJ
,
Thompson
IM
. 
Alcohol consumption, finasteride, and prostate cancer risk: results from the Prostate Cancer Prevention Trial
.
Cancer
2009
;
115
:
3661
9
.
40.
Bagnardi
V
,
Randi
G
,
Lubin
J
,
Consonni
D
,
Lam
TK
,
Subar
AF
, et al
Alcohol consumption and lung cancer risk in the Environment and Genetics in Lung Cancer Etiology (EAGLE) study
.
Am J Epidemiol
2010
;
171
:
36
44
.

Supplementary data