Fabella prevalence rate increases over 150 years, and rates of other sesamoid bones remain constant: a systematic review

Abstract The fabella is a sesamoid bone located behind the lateral femoral condyle. It is common in non‐human mammals, but the prevalence rates in humans vary from 3 to 87%. Here, we calculate the prevalence of the fabella in a Korean population and investigate possible temporal shifts in prevalence rate. A total of 52.83% of our individuals and 44.34% of our knees had fabellae detectable by computed tomography scanning. Men and women were equally likely to have a fabella, and bilateral cases (67.86%) were more common than unilateral ones (32.14%). Fabella presence was not correlated with height or age, although our sample did not include skeletally immature individuals. Our systematic review yielded 58 studies on fabella prevalence rate from 1875–2018 which met our inclusion criteria, one of which was an outlier. Intriguingly, a Bayesian mixed effects generalized linear model revealed a temporal shift in prevalence rates, with the median prevalence rate in 2000 (31.00%) being ~ 3.5 times higher than that in 1900 (7.64%). In all four countries with studies before and after 1960, higher rates were always found after 1960. Using data from two other systematic reviews, we found no increase in prevalence rates of 10 other sesamoid bones in the human body, indicating that the increase in fabella prevalence rate is unique. Fabella presence/absence is due to a combination of genetic and environmental factors: as the prevalence rates of other sesamoid bones have not changed in the last 100 years, we postulate the increase in fabella prevalence rate is due to an environmental factor. Namely, the global increase in human height and weight (due to improved nutrition) may have increased human tibial length and muscle mass. Increases in tibial length could lead to a larger moment arm acting on the knee and on the tendons crossing it. Coupled with the increased force from a larger gastrocnemius, this could produce the mechanical stimuli necessary to initiate fabella formation and/or ossification.


Introduction
The fabella (Latin for 'little bean') is a sesamoid bone located in the knee joint behind the lateral femoral condyle. Embedded in the tendon of the lateral head of the gastrocnemius muscle, it is stabilized by the fabellofibular ligament, connecting the distal insertion of the fabella to the fibular head (Minowa et al. 2004;Piyawinijwong et al. 2012;Driessen et al. 2014;Hauser et al. 2015;Kurto glu et al. 2015) and the posterior capsule of the knee. In rare instances, it serves as an additional origin for a muscle bundle of the popliteal muscle (Duc et al. 2004). Fabella prevalence in humans ranges from 3 to 87% (Silva et al. 2010;Zeng et al. 2012), making it a normal variant in human anatomy. The highest rates reported are in Asians and Australians, and the lowest rates in Europeans and South Americans (Minowa et al. 2004;Silva et al. 2010;Zeng et al. 2012;Hauser et al. 2015). Although its exact function is unknown, the fabella is more common in non-human mammals (Pearson & Davin, 1921;Sarin et al. 1999), which has prompted functional and evolutionary debates about the role of the fabella in locomotion (Sarin et al. 1999;Jin et al. 2017).
Most studies reporting on prevalence rates in humans have determined the presence of the fabella through surgeries/dissections (Agathangelidis et al. 2016), X-rays (Pancoast, 1909), computed tomography (CT) scans (Hauser et al. 2015), and magnetic resonance imaging (MRI) scans (Hedderwick et al. 2017). Ultrasound (Sekiya et al. 2002) and PET-CT (Usmani et al. 2017) have been used to examine the fabella, but no studies have employed these methods to calculate fabella prevalence rate. Problems in calculating prevalence rates can occur depending on the method employed. For example, fabellae are sometimes so small they are difficult to detect on MRI scans, particularly if the knee is not positioned correctly (Yu et al. 1996;Ehara, 2014). Additionally, as the fabella is occasionally cartilaginous (Jin et al. 2017), its presence may not always be detected by X-rays or CT scans. For example, a recent study on a Chinese population reported that 57.9% of the cartilaginous fabella were not visible on radiographs (Zeng et al. 2012). This highlights an issue with comparing prevalence rates between studies, as some consider only osseous fabellae, whereas others also consider cartilaginous ones. Comparing prevalence rates is further complicated as several studies do not specify whether the prevalence rates included cartilaginous with more recent ones (Hessen, 1946). This is true regardless of whether the more recent studies calculate their prevalence rates using bony or bony and cartilaginous fabellae.
Finally, the fabella can be problematic in cases of total knee arthroplasty (Larson & Becker, 1993;Wang, 1995;Erichsen, 1997;Segal et al. 2004;Theodorou et al. 2005;Jung et al. 2007;Hou, 2016;Kwee et al. 2016;Okano et al. 2016). The absence of an articulating groove in the back of the lateral femoral condyle, which serves to stabilize the fabella and is present in some anatomical variants (e.g. Chew et al. 2014), leads to a fabella medio-lateral instability, causing it to painfully 'snap' over the replacement condyle. The reason for this pain is not known. Hou (2016) recently investigated the effects of the fabella on posterolateral pain and palsy of common peroneal nerve following total knee arthroplasty. During trials, fabellae were excised from some patients but left in others. Post-surgery, posterolateral pain and palsy of common peroneal nerve were only observed in patients who still had fabellae. Accordingly, Hou recommended removing the fabella when knee replacement surgery is performed.
Here, we present the prevalence rate of the fabella in a population of Koreans using a randomized previously gathered dataset. As factors related to sex and length/speed of growth and development are correlated to bone formation (i.e. men are generally taller, and tall people have longer bones that are generally mechanically loaded more heavily), we investigate the effects of sex, age, and height on fabella prevalence rate. In addition, as other studies have reported higher rates for bilateral fabellae than for unilateral ones (Phukubye & Oyedele, 2011;Piyawinijwong et al. 2012;Egerci et al. 2017), we investigate whether bilateral or unilateral fabellae are more common.
To contextualize our prevalence rate results, we performed a systematic review to determine how Koreans compare with other populations, investigated possible changes in prevalence rate through time, and compared this with prevalence rates of other sesamoid bones.

Prevalence rate
Sample A randomized sample of previously collected CT scans, totalling 212 knees from 106 individuals (f = 55, m = 51), were investigated for the presence of the fabella (Dai et al. 2012). Scans were gathered as part of a larger project to examine human anatomy, and represent a randomized sample of Koreans. Ages of the individuals ranged from 21 to 60 years (mean/median = 52.45/55 years) and heights from 146 to 178 cm (mean/median = 160.65/160 cm; Table 1). The resolution of the scans ranged from 0.8220 9 0.8220 mm 2 to 0.9626 9 0.9626 mm 2 with a slice thickness of 1.0000 mm.
CT scans prohibit the distinction between highly dense, cartilaginous and ossified fabellae, and detection of lower density, cartilaginous fabellae. Accordingly, we made no distinction between cartilaginous and bony fabellae. As it is likely that many cartilaginous fabellae are missed by CT scans, this reported prevalence rate represents a minimum rate for this sample.

Data collection
We recorded the presence/absence of the fabella on both right and left knees. Although the fabella is located behind the lateral condyle of the femur, the rest of the knee was inspected for sesamoid bones as (1) fabella presence is often correlated with the presence of other sesamoid bones (Sarin et al. 1999) and (2) some studies have reported fabellae in the medial head of the gastrocnemius (Kawashima et al. 2007;Zeng et al. 2012). Due to the resolution of the CT scans and the miniscule size of some of the fabellae (Fig. 1), fabella dimensions were not measured.

Data sources
To complete a comprehensive literature review, the following search strategies were used for the systematic review: (1) computer search of databases and (2) review of bibliographies of all articles retrieved. Textbooks were not utilized unless they specifically came up in the computer search or bibliographies. This strategy is in accordance with Stroup et al. (2000).

Computer search
We searched google.scholar.co.uk for articles pertaining to the fabella in April 2018 and updated our results in October 2018. The search term fabella yielded 9140 results, many of which were not relevant to this study. To narrow the results, the following search terms were employed: fabella sesamoid, fabellae sesamoid -fabella, fabella knee -sesamoid, fabellae knee -sesamoid -fabella, cyamella -fabella -fabellae, fabella incidence rate -sesamoid -knee, fabellae incidence rate -sesamoid -knee -fabella, fabella prevalence ratesesamoid -knee -incidence, and fabellae prevalence rate -sesamoidknee -fabella -incidence. A hyphen before a word indicates the following word was excluded from that search, preventing the same article/citation from appearing in multiple searches.
Abstracts were reviewed first by M.A.B., and later by E.D.F. if necessary, and selected for further review if they met the following criteria: (1) the studies were on humans, (2) the studies were anatomical or medical in nature, (3) not case studies, and (4) a link was provided through which the article could be accessed. Full texts were reviewed by M.A.B.. Studies were excluded if they (1) did not report on prevalence or incidence rates based on data gathered in that study, (2) calculated rates with samples <12 knees, (3) did not report on the number of knees analysed in the study or (4) did not use a randomized sample (e.g. studies on fabella syndrome). If studies were not written in English, they were translated either by people fluent in those languages or using google translate. While imperfect, google translate worked well enough to extract the necessary data.

Review of bibliographies
If papers referenced other studies on prevalence rates, full texts of those studies were obtained through scholar.google.co.uk or interlibrary loan. If the original studies could not be located (as was the case with several older studies), data were extracted from the paper that referenced the original study, when possible. If not possible, the original study was excluded.

Statistical analysis
Korean dataset R and RSTUDIO were used for statistical analyses (R Team, 2015; R Core Team, 2018). Prevalence rates for the Korean population were calculated as the percentage of knees with fabellae and individuals with fabellae. For those with fabellae, the percentage of bilateral and unilateral cases was calculated. Pearson's chi-square tests were performed using the chisq.test function (simulate P-value = TRUE, B = 10 000) to investigate the correlation between sex and prevalence rates. The simulate-P-value simulates datasets using Monte Carlo simulations to estimate P-values for chi-square tests. A Pearson's chi-square test was performed with the unilateral data to investigate bilateral asymmetry. Point biserial correlations with

Systematic review
Published studies report on both knee and individual prevalence rates. We chose to transform all prevalence rates to knee prevalence rates for two reasons. First, several studies prior to 1950, reported knee and not individual prevalence rates, and it is not possible to know how many individuals were in the sample, especially considering some samples had an odd number of knees. Secondly, some studies were carried out on single legs and on whole individuals (e.g. X-rays taken of just one knee, or only one knee was available for dissection). Studies in which this transformation could not be performed were excluded.
A Bayesian mixed effects generalized linear model was executed to investigate temporal changes in prevalence rate while accounting for the random effects of country and method for data collection using the rethinking package in R (McElreath, 2016). A logistic regression was utilized to ensure prevalence rates were between 0 and 1. The model predicts the number of fabellae present for a given sample size, allowing the regression to take study sample size into account, which varied greatly, from 12 to 2340 knees. 'Country' was the country in which the study was conducted, unless the study specified the race of the sample used. If more than one race was specified, each race was treated as an individual data point (Mia skiewicz & Partyka, 1984). 'Method' was the method for data collection, either X-ray, CT scans, MRI or anatomical dissection. If more than one method was used, each method was treated as an individual data point (Chew et al. 2014;Hedderwick et al. 2017).
The Bayesian model followed protocol set out by McElreath (2016). The map2stan function was used to create a binomial distribution using the number of knees and fabellae in the published studies. The probability that a fabella would be present was defined as follows: where broad, weakly regularizing priors were used for the fixed (a, b) and random (a country , a method ) effects (see Data S1 for further details). Markov-chain Monte Carlo (MCMC) estimation was used to estimate the posterior probability distribution (4 chains, 10 000 iterations, 1000 iterations warmup).
To determine whether there were any outliers, a Pearson's linear regression was run between the natural log of knee and fabella count. A log transformation was used, as the density plot of study sample sizes were non-normally distributed. If any study had an unusually high or low variance (i.e. an unusually high or low number of fabellae for that sample size), it was considered an outlier and removed from further analyses. After outliers were removed, the missing data concerning 'method for data collection' was imputed using the mice package in R (van Buuren & Groothuis-Oudshoorn, 2011) . Two methods were used to impute the data, 'pmm' and 'lmer,' to create 20 imputed datasets (10/method). The consensus results were used for further analysis.
One study investigating prevalence rates in a Japanese population identified a correlation between fabella prevalence rate and age, finding a lower prevalence rate in younger (< 50 years, 31%) than older individuals (> 50 years, 47%) (Kato et al. 2012). In our dataset, individuals < 50 years old were no more or less likely to have a fabella than were individuals > 50 years old (younger = 23/94, older = 35/118, v 2 = 0.7099, P = 4448). In this study, fabellae ranged in size from small (just a few pixels) to large (Fig. 1). In general, fabellae did not appear to articulate with the lateral femoral condyle. However, the CT scans were acquired postmortem, and soft tissues were severely deformed in most individuals, making it possible that some fabellae would have articulated with the condyle in life but were separated in death. Some large fabellae were still articulated with the posterior surface of the lateral femoral condyle, the most drastic of which was observed in female 005 (Fig. 2), which shows a large articulating surface in the femur.

Systematic review
Our searches revealed 2631 abstracts on fabella prevalence rates between 1875 and 2018, written in seven languages (English, German, French, Spanish, Italian, Japanese, and Therefore, there are 97 fabellae in total. j 11/31 individuals had fabellae. k Estimated 18 fabellae based on an prevalence rate of 27%. l Reports on fabellae in medial headignored here, as it is unusually high, particularly given the lack of medial fabellae in other studies. m Reports a couple of medial fabellaenot possible to tease them out, prevalence rate may be too high. n Prevalence rate of 31.25% (25/80) for individual. Unknown if one or two knees were inspected per individual.
A total of 185 full-text articles/conference proceedings were reviewed, 66 of which reported on fabella prevalence rates. Of the 66 studies, five were discarded from further analysis as they did not fit the inclusion criteria. Pancoast (1909) and Hukuda et al. (1983) reported that 67/529 and 11/31 individuals from the USA and Japan, respectively, had fabellae, but these could not be transformed into a knee prevalence rate (Pancoast, 1909;Hukuda et al. 1983). Chew et al. (2014 reported on a prevalence rate of 31.25% (25/ 80) in 'Asians', but we could not determine whether this was an individual or knee rate (Chew et al. 2014). Siina (1931, taken from Table 1 (tabelle I) in Hessen (1946), and Munshi et al. (2003), had a sample sizes of 10 and 8 knees, respectively (Munshi et al. 2003). Finally, three studies claimed to have data on fabella presence/absence, but the data were not present, at least not in the versions of the papers we had access to (Nishimura & Shimizu, 1963;Orzincolo et al. 1987;Osti et al. 2013). Our final analysis included 21 676 knees and represented studies done in 27 countries. It should be noted that Taiwan was part of Japan from 1895 to 1945, at the time of studies of Kitahara (1935) andHanamuro (1927). According to Hessen (1946), Kitahara's (1935) sample was 'Formosawilde', indicating it consisted of There were no sex-based differences. Of the 56 individual cases, bilateral cases were significantly more prevalent than unilateral ones. Bilateral cases were more prevalent than unilateral in males (n = 24), but there was no difference in females (n = 32). Within unilateral cases, fabellae were equally likely to be present in the right or left knee. There were no differences between the sexes (see text for test statistics and P-values).  the natives of Taiwan. As such, we have classified this sample as being from Taiwan, even though no such political entity existed at the time. According to Hessen (1946), Hanamuro (1927 included individuals from Formosa as well, but classified them as 'Formosa-Chinesen', indicating they were immigrants from mainland China into Taiwan. As such, we classified their sample as being from China. A summary of prevalence rates reported in the literature can be found in Table 2. We identified one outlier in our dataset (Fig. 3), as the number of fabellae (n = 2) was exceptionally low for that number of knees (n = 62). This is not to say the data are incorrect, only that it is an outlier from the other 56 studies, and thus was excluded from further analyses.
There were five studies for which the method remained 'unknown', either because the method was not mentioned in the study or we were not able to obtain the original study and identify the method. We assumed Parsons & Keith (1897) used anatomical dissections, as the X-ray was invented in 1895, making it unlikely they used X-rays to collect their data. For the four other studies, all imputed datasets yielded consistent results for Sugiyama (1914), Ooi/Oi (1930), andMikami (1932), classifying the first two as anatomical dissections and the third as X-ray. According to the imputed data, Pichler (1918) was categorized as X-ray 15/20 times, MRI 3/20 times, and CT 2/20 times. As MRI and CT scanners were not invented in 1918, we assume Pichler used X-rays to collect their data.
The logistic regression revealed a strong increase in prevalence rates through time (P slope < 0.01, P intercept < 0.01; Fig. 4). The r code and raw data used to conduct the analysis are available in the Data S1 and Table S1. Assuming median random and fixed effects, the results show that: logitðPrevalenceÞ ¼ À33:3390 þ ð1:6314 Ã 10 À2 Þ Ã Year Interestingly, recent studies show a higher variance in prevalence rates compared with older studies. This is because there is an increase in maximum prevalence rates, with no real increase in minimum prevalence rates, causing a larger spread of the data. Although different populations were examined before and after 1960, and a genetic component may be involved in population-related fabella prevalence rates (Sarin et al. 1999), the authors are confident that the observed increase in fabella prevalence rates is not affected by these factors, as described below.
Prevalence rates were reported in four countries both before and after 1960: China, Japan, Korea, and USA. For China and Korea, there was one study before and one study after 1960; in both countries, the more recent study had a higher prevalence rate (Fig. 5). For USA and Japan, there were several studies both before and after 1960, and Pearson's linear regressions revealed positive relationships between prevalence rate and time in both countries. As there were relatively few studies in each country, we chose simpler Pearson's linear regressions in lieu of binomial mixed effect models to provide a visualization of the average change in prevalence rate over time. As random effects were ignored, little faith should be put in the regression equations and their P-values (Fig. 5). Although it is not possible to hold genetics constant between the older and newer studies, particularly in countries that have large levels of genetic diversity, such as USA, this evidence supports the idea that the increase in prevalence rates is not a by-product of different populations being used in studies before and after 1960.
Why would there be an increase in fabella prevalence rate over time? Skeletal phenotypes result from a combination of genetic and environmental factors. Although fabella formation appears to have a genetic component, it is improbable a genetic mutation is responsible for the worldwide increase in prevalence rates; the probability of a mutation occurring in Homo sapiens and spreading throughout the entire species in the past 100 years is an unprecedented and unlikely scenario.
Environmentally, it is possible that the increase in prevalence rates could be due to a hormonal or epigenetic shift. Fig. 3 Plot of the natural log of sample size (number of knees) and number of fabellas for the 57 studies considered for this analysis. A Pearson's correlation revealed a statistically significant relationship between the two variables (y = 0.82350 * x À0.60879; tvalue = 11.149, P = 2.96e-16), with an intercept that is not statistically different from zero (t-value = À1.541, P = 0.129). The data for Brazil (Silva et al., 2010) represent an outlier for this dataset.
Since the mid-20th century, there has been a marked increase in plastic usage (Zalasiewicz et al. 2016), and plastics are known to have deleterious effects on growth and development. For example, several chemicals found in plastics are known to disrupt hormonal pathways in vertebrates and other animals. It is therefore possible that plastics could have affected human skeletal growth and development, and be responsible for the increase in fabella prevalence rates. If a hormonal or epigenetic pathway were responsible, it is reasonable to assume the effects would be systematic, influencing all the sesamoid bones in the human body.
To test this idea, we investigated temporal changes in prevalence rates in other sesamoid bones in the human body. We identified two systematic reviews investigating sesamoid bone prevalence rates in the human hand (Yammine, 2014) and foot (Yammine, 2015) with data from 1892 onwards. Using these reviews, we investigated temporal changes in prevalence rate in six sesamoid bones in the hand and four sesamoid bones in the foot.
Due to the low number of studies investigating prevalence rates for these bones (16 across 120 years for the hand and 16 across 121 years for the foot), we ran binomial regressions without random effects using the glm function in R to investigate possible temporal changes. Our analyses revealed there were no temporal changes in sesamoid bone prevalence rates in either the hand or the foot (Tables 5  and 6; Figs 6 and 7). These results imply the increase in fabella prevalence rate does not have a hormonal or epigenetic origin, and the increase in fabella prevalence rate is unique.
Sesamoid bones form in areas of high mechanical stimuli, such as pressure, friction or stress (Sarin & Carter, 2000), and act to modify/reduce pressure, friction or stress. It is therefore possible that some change in mechanical loading could have caused an increase in fabella prevalence rate. Differences in loading could be due to differences in kinematics or muscle mass/bone lengths. We do not believe the differences are due to kinematics for the following reasons. First, it is unlikely that all humans, worldwide, have begun to move their lower limbs in a consistently different manner in the last 100 years. Secondly, there appears to be no correlation between magnitude of mechanical loading over one's lifetime and fabella presence in people today, with fabellae being found in both active individuals, such as non-professional (Dashefsky, 1977;Kuur, 1986) and Olympic level athletes (Zenteno et al. 2010), and inactive individuals, such as foetuses (Minowa et al. 2005;Jin et al. 2017) and the elderly (Laird, 1991;Ando et al. 2017). Finally, unlike in other mammals, the fabella likely offers no significant mechanical advantage in humans, as when excised (common practice to address fabella syndrome), no ill mechanical effects are observed (Weiner & Macnab, 1982;Zenteno et al. 2010;Agathangelidis et al. 2016;Okano et al. 2016). This implies there may be no significant mechanical, evolutionary advantage to having a fabella (Sarin et al. 1999). It is, however, possible global changes in muscle mass/ bone lengths could be responsible. Worldwide, there has been a general increase in dietary quality and nutrition over the last 100 years, which has allowed humans to come much closer to achieving their genetic potential. 1 This Fig. 4 There is a statistically significant relationship between prevalence rate and time, with people being, on average, nearly 3.5 times more likely to have a fabella in 2018 than in 1918. The confidence intervals are, from widest to narrowest, 99, 95, 75, and 50%. The raw data used to create this figure are available in the Table S2.  5 Four countries (China, Japan, Korea, and USA) had prevalence rates reported both before and after 1960. For China and Korea, there was only one study before and one study after 1960, and the lines connect these studies. For the USA and Japan, there were several, and Pearson's linear regressions were run. There is no statistically significant relationship in the USA (P = 0.0793), but there is a significant relationship in Japan (prevalence rates = 0.5064 * year À947.9; P = 2.25e-4).

1
The term 'genetic potential' often refers to the idea that humans have a genetically determined upper limit to their adult stature and anthropometric dimensions (Bogin, 2006). Although this comes dangerously close to supporting the concept of genetic determinism, we are not using it in that manner here. means people are taller, weigh more, and have bigger muscles today than they did 100 years ago. Increases in tibial length could lead to a larger moment arm acting on the knee and on the tendons crossing it. Coupled with the increased force from a larger gastrocnemius, this could produce the mechanical stimuli necessary to initiate fabella formation and/or ossification. However, these factors do not explain the high prevalence of cartilaginous fabellae in foetuses, or why there was no relationship between presence and height in our sample.
Lastly, it is possible there is no shift in fabella prevalence rate, but the increase in prevalence rates is due to a change in fabella identification, where fabellae that were being previously ignored are now being identified. We believe this is highly unlikely for two reasons. First, there were no other changes in the prevalence of sesamoid bones in the hand or foot, and if there was a change in sesamoid bone identification protocol, it would likely not be isolated to the fabella. Secondly, the inclusion of X-ray and CT scans to determine prevalence rates in recent studies should lead to a decrease, not an increase, in prevalence rates through time, as cartilaginous fabellae, which may or may not have been included in previous studies, cannot be detected by Xrays and CT scans.
In this study, we investigated the prevalence rate of the fabella in a Korean population using published CT scans. Our prevalence rate of 52.83 and 44.34% for individuals and knees, respectively, falls within the range of those reported in the literature and shows an increase in fabella prevalence in Koreans over the past 80 years. In addition, we found bilateral fabellae to be more common than unilateral ones, there were no sex differences in prevalence Data taken from Table 2 in Yammine (2014). Although Yammine (2014) reported differences in prevalence due to sex and race, all data were pooled here, as there were only 16 studies stretching over 120 years. Prevalence rates were given per hand. In cases where ulnar and radial sesamoid bones were reported separately, the higher value was used, as it was not possible to determine whether the sesamoid bones were always from the same or different individuals. Z-value = test statistic. A Bonferroni-corrected P-value of 0.00833 (P = 0.05/6) shows a lack of any statistically significant trends. Data taken from Table 6 in Yammine (2015): data on the hallux (Table 2) were not analysed because they were highly mixed. Similar to the data with the sesamoid bones in the data, all data were pooled here, as there were only 16 studies stretching over 121 years. Prevalence rates were given per foot. In cases where tibial and ulnar sesamoid bones were reported separately, the higher value was used, as it was not possible to determine if the sesamoid bones were always from the same or different individuals. Z-value = test statistic. A Bonferroni-corrected Pvalue of 0.0125 (P = 0.05/4) shows a lack of any statistically significant trends.  (Table 5).  (2015) (n = 16 studies). Similar to the sesamoid bones in the hand, there was no correlation between foot sesamoid bone prevalence and time (Table 6).
rates, and presence of a fabella was uncorrelated with height and age. We also found a significant increase in fabella prevalence rates through time, but we are unsure why this has occurred and why there has not been an increase in other sesamoid bones in the human body during the same time span.