Article Text

Download PDFPDF

A population-based cross-sectional study of the association between facial morphology and cardiometabolic risk factors in adolescence
  1. Jelena Djordjevic1,
  2. Debbie A Lawlor2,3,
  3. Alexei I Zhurov1,
  4. Arshed M Toma1,
  5. Rebecca Playle1,
  6. Stephen Richmond1
  1. 1Department of Applied Clinical Research and Public Health, School of Dentistry, Cardiff University, Cardiff, UK
  2. 2MRC Centre for Causal Analyses in Translational Epidemiology, University of Bristol, Bristol, UK
  3. 3School of Social and Community Medicine, University of Bristol, Bristol, UK
  1. Correspondence to Dr Jelena Djordjevic; alexmjz{at}


Objective To determine whether facial morphology is associated with fasting insulin, glucose and lipids independent of body mass index (BMI) in adolescents.

Design Population-based cross-sectional study.

Setting Avon Longitudinal Study of Parents and Children (ALSPAC), South West of England.

Participants From the ALSPAC database of 4747 three-dimensional facial laser scans, collected during a follow-up clinic at the age of 15, 2348 white British adolescents (1127 males and 1221 females) were selected on the basis of complete data on cardiometabolic parameters, BMI and Tanner's pubertal stage.

Main outcome measures Fasting insulin, glucose and lipids (triglycerides, high-density lipoprotein cholesterol (HDLc) and low-density lipoprotein cholesterol (LDLc)).

Results On the basis of the collection of 63 x, y and z coordinates of 21 anthropometric landmarks, 14 facial principal components (PCs) were identified. These components explained 82% of the variation in facial morphology and were used as exposure variables. With adjustment for age, gender and pubertal stage, seven PCs were associated with fasting insulin, none with glucose, three with triglycerides, three with HDLc and four with LDLc. After additional adjustment for BMI, four PCs remained associated with fasting insulin, one with triglycerides and two with LDLc. None of these associations withstood adjustment for multiple comparisons.

Conclusions These initial hypotheses generating analyses provide no evidence that facial morphology is importantly related to cardiometabolic outcomes. Further examination might be warranted. Facial morphology assessment may have value in identifying other areas of disease risk.

  • Epidemiology
  • Preventive Medicine
  • Statistics & Research Methods

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: and

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Article summary

Article focus

  • Three-dimensional imaging opens up a new chapter in investigations of facial morphology. Previous research revealed associations of facial morphology with obesity in adolescents, but whether facial morphology can be used to identify those at future risk of adverse cardiometabolic outcomes is unknown.

Key messages

  • Our results suggest that facial morphology is not strongly or consistently associated with fasting insulin, glucose or lipids, particularly after adjustment for body mass index, in white British adolescents. Facial morphology is therefore unlikely to be useful in identifying white British adolescents at future risk of adverse cardiometabolic outcomes.

  • A suggested methodology can be used in future studies to explore the associations between facial parameters and other health outcomes. It might provide valuable insights into how facial morphology can be indicative of health.

Strengths and limitations of this study

  • The strengths of this study are a large sample size and the homogeneity of the sample: all participants were of white origin, born and brought up in the same region of the UK. A non-invasive, accurate and reliable method was used for capturing details of facial soft tissue morphology. A comprehensive statistical analysis was undertaken to extract principal components of facial morphology.

  • The study has some limitations. First of all, the study is ethnic-specific. Second, a face could not be easily represented as a single exposure due to the complexity of its morphology and the vast amount of data captured by the laser scanning system. Therefore, some data reduction was necessary prior to the analysis. Furthermore, it was not possible to control all the confounding factors in a cross-sectional study design. Since faces of adolescents are still developing, changing their shape and size, future studies might have to investigate the relationship between these changes and cardiometabolic characteristics over time.


Recent technological advancements in imaging methods marked a transition from a two-dimensional to three-dimensional (3D) approach in craniofacial research, thus opening a new era. A special emphasis has been placed on the development and application of non-invasive methods to capture the human face accurately and reliably.1 ,2 Among these, laser surface scanning and stereophotogrammetry have gained wide acceptance of the research community.3 So far, a large spectrum of medical disciplines have utilised these methods in the investigations of facial growth, facial dysmorphology, craniofacial identification, as well as the influence of different medical conditions on facial phenotype.4–12 Therefore, an exciting opportunity has occurred to explore whether facial characteristics can serve as new diagnostic measures of illnesses.

Childhood obesity is becoming an epidemic health problem.13 It is evident from many studies that it is associated with an increased risk of type 2 diabetes and cardiovascular disease later in life.14–16 Despite this fact, the relation between obesity and craniofacial development has been rarely investigated. Bimaxillary prognathism (overdeveloped jaws in the sagittal direction) and increased transverse facial dimensions seem to indicate the difference between obese adolescents and their normal-weighted peers.17–19 However, the association between metabolic phenotype and facial form has not been addressed previously.

In order to investigate this problem, a large sample and a comprehensive 3D approach to facial measurements are needed. In this cross-sectional study, which can be considered to be hypothesis generating, we examined the associations of facial soft tissue morphology with metabolic phenotype (fasting insulin, glucose, triglycerides, high-density lipoprotein cholesterol (HDLc) and low-density lipoprotein cholesterol (LDLc)) in a large general population cohort of adolescents using an existing database of 3D facial laser scans.

Material and methods


We used the data from the Avon Longitudinal Study of Parents and Children (ALSPAC), a UK-based longitudinal birth cohort study designed to explore genetic and environmental influences on health and well-being.20 ,21 All pregnant women were eligible to participate in ALSPAC if their estimated delivery date fell between 1 April 1991 and 31 December 1992, inclusive. In total, 14 541 pregnant women were recruited, and from these women there were 14 676 live-born infants. Since age 7, surviving offspring have been invited to regular follow-up clinics.

The current study was approved by the ALSPAC Law and Ethics Committee and the Local Research Ethics Committee and informed consent was obtained from children and their parents or guardians. The data collected during an annual follow-up clinic at the age of 15, which was attended by 5235 adolescents, were examined. On that occasion, facial laser scanning was performed, and after a dropout of 488 individuals due to the low quality of the scans or some sort of facial dysmorphology, a database of 4747 individuals (2233 males and 2514 females) was formed.22 Out of these, we selected 2348 white adolescents (1127 males and 1221 females), with complete data related to the outcome and confounding variables (see below), as facial laser scans were used to derive exposure variables. The flow chart diagram (figure 1) shows a gradual selection of individuals who comprised the final sample. In order to make sure that there was no selection bias, we first compared facial principal components (PCs; ie, exposure variables; see below) of the study sample (2348 adolescents) with those of 4747 adolescents forming a 3D facial database22 and concluded that there was no reason to believe that selected individuals were significantly different in terms of facial morphology. Second, we compared observed values for outcomes and confounding variables in the study sample (2348 adolescents) with imputed variables in the eligible sample of the follow-up clinic (5235 adolescents), which were published as supplementary online material of the previous study.23 The distributions in imputed datasets were very similar to those observed, providing some evidence that the missing data were missing at random.

Figure 1

Flow chart showing the selection of a study sample from a 15+ year follow-up clinic of the Avon Longitudinal Study of Parents and Children. All analyses presented in this paper are based on 2348 participants with complete data on facial soft tissue morphology (exposure), blood-based indicators of insulin resistance and associated cardiometabolic risk factors (outcomes) and body mass index and pubertal stage (covariables).


Exposure variables

Facial laser scans were used to derive PCs of facial morphology, which served as the exposure variables. This is explained in detail in the section on statistical analysis. Prior to this, it was necessary to perform three steps, which will be described here. First of all, facial scans were processed. The validity and reliability of the laser scanning procedure, as well as the processing stages of the scans, have been previously investigated.24–27 Second, 21 anthropometric landmarks were manually identified on facial scans by one experienced examiner (figure 2), according to their respective definitions by Farkas,28 and their x, y and z coordinates were saved for the subsequent analysis. Previous research showed that these landmarks are clinically reliable.29 ,30 Finally, facial scans were initially normalised according to the natural head position, with the origin of the coordinate system set at the point halfway between the inner corners of the eyes (mid-endocanthion). The x-axis was pointing left, from the right eye to the left eye, the y-axis was pointing vertically upwards from the chin to the forehead and the z-axis was pointing outwards, in the nose direction. The coronal, sagittal and transverse planes were taken as the xy, yz and xz planes, respectively.1 ,2 ,8 ,22 ,29

Figure 2

Twenty-one anthropometric landmarks which were identified on facial laser scans of participants. (1) Glabella (g); (2) Nasion (n); (3) Endocanthion left (enl); (4) Endocanthion right (enr); (5) Exocanthion left (exl); (6) Exocanthion right (exr); (7) Palpebrale superius left (psl); (8) Palpebrale superius right (psr); (9) Palpebrale inferius left (pil); (10) Palpebrale inferius right (pir); (11) Pronasale (prn); (12) Subnasale (sn); (13) Alare left (all); (14) Alare right (alr); (15) Labiale superius (ls); (16) Crista philtri left (cphl); (17) Crista philtri right (cphr); (18) Labiale inferius (li); (19) Cheilion left (chl); (20) Cheilion right (chr); (21) Pogonion (pg). Definitions by Farkas28 were used. Reprinted from the author's previous publication with permission from ‘John Wiley and Sons’.

Outcome variables

Fasting insulin, fasting glucose, triglycerides, HDLc and LDLc were taken as the outcome variables. Full details of their assessment have been previously reported.23

Confounding variables

Since this study is exploratory (being the first to examine these associations) and our main motivation was to understand whether facial morphology might be able to predict those at risk of cardiometabolic disease over and above simple measurement of adiposity, we did not adjust for a wide range of confounding variables. However, we adjusted for age, pubertal stage and body mass index (BMI), as these are potentially important predictors of cardiometabolic risk and we would want to be clear that facial morphology predicted outcome over and above these. The age of the participants was recorded in months as they arrived at the clinic. Pubertal status was assessed on participants’ self report with Tanner’s questionnaires.

Statistical analyses

Participant characteristics were summarised with means (SD) for continuous, approximately normally distributed variables, median (IQR) for continuous right skewed variables, and the number (%) of categorical variables.

Generalised Procrustes Analysis (GPA) was performed on landmark configurations (each consisting of 63 x, y and z coordinates of 21 facial landmarks) in order to remove differences in landmark positions attributable to translation and rotation.31–33 Scaling was not performed in order to preserve facial size. Principal Component Analysis (PCA) was used to reduce the set of 63 coordinates into a smaller number of independent components of facial morphology. According to the ‘Kaiser-Guttman criterion’, PCs with eigenvalues greater than the average eigenvalue value were retained34–36 and saved as new exposure variables. The rotation method used for PCA was varimax with Kaiser normalisation.37 GPA was performed in the open source software R project and PCA in SPSS V.17.0 (SPSS Inc, Chicago, Illinois, USA).

While this is a cross-sectional study, in all our analyses we examined the association of PCs of facial morphology (as exposures) with fasting insulin, glucose, triglycerides, HDLc and LDLc (as outcomes) using multivariable linear regression models. No evidence was found for any gender interactions (all p values ≥0.1), and therefore analyses are presented with both genders combined. In the first model, we adjusted for age, gender and pubertal stage. In the second model, we adjusted for age, gender, pubertal stage and BMI and examined how much this reduced any associations of facial PCs with the outcomes. Fasting insulin and triglyceride levels were right (positively) skewed and their logged values were used in the linear regression models, which ensured that the model residuals were approximately normally distributed. The resultant regression coefficients with 95%CI are presented.

In these multivariable analyses, 140 comparisons were made (14 exposures with 5 outcomes and 2 models). In initial analyses, we considered the conventional 0.05 level of statistical significance. We then adjusted for multiple comparisons using a Bonferroni correction by dividing 0.05 by 140; thus, for these corrected analyses, a p value of 0.0004 would be considered to be statistically significant at the 0.05 level. All statistical analyses were performed in SPSS V.17.0 (SPSS Inc).


Table 1 shows the characteristics of the study participants. The PCA identified 14 PCs of facial morphology (table 2). Each PC consisted of a number of coordinates of anthropometric landmarks. For example, the first PC (PC 1) comprised 17 y coordinates of landmarks located in the upper and lower thirds of the face. These coordinates represented facial height (size). In order to facilitate the understanding and interpretation of individual PCs, they are presented graphically in figure 3. The first three PCs (facial size, inter-eye distance and prominence of the nose and lower lip) accounted for almost half of the total variation (45.7%). The other 11 PCs contributed to facial variation to a much lesser extent (between 1.6 and 5%), but marked those subtle features which make the faces unique.

Table 1

Characteristics of the study sample

Table 2

The results of the principal component analysis showing partial correlation coefficients between coordinates of anthropometric landmarks and facial principal components

Figure 3

Facial principal components (PCs). Numbers indicate percentages of normal facial variation explained by the given principal component. Coordinates which constitute each principal component are marked on the face (refer to table 2), and arrows indicate the x, y and z directions.

The multivariable associations of the 14 PCs with cardiometabolic outcomes are shown in tables 37. With adjustment for age, gender and pubertal stage (model 1), seven PCs were associated with fasting insulin, none with fasting glucose, three PCs with triglycerides and HDLc and four PCs with LDLc. After additional adjustment for BMI (model 2), four PCs remained associated with fasting insulin, none with glucose, one with triglycerides, none with HDLc and two with LDLc. However, none of these associations withstood adjustment for multiple comparisons.

Table 3

Multivariable association of 14 facial principal components (exposures) with fasting insulin as an outcome

Table 4

Multivariable association of 14 facial principal components (exposures) with fasting glucose as an outcome

Table 5

Multivariable association of 14 facial principal components (exposures) with triglycerides as an outcome

Table 6

Multivariable association of 14 facial principal components (exposures) with high-density lipoprotein cholesterol as an outcome

Table 7

Multivariable association of 14 facial principal components (exposures) with low-density lipoprotein cholesterol as an outcome


Laser surface scanning is a non-invasive technology which enables the accurate and precise analysis of facial morphology.1 ,2 ,24–26 Owing to its portability, easy application and relatively low cost, this technique is very suitable for epidemiological field studies. The vast amount of data captured by the system (more than 40 000 points, each consisting of x, y and z coordinates) is a testimony to the complexity of facial surfaces. For this reason, the face cannot be easily represented as a single exposure.

Therefore, it was necessary to make some facial data-reduction prior to its meaningful use. First of all, GPA (a widely established method in statistical shape analysis) was used to place landmark coordinates in the same space reducing confounding errors (rotation and translation). Second, PCA was applied on the set of coordinates and 14 facial PCs were identified, which accounted for almost 82% of the total variation in normal facial form, consisting of size and shape. Normal facial variation was recently analysed on a complete sample of 4747 faces from the ALSPAC database and the same number of PCs was extracted, with an almost identical order of individual PCs and very similar percentages of variation.22

The application of this statistical technique is not new. Previously, PCA was performed on two-dimensional data sets, obtained from either lateral skull radiographs or photographs of both children and adults.38–41 The resultant number of PCs in these studies was between 6 and 8, and these explained up to 90% of the total variance in facial profile, based on linear measurements between anthropometric landmarks, or their coordinates. However, with the introduction of sophisticated 3D imaging techniques, the amount of data entering PCA increased significantly. Therefore, the number of PCs which represent facial variation also increased: between 14 and 16 PCs have been reported to account for between 86% and 92% of the total variation.10 ,11 ,41 ,42

Although the first three components in the current study explain almost half of the total variation, other components are also important, since they represent subtle changes that make the face unique. Therefore, a decision was made to keep all of them in the subsequent multivariable analyses. Following adjustment for BMI and taking account of multiple statistical testing, we did not find that any of these PCs were associated with fasting insulin or associated cardiometabolic risk factors, suggesting that facial morphology is unlikely to be a reliable way of predicting young people at future risk of type 2 diabetes or cardiovascular disease. Consistent with other large epidemiological studies conducted in healthy general population samples, we were not able to directly measure insulin resistance using the gold standard euglycaemic hyperinsulinaemic clamp. Fasting insulin has been shown to have modest-to-strong correlations with clamp-assessed insulin resistance (correlation coefficients 0.5–0.9) in children and adolescents.43 ,44 Any measurement error is likely to be non-differential and therefore would be expected to bias results towards the null. Since strong associations of these outcomes with BMI have been shown in ALSPAC,23 any associations with a better measure of insulin resistance are unlikely to be stronger than those of BMI.

The study has some limitations. First of all, it is ethnic-specific, and therefore future studies will have to address the research question in different ethnic groups. Second, facial variation can be affected by many different factors. While it is possible to control the age, gender and ethnicity of the sample, environmental factors present a greater challenge, even with a good research strategy, as many of them can be unknown at the time of the study. The face changes throughout life, increasing in size and changing in shape.1 ,2 ,7 This holds true for the present sample consisting of 15-year-old adolescents. The cross-sectional design of the study did not allow us to track these changes and analyse their relationship with metabolic phenotype through time. That may be more important than the assessment of variation among individuals and thus should be considered in future studies.


Our results do not provide strong evidence that facial morphology is robustly and importantly associated with cardiometabolic risk factors. The associations identified were not consistent across outcomes, were weak in magnitude, attenuated with adjustment for BMI, and did not withstand correction for multiple statistical testing. Further study of facial parameters with cardiometabolic and/or other health outcomes might provide valuable insights into how facial morphology can be indicative of health.


We are extremely grateful to all the mothers, their parents and the children who were taking part and to the midwives for their cooperation and help in recruitment. The whole ALSPAC study team comprises interviewers, computer technicians, laboratory technicians, clerical workers, research scientists, volunteers and managers who continue to make the study possible. This study could not have been undertaken without the financial support of the Medical Research Council, the Wellcome Trust, the UK Department of Health, the Department of the Environment, and DfEE, the National Institutes of Health, a variety of medical research charities and commercial companies. The ALSPAC study is part of the WHO initiated European Longitudinal Study of Pregnancy & Childhood.



  • Contributors JD, DAL and SR were responsible for the conception and the design of the study. JD and SR initiated the study. JD and AMT collected the data. JD, DAL and RP were responsible for statistical analyses. JD, AMT, AIZ and SR analysed and interpreted the data on facial parameters. JD, DAL and RP analysed and interpreted the data on metabolic parameters. JD and DAL wrote the first draft of the paper. All authors contributed to and approved the final version of the paper.

  • Funding The research leading to these results has received funding from the Cardiff University for the assessment of facial morphology and from USA National Institute of Health (NIH): National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) grant number (R01 DK077659) for the measurement of fasting glucose and insulin. The UK Medical Research Council and Wellcome Trust (092731), together with the University of Bristol, provide core support for the ALSPAC study. DAL works in a centre that receives infrastructure funding from the UK Medical Research Council (G0600705). None of the funders influenced the data analyses or the interpretation of results. The comments made in this paper are those of the authors and not necessarily those of any funders.

  • Competing interests None.

  • Ethics approval The study was approved by the ALSPAC Law and Ethics Committee and the Local Research Ethics Committee.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.