Next Article in Journal
On a Method For Reconstructing Computed Tomography Datasets from an Unstable Source
Next Article in Special Issue
Comparative Study of Contact Repulsion in Control and Mutant Macrophages Using a Novel Interaction Detection
Previous Article in Journal
Subpixel Localization of Isolated Edges and Streaks in Digital Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multilevel Analysis of the Influence of Maternal Smoking and Alcohol Consumption on the Facial Shape of English Adolescents †

School of Dentistry, Cardiff University, Heath Park, Cardiff CF14 4XY, UK
*
Author to whom correspondence should be addressed.
This article is an extended version of our paper published in Farnell, D.J.J.; Galloway, J.; Zhurov, A.I.; Richmond, S. Multilevel models of age-related changes in facial shape in adolescents. In Proceedings of the Medical Image Understanding and Analysis, MIUA 2019, Liverpool, UK, 24–26 July 2019; Zheng, Y., Williams, B., Chen, K., Eds.; Springer: Cham, Switzerland, 2019; Volume 1065, pp. 101–113.
Submission received: 26 March 2020 / Revised: 10 May 2020 / Accepted: 15 May 2020 / Published: 18 May 2020
(This article belongs to the Special Issue MIUA2019)

Abstract

:
This cross-sectional study aims to assess the influence of maternal smoking and alcohol consumption during pregnancy on the facial shape of non-syndromic English adolescents and demonstrate the potential benefits of using multilevel principal component analysis (mPCA). A cohort of 3755 non-syndromic 15-year-olds from the Avon Longitudinal Study of Parents and Children (ALSPAC), England, were included. Maternal smoking and alcohol consumption during the 1st and 2nd trimesters of pregnancy were determined via questionnaire at 18 weeks gestation. 21 facial landmarks, used as a proxy for the main facial features, were manually plotted onto 3D facial scans of the participants. The effect of maternal smoking and maternal alcohol consumption (average 1–2 glasses per week) was minimal, with 0.66% and 0.48% of the variation in the 21 landmarks of non-syndromic offspring explained, respectively. This study provides a further example of mPCA being used effectively as a descriptive analysis in facial shape research. This is the first example of mPCA being extended to four levels to assess the influence of environmental factors. Further work on the influence of high/low levels of smoking and alcohol and providing inferential evidence is required.

1. Introduction

In facial shape research, we are faced with a large number of outcome variables in the form of facial landmarks. Principal component analysis (PCA) is a well-established technique for dimension reduction and can be used to help make analyses more manageable. However, conventional PCA calculates the principal components (PCs) and their corresponding eigenvalues using the dataset as whole, which ignores any groups or clusters that occur naturally in the population of subjects or shapes [1]. Such clusters often occur due to factors/covariates such as sex, age, or ethnicity. Furthermore, we cannot prevent the effects of different factors/covariates becoming “mixed together” in conventional PCA because PCs are orthogonal to each other. Further, the effects of some covariates on facial shape might be small and conventional PCA therefore limits our ability to assess subtle differences in facial shape because this might be “lost in the tail” of the eigenvalues. Halazonetis [2] points out yet another potential limitation of conventional PCA, stating: “A significant disadvantage of PCA is that the resulting components, because they are derived statistically, do not necessarily have a clear biological interpretation. Usually, only the first, or the first few, can be described satisfactorily...” (p. 580).
We believe taking a multilevel approach to PCA is an elegant way to overcome these issues and can be seen as a non-trivial extension to PCA. Multilevel PCA (mPCA) addresses these problems directly by representing (and so isolating to some extent) specific effects at separate levels of the model. Furthermore, although PCs are still orthogonal within each level using mPCA, they do not have to be orthogonal between levels. This process should, in principle, reduce the amount of mixing due to orthogonality in conventional PCA. Therefore, mPCA allows dimension reduction to be carried out in a supervised manner, thus increasing our confidence that the information we retain within the level-specific PCs describes our covariates of interest. mPCA has been described by Lecron et al. [3] and used in the segmentation of spine radiographs. Subsequently, mPCA has been used to assess dental radiographs [4] and has been successful in assessing the influence of various covariates on facial shape, including ethnicity, sex, smiling and age [5,6,7,8,9]. These papers document the progression of the mPCA technique and provide the mathematical background.
mPCA gives a number of benefits. Firstly, we are able to calculate the percentage of the total variation due to each covariate. This provides useful information and allows us to determine the importance of certain covariates in facial shape research. Secondly, mPCA allows us to be more confident that differences between groups are due to our covariate of interest and allows the separation of the influence of these from variation not explained by the other levels of the model by utilising a subject level. Finally, mPCA allows us to assess covariates that may have a subtle influence on facial shape, as the variation due to these covariates is less likely to be “hidden” during the dimension reduction process. An explanation of the relationship between mPCA [10,11,12] and multilevel methods implemented, e.g., using structural equation models (SEMs), is given by Timmerman [10]: “traditionally, SEM aims at inference, whereas CA is used as a descriptive technique”. Indeed, we use mPCA here as a purely descriptive approach for our data set.
Maternal smoking and alcohol consumption are two covariates that may have a subtle influence on the facial shape of non-syndromic adolescents. Maternal smoking and alcohol consumption are potential teratogens (agents that are capable of disrupting normal development during pregnancy). Maternal smoking may influence facial shape by inducing hypoxic states in the foetus [13] or through epigenetics [14,15]. Maternal smoking has been linked to craniosynostosis [16,17], head circumference [18,19,20], facial asymmetry due to hemifacial microsomia [21] and cleft lip and palate [16,22]. Maternal alcohol consumption during pregnancy has been linked to foetal alcohol syndrome (FAS). Common facial features associated with FAS include short palpebral fissures, reduced upper lip height, maxillary hypoplasia and a short nose [23], alongside cognitive difficulties [24]. Historically, FAS has been attributed to high levels of alcohol consumption during pregnancy. However, it could be hypothesised that the features affected in FAS may be influenced to a lesser extent when lower alcohol levels are consumed by mothers during pregnancy.
It is clearly advantageous that expectant mothers reduce the exposure of their child to teratogens. Unfortunately, despite these risks, expectant mothers still smoke and consume alcohol during pregnancy. With the advent of social media, facial appearance is becoming an increasing priority to the general public. Therefore, demonstrating the effects of smoking and alcohol on normal facial development, to supplement the risk of syndrome development, may act as a further deterrent for expectant mothers.
Previous studies have suggested that even low levels of alcohol and smoking during pregnancy may have an influence on the facial shape of one-year-old children [25]. However, others suggest that maternal alcohol consumption has a minimal effect on the facial shape of non-syndromic children at 15 years old [26]. Given that maternal smoking and alcohol consumption may have a subtle influence on the facial shape of adolescents, they provide excellent examples for demonstrating where mPCA can be useful.
The aim of this study is therefore to investigate the influence of maternal smoking and alcohol consumption on the facial shape of their offspring at 15 years old. We demonstrate the advantages of using mPCA by using an extended mPCA model of four levels to assess the influence of (1) maternal smoking, (2) maternal alcohol consumption, (3) biological sex and (4) subject on the facial shape of non-syndromic English adolescents who presented without any obvious craniofacial defects. We also aim to demonstrate the advantages of mPCA, in particular the descriptive visualisations possible with mPCA. We provide the equivalent visualisations for conventional PCA as a comparison. It is hoped readers could extrapolate the potential advantages to their field of study.
Maternal smoking and alcohol consumption were categorised according to exposure of the foetus during the 1st and/or 2nd trimesters of pregnancy. The effect of maternal smoking and maternal alcohol consumption (average 1–2 glasses per week) was minimal, with 0.66% and 0.48% of the variation in the 21 landmarks of non-syndromic offspring explained, respectively.

2. Materials and Methods

This is a cross-sectional study, utilising the data available from the Avon Longitudinal Study of Parents and Children (ALSPAC). The ALSPAC Ethics and Law Committee, and the Local Research Ethics Committees, provided approval for the study (B3166). Informed consent for the use of data collected via questionnaires and clinics was obtained from participants following the recommendations of the ALSPAC Ethics and Law Committee at the time.
The ALSPAC cohort includes the children and partners of women who were expected to give birth between 1 April 1991 and 31 December 1992. Further cohorts of children, who were missed from the original subject recruitment phase were subsequently recruited, with a total of 15,589 foetuses included. Of these, 14,901 were alive at 1 years old, with 3D facial scans available for 4747 children at 15 years old. The mean age of the subjects was 15.43 (+/− 0.28 SD) years old. Maternal smoking during the 1st and 2nd trimesters was determined via questionnaire at 18 weeks gestation. Maternal alcohol consumption during the 1st trimester and when the baby was first felt to have moved (as a proxy for the 2nd trimester) was also determined via questionnaire. These were added to the statistical model as separate levels. Biological sex was added as a confounder as a separate level in the model. The ALSPAC website contains details of all the data that is available through a fully searchable data dictionary and variable search tool (http://www.bristol.ac.uk/alspac/researchers/our-data/). Further details on the ALSPAC can be found at Boyd et al. [27] and Fraser et al. [28].
The 3D laser facial scans were collected in natural head position using Konica Minolta VI-900 laser scanners. Those with craniofacial abnormalities and twins/triplets were excluded from this study. The methodology for the laser scans can be found in Kau et al. [29]. To facilitate landmarking, the faces were standardised by aligning on a vertical cylinder and using mid-endocanthion as the origin [30]. The landmark configurations were all superimposed using generalised Procrustes analysis. The scaling of the faces attempts to isolate differences in shape. Further information on the image processing procedure can be found at Toma et al. [31].
The outcome measure used was differences in 21 facial landmarks (used as a proxy for the main facial features), as described by Farkas [32]. These were plotted manually onto 3D facial scans of the subjects, using Rapidform 2006 (Figure 1). The researcher plotting the landmarks underwent a calibration exercise and was blinded to the smoking and alcohol consumption of the subjects’ mothers during pregnancy. The facial images were plotted in grey scale to remove any influence of skin tone or makeup. The landmarks chosen are accurate and reproducible [33].
mPCA was used for the analysis. This method requires categorisation of the data. The sample sizes per grouping are presented in Table 1 and the flow of inclusion of subjects in Figure 2. A sample size calculation was not possible due to the novelty of the analysis. The mPCA technique does not currently manage missing data. Analyses were restricted in the following ways:
  • It was not possible to assess the influence of smoking and alcohol levels throughout pregnancy (e.g., low/high) due to sample sizes in each grouping, and this is the subject of further work.
  • Subjects who smoked or consumed alcohol in the 2nd trimester only were excluded due to reduced sample sizes in some of the groupings.
mPCA was used as a supervised dimension reduction technique and was carried out in MatLab 2017b. Maternal smoking, maternal alcohol consumption and biological sex of subject were explicitly modelled here. An attempt is made to control for “variation not explained by levels 1, 2 or 3” by adding a “subject” level to the model. For each of the following four levels, the eigenvalues and eigenvectors were found by determining covariance matrices for each level individually:
  • Level one: Maternal smoking during pregnancy;
  • Level two: Maternal alcohol consumption during pregnancy;
  • Level three: Biological sex of the subject;
  • Level four: Variations that are consistent across levels 1 to 3, i.e., the variation not explained by levels 1 to 3.
The number of eigenvalues retained in the mPCA analysis was determined by visualisation of the eigenvalue plot (Figure 3). The maximum number of eigenvalues that can be retained at each level is the number of groups minus 1, which is because the ranks of covariance matrices are constrained by the number of groups or subjects used in obtaining them. As sex only has two possible groupings (male/female), only one eigenvalue was retained at this level.
The mPCA model is fitted to the new set of landmark points by using all four levels in the model at the same time. This is written as: new shape = grand mean shape (over all subjects, all smoking groups, all alcohol groups and both sexes) + weighted modes/components due to level 1 (maternal smoking) + weighted modes/components due to level 2 (maternal alcohol) + weighted modes/components due to level 3 (biological sex of subject) + weighted modes/components due to level 4 (subjective variations that are consistent across all groups). These weights are also called component "scores". The model fit was found by using a global optimisation procedure of an overall point-to-point cost function using MatLab.
Farnell et al. [5,6,7,8,9] describe the mathematical concepts and various applications of mPCA, including nested, non-nested and mixed approaches. A non-nested approach was used here, with median averaging of the covariance matrices used to attempt to deal with outliers. The terms “nested” and “non-nested” refer to the population of subjects, shapes or images being studied, and the analysis follows the structure of the data. Fully “nested” cases are those where shapes, subjects, or groups belong to exactly and only one group in the level above it and (importantly) at all levels. A classic example here are clusters by school or class for some arbitrary “outcome” (e.g., exam results or here, facial shape). Each child belongs to one class only and each class belongs to only one school. Thus, this design is “fully nested” at all levels. Fully “non-nested” cases have groups at a given level that can belong to more than one group in the levels above it. Our dataset is a good example of a non-nested case. We represent maternal smoking, maternal alcohol consumption and biological sex at different levels (e.g., 1, 2 and 3) of the model and subjects at the bottom level (4). Note that one can have males and females in any of the smoking or alcohol groups. Note that there is no actual order as such for the smoking, alcohol and biological sex levels (1, 2 and 3) in this example, although we might impose an order on the data artificially. In practice, all that this would mean for our analysis using mPCA is that covariance matrices are found in a slightly different manner. Finally, some multilevel cases contain both fully nested and non-nested elements, which one might refer to as “mixed”. “Mixed” approaches are appropriate, for example, when assessing multiple images from the same subject. The use of mPCA in facial shape research can be seen as a complimentary descriptive technique to traditional statistical methods in this field.
Estimates of the (standard) errors associated with the eigenvalues were found using Monte Carlo simulation (number of MC replicates = 1000). A simple Monte Carlo method was used here only in order to estimate the (standard) errors of eigenvalues. The mean shape vector (mu) and the covariance matrix (sigma) were found for the original data set of sample size n. Note that the x-, y-, and z-components of the landmark points in the original data set were found to be normally distributed. The MATLAB command “mvnrnd(mu,sigma,n)” was used to create Monte Carlo replicates of the data set, assuming an underlying (correlated) multivariate normal distribution. The 1000 MC replicates of the data were thus created, and eigenvalues were determined for each replicate. In this manner, standard deviations for the eigenvalues over the 1000 MC replicates provided broad estimates of standard errors. A limitation of this simple approach is that some of the complexity of the original data might be lost; future studies will consider more complicated MC schemes that retain this complexity.
In summary, the following analysis was conducted:
  • Non-nested, 4-level analysis;
  • Levels: 1) smoking, 2) alcohol, 3) sex, and 4) subject;
  • Eigenvalues retained: smoking: 2, alcohol: 2, sex: 1, and subject: 35.
Conventional PCA was performed using the PCA function in MatLab 2017b in order to generate component scores for a comparison with mPCA. The sign attributed to the scores is reversed in comparison to those generated using the inhouse algorithm associated with mPCA. However, the patterns visualised in the data are the same.
The total percentage of variation explained by each mPCA level was calculated by dividing the retained eigenvalue magnitude for the respective levels by the total eigenvalue magnitude. The component scores for conventional PCA and each level of mPCA were standardised by dividing by the square root of the respective eigenvalue. Where appropriate, the differences in each of the 21 landmarks were determined by finding the difference between the landmarks of the average face plus 2*square root of the eigenvalue magnitude and average face minus 2*square root of the eigenvalue magnitude.

3. Results

3.1. Demographics

In total, 1775 males and 1980 females were included in this study. Smoking and alcohol consumption levels of their mothers during pregnancy are presented in Table 1, with more detailed information available as supplementary material (Supplementary Tables S2–S5). Smoking during pregnancy was rarer than alcohol consumption. In total, 593 (15.8%) mothers smoked during the 1st trimester, with 159 of the mothers included in the study stopping smoking before the 2nd trimester. The smoking levels of mothers in the 1st and 2nd trimester varied, with a mean ranging from 5 to 11 cigarettes per day in each group. In total, 2378 mothers (63.3%) consumed alcohol within the 1st trimester of pregnancy and 1798 mothers (47.9%) during the 1st and 2nd trimesters. Alcohol consumption during pregnancy was therefore relatively common for the ALSPAC mothers in 1990. The alcohol consumption levels of mothers in the 1st and 2nd trimesters varied from <1 glass to 2 glasses per week on average. Very high levels of alcohol consumption were rare, with five mothers reporting drinking 3–9 glasses/day and two mothers reporting 10+ glasses/day (Supplementary Tables S4 and S5). The mothers were asked to take “glass” to mean “a pub measure of spirits, half a pint of lager/cider, a glass of wine etc.”.

3.2. Maternal Smoking Level

One of the main benefits of mPCA is that we can easily calculate the percentage of the total variation due to maternal smoking. Analysis of the eigenvalues reveals that the maternal smoking level had a minimal influence (0.66%) on the total variation in the 21 facial landmarks (Table 2, Figure 3). Estimates of the (standard) errors associated with the eigenvalues found using Monte Carlo simulation are shown in Supplementary Table S1. Given that a small amount of variation is attributed to maternal smoking, it follows that the variance associated with this covariate is small. Therefore, in conventional PCA, the true influence of maternal smoking may be “hidden” in multiple PCs, particularly if a large amount of the variance associated with these PCs is explained by another covariate.
The standardised component scores for conventional PCA PC1 and PC2 do not show any obvious patterns in the separation of the group means (Figure 4a). Inspection of the component scores up to PC10 suggests that there may be some subtle separation of the group means due to maternal smoking at PC7/8 (Figure 4b, Supplementary Figure S1). However, as discussed previously, when using conventional PCA, it is difficult to assess this with certainty from a visual plot, as each PC cannot be solely attributed to one covariate. This is because the PCs are orthogonal and were generated in an unsupervised manner, not taking into account any of the covariates at the beginning of the analysis. The variation explained by PC7/8 may include the influence of maternal smoking but will be “mixed in” with variation due to other covariates.
Multilevel PCA allows us to be more certain about the differences we visualise, as the PCs at level one were generated taking maternal smoking into account. Indeed, we see a similar pattern in the separation of the group means at PC2 of the mPCA model as we see with conventional PCA PC7/8. The offspring of mothers who smoked in the 1st trimester only and smoked in both the 1st and 2nd trimesters are possibly differentiable from those whose mothers did not smoke during pregnancy (Figure 4c). The ease at which any potential patterns can be visualised demonstrates one of the strengths of mPCA. Eight PCs needed to be visualised when using conventional PCA to gain a suggestion of the influence of maternal smoking, whereas only two PCs needed to be visualised with mPCA. This is more efficient.
Our analyses suggest that the nasal bridge shape/length and orbital position/shape may be most influenced by the maternal smoking level. Differences were also seen in the prominence of the brow-ridge, lower lip and chin, and the position of the corners of the mouth (Table 3). These differences may suggest that smoking could influence shared facial genes. However, as there is an outlying group of females who smoked in the 1st trimester and consumed no alcohol (Figure 4c), only a suggestion of a pattern in the separation of the group means, differential sample sizes in the groups and a large range of component scores in each grouping (Table 4), any differences should be interpreted with caution at this stage.

3.3. Maternal Alcohol Consumption Level

The maternal alcohol consumption level, in this population, influenced the variation in the 21 facial landmarks by 0.48%, and thus had a minimal effect (Table 2, Figure 3). Estimates of the (standard) errors associated with the eigenvalues found using Monte Carlo simulation are shown in Supplementary Table S1. The scatter plots of the population centroids suggest that there is little pattern in the separation of the group means at this level in mPCA (Figure 4d) and when visualising the first 10 PCs of conventional PCA (Figure 4c, Supplementary Figure S1). This further confirms that alcohol consumption in this population does not appear to have had a long-term influence on the children’s main facial features. However, as the maternal alcohol consumption in this population was generally low at 1–2 glasses per week, and only non-syndromic children without obvious craniofacial abnormality were assessed here, this finding does not condone the consumption of alcohol during pregnancy.

3.4. Biological Sex Level

Multilevel PCA suggests that the biological sex level explains 8.88% of the variation in the 21 facial landmarks (Table 2, Figure 3). Estimates of the (standard) errors associated with the eigenvalues found using Monte Carlo simulation are shown in Supplementary Table S1. Conventional PCA shows that there is a clear pattern in the separation of the group means of the subjects due to biological sex when visualising the component scores at PC1 and PC2 (Figure 5a). However, as discussed previously, it is very likely that other covariates are contributing to this effect. The effects of the unknown covariates may enhance or reduce the separation that we see. Clear separation of the group means by biological sex is also visible at mPCA PC1 (Figure 5b). Only one PC is possible in mPCA due to the rank of the matrices.

3.5. Subject Variation (Variation not Explained by Levels 1–3)

Subject variation is not explicitly modelled in conventional PCA. Instead the effects of other covariates are “mixed in” with the effects of maternal smoking, maternal alcohol consumption and biological sex. Being able to separate these effects in mPCA is a distinct advantage. Multilevel PCA suggests that 85.85% of the variation in the 21 landmarks is not due to smoking, alcohol or biological sex (Table 2). Estimates of the (standard) errors associated with the eigenvalues found using Monte Carlo simulation are shown in Supplementary Table S1. It follows that in conventional PCA, the effect of other covariates (not smoking, alcohol or biological sex) on the PCs, and thus the component scores, is high. This can be visualised in the eigenvalue plot, where the eigenvalues for conventional PCA follow a similar pattern as the subject level eigenvalues in mPCA (Figure 3).
When visualising the standardised component scores from the subject level of mPCA, the population centroids by any of the groupings (i.e., maternal smoking, maternal alcohol level or biological sex) are generally centred at the origin at the subject level (Figure 5c), as anticipated. There are a few groups that are more distant to the origin, which may be explained by the small sample sizes in some of the groups.

4. Discussion

4.1. Smoking

The influence of maternal smoking during pregnancy on non-syndromic facial shape appears to be subtle (0.66% of the total variation in the 21 landmarks). The differences in the facial features in this population appear to be most evident when comparing subjects whose mothers did not smoke in pregnancy with those whose mothers smoked in the 1st trimester or the 1st and 2nd trimesters. A previous study on a one-year-old population suggested in their supplementary material that some minor changes, in the region of 0.4–0.5 mm, may be seen in the chin, lower lip prominence, forehead prominence and nasal bridge shape [25]. Our findings agree with most of these findings, with the addition of possible differences in the orbits. However, given that the differences explain less than 1% of the total variation in the 21 landmarks, our overall findings also agree with Koziel et al. [34], who suggested that, after 10 years old, the influence of maternal smoking may be minimal, particularly in females. It should also be noted that our findings are based on individuals who have no history of craniofacial abnormality, do not separate high/low levels of exposure and assessed 21 facial landmarks only. This study primarily aims to provide an example of mPCA being used for the exploration of environment factors on facial shape at this stage and thus a more detailed analysis is the subject of further work.

4.2. Alcohol

There have been conflicting results on the impact of low levels of maternal alcohol consumption on facial shape in the literature. Muggli et al. [25] suggested that alcohol during the 1st trimester could have an impact on facial shape at 1 years old. They suggest that low levels of alcohol consumption could influence forehead shape, moderate to high alcohol consumption could affect orbital shape, midface and chin prominence, and binge levels of alcohol could affect the chin. These differences appear to be in the range of 1 mm. Conversely, Howe et al. [26], who also looked at the influence of low alcohol consumption on the ALSPAC cohort, found that alcohol had minimal effect on facial shape at 15 years old, using alternative methodology to mPCA. It could therefore be hypothesised that any change in facial shape due to maternal alcohol (without a FAS diagnosis) is overshadowed by the influence of other factors as the child develops into adolescence. However, these differences may also be explained by the fact that different populations and different methodologies were used.

4.3. Strengths and Limitations of mPCA

This study is a starting point, using mPCA methodology in the field of facial shape research. mPCA methodology is a flexible dimension reduction technique. Separating the covariates into levels allows us to be more confident that subtle variation due to smoking and alcohol exposure during pregnancy was retained and increases confidence that the differences visualised between the groups using scatter plots were due to the covariate assessed. The addition of a “subject” level has also attempted to isolate the effect of variation not explained by smoking, alcohol or sex. As many of the covariates in facial shape research are likely to be unknown, studies which have utilised alternative methodologies may not correct for all of the necessary covariates. The inclusion of a “subject” level is therefore a step towards modelling this unknown variation.
As mPCA requires categorisation of the covariates, it can be difficult to assess naturally continuous variables. Furthermore, as the number of eigenvalues retained is limited by the rank of the matrices, the number of eigenvalues available for assessment is reduced. The addition of inferential evidence is clearly the next step for this analysis and work is ongoing to develop a bootstrap approach to this.

4.4. Strengths and Limitations of Study

This study utilises a large overall sample and begins to investigate the long-term impact of maternal smoking and alcohol consumption on the facial shape of offspring. Previous studies have assessed the influence of possible dose-dependent relationships and the influence of the individuals’ alcohol susceptibility [25]. It was not possible to assess these factors in this population, using this methodology, due to sample size restrictions. This, however, could be the subject of further work. Further, smoking and alcohol exposures were determined via self-questionnaire, which have conflicting levels of reliability in the literature [35,36]. Furthermore, the image data was reduced to 21 facial landmarks in this study. Further work will focus on the use of many more landmarks through the form of dense facial meshes, such as those used in Muggli et al. [25]. Scaled faces were used here in order to isolate the influence on facial shape, as is common in facial shape research, but this removes the influence on size. Investigating the influence on facial shape and size may be interesting in future. Finally, the influence of protective factors such as folic acid [37] and multivitamins [38] alongside the influence of nicotine replacement therapy, paternal smoking during pregnancy, and smoking exposure after pregnancy would all also be interesting to investigate in future.

5. Conclusions

  • mPCA is a useful descriptive tool for assessing the influence of environmental factors on facial shape.
  • The effect of maternal smoking and maternal alcohol consumption (average 1–2 glasses per week) was minimal, with 0.66% and 0.48% of the variation in the 21 landmarks of non-syndromic offspring explained, respectively.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2313-433X/6/5/34/s1, Figure S1: Scatter plots of the standardized component scores for conventional PCA up to PC10. Group centroids (18 groups) are represented as smaller landmarks, and global centroids are represented as larger landmarks. A subtle pattern in the separation of the group means is visible at PC7/8, with the global centroid of individuals whose mothers did not smoke during pregnancy distant to the global centroid of those whose mothers smoked during pregnancy. No other obvious patterns in the separation of the group means is visualised. Table S1: First 10 eigenvalues for conventional PCA and mPCA (four levels). Estimates of the (standard) errors (SEs) associated with the eigenvalues can be found using Monte Carlo simulation (number of MC replicates = 1000), assuming a multivariate normal distribution for the landmark point data. Table S2: Frequencies and mean smoking levels (smoking, 1st trimester); Table S3: Frequencies and mean smoking levels (smoking, 2nd trimester); Table S4: Frequencies and mean alcohol levels (alcohol, 1st trimester); Table S5: Frequencies and mean alcohol levels (alcohol, 2nd trimester).

Author Contributions

Conceptualization, J.G., D.J.J.F., S.R. and A.I.Z.; Methodology, J.G., D.J.J.F., S.R. and A.I.Z.; Software, D.J.J.F. and A.I.Z.; Validation, J.G., D.J.J.F. and A.I.Z.; Formal Analysis, J.G.; Investigation, J.G.; Resources, D.J.J.F., S.R. and A.I.Z.; Data Curation, J.G.; Writing—Original Draft Preparation, J.G. Writing—Review and Editing, J.G., D.J.J.F., S.R. and A.Z.; Visualization, J.G. Supervision, D.J.J.F., S.R. and A.I.Z.; Project Administration, J.G., D.J.J.F., and S.R.; Funding Acquisition, S.R. All authors have read and agreed to the published version of the manuscript.

Funding

The UK Medical Research Council and Wellcome (grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for the ALSPAC. This publication is the work of the authors, who will serve as guarantors for the contents of this paper. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). Gathering of the facial scans used in this research was funded by Cardiff University. Gathering of the smoking level information was funded by MRC (MR/M006727/1). Gathering of information on maternal alcohol consumption was funded by MRC (MR/L022206/1).

Acknowledgments

We are extremely grateful to all the families who took part in the ALSPAC, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gelman, A.; Hill, J. Data Analysis Using Regression and Multilevel/Hierarchical Models; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  2. Halazonetis, D.J. Morphometrics for cephalometric diagnosis. Am. J. Orthod. Dentofac. 2004, 125, 571–581. [Google Scholar] [CrossRef] [PubMed]
  3. Lecron, F.; Boisvert, J.; Benjelloun, M.; Labelle, H.; Mahmoudi, S. Multilevel statistical shape models: A new framework for modeling hierarchical structures. In Proceedings of the 9th IEEE International Symposium on Biomedical Imaging (ISBI), Barcelona, Spain, 2–5 May 2012; pp. 1284–1287. [Google Scholar]
  4. Farnell, D.J.J.; Popat, H.; Richmond, S. Multilevel principal component analysis (mPCA) in shape analysis: A feasibility study in medical and dental imaging. Comput. Method Prog. Biomed. 2016, 129, 149–159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Farnell, D.J.J.; Galloway, J.; Zhukov, A.; Richmond, S.; Perttiniemi, P.; Katic, V. Initial results of multilevel principal component analysis of facial shape. In Proceedings of the Annual Conference on Medical Image Understanding and Analysis, MIUA 2017, Edinburgh, UK, 11–13 July 2017; Vales, M., Gonzalez-Castro, V., Eds.; Springer: Cham, Switzerland, 2017; Volume 723, pp. 1–12. [Google Scholar]
  6. Farnell, D.J.J.; Galloway, J.; Zhurov, A.I.; Richmond, S.; Perttiniemi, P.; Lähdesmäki, R. What’s in a Smile? Initial Results of Multilevel Principal Components Analysis of Facial Shape and Image Texture. In Proceedings of the Medical Image Understanding and Analysis, MIUA 2018, Southampton, UK, 9–11 July 2018; Springer: Cham, Switzerland, 2018; Volume 894, pp. 177–188. [Google Scholar]
  7. Farnell, D.J.J.; Galloway, J.; Zhurov, A.I.; Richmond, S.; Marshall, P.L.; Rosin, K.; Al-Meyah, K.; Pirttiniemi, P.; Lahdesmaki, R. What’s in a smile? Initial analyses of dynamic changes in facial shape and appearance. J. Imaging 2019, 5, 2. [Google Scholar] [CrossRef] [Green Version]
  8. Farnell, D.J.J.; Galloway, J.; Zhurov, A.I.; Richmond, S. Multilevel models of age-related changes in facial shape in adolescents. In Proceedings of the Medical Image Understanding and Analysis, MIUA 2019, Liverpool, UK, 24–26 July 2019; Zheng, Y., Williams, B., Chen, K., Eds.; Springer: Cham, Switzerland, 2019; Volume 1065, pp. 101–113. [Google Scholar]
  9. Farnell, D.J.J.; Richmond, S.; Galloway, J.; Zhurov, A.I.; Pirttiniemi, P.; Heikkinen, T.; Harila, V.; Matthews, H.; Claes, P. Multilevel principal components analysis of three-dimensional facial growth in adolescents. Comput. Method Prog. Biomed. 2020, 188, 105272. [Google Scholar] [CrossRef]
  10. Timmerman, M.E. Multilevel component analysis. Brit. J. Math. Stat. Psychol. 2006, 59, 301–320. [Google Scholar] [CrossRef]
  11. Harnquist, K. Primary mental abilities at collective and individual levels. J. Educ. Psychol. 1978, 70, 706–716. [Google Scholar]
  12. STATA FAQ: How Can I Do Multilevel Principal Components Analysis? Available online: https://stats.idre.ucla.edu/stata/faq/how-can-i-do-multilevel-principal-components-analysis/ (accessed on 14 April 2020).
  13. Socol, M.L.; Manning, F.A.; Muata, Y.; Druzin, M.L. Maternal smoking causes fetal hypoxia: Experimental evidence. AJOG 1982, 142, 214–218. [Google Scholar] [CrossRef]
  14. Richmond, R.C.; Simpkin, A.J.; Woodward, G.; Gaunt, T.R.; Lyttleton, O.; McArdle, W.L.; Ring, S.M.; Smith, A.D.A.C.; Timpson, N.J.; Tilling, K.; et al. Prenatal exposure to maternal smoking and offspring DNA methylation across the lifecourse: Findings from the Avon Longitudinal Study of Parents and Children. Hum. Mol. 2015, 24, 2201–2217. [Google Scholar] [CrossRef] [Green Version]
  15. Rauschert, S.; Melton., P.E.; Burdge, G.; Craig, J.M.; Godfrey, K.M.; Holbrook, J.D.; Lillycrop, K.; Mori, T.A.; Beilin, L.J.; Oddy, W.H.; et al. Maternal smoking during pregnancy induces persistent epigenetic changes into adolescence, independent of postnatal smoke exposure and is associated with cardiometabolic risk. Front. Genet. 2019, 10, 770. [Google Scholar] [CrossRef]
  16. Hackshaw, A.; Rodeck, C.; Boniface, S. Maternal smoking in pregnancy and birth defects: A systematic review based on 173 687 malformed cases and 11.7 million controls. Hum. Reprod. Update 2011, 17, 589–604. [Google Scholar] [CrossRef] [Green Version]
  17. Kallen, K. Maternal smoking and craniosynostosis. Teratology 1999, 60, 146–150. [Google Scholar] [CrossRef]
  18. Kallen, K. Maternal smoking during pregnancy and infant head circumference at birth. Early Hum. Dev. 2000, 28, 197–204. [Google Scholar] [CrossRef]
  19. Wang, X.; Tager, I.B.; Van Vunakis, H.; Speizer, F.E.; Hanrahan, J.P. Maternal smoking during pregnancy, urine cotinine concentrations, and birth outcomes. A prospective cohort study. Int. J. Epidemiol. 1997, 26, 978–988. [Google Scholar] [CrossRef] [PubMed]
  20. Lindley, A.A.; Becker, S.; Gray, R.H.; Herman, A.A. Effect of continuing or stopping smoking during pregnancy on infant birth weight, crown-heel length, head circumference, pinderal index, and brain:body weight ratio. Am. J. Epidemiol. 2000, 152, 219–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Werler, M.M.; Starr, J.R.; Cloonan, Y.K.; Speltz, M.L. Hemifacial microsomia: From gestation to childhood. J. Craniofac. Surg. 2009, 20, 664–669. [Google Scholar] [CrossRef]
  22. Little, J.; Cardy, A.; Munger, R.G. Tobacco smoking and oral clefts: A meta-analysis. Bull. World Health Organ. 2004, 82, 213–218. [Google Scholar]
  23. Astley, S.J.; Sterling, K.; Clarren, M.D.; Little, R.E.; Sampson, P.D.; Daling, J.R. Analysis of facial shape in children gestationally exposed to marijuana, alcohol, and/or cocaine. Pediatrics 1992, 89, 67–77. [Google Scholar]
  24. Jones, K.L.; Smith, D.W. Recognition of the fetal alcohol syndrome in early infancy. Lancet 1973, 302, 999–1001. [Google Scholar] [CrossRef]
  25. Muggli, E.; Matthews, H.; Penington, A.; Claes, P.; O’Leary, C.; Forster, D.; Donath, S.; Anderson, P.J.; Lewis, S.; Nagle, C.; et al. Association between prenatal alcohol exposure and craniofacial shape of children at 12 months of age. JAMA Pediatr. 2017, 171, 771–780. [Google Scholar] [CrossRef] [Green Version]
  26. Howe, J.L.; Sharp, G.C.; Hemani, G.; Zuccolo, L.; Richmond, S.; Lewis, S.J. Prenatal alcohol exposure and facial morphology in a UK cohort. Drug Alcohol Depend. 2019, 197, 42–47. [Google Scholar] [CrossRef] [Green Version]
  27. Boyd, A.; Golding, J.; Macleod, J.; Lawlor, D.A.; Fraser, A.; Henderson, J.; Molloy, L.; Ness, A.; Ring, S.; Davey Smith, G. Cohort profile: The ‘Children of the 90s’; the index offspring of The Avon Longitudinal Study of Parents and Children (ALSPAC). Int. J. Epidemiol. 2013, 42, 111–127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Fraser, A.; Macdonald-Wallis, C.; Tilling, K.; Boyd, A.; Golding, J.; Davey Smith, G.; Henderson, J.; Macleod, J.; Molloy, L.; Ness, A.; et al. Cohort profile: The Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. Int. J. Epidemiol. 2013, 42, 97–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kau, C.H.; Richmond, S.; Zhurov, A.I.; Knox, J.; Chestnutt, I.; Hartles, F.; Playle, R. Reliability of measuring facial morphology with a 3-dimensional laser scanning system. Am. J. Orthod. Dentofac. Orthop. 2005, 128, 424–430. [Google Scholar] [CrossRef] [PubMed]
  30. Zhurov, A.; Richmond, S.; Kau, C.H.; Toma, A. Averaging facial images. In Three-Dimensional Imaging for Orthodontics and Maxillofacial Surgery; Kau, C.H., Richmond, S., Eds.; Wiley-Blackwell: Hoboken, NJ, USA, 2010; pp. 126–144. [Google Scholar]
  31. Toma, A.M.; Zhurov, A.I.; Playle, R.; Marshall, D.; Rosin, P.L.; Richmond, S. The assessment of facial variation in 4747 British school children. Eur. J. Orthod. 2012, 34, 655–664. [Google Scholar] [CrossRef] [PubMed]
  32. Farkas, L.G. Anthropometry of the head and face in clinical practice. In Anthropometry of the Head and Face; Farkas, L.G., Ed.; Raven Press: New York, NY, USA, 1994. [Google Scholar]
  33. Toma, A.M.; Zhurov, A.; Playle, R.; Ong, E.; Richmond, S. Reproducibility of facial soft tissue landmarks on 3D laser-scanned facial images. Orthod. Craniofac. Res. 2009, 12, 33–42. [Google Scholar] [CrossRef]
  34. Koziel, S.; Zadzinska, E.; Gomula, A. Parental smoking during pregnancy and head shape and size in school children. Ann. Hum. Biol. 2018, 45, 401–405. [Google Scholar] [CrossRef]
  35. Kvalvik, L.G.; Nilsen, R.M.; Skjaerven, R.; Vollset, S.E.; Midttun, O.; Ueland, P.M.; Haug, K. Self-reported smoking status and plasma cotinine concentrations among pregnant women in the Norwegian Mother and Child Cohort Study. Pediatr. Res. 2012, 72, 101–107. [Google Scholar] [CrossRef]
  36. Shipton, D.; Tappin, D.M.; Vadiveloo, T.; Crosslet, J.A.; Aitken, D.A.; Chalmers, J. Reliability of self reported smoking status by pregnant women for estimating smoking prevalence: A retrospective, cross sectional study. BMJ 2009, 339, b4347. [Google Scholar] [CrossRef] [Green Version]
  37. Carmichael, S.L.; Ma, C.; Rasmussen, S.A.; Honein, M.A.; Lammer, E.J.; Shaw, G.M. The National Birth Defects Prevention Study Craniosynostosis and maternal smoking. Birth Defects Res. (Part A) 2008, 82, 78–85. [Google Scholar] [CrossRef]
  38. Shaw, G.M.; Nelson, V.; Carmichael, S.L.; Lammer, E.J.; Finnell, R.H.; Rosenquist, T.H. Maternal periconceptional vitamins: Interactions with selected factors and congenital anomalies? Epidemiology 2002, 13, 625–630. [Google Scholar] [CrossRef]
Figure 1. The 21 landmarks as described by Farkas [32]: (1) Glabella, (2) Nasion, (3) Endocanthion (Left), (4) Endocanthion (Right), (5) Exocanthion (Left), (6) Exocanthion (Right), (7) Palpebrale Superius (Left), (8) Palpebrale Superius (Right), (9) Palpebrale Inferius (Left), (10) Palpebrale Inferius (Right), (11) Pronasale, (12) Subnasale, (13) Alare (Left), (14) Alare (Right), (15) Labiale Superius, (16) Crista Phitri (Left), (17) Crista Philtri (Right), (18) Labiale Inferius, (19) Cheilion (Left), (20) Cheilion (Right), and (21) Pogonion.
Figure 1. The 21 landmarks as described by Farkas [32]: (1) Glabella, (2) Nasion, (3) Endocanthion (Left), (4) Endocanthion (Right), (5) Exocanthion (Left), (6) Exocanthion (Right), (7) Palpebrale Superius (Left), (8) Palpebrale Superius (Right), (9) Palpebrale Inferius (Left), (10) Palpebrale Inferius (Right), (11) Pronasale, (12) Subnasale, (13) Alare (Left), (14) Alare (Right), (15) Labiale Superius, (16) Crista Phitri (Left), (17) Crista Philtri (Right), (18) Labiale Inferius, (19) Cheilion (Left), (20) Cheilion (Right), and (21) Pogonion.
Jimaging 06 00034 g001
Figure 2. A flow chart of participant inclusion.
Figure 2. A flow chart of participant inclusion.
Jimaging 06 00034 g002
Figure 3. Eigenvalue plot allowing visualisation of the eigenvalue magnitude for each principal component, at each level of the model.
Figure 3. Eigenvalue plot allowing visualisation of the eigenvalue magnitude for each principal component, at each level of the model.
Jimaging 06 00034 g003
Figure 4. Scatter plots of the standardised component scores at PC1 and PC2 for conventional PCA (a,b,d) and mPCA (c,e). No obvious pattern in the separation of the group means is evident in conventional PCA PC1/2 (a). There is a suggestion of a pattern in conventional PCA PC7/8 (b) and mPCA smoking level PC2 (c). Subjects whose mothers did not smoke during pregnancy are possibly separated from those whose mothers smoked during the 1st trimester or both the 1st and 2nd trimesters. There is no obvious pattern in the separation of the group means with alcohol consumption during pregnancy (d,e). Scatter plots for conventional PCA PC3–10 are available as supplementary material.
Figure 4. Scatter plots of the standardised component scores at PC1 and PC2 for conventional PCA (a,b,d) and mPCA (c,e). No obvious pattern in the separation of the group means is evident in conventional PCA PC1/2 (a). There is a suggestion of a pattern in conventional PCA PC7/8 (b) and mPCA smoking level PC2 (c). Subjects whose mothers did not smoke during pregnancy are possibly separated from those whose mothers smoked during the 1st trimester or both the 1st and 2nd trimesters. There is no obvious pattern in the separation of the group means with alcohol consumption during pregnancy (d,e). Scatter plots for conventional PCA PC3–10 are available as supplementary material.
Jimaging 06 00034 g004
Figure 5. Scatter plots of the standardised component scores at PC1 and PC2 for conventional PCA (a), PC1 for sex level mPCA (b) and PC1 and PC2 for subject level mPCA (c). Clear separation of the group means of the biological sexes is evident for both conventional PCA and mPCA. At a subject level, the centroids seem to be centred around the origin, although there is some deviation from this, perhaps due to group sample sizes.
Figure 5. Scatter plots of the standardised component scores at PC1 and PC2 for conventional PCA (a), PC1 for sex level mPCA (b) and PC1 and PC2 for subject level mPCA (c). Clear separation of the group means of the biological sexes is evident for both conventional PCA and mPCA. At a subject level, the centroids seem to be centred around the origin, although there is some deviation from this, perhaps due to group sample sizes.
Jimaging 06 00034 g005
Table 1. Subject demographics 1.
Table 1. Subject demographics 1.
GroupingSexSample SizeAgeSmoking Levels
Mean +/- SD (Cigarettes/Day)
Alcohol Levels
Mean +/- SD (Glasses/Week) 2
1st Trimester2nd Trimester1st Trimester2nd Trimester
No smoking
No alcohol
Male57415.41 +/- 0.25----
Female62215.42 +/- 0.28----
No smoking
Alcohol 1st trimester only
Male21415.42 +/- 0.26--0.82 +/- 2.91-
Female27515.40 +/- 0.25--0.68 +/- 0.90-
No smoking
Alcohol 1st and 2nd trimesters
Male70215.42 +/- 0.27--1.17 +/- 3.401.22 +/- 3.42
Female77515.44 +/- 0.30--1.02 +/- 2.041.10 +/- 2.87
Smoking 1st trimester only
No alcohol
Male2315.51 +/- 0.307.10 +/- 7.16---
Female1815.60 +/- 0.347.19 +/- 5.35---
Smoking 1st trimester only
Alcohol 1st trimester only
Male1515.33 +/- 0.145.50 +/- 4.27-2.07 +/- 3.46-
Female1815.56 +/- 0.459.26 +/- 6.28-1.24 +/- 2.22-
Smoking 1st trimester only
Alcohol 1st and 2nd trimester
Male4115.40 +/- 0.257.24 +/- 6.30-1.16 +/- 1.560.84 +/- 0.47
Female4415.53 +/- 0.427.77 +/- 6.74-2.26 +/- 6.401.25 +/- 2.07
Smoking 1st and 2nd trimesters
No alcohol
Male6815.44 +/- 0.2810.48 +/- 6.2710.47 +/- 7.40--
Female7215.47 +/- 0.3610.88 +/- 6.9011.47 +/- 6.34--
Smoking 1st and 2nd trimesters
Alcohol 1st trimester only
Male2315.41 +/- 0.2810.70 +/- 7.068.48 +/- 5.940.67 +/- 0.38-
Female3515.44 +/- 0.3211.19 +/- 5.739.56 +/- 6.330.70 +/- 0.40-
Smoking 1st and 2nd trimesters
Alcohol 1st and 2nd trimesters
Male11515.40 +/- 0.2110.54 +/- 6.819.66 +/- 6.561.97 +/- 5.631.20 +/- 1.83
Female12115.42 +/- 0.2511.26 +/- 7.6410.81 +/- 7.462.29 +/- 5.811.28 +/- 2.15
1 The raw data received from the Avon Longitudinal Study of Parents and Children (ALSPAC) had been previously categorised. The sample means were therefore calculated using: mean = sum(frequency x mid interval point of category)/sum frequency. The frequency tables are available as supplementary material (Supplementary Tables S2–S5). 2 The mothers were asked to take “glass” to mean “a pub measure of spirits, half a pint of lager or cider, a glass of wine etc.”.
Table 2. Percentage of the total variation in the 21 facial landmarks as explained by each level.
Table 2. Percentage of the total variation in the 21 facial landmarks as explained by each level.
Level% of Total Variation in 21 Landmarks
PC1PC2All Retained PCs
Level 1: Maternal Smoking0.38%0.28%0.66%
Level 2: Maternal Alcohol Consumption0.31%0.17%0.48%
Level 3: Biological Sex8.88%N/A8.88%
Level 4: Subject16.75%9.99%85.85%
Table 3. Differences in each landmark due to maternal smoking level PC2. These results should be interpreted with caution at this stage. They should be taken as a suggestion of which landmarks are likely to vary due to maternal smoking rather than providing exact values.
Table 3. Differences in each landmark due to maternal smoking level PC2. These results should be interpreted with caution at this stage. They should be taken as a suggestion of which landmarks are likely to vary due to maternal smoking rather than providing exact values.
LandmarkDifference in mm between
Average Face + 2* Square Root Eigenvalue and Average Face − 2* Square Root Eigenvalue
xyz
Glabella−0.08−0.24−0.53
Nasion−0.02−0.84−0.22
Endocanthion (Left)−0.76−0.110.20
Endocanthion (Right)0.880.040.05
Exocanthion (Left)0.410.350.27
Exocanthion (Right)−0.070.320.02
Palpebrale Superius (Left)−0.510.45−0.15
Palpebrale Superius (Right)0.500.30−0.23
Palpebrale Inferius (Left)−0.470.15−0.09
Palpebrale Inferius (Right)0.540.110.04
Pronasale−0.110.860.70
Subnasale−0.140.600.08
Alare (Left)−0.490.070.18
Alare (Right)0.290.200.12
Labiale Superius−0.07−0.30−0.24
Labiale Inferius−0.07−0.200.56
Crista Phitri (Left)−0.13−0.42−0.27
Crista Phitri (Right)−0.07−0.33−0.20
Cheilion (Left)0.25−0.62−0.65
Cheilion (Right)−0.04−0.65−0.37
Pogonion0.180.250.50
Table 4. Mean standardised component score and standard deviation for each grouping.
Table 4. Mean standardised component score and standard deviation for each grouping.
GroupingSexSmoking Level mPCAAlcohol Level mPCA
PC1SDPC2SDPC1SDPC2SD
No smoking
No alcohol
Male−0.3595.9970.1955.104−0.2383.302−0.1283.720
Female−0.1845.9180.2774.7340.0123.360−0.3823.381
No smoking
Alcohol 1st trimester only
Male0.5915.7011.2554.758−0.2523.365−0.3423.606
Female0.3396.2220.7524.324−0.1013.658−0.1173.607
No smoking
Alcohol 1st and 2nd trimesters
Male0.3465.8520.2274.968−0.1563.605−0.0993.709
Female−0.1675.8560.3484.894−0.2363.689−0.4443.396
Smoking 1st trimester only
No alcohol
Male1.4334.6910.0763.5100.0162.796−0.4333.262
Female−0.7585.2982.2914.108−0.6562.8970.5653.842
Smoking 1st trimester only
Alcohol 1st trimester only
Male0.8734.711−0.9923.2650.6972.7420.8874.030
Female−0.4616.934−1.1993.381−1.0513.456−0.7623.747
Smoking 1st trimester only
Alcohol 1st and 2nd trimester
Male−0.5876.3010.3755.091−0.6772.9400.2853.411
Female1.2315.813−0.4315.318−0.0793.1111.2612.944
Smoking 1st and 2nd trimesters
No alcohol
Male−0.3324.674−0.2494.717−0.2473.411−0.6783.841
Female1.2916.3940.2834.985−0.5113.141−0.3262.925
Smoking 1st and 2nd trimesters
Alcohol 1st trimester only
Male−0.4726.068−0.5015.5690.9254.009−0.1923.597
Female−0.9155.861−0.3925.124−0.5193.733−0.4413.825
Smoking 1st and 2nd trimesters
Alcohol 1st and 2nd trimesters
Male0.6095.7880.7744.641−0.2233.489−0.5753.266
Female−1.1955.121−0.1464.6630.0823.534−0.8433.265

Share and Cite

MDPI and ACS Style

Galloway, J.; Farnell, D.J.J.; Richmond, S.; Zhurov, A.I. Multilevel Analysis of the Influence of Maternal Smoking and Alcohol Consumption on the Facial Shape of English Adolescents. J. Imaging 2020, 6, 34. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging6050034

AMA Style

Galloway J, Farnell DJJ, Richmond S, Zhurov AI. Multilevel Analysis of the Influence of Maternal Smoking and Alcohol Consumption on the Facial Shape of English Adolescents. Journal of Imaging. 2020; 6(5):34. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging6050034

Chicago/Turabian Style

Galloway, Jennifer, Damian J.J. Farnell, Stephen Richmond, and Alexei I. Zhurov. 2020. "Multilevel Analysis of the Influence of Maternal Smoking and Alcohol Consumption on the Facial Shape of English Adolescents" Journal of Imaging 6, no. 5: 34. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging6050034

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop