Research Paper Volume 13, Issue 23 pp 24989—25003

Longitudinal trajectories and determinants of human fungiform papillae density

Ajoy C. Karikkineth1, *, , Eric Y. Tang2, *, , Pei-lun Kuo3, , Luigi Ferrucci3, , Josephine M. Egan2, , Chee W. Chia2, ,

  • 1 Clinical Core Laboratory and Biorepository, Intramural Research Program, National Institute on Aging, National Institutes of Health, Baltimore, MD 21224, USA
  • 2 Laboratory of Clinical Investigation, National Institute on Aging, National Institutes of Health, Baltimore, MD 21224, USA
  • 3 Translational Gerontology Branch, National Institute on Aging, National Institutes of Health, Baltimore, MD 21224, USA
* Co-first authors

Received: July 29, 2021       Accepted: November 23, 2021       Published: December 2, 2021      

https://doi.org/10.18632/aging.203741
How to Cite

Copyright: © 2021 Karikkineth et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Tongue fungiform papillae contain taste buds crucial for taste and hormone-producing taste receptor cells; therefore, they may be considered as endocrine organs and have important age-associated physiological implications. We examine the cross-sectional and longitudinal trajectories of fungiform papillae density in 1084 participants from the Baltimore Longitudinal Study of Aging using linear regression models and mixed effects models. At baseline, the mean age was 67.86 ± 14.20 years, with a mean follow-up time among those with repeat visits of 4.24 ± 1.70 years. Women (53%) were younger (66.85 ± 13.78 vs. 69.04 ± 14.61 years, p < 0.001) and had a higher fungiform papillae density than men (16.14 ± 9.54 vs. 13.77 ± 8.61 papillae/cm2, p < 0.001). Whites (67%) had a lower fungiform papillae density than non-Whites after adjusting for age and sex. Factors cross-sectionally associated with a lower fungiform papillae density included a higher waist-hip ratio (β = −8.525, p = 0.029), current smoking status (β = −5.133, p = 0.014), and alcohol use within the past 12 months (β = −1.571, p = 0.025). Longitudinally, fungiform papillae density decreased linearly with follow-up time (β = −0.646, p < 0.001). The rate of decline was not affected by sex, race, BMI, waist-hip ratio, smoking, or alcohol use. The longitudinal decline of fungiform papillae density over time needs to be explored further in order to identify other possible age-associated physiological determinants.

Introduction

Taste is one of the five fundamental primary senses and plays a crucial role in our interactions with the environment. Increased interest in research related to both taste and smell has now been fueled by the present SARS-CoV-2 pandemic because as many as 80% of people infected report taste and smell loss and/or changes [1]. An unknown but increasing percentage of people with long-Covid syndrome also continue to report taste and smell changes.

Taste signaling begins in specialized chemosensory taste receptor cells within taste buds that are located within specialized structures called taste papillae on our tongues [2]. Fungiform papillae, one type of taste papillae, are densely packed on the anterior tongue (Figure 1A). There is a wide variation in the number of fungiform papillae across individuals, with reported estimates ranging from less than 10/cm2 to greater than 200/cm2 [3, 4]. The underlying reason for this wide variation is not clear. Studies in different populations have suggested that fungiform papillae density declines with aging. However, most of the published studies on the association of fungiform papillae density with age to date have been cross-sectional [4], and mostly in small study populations of less than a hundred individuals [57]. These studies are unable to distinguish between age related changes and cohort effects. One study which looked at longitudinal changes of fungiform papillae density examined a narrow age cohort, and was mainly designed to look at changes with obesity rather than longitudinal changes with age [8].

(A) Method to assess fungiform papillae density. Blue food coloring is used to provide optimal contrast between fungiform papillae (do not take up blue dye and appear pink) and other tongue structures (coated blue). Two clear plastic hole reinforcement templates (7 mm in diameter) are placed posterior to the apex of the tongue on each side of the median sulcus. Tongue images containing the two templates were taken using a digital camera. The fungiform papillae present within the two 7 mm holes were then counted and normalized to the area of the holes and expressed as fungiform papillae density (number of fungiform papillae/cm2). (B) Representative tongue images from 16 participants with age spans from the 20s to the 90s. As shown, fungiform papillae density varies widely among individuals and across lifespan. The top panel shows individuals with higher fungiform papillae density versus those individuals with lower fungiform papillae density in the lower panel.

Figure 1. (A) Method to assess fungiform papillae density. Blue food coloring is used to provide optimal contrast between fungiform papillae (do not take up blue dye and appear pink) and other tongue structures (coated blue). Two clear plastic hole reinforcement templates (7 mm in diameter) are placed posterior to the apex of the tongue on each side of the median sulcus. Tongue images containing the two templates were taken using a digital camera. The fungiform papillae present within the two 7 mm holes were then counted and normalized to the area of the holes and expressed as fungiform papillae density (number of fungiform papillae/cm2). (B) Representative tongue images from 16 participants with age spans from the 20s to the 90s. As shown, fungiform papillae density varies widely among individuals and across lifespan. The top panel shows individuals with higher fungiform papillae density versus those individuals with lower fungiform papillae density in the lower panel.

Understanding whether fungiform papillae density decreases with age is important because beyond their role in taste perception, which by itself can affect eating behavior, fungiform papillae play other important and still not elucidated physiological roles. For example, fungiform papillae provide a source of undifferentiated stem cells to replenish taste receptor cells as they are lost [9]. Therefore, physiologically it would seem very likely that a lower fungiform papillae density would lead to lower taste perception. A lower fungiform papillae density has been associated with lower taste intensity perception in both healthy individuals [3, 10] and in those with disease conditions [11]. The lack of association between fungiform papillae density and perceived taste perception in a recent study [4] is surprising, and may be due to methodological issues in assessing taste perception. Electrogustatory thresholds, which are more objective measures of taste perception than perceived taste intensity, are associated not only with the number of fungiform papillae but also with their structure and vascularization [12, 13]. A lower taste perception is an important issue because it impairs quality of life, and it has also been associated with poor dietary choices [14], obesity [15], increased rates of smoking [16], and alcohol use [17].

As our understanding of the physiological roles of fungiform papillae and taste receptors improve, it is becoming increasingly clear that the function of fungiform papillae, taste receptor cells, and taste receptors are intertwined with other physiological systems in the body. Taste receptors are expressed in extra-gustatory tissues where they serve various functions. The bitter taste receptor TAS2R38 has been shown to regulate differentiation and de-lipidation in adipocytes [18], and has been associated with Parkinson’s disease [19] and gastric cancer [20]. Additionally, hormone receptors normally associated with extra-gustatory tissues, and hormones such as CCK [21], glucagon [22], glucagon-like peptide-1 (GLP-1) [23, 24], ghrelin [25], VIP [26], and insulin [27], have been shown to be produced in taste receptor cells. Indeed, emerging research suggests that fungiform papillae may be considered as endocrine organs.

Few studies have examined factors that may affect the rate of longitudinal changes in fungiform papillae density with aging. For example, some studies have reported that male sex [4, 28], smoking [29], and alcohol use [4] are associated with a lower fungiform papillae density, but we do not know what role these factors play in the change in fungiform papillae density over time, or the direction of their associations. Racial differences have been reported in taste sensitivity [17, 30], but as far as we know, there are no studies reporting racial differences in fungiform papillae density either cross-sectionally or longitudinally. The associations between obesity and fungiform papillae density are also complicated. Both higher taste thresholds [16, 31] and increased taste sensitivity [32] have been shown to be associated with obesity. However, it has also been reported that mice and men who were obese had a loss of fungiform papillae, possibly mediated through inflammation [8]. Therefore, the exact direction of the associations between obesity, taste, and fungiform papillae density remains unclear. With fungiform papillae being considered as possible endocrine organs, the understanding of their age-associated changes can have important implications.

We studied fungiform papillae density in men and women from the Baltimore Longitudinal Study of Aging (BLSA) with the aims to: (1) delineate longitudinal patterns of change of fungiform papillae density; (2) explore sex and racial differences in the rate of change; and (3) identify other possible determinants of the rate of longitudinal changes in fungiform papillae density.

Results

Our study population included 1084 participants, 53 percent of whom were women (Table 1). The mean age at baseline was 67.86 ± 14.20 years, with a mean follow-up time among those with repeat visits of 4.24 ± 1.70 years (range 0.90 to 8.21 years) over 3.04 ± 1.18 visits (range 2 to 9 visits); 66.97% of the study population were Whites. Women were younger than men (66.85 ± 13.78 years vs. 69.04 ± 14.61 years, p < 0.001), and fungiform papillae density at baseline was higher among women than men (16.14 ± 9.54 papillae/cm2 vs. 13.77 ± 8.61 papillae/cm2, p < 0.001).

Table 1. Characteristics of study population at baseline visit classified by sex.

WomenMenTotal
n = 579n = 505n = 1084
Fungiform papillae density, number/cm216.14 ± 9.5413.77 ± 8.6115.04 ± 9.19
Demographics
 Age, years66.85 ± 13.7869.04 ± 14.6167.86 ± 14.20
 Follow-up time, years*4.28 ± 1.694.19 ± 1.724.24 ± 1.70
 Number of visits*3.02 ± 1.143.07 ± 1.233.04 ± 1.18
 Whites, %62.5272.0866.97
Anthropometric
 BMIc, kg/m226.93 ± 5.1727.56 ± 4.127.23 ± 4.71
 % Obese (BMI >30kg/m2)25.0423.7624.45
 Waist hip ratio0.80 ± 0.070.93 ± 0.070.86 ± 0.10
Smoking
 Smokers, %
  -Never smokers68.3562.8265.77
  -Quit more than 10 years ago27.483330.06
  -Quit less than 10 years ago3.131.792.5
  -Current smoker1.042.391.67
 Pack yearsa4.08 ± 9.848.26 ± 17.966.03 ± 14.36
Alcohol
 Used alcohol in past 12 months, %79.4885.1782.12
 Alcoholic drinks in a typical week, continuousb0.96 ± 1.141.43 ± 1.351.18 ± 1.26
Alcoholic drinks in a typical week, categoricalb
 0 to 171.5554.5163.62
 2 to 326.1839.0832.18
 4 or more2.276.414.2
Values reported are mean ± SD or percentages. *Among those with repeat visits. aAmong smokers. bAmong those who drink alcohol. cBMI: body mass index.

In cross-sectional analysis (Table 2), fungiform papillae density was lower by 1.81 papillae/cm2 for every decade of baseline age (p < 0.001). Men had a lower fungiform papillae density than women independent of age (β = −1.840, p = 0.001). Age and sex adjusted cross-sectional comparisons at baseline showed that Whites had a lower fungiform papillae density (β = −1.435, p = 0.012) than non-Whites. A lower fungiform papillae density was associated with a higher waist hip ratio (β = −8.525, p = 0.029), a smoking status of being a current smoker (β = −5.133, p = 0.014) or having quit within 10 years (β = −3.693, p = 0.031) compared to never smokers, and alcohol use within the past 12 months (β = −1.571, p = 0.025). BMI expressed either as continuous or categorical, pack years of cigarette use, and number of drinks in a week expressed either as continuous or categorical were not associated with fungiform papillae density at baseline (all p > 0.05).

Table 2. Cross-sectional associations with fungiform papillae density at baseline.

βP value
Demographics
 Age–0.181<0.001
 Sex (ref = women)–1.8400.001
 Race (ref = non-White)–1.4350.012
Anthropometric measuresa
 Body mass indexb0.0760.178
 Obesity (BMI > 30kg/m2)0.6530.211
 Waist hip ratio–8.5250.029
Smoking measures
 Smoking statusc
  Never smokersRefRef
  Quit more than 10 years ago–0.6750.258
  Quit less than 10 years ago–3.6930.031
  Current smoker–5.1330.014
  Pack years of smoking–0.0230.220
Alcohol use measures
 Used alcohol in past 12 months–1.5710.025
 Alcoholic drinks in the past 12 months, continuous-0.0930.675
Alcohol drinks in the past 12 months, categorical
 0 to 1RefRef
 2 to 3–0.6890.251
 4 or more0.5250.699
aAll models were adjusted for age and sex. bSeparate age and sex adjusted models were created for each of the variables listed above. cSmoking status was analyzed as an indicator variable, with never smokers as the baseline category.

In longitudinal models (Table 3, Model 1), fungiform papillae density decreased with time of follow-up at a rate of approximately 6.46 papillae/cm2 per decade (p < 0.001). There was no interaction between the baseline age and time variable, indicating that the rate of decline was similar at all ages. The rate of decline in men and women over follow-up time was similar (Figure 2), indicated by the non-significant Sex X Time interaction term (p > 0.05). Similarly, interaction of race with time was not significant (p > 0.05), indicating that the decline in fungiform papillae density over time did not differ by race.

Table 3. Mixed effects models showing longitudinal associations of anthropometric measures.

Model 1Model 2aModel 2bModel 2c
Base modelBMIObesityWaist hip ratio
βP valueβP valueβP valueβP value
Entry age–0.181<0.001–0.166<0.001–0.180<0.001–0.176<0.001
Time–0.646<0.001–0.502<0.001–0.645<0.001–1.5960.001
Sex (ref = women)–1.3420.002–5.0770.006–1.3470.002–0.7940.107
Race (ref = non-White)–1.869<0.001–1.739<0.001–1.798<0.001–1.857<0.001
Anthropometric measures X Time
 BMI–0.7500.453
 Obesity (BMI >30 kg/m2)–0.0110.933
 Waist hip ratio1.0570.057
Constant29.207<0.00127.777<0.00128.933<0.00134.505<0.001
Entry age is the age at baseline. BMI: body mass index. ‘X’ indicates an interaction term.
Longitudinal changes in fungiform papillae density in men and women. The figure shows that men and women have decreasing fungiform papillae density over age and time, and the longitudinal rate of decline in men and women is the same.

Figure 2. Longitudinal changes in fungiform papillae density in men and women. The figure shows that men and women have decreasing fungiform papillae density over age and time, and the longitudinal rate of decline in men and women is the same.

None of the anthropometric measures had a significant interaction with follow-up time (all p > 0.05), which means that they did not affect the fungiform papillae density rate of decline. Similarly, none of the smoking variables and none of the alcohol use variables affected the rate of decline, either tested in the model one at a time or together in a multivariate model (Tables 46).

Table 4. Mixed effects models showing longitudinal associations of smoking measures.

Model 1Model 3aModel 3b
Base modelSmoking statusPack years
βP valueβP valueβP value
Entry age–0.181<0.001–0.183<0.001–0.179<0.001
Time–0.646<0.001–0.656<0.001–0.670<0.001
Sex (ref = women)–1.3420.002–1.3370.002–1.3040.003
Race (ref = non-White)–1.869<0.001–1.908<0.001–1.866<0.001
Smoking measures X Time
 Smoking status
  NeverRefRef
  Quit > 10 years ago–0.0280.822
  Quit < 10 years ago0.5650.335
  Current smoker–0.0780.893
Pack years0.0040.350
Constant29.207<0.00129.571<0.00129.142<0.001
Entry age is the age at baseline. ‘X’ indicates an interaction term. Since we were interested mainly in the longitudinal effects over time, only the time interaction terms for the smoking measures are reported in the table above. Separate models were created for each of the smoking measures – one model with the categorical smoking status, and one model with the continuous variable Pack years.

Table 5. Mixed effects models showing longitudinal associations of alcohol use measures.

Model 1Model 4aModel 4bModel 4c
Base modelAlcohol use or notDrinks in a week (continuous)Drinks in a week (categorized)
βP valueβP valueβP valueβP value
Entry age–0.181<0.001–0.182<0.001–0.168<0.001–0.168<0.001
Time–0.646<0.001–0.550<0.001–0.663<0.001–0.654<0.001
Sex (ref = women)–1.3420.002–1.3440.002–1.4810.033–1.4950.002
Race (ref = non-White)–1.869<0.001–1.7820.006–1.4630.002–1.4650.005
Alcohol use measures X Time
 Alcohol use, yes or no–0.1200.414
 Drinks in a typical week, continuous0.0040.940
 Drinks in a week, categorical
  0 or 1RefRef
  2 or 30.0490.695
  4 or more–0.0720.797
Constant29.207<0.00127.933<0.00123.151<0.00129.13911<0.001
Entry age is the age at baseline. ‘X’ indicates an interaction term. Since we were interested mainly in the longitudinal effects over time, only the time interaction terms for the alcohol use measures are reported in the table above. Separate models were created for each of the alcohol use measures – one model for “Alcohol use, yes or no”, one model for “Drinks in a typical week, continuous”, and one model for “Drinks in a week, categorical”.

Table 6. Mixed effects model showing multivariate longitudinal associations.

Multivariate model
βP value
Entry age–0.183<0.001
Time–0.2890.396
Sex (ref = women)–1.3570.002
Race (ref = non-White)–1.776<0.001
BMI0.0550.305
BMI x Time–0.0100.371
Smoke history
 Former smokerRefRef
 Quit > 10 years ago–0.4190.432
 Quit < 10 years ago–1.9250.220
 Current smoker–4.0990.037
Smoke history X Time
 Former smokerRefRef
 Quit > 10 years ago0.0040.973
 Quit < 10 years ago0.5960.310
 Current smoker–0.0730.899
Alcohol use–0.1820.759
Alcohol use X Time–0.1230.410
Constant28.121<0.001
BMI: body mass index.

Discussion

We provide the first demonstration that fungiform papillae density in humans decreases with age both cross-sectionally and longitudinally in a large, well characterized community-based cohort. Cross-sectionally, fungiform papillae density is lower in men, those with a higher waist hip ratio, and those who smoke and/or drink alcohol. Longitudinally, fungiform papillae density decreases over time. The rate of decline of fungiform papillae density is similar among men and women, and similar among Whites and non-Whites. The longitudinal decline in fungiform papillae density is not affected by obesity, smoking, or alcohol use. Therefore, we conclude that decline in fungiform papillae density is an age-related phenomenon.

Fungiform papillae density

To our knowledge, the Beaver Dam Offspring Study is the only large cohort study that has investigated the cross-sectional association of age with fungiform papillae density. The Beaver Dam Offspring Study showed a decrease in fungiform papillae density with age in 2371 participants (4). Prior studies on fungiform papillae involved smaller number of subjects (number ranged from 12 to 182) with a narrow age-range (18 to 35 years of age) [3336]. Since fungiform papillae density varies widely among individuals (Figure 1B), as noted in other studies [4, 3337], these small studies were unable to reliably estimate the effect of aging on fungiform papillae density in the general population. The mean fungiform papillae density for BLSA cohort is 15.04 ± 9.19 papillae/cm2 with a range of 0–67 papillae/cm2. Our fungiform papillae counts, obtained using the position we selected for counting the fungiform papillae – approximately 7mm posterior to the apex of the tongue and 7mm lateral to the median sulcus (Figure 1A), do match the counts obtained from autopsy tongue tissues. Specifically, Cheng and colleague sectioned the tongue into 1cm segments; the fungiform papillae density in the first two segments starting from the apex of the tongue, similar in position to where we placed our templates, is approximately 17 papillae/cm2 [38]. Miller also showed a decrease in fungiform papillae density from apex to mid-region of the tongue [3]. The wide discrepancy in fungiform papillae density among different studies is most likely due to the location that was used to quantify the fungiform papillae. Studies using area at the apex of the tongue, such as the Beaver Dam Offspring study, yielded higher counts, and studies using tongue region more posterior to the apex yielded lower counts. In addition, studies with a younger cohort would likely yield a higher count as fungiform papillae density has been shown to decrease with age [4]. Our BLSA study cohort, with its wide range of participant age, can provide a more accurate estimate of the aging effect.

Decline with age

To our knowledge, this is the first report of longitudinal trajectories of fungiform papillae density with age. While other studies have reported cross-sectional changes of fungiform papillae density across different ages, they are not able to provide information on the rate of change over time and are subjected to bias from time-invariant unmeasured confounders and cohort effect. In contrast, longitudinal data can provide us with the precious opportunity to distinguish effects of aging from cohort effects. Reports on the association of fungiform papillae density with age have generally shown a lower fungiform papillae density at older ages [3, 4, 12]. We also report a lower fungiform papillae density at older ages, and additionally demonstrate decline of fungiform papillae density over time.

Many reasons can explain the longitudinal decline of fungiform papillae density. Fungiform papillae are probably subject to the same aging-related pathologies observed in other tissues, such as on-going low-grade inflammation, deficient or defective repair after cellular insults, and decreased and defective replacement of cells from a stem cell population. The fact that fungiform papillae have the potential to be directly exposed to noxious stimuli such as hot, spicy, or pungent food may increase the likelihood of such insults. Fungiform papillae are innervated by the chorda tympani nerves, and have a rich blood supply; declines in the health of these systems with aging could contribute to the aging associated decline in fungiform papillae density [2]. Subjects with taste dysfunction were found to have flat and irregular fungiform papillae with poor blood flow [39]. Transection of the chorda tympani nerve is associated with atrophy of fungiform papillae, which begin to resemble filiform papillae [40], as well as lower fungiform papillae density [6, 11]. Loss of amiloride sensitive neural responses in Skn-1a-deficient mice was correlated with the disappearance of taste receptor cells [41]. Diseases with abnormal neurotrophic support, such as Stüve-Wiedemann syndrome and Machado-Joseph Disease, exhibit taste dysfunction and lower numbers or complete absence of fungiform papillae [42]. By observing the fungiform papillae of all these BLSA participants, the authors noted that the morphology and size are highly variable between individuals, and the distinction between different types of papillae, fungiform versus filiform for example, become less obvious with older age.

Therefore, with fungiform papillae possibly being endocrine organs, decrease in the density or number of fungiform papillae might be an indicator of health for multiple cellular systems at the local level, as well as body systems such as the cardiovascular, peripheral nervous system, endocrine system, or metabolic system; this raises the question of whether fungiform papillae density and taste function could be a marker of aging. Interestingly, certain haplotypes of the TAS2R38 receptor have recently been found to be associated with longevity [43].

Sex differences

Women had a higher fungiform papillae density than men. We also report for the first time that the rates of decline of fungiform papillae density over time is similar in men and women. Previous reports on the associations between sex, taste, and fungiform papillae density have shown conflicting findings. Some studies have reported no sex related differences in fungiform papillae density [29] or taste [44], while others did find such differences [4, 6, 28]. Women have been shown to have higher taste sensitivity, and supertasters may be more common among women [34]. EGM thresholds have been shown to be lower in young women [6].

The etiology of the sex differences is not clear. It could be related to the sex hormones estrogen and progesterone, as demonstrated by changes in taste intensity and hedonic preferences in pregnancy, such as an increase in bitter intensity [45]. Sex steroid hormones have been shown to affect gustatory processing at the levels of the taste receptor, peripheral nerve, and the central nervous system, with receptors for sex hormones prominently present in several nuclei associated with central gustatory pathways [46].

Racial differences

Whites had a lower fungiform papillae density than non-Whites, but both groups have a similar rate of decline over time. To our knowledge, racial differences in fungiform papillae density have never been reported previously. However, there have been reports of racial differences in taste sensitivity [47]. Non-Hispanic Blacks have been noted to have a higher prevalence of taste impairment [17], but have higher-rated taste sensations [30]. Taste impairment may be associated with a compensatory increase in fungiform papillae density to help with taste sensation. Proposed mechanisms for these differences include genetic differences such as differences in the frequency distribution and functional variants of TAS2R16 and TAS2R38 haplotypes [30, 48], as well as higher rates of nerve impairment among African Americans due to various factors such as trauma or upper respiratory infections [30]. It is possible that at least some of the effects of these factors on taste perception are mediated via fungiform papillae density, or that these factors could affect fungiform papillae density directly.

BMI and waist hip ratio

In our study, BMI and obesity (defined as BMI > 30 kg/m2) were not associated cross-sectionally or longitudinally with fungiform papillae density. Central obesity measured by waist hip ratio was cross-sectionally associated with lower fungiform papillae density, but had no longitudinal associations with fungiform papillae density. The discrepancy between BMI and waist hip ratio might be due to the fact that waist hip ratio is a better indicator of central or visceral obesity and overall risk than BMI [49]. Additionally, among older adults, a lower BMI might be an indicator of poor nutritional status rather than being a marker of obesity [50]. Previous reports exploring the association between obesity and fungiform papillae density have been conflicted. Proserpio et al. reported a lower fungiform papillae density in those who have a BMI ≥ 30 kg/m2, but the study was cross-sectional, with a small sample size and narrow age range [31]. Additionally, no adjustments were made for age differences among the subjects. Fischer et al. did not find an association between obesity and fungiform papillae density in a large population-based cross-sectional study that also used BMI ≥ 30 kg/m2 as the measure of obesity [4].

Kaufman et al. found that mice fed with a high fat diet became obese and had a lower fungiform papillae density [51]; this association appeared to be mediated via inflammation, specifically TNF-alpha. However, studies showed mixed results when conducted among humans. A higher BMI was not associated with a lower fungiform papillae density. On the other hand, a study using a cohort of college-aged students found that an increase in neck circumference, taken as surrogate biomarker of obesity, was associated with a longitudinal decrease in fungiform papillae density [8]. However, neck circumference could be affected by other factors such as alcohol use and could be an indicator of other disease conditions such as sleep apnea or cardiovascular disease, which may also affect fungiform papillae density. In addition, the age range of the college-aged cohort studied was narrow. Finally, there are significant differences between the biological systems for taste in mice and men. In our study, we examined a large number of participants from a broad age range, adjusting for age, sex, and race. We also examined obesity in a variety of ways, using indicators of overall obesity as well as central adiposity, and using continuous as well as categorical measures. Our participants also had a broad range of BMI and waist hip ratios.

While obesity may not be related to fungiform papillae density over time, it is still possible that the reverse is true i.e., that fungiform papillae density may be predictive of weight gain over time, and this needs to be examined in future studies. fungiform papillae density has been shown to be associated with taste intensity [3, 10], and it seems plausible that taste and resulting dietary preferences could be associated with weight gain and obesity over time. While choosing a healthy diet high in fruits and vegetables is associated with smaller gains in BMI and waist circumference [52], those prone to obesity may have a heightened hedonic response to unhealthy foods such as sweet and salty food [16, 53, 54], or energy dense foods in general [31].

Smoking

Current smokers and those who quit less than 10 years ago had a lower fungiform papillae density cross-sectionally in our study. This is in contrast to some previous studies [55], but consistent with more recent studies [4, 29]. Smokers have also been shown to have higher electrogustatory thresholds and lower taste sensitivity than non-smokers [16].

However, none of the various smoking measures that we used were associated with longitudinal changes in fungiform papillae density. Nevertheless, these results should be considered with caution because the BLSA population tend to be healthy, with a low prevalence of smokers and especially, of heavy smokers. Therefore, the lack of significant longitudinal association could be due to the small proportion of smokers that we had, although the association we observed cross-sectionally with categories of smoking status would argue against that. It is also possible that the reverse is true in that fungiform papillae density, through its associations with taste intensity, could be associated with changes in smoking preferences over time.

Alcohol use

Alcohol consumption within the past 12 months was associated with lower fungiform papillae density cross-sectionally in our study. The number of drinks per day was not associated with fungiform papillae density cross-sectionally. Additionally, none of the alcohol use variables were associated with longitudinal changes in fungiform papillae density. Again, this could be due to the small proportion of heavy drinkers in our study. In the study by Fischer et al., fungiform papillae density was lower among those consuming 4 or more drinks per day [4], but in our study, those consuming 4 or more drinks were less than 2 percent of the study population.

Strengths and limitations

The strengths of our study include a well characterized cohort of over 1000 participants with wide age-range and multiple follow-up visits. More importantly, this is the first study to follow fungiform papillae density longitudinally and the first to provide a precious opportunity to distinguish effects of aging from cohort effects. However, our study has limitations. The morphology of the fungiform papillae varies among individuals and distinction between fungiform papillae and other papillae becomes less obvious with age, which may introduce variability into the counting of fungiform papillae.

Conclusions

Fungiform papillae density decreases longitudinally over time. Further studies are needed to explore other possible determinants of the longitudinal decline in fungiform papillae density. It is also important to know whether the observed decline in fungiform papillae density leads to effects such as modifications in dietary patterns, smoking or alcohol use, and morbidity. As potential endocrine organs, fungiform papillae may provide us with a new window into human physiology that is easily accessible and observable, one that could be utilized in biopsy for further testing and examination. Along those lines, we have recently reported that SARS-CoV-2 infect taste receptor cells in taste buds of biopsied fungiform papillae [56]. With all these relevant areas of research in mind, we have developed an in vivo model and mechanistic studies are currently underway to study age-related and SARS-CoV-2 related taste loss or changes. In summary, fungiform papillae density may potentially serve as an indicator of aging and an indicator for endocrine- and disease-related processes.

Methods

Study population

Participants were selected from the BLSA, a study of normative human aging established in 1958, with a well-characterized cohort, and currently conducted by the Intramural Research Program of NIA. Participants are followed for life every 1 to 4 years, depending on age. The BLSA protocol was approved by the Institutional Review Board of the National Institutes of Health and written informed consent was obtained from all participants at each visit. Eligible participants for this study had to have measures of fungiform papillae density, and have anthropometric measures detailed below.

Fungiform papillae density

In the BLSA, we began to quantify fungiform papillae density in 2011 using modified procedure originally developed by Miller and Reedy [33]. Briefly, blue food coloring is used to provide optimal contrast between fungiform papillae (which do not take up blue dye and appear pink) and other tongue structures (coated blue). Two clear plastic hole reinforcement templates (7 mm inner diameter) are placed posterior to the apex of the tongue with one template on each side of the median sulcus (Figure 1A). Tongue images containing the two templates were taken using a digital camera. The fungiform papillae present within the two 7 mm holes were then counted and normalized to the area of the holes and expressed as fungiform papillae density (number of fungiform papillae/cm2).

BMI and waist hip ratio

Weight was measured in kilograms and height in centimeters. Body mass index (BMI) was calculated as weight in kilograms divided by the square of height in meters. BMI was used both as a continuous variable, as well as a categorical variable, defined as either obese (BMI > 30 kg/m2) or not obese (BMI ≤ 30 kg/m2). Waist circumference was measured as the minimal abdominal circumference between the lower edge of the rib cage and the iliac crests. The hip was defined as the maximal circumference around the gluteal muscles below the iliac crests. Waist hip ratio was then calculated, as a continuous measure, by dividing waist circumference by hip circumference.

Smoking status

Smoking status was examined two different ways. One method is to defined smoking using a categorical variable – “never smoker”, “quit more than 10 years ago”, “quit less than 10 years ago”, or “current smoker”. We also coded smoking as pack years of smoking, calculated as average number of packs of cigarettes smoked multiplied by years of use. Those who had never smoked were coded as 0 pack years.

Alcohol use measures

At each visit, participants were asked if they had used alcohol within the past 12 months. This was used as a binary variable consisting of “Yes” or “No” responses. If they had used alcohol, they were asked how many drinks they had on average each day, and this was then categorized in to 3 categories – “0 or 1”, “2 or 3”, and “4 or more”. Those who had not used alcohol within the past 12 months were coded as 0 drinks per day.

Statistical analysis

First, we compared baseline characteristics between men and women using Student’s t-tests and chi-squared tests for continuous variables and categorical variables, respectively. Then we examined the baseline cross-sectional associations of fungiform papillae density with various covariates using age- and sex-adjusted linear regression models. Finally, we tested for longitudinal changes and determinants of fungiform papillae density using linear mixed effects models, which can handle unbalanced, unequally spaced observations that are typical of a longitudinal study such as the BLSA.

Initial longitudinal models tested for the longitudinal changes of fungiform papillae density with baseline age, follow-up time, as well as with sex and race and their interaction terms with follow-up time. Longitudinal trajectories for fungiform papillae density for men and women were graphed using the base model. Subsequently, we tested, one at a time, for any longitudinal associations of the different measures with fungiform papillae density by testing for the significance of their interaction term with time. Finally, we tested the different covariates together in a multivariate model. All analysis was done by Stata Statistical Software (Release 16. College Station, TX: StataCorp LLC. StataCorp. 2017).

Author Contributions

CWC, JME conceptualized and designed the study. CWC, JME, LF conducted the clinical study. CWC, JME, LF, ACK, PK, EYT interpreted the data. EYT, ACK, PK performed statistical analyses. CWC, ACK wrote the original draft of the manuscript. All authors contributed to review and editing of the manuscript.

Acknowledgments

The authors would like to thank all the personnel, staff, and participants who were involved in the study. The authors are especially grateful to Lauren Brick of the NIDA/NIA Visual Media Section for illustrating the graphic figure in Figure 1A.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

Funding

This work was supported by the Intramural Research Program of the National Institutes of Health, National Institute on Aging.

References

  • 1. Lechien JR, Chiesa-Estomba CM, Beckers E, Mustin V, Ducarme M, Journe F, Marchant A, Jouffe L, Barillari MR, Cammaroto G, Circiu MP, Hans S, Saussez S. Prevalence and 6-month recovery of olfactory dysfunction: a multicentre study of 1363 COVID-19 patients. J Intern Med. 2021; 290:451–61. https://doi.org/10.1111/joim.13209 [PubMed]
  • 2. Stamps JJ. Chemosensory Function during Neurologically Healthy Aging. In: Heilman KM, Nadeau SE, eds. Cognitive Changes and the Aging Brain. Cambridge: Cambridge University Press. 2019; 68–94. https://doi.org/10.1017/9781108554350.006
  • 3. Miller IJ Jr, Reedy FE Jr. Variations in human taste bud density and taste intensity perception. Physiol Behav. 1990; 47:1213–19. https://doi.org/10.1016/0031-9384(90)90374-D [PubMed]
  • 4. Fischer ME, Cruickshanks KJ, Schubert CR, Pinto A, Klein R, Pankratz N, Pankow JS, Huang GH. Factors related to fungiform papillae density: the beaver dam offspring study. Chem Senses. 2013; 38:669–77. https://doi.org/10.1093/chemse/bjt033 [PubMed]
  • 5. Segovia C, Hutchinson I, Laing DG, Jinks AL. A quantitative study of fungiform papillae and taste pore density in adults and children. Brain Res Dev Brain Res. 2002; 138:135–46. https://doi.org/10.1016/s0165-3806(02)00463-7 [PubMed]
  • 6. Just T, Pau HW, Witt M, Hummel T. Contact endoscopic comparison of morphology of human fungiform papillae of healthy subjects and patients with transected chorda tympani nerve. Laryngoscope. 2006; 116:1216–22. https://doi.org/10.1097/01.mlg.0000224509.61099.29 [PubMed]
  • 7. Correa M, Hutchinson I, Laing DG, Jinks AL. Changes in fungiform papillae density during development in humans. Chem Senses. 2013; 38:519–27. https://doi.org/10.1093/chemse/bjt022 [PubMed]
  • 8. Kaufman A, Kim J, Noel C, Dando R. Taste loss with obesity in mice and men. Int J Obes (Lond). 2020; 44:739–43. https://doi.org/10.1038/s41366-019-0429-6 [PubMed]
  • 9. Feng P, Huang L, Wang H. Taste bud homeostasis in health, disease, and aging. Chem Senses. 2014; 39:3–16. https://doi.org/10.1093/chemse/bjt059 [PubMed]
  • 10. Delwiche JF, Buletic Z, Breslin PA. Relationship of papillae number to bitter intensity of quinine and PROP within and between individuals. Physiol Behav. 2001; 74:329–37. https://doi.org/10.1016/s0031-9384(01)00568-6 [PubMed]
  • 11. Saito T, Ito T, Ito Y, Manabe Y, Sano K. Relationship between gustatory function and average number of taste buds per fungiform papilla measured by confocal laser scanning microscopy in humans. Eur J Oral Sci. 2017; 125:44–48. https://doi.org/10.1111/eos.12329 [PubMed]
  • 12. Pavlidis P, Gouveris H, Anogeianaki A, Koutsonikolas D, Anogianakis G, Kekes G. Age-related changes in electrogustometry thresholds, tongue tip vascularization, density, and form of the fungiform papillae in humans. Chem Senses. 2013; 38:35–43. https://doi.org/10.1093/chemse/bjs076 [PubMed]
  • 13. Pavlidis P, Gouveris H, Kekes G, Maurer J. Electrogustometry thresholds, tongue tip vascularization, and density and morphology of the fungiform papillae in diabetes. B-ENT. 2014; 10:271–78. [PubMed]
  • 14. Li Q, Jin R, Yu H, Lang H, Cui Y, Xiong S, Sun F, He C, Liu D, Jia H, Chen X, Chen S, Zhu Z. Enhancement of Neural Salty Preference in Obesity. Cell Physiol Biochem. 2017; 43:1987–2000. https://doi.org/10.1159/000484122 [PubMed]
  • 15. Coltell O, Sorlí JV, Asensio EM, Fernández-Carrión R, Barragán R, Ortega-Azorín C, Estruch R, González JI, Salas-Salvadó J, Lamon-Fava S, Lichtenstein AH, Corella D. Association between taste perception and adiposity in overweight or obese older subjects with metabolic syndrome and identification of novel taste-related genes. Am J Clin Nutr. 2019; 109:1709–23. https://doi.org/10.1093/ajcn/nqz038 [PubMed]
  • 16. Park DC, Yeo JH, Ryu IY, Kim SH, Jung J, Yeo SG. Differences in taste detection thresholds between normal-weight and obese young adults. Acta Otolaryngol. 2015; 135:478–83. https://doi.org/10.3109/00016489.2014.975370 [PubMed]
  • 17. Liu G, Zong G, Doty RL, Sun Q. Prevalence and risk factors of taste and smell impairment in a nationwide representative sample of the US population: a cross-sectional study. BMJ Open. 2016; 6:e013246. https://doi.org/10.1136/bmjopen-2016-013246 [PubMed]
  • 18. Cancello R, Micheletto G, Meta D, Lavagno R, Bevilacqua E, Panizzo V, Invitti C. Expanding the role of bitter taste receptor in extra oral tissues: TAS2R38 is expressed in human adipocytes. Adipocyte. 2020; 9:7–15. https://doi.org/10.1080/21623945.2019.1709253 [PubMed]
  • 19. Cossu G, Melis M, Sarchioto M, Melis M, Melis M, Morelli M, Tomassini Barbarossa I. 6-n-propylthiouracil taste disruption and TAS2R38 nontasting form in Parkinson's disease. Mov Disord. 2018; 33:1331–39. https://doi.org/10.1002/mds.27391 [PubMed]
  • 20. Choi JH, Lee J, Choi IJ, Kim YW, Ryu KW, Kim J. Genetic Variation in the TAS2R38 Bitter Taste Receptor and Gastric Cancer Risk in Koreans. Sci Rep. 2016; 6:26904. https://doi.org/10.1038/srep26904 [PubMed]
  • 21. Herness S, Zhao FL. The neuropeptides CCK and NPY and the changing view of cell-to-cell communication in the taste bud. Physiol Behav. 2009; 97:581–91. https://doi.org/10.1016/j.physbeh.2009.02.043 [PubMed]
  • 22. Elson AE, Dotson CD, Egan JM, Munger SD. Glucagon signaling modulates sweet taste responsiveness. FASEB J. 2010; 24:3960–69. https://doi.org/10.1096/fj.10-158105 [PubMed]
  • 23. Shin YK, Martin B, Golden E, Dotson CD, Maudsley S, Kim W, Jang HJ, Mattson MP, Drucker DJ, Egan JM, Munger SD. Modulation of taste sensitivity by GLP-1 signaling. J Neurochem. 2008; 106:455–63. https://doi.org/10.1111/j.1471-4159.2008.05397.x [PubMed]
  • 24. Calvo SS, Egan JM. The endocrinology of taste receptors. Nat Rev Endocrinol. 2015; 11:213–27. https://doi.org/10.1038/nrendo.2015.7 [PubMed]
  • 25. Shin YK, Martin B, Kim W, White CM, Ji S, Sun Y, Smith RG, Sévigny J, Tschöp MH, Maudsley S, Egan JM. Ghrelin is produced in taste cells and ghrelin receptor null mice show reduced taste responsivity to salty (NaCl) and sour (citric acid) tastants. PLoS One. 2010; 5:e12729. https://doi.org/10.1371/journal.pone.0012729 [PubMed]
  • 26. Martin B, Shin YK, White CM, Ji S, Kim W, Carlson OD, Napora JK, Chadwick W, Chapter M, Waschek JA, Mattson MP, Maudsley S, Egan JM. Vasoactive intestinal peptide-null mice demonstrate enhanced sweet taste preference, dysglycemia, and reduced taste bud leptin receptor expression. Diabetes. 2010; 59:1143–52. https://doi.org/10.2337/db09-0807 [PubMed]
  • 27. Doyle ME, Fiori JL, Gonzalez Mariscal I, Liu QR, Goodstein E, Yang H, Shin YK, Santa-Cruz Calvo S, Indig FE, Egan JM. Insulin Is Transcribed and Translated in Mammalian Taste Bud Cells. Endocrinology. 2018; 159:3331–39. https://doi.org/10.1210/en.2018-00534 [PubMed]
  • 28. Shen Y, Kennedy OB, Methven L. Exploring the effects of genotypical and phenotypical variations in bitter taste sensitivity on perception, liking and intake of brassica vegetables in the UK. Food Qual Prefer. 2016; 50:71–81. https://doi.org/10.1016/j.foodqual.2016.01.005
  • 29. Pavlidis P, Gouveris C, Kekes G, Maurer J. Changes in electrogustometry thresholds, tongue tip vascularization, density and form of the fungiform papillae in smokers. Eur Arch Otorhinolaryngol. 2014; 271:2325–31. https://doi.org/10.1007/s00405-014-3003-9 [PubMed]
  • 30. Williams JA, Bartoshuk LM, Fillingim RB, Dotson CD. Exploring Ethnic Differences in Taste Perception. Chem Senses. 2016; 41:449–56. https://doi.org/10.1093/chemse/bjw021 [PubMed]
  • 31. Proserpio C, Laureati M, Bertoli S, Battezzati A, Pagliarini E. Determinants of Obesity in Italian Adults: The Role of Taste Sensitivity, Food Liking, and Food Neophobia. Chem Senses. 2016; 41:169–76. https://doi.org/10.1093/chemse/bjv072 [PubMed]
  • 32. Fischer ME, Cruickshanks KJ, Pankow JS, Pankratz N, Schubert CR, Huang GH, Klein BE, Klein R, Pinto A. The associations between 6-n-propylthiouracil (PROP) intensity and taste intensities differ by TAS2R38 haplotype. J Nutrigenet Nutrigenomics. 2014; 7:143–52. https://doi.org/10.1159/000371552 [PubMed]
  • 33. Miller IJ Jr, Reedy FE Jr. Quantification of fungiform papillae and taste pores in living human subjects. Chem Senses. 1990; 15:281–94. https://doi.org/10.1093/chemse/15.3.281
  • 34. Bartoshuk LM, Duffy VB, Miller IJ. PTC/PROP tasting: anatomy, psychophysics, and sex effects. Physiol Behav. 1994; 56:1165–71. https://doi.org/10.1016/0031-9384(94)90361-1 [PubMed]
  • 35. Essick GK, Chopra A, Guest S, McGlone F. Lingual tactile acuity, taste perception, and the density and diameter of fungiform papillae in female subjects. Physiol Behav. 2003; 80:289–302. https://doi.org/10.1016/j.physbeh.2003.08.007 [PubMed]
  • 36. Zhang GH, Zhang HY, Wang XF, Zhan YH, Deng SP, Qin YM. The relationship between fungiform papillae density and detection threshold for sucrose in the young males. Chem Senses. 2009; 34:93–99. https://doi.org/10.1093/chemse/bjn059 [PubMed]
  • 37. Miller IJ Jr. Human taste bud density across adult age groups. J Gerontol. 1988; 43:B26–30. https://doi.org/10.1093/geronj/43.1.m26 [PubMed]
  • 38. Cheng LHH, Robinson PP. The distribution of fungiform papillae and taste buds on the human tongue. Arch Oral Biol. 1991; 36:583–89. https://doi.org/10.1016/0003-9969(91)90108-7 [PubMed]
  • 39. Negoro A, Umemoto M, Fukazawa K, Terada T, Sakagami M. Observation of tongue papillae by video microscopy and contact endoscopy to investigate their correlation with taste function. Auris Nasus Larynx. 2004; 31:255–59. https://doi.org/10.1016/j.anl.2004.01.009 [PubMed]
  • 40. Nagato T, Matsumoto K, Tanioka H, Kodama J, Toh H. Effect of denervation on morphogenesis of the rat fungiform papilla. Acta Anat (Basel). 1995; 153:301–309. https://doi.org/10.1159/000147739 [PubMed]
  • 41. Ohmoto M, Jyotaki M, Foskett JK, Matsumoto I. Sodium-Taste Cells Require Skn-1a for Generation and Share Molecular Features with Sweet, Umami, and Bitter Taste Cells. eNeuro. 2020; 7:ENEURO.0385-20.2020. https://doi.org/10.1523/ENEURO.0385-20.2020 [PubMed]
  • 42. Gardiner J, Barton D, Vanslambrouck JM, Braet F, Hall D, Marc J, Overall R. Defects in tongue papillae and taste sensation indicate a problem with neurotrophic support in various neurological diseases. Neuroscientist. 2008; 14:240–50. https://doi.org/10.1177/1073858407312382 [PubMed]
  • 43. Melis M, Errigo A, Crnjar R, Pes GM, Tomassini Barbarossa I. TAS2R38 bitter taste receptor and attainment of exceptional longevity. Sci Rep. 2019; 9:18047. https://doi.org/10.1038/s41598-019-54604-1 [PubMed]
  • 44. Wang JJ, Liang KL, Lin WJ, Chen CY, Jiang RS. Influence of age and sex on taste function of healthy subjects. PLoS One. 2020; 15:e0227014. https://doi.org/10.1371/journal.pone.0227014 [PubMed]
  • 45. Duffy VB, Bartoshuk LM, Striegel-Moore R, Rodin J. Taste changes across pregnancy. Ann N Y Acad Sci. 1998; 855:805–809. https://doi.org/10.1111/j.1749-6632.1998.tb10663.x [PubMed]
  • 46. Martin LJ, Sollars SI. Contributory role of sex differences in the variations of gustatory function. J Neurosci Res. 2017; 95:594–603. https://doi.org/10.1002/jnr.23819 [PubMed]
  • 47. Doty RL, Applebaum S, Zusho H, Settle RG. Sex differences in odor identification ability: a cross-cultural analysis. Neuropsychologia. 1985; 23:667–72. https://doi.org/10.1016/0028-3932(85)90067-3 [PubMed]
  • 48. Wang JC, Hinrichs AL, Bertelsen S, Stock H, Budde JP, Dick DM, Bucholz KK, Rice J, Saccone N, Edenberg HJ, Hesselbrock V, Kuperman S, Schuckit MA, et al. Functional variants in TAS2R38 and TAS2R16 influence alcohol consumption in high-risk families of African-American origin. Alcohol Clin Exp Res. 2007; 31:209–15. https://doi.org/10.1111/j.1530-0277.2006.00297.x [PubMed]
  • 49. Huxley R, Mendis S, Zheleznyakov E, Reddy S, Chan J. Body mass index, waist circumference and waist:hip ratio as predictors of cardiovascular risk--a review of the literature. Eur J Clin Nutr. 2010; 64:16–22. https://doi.org/10.1038/ejcn.2009.68 [PubMed]
  • 50. Cook Z, Kirk S, Lawrenson S, Sandford S. Use of BMI in the assessment of undernutrition in older subjects: reflecting on practice. Proc Nutr Soc. 2005; 64:313–17. https://doi.org/10.1079/pns2005437 [PubMed]
  • 51. Kaufman A, Choo E, Koh A, Dando R. Inflammation arising from obesity reduces taste bud abundance and inhibits renewal. PLoS Biol. 2018; 16:e2001959. https://doi.org/10.1371/journal.pbio.2001959 [PubMed]
  • 52. Newby PK, Muller D, Hallfrisch J, Qiao N, Andres R, Tucker KL. Dietary patterns and changes in body mass index and waist circumference in adults. Am J Clin Nutr. 2003; 77:1417–25. https://doi.org/10.1093/ajcn/77.6.1417 [PubMed]
  • 53. Salbe AD, DelParigi A, Pratley RE, Drewnowski A, Tataranni PA. Taste preferences and body weight changes in an obesity-prone population. Am J Clin Nutr. 2004; 79:372–78. https://doi.org/10.1093/ajcn/79.3.372 [PubMed]
  • 54. Fischer ME, Cruickshanks KJ, Schubert CR, Pinto A, Huang GH, Klein BE, Klein R, Pankow JS. The association of taste with change in adiposity-related health measures. J Acad Nutr Diet. 2014; 114:1195–202. https://doi.org/10.1016/j.jand.2014.04.013 [PubMed]
  • 55. Konstantinidis I, Chatziavramidis A, Printza A, Metaxas S, Constantinidis J. Effects of smoking on taste: assessment with contact endoscopy and taste strips. Laryngoscope. 2010; 120:1958–63. https://doi.org/10.1002/lary.21098 [PubMed]
  • 56. Doyle ME, Appleton A, Liu QR, Yao Q, Mazucanti CH, Egan JM. Human Type II Taste Cells Express Angiotensin-Converting Enzyme 2 and Are Infected by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Am J Pathol. 2021; 191:1511–19. https://doi.org/10.1016/j.ajpath.2021.05.010 [PubMed]